[gradsusr] Median for ensembles.

[Ryglicki, David CIV FNMOC]

Hello, everyone.

The frame for this question is plotting a median line (as opposed to a mean) for a box-whisker plot. The logic surrounding the plot, I get, but my question is specifically about the middle of my script.

The way I have it now, I set lat/lon/height, loop the ensemble members in a while loop, and store the value at that point for each ensemble in a dummy variable (emem.ee). I then run that through a bubble sort, figure out if there's an odd or even number of ensemble members, do the math, and boom, we have a median - med.tt, with tt being the time index.

My question is how do I get the variables in med.tt into a displayable grads variable - say, "median" - so that when I'm finally getting ready to plot, I can just type in "d median?" Is this even possible via some combination of tloop and define? Should I remove this process from GrADS entirely and compute it in a Fortran program?

Additionally, if anyone out there in gradsusr land has a script for computing and plotting median in a line graph, I would not say no to that. I couldn't find one in several attempted Google searches.

David Ryglicki, Ph.D.
N3 Models
Fleet Numerical Meteorology and Oceanography Center
Bldg 702, Rm 146
Office: 831|656|4519


-------------- next part --------------
A non-text attachment was scrubbed...
Name: smime.p7s
Type: application/x-pkcs7-signature
Size: 5686 bytes
Desc: not available
Url : http://gradsusr.org/pipermail/gradsusr/attachments/20120507/a97f4fb2/attachment-0003.bin