<div dir="ltr"><div>Hello,<br>I have some netCDF files that I created using my own FORTRAN code. I wish to display some arrays in those files that have some non-world-coordinate dimensions. For example, my data contain probabilities of precipitation from three different ensembles using five different precipitation accumulation thresholds. Therefore, I have two non-world-coordinate dimensions. I see in the documentation how to handle one (see below):<br><p class=""><span style="font-family:georgia,serif"><i>If your data file contains a variable that also
varies in a non-world-coordinate dimension (e.g. histogram interval,
spectral band, ensemble number) then you can put a non-negative
integer in the list of varying dimensions that will become the array
index of the extra dimension. For example:
</i></span></p><ul><p class=""><span style="font-family:georgia,serif"><i> VAR=>hist0 0 0,y,x
First historgram interval for VAR<br>
VAR=>hist1 0 1,y,x Second
historgram interval for VAR<br>
VAR=>hist2 0 2,y,x Third
histogram interval for VAR </i></span></p></ul><span style="font-family:georgia,serif"><i>
</i></span><p class=""><span style="font-family:georgia,serif"><i>Another option in this example would be to fill
the unused Z axis with the histogram intervals: </i></span></p><span style="font-family:georgia,serif"><i>
</i></span><ul><p class=""><span style="font-family:georgia,serif"><i> zdef 3 linear 1 1<br>
... <br>
VAR=>hist 3 z,y,x VAR Histogram</i></span></p></ul><span style="font-family:georgia,serif"><i>
</i></span><p class=""><span style="font-family:georgia,serif"><em>In this case, it would appear to GrADS that variable
'hist' varies in Z, but the user would have to remember that the
Z levels correspond to histogram intervals. The latter technique
makes it easier to slice through the data, but is not the most accurate
representation. And if you don't have an unsued world-coordinate
axis available, then you still have a way to access your data.</em></span></p>However, I am unsure of how to handle two. Let alone, if you have no "free" world coordinates, how do you reference these other dimensions when displaying the grids in Grads?<br><br></div>Here is an ncdump of one of my files:<br><br>netcdf \0527_neighborhood_PoP_007200 {<br>dimensions:<br> X = 803 ;<br> Y = 603 ;<br> threshold = 5 ;<br> ensemble = 3 ;<br>variables:<br> float forecast_PoP(ensemble, threshold, Y, X) ;<br> float observed_PoP(threshold, Y, X) ;<br>}<br><br clear="all"><div><div><div><div>Jeff Duda<br></div><div>-- <br>Jeff Duda<br>Graduate research assistant<br>University of Oklahoma School of Meteorology<br>Center for Analysis and Prediction of Storms<br>
</div></div></div></div></div>