Apply function to Dataframe based on Multi-level Grouping in R

Question

I am trying to apply a function to a dataframe to add a column which calculates the percentile rank for each record based on Weather Station ID (WSID) and Season Grouping.

<