<div dir="ltr"><div>Jennifer,<br>Thank you for the reply. I had tried something very similar to that, but it didn't work. Now I can't remember what that was. Regardless, it works now.<br><br></div>Thanks!<br><br>Jeff<br></div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Sep 24, 2014 at 6:56 AM, Jennifer Adams <span dir="ltr"><<a href="mailto:jma@cola.iges.org" target="_blank">jma@cola.iges.org</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div style="word-wrap:break-word">Hi, Jeff — <div>Your data file has 4 dimensions, so you can map your data’s dimensions to the world coordinate grid dimensions and not worry about putting array indices in your variable declarations. May I suggest the following:</div><div><br></div><div>xdef 803 linear 1 1</div><div>ydef 603 linear 1 1</div><div>zdef 5 linear 1 1</div><div>edef 3 names 1 2 3</div><div>tdef 1 linear 01jan0001 1mo</div><div>vars 2</div><div>forecast_PoP=>fcst 5 e,z,y,x forecast PoP</div><div><div dir="ltr">observed_PoP=>obs 5 z,y,x observed PoP</div></div><div dir="ltr">endvars</div><div><br></div><div>It is more confusing to use T for a dimension that isn’t really a time coordinate, so I have left that one unused, with a default time and a size of 1. </div><div>—Jennifer</div><div><br></div><div><br></div><div><br><div><div><div class="h5"><div>On Sep 23, 2014, at 11:42 AM, Jeff Duda <<a href="mailto:jeffduda319@gmail.com" target="_blank">jeffduda319@gmail.com</a>> wrote:</div><br></div></div><blockquote type="cite"><div><div class="h5"><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><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><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><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><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><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>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>
_______________________________________________<br>gradsusr mailing list<br><a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a><br><a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br></blockquote></div><br><div>
<span style="border-collapse:separate;border-spacing:0px"><span style="border-collapse:separate;color:rgb(0,0,0);font-family:Helvetica;font-size:12px;font-style:normal;font-variant:normal;font-weight:normal;letter-spacing:normal;line-height:normal;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px"><div>--</div><div>Jennifer M. Adams<br>Center for Ocean-Land-Atmosphere Studies (COLA)<br>111 Research Hall, Mail Stop 2B3<br>George Mason University<br>4400 University Drive<br>Fairfax, VA 22030 <br><br></div><div><br></div><br></span></span><br>
</div>
<br></div></div><br>_______________________________________________<br>
gradsusr mailing list<br>
<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
<a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
<br></blockquote></div><br><br clear="all"><br>-- <br>Jeff Duda<br>Graduate research assistant<br>University of Oklahoma School of Meteorology<br>Center for Analysis and Prediction of Storms<br>
</div>