<div dir="ltr">Hello,<div><br></div><div>I actually ran it on a westgrid system as a batch job with the command:</div><div><br></div><div><p style="margin:0px;font-size:11px;font-family:Menlo">qsub -l procs=90,pmem=2000mb,walltime=01:00:00 run_MOM_solo.pbs</p>
<p style="margin:0px;font-size:11px;font-family:Menlo"><br></p><p style="margin:0px"><font face="arial, helvetica, sans-serif">The file run_MOM_solo.pbs is attached as well as the MOM_run.csh file that it calls. </font></p>
<p style="margin:0px"><font face="arial, helvetica, sans-serif"><br></font></p><p style="margin:0px"><font face="arial, helvetica, sans-serif">Thanks,</font></p><p style="margin:0px"><font face="arial, helvetica, sans-serif"><br>
</font></p><p style="margin:0px"><font face="arial, helvetica, sans-serif">David</font></p></div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Mon, Jul 28, 2014 at 1:25 PM, Charles Seman - NOAA Federal <span dir="ltr">&lt;<a href="mailto:charles.seman@noaa.gov" target="_blank">charles.seman@noaa.gov</a>&gt;</span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">David,<br>
<br>
How did you get the data?  Did you run the GFDL model locally?  If so,<br>
did you use &quot;frerun&quot; with an xml file?  If so, there is a companion<br>
&quot;frepp&quot; command which can be used to post-process the history data.<br>
Could you provide more details on the source of your file?<br>
<br>
Thanks,<br>
Chuck<br>
<div class=""><br>
On 07/26/2014 07:38 PM, David Stanfel wrote:<br>
&gt; It is GFDL Model data using the Molecular Ocean Model MOM_solo type.<br>
&gt; This is for the box1 test case. I did try sdfopen and it came out with<br>
&gt; the following error message:<br>
&gt;<br>
&gt; Scanning self-describing file: <a href="http://19800101.ocean_neutral.nc" target="_blank">19800101.ocean_neutral.nc</a><br>
</div>&gt; &lt;<a href="http://19800101.ocean_neutral.nc" target="_blank">http://19800101.ocean_neutral.nc</a>&gt;<br>
<div class="">&gt;<br>
&gt; SDF Error: 365 day calendars are no longer supported by sdfopen.<br>
&gt;<br>
&gt;    To open this file with GrADS, use a descriptor file with<br>
&gt;<br>
&gt;    a complete TDEF entry and OPTIONS 365_day_calendar.<br>
&gt;<br>
&gt;<br>
&gt; I do not have any pp data and I don&#39;t know how to go about creating the<br>
&gt; grid data. How could I do that, putting it into lat, long coordinates?<br>
&gt;<br>
&gt; Thanks,<br>
&gt;<br>
&gt; David<br>
&gt;<br>
&gt;<br>
&gt; On Fri, Jul 25, 2014 at 2:39 PM, Charles Seman - NOAA Federal<br>
</div><div class="">&gt; &lt;<a href="mailto:charles.seman@noaa.gov">charles.seman@noaa.gov</a> &lt;mailto:<a href="mailto:charles.seman@noaa.gov">charles.seman@noaa.gov</a>&gt;&gt; wrote:<br>
&gt;<br>
&gt;     David,<br>
&gt;<br>
&gt;     This file looks like it&#39;s from a GFDL model?  Based on the names of the<br>
&gt;     coordinates, and the global attributes &quot;grid_type&quot; and &quot;grid_tile&quot;, it<br>
&gt;     looks like it is a &quot;history&quot; file which is on the model grid and has yet<br>
&gt;     to be regridded to a lat/lon grid.  I&#39;m not a user of the ocean data,<br>
&gt;     but am working with a GFDL history data file from the atmospheric model<br>
&gt;     on the cubed-sphere model grid tile4.  Be very careful!  The grid<br>
&gt;     coordinates xt_ocean and yt_ocean may need to be mapped into actual<br>
&gt;     global coordinates.  To get a view of the data, have you tried opening<br>
&gt;     the file in GrADS without using a GrADS ctl file?  GrADS is able to open<br>
&gt;     and view the tile4 file I&#39;m working with, but one must be aware that the<br>
&gt;     map background that GrADS assigns is not physically correct.  Can you<br>
&gt;     use post-processed &quot;pp&quot; data (if you have this, it should be in a &quot;pp&quot;<br>
&gt;     subdirectory)?  These &quot;pp&quot; data have been put on a physical grid.<br>
&gt;     I talked with someone here at GFDL who works with ocean data, and he<br>
&gt;     says Ferret is able to open the datasets he works with...<br>
&gt;<br>
&gt;     Here&#39;s an example of what&#39;s in a pp/ocean/<a href="http://ocean.static.nc" target="_blank">ocean.static.nc</a><br>
</div>&gt;     &lt;<a href="http://ocean.static.nc" target="_blank">http://ocean.static.nc</a>&gt; file.  Note<br>
<div><div class="h5">&gt;     how the tracer latitude and longitudes are functions of<br>
&gt;     yt_ocean,xt_ocean: geolat_t(yt_ocean,xt_ocean) and<br>
&gt;     geolon_t(yt_ocean,xt_ocean).  So the lat and lon are themselves<br>
&gt;     functions of grid coordinate variables. Here&#39;s an example of<br>
&gt;     two-dimensional lat, lon coordinates:<br>
&gt;     <a href="http://cfconventions.org/Data/cf-convetions/cf-conventions-1.6/build/cf-conventions.html#idp5559280" target="_blank">http://cfconventions.org/Data/cf-convetions/cf-conventions-1.6/build/cf-conventions.html#idp5559280</a><br>

&gt;<br>
&gt;     netcdf ocean.static {<br>
&gt;<br>
&gt;     dimensions:<br>
&gt;<br>
&gt;               xu_ocean = 360 ;<br>
&gt;<br>
&gt;               yu_ocean = 200 ;<br>
&gt;<br>
&gt;               xt_ocean = 360 ;<br>
&gt;<br>
&gt;               yt_ocean = 200 ;<br>
&gt;<br>
&gt;     variables:<br>
&gt;<br>
&gt;               double xt_ocean(xt_ocean) ;<br>
&gt;<br>
&gt;                       xt_ocean:long_name = &quot;tcell longitude&quot; ;<br>
&gt;<br>
&gt;                       xt_ocean:units = &quot;degrees_E&quot; ;<br>
&gt;<br>
&gt;                       xt_ocean:cartesian_axis = &quot;X&quot; ;<br>
&gt;<br>
&gt;               double yt_ocean(yt_ocean) ;<br>
&gt;<br>
&gt;                       yt_ocean:long_name = &quot;tcell latitude&quot; ;<br>
&gt;<br>
&gt;                       yt_ocean:units = &quot;degrees_N&quot; ;<br>
&gt;<br>
&gt;                       yt_ocean:cartesian_axis = &quot;Y&quot; ;<br>
&gt;<br>
&gt;               double xu_ocean(xu_ocean) ;<br>
&gt;<br>
&gt;                       xu_ocean:long_name = &quot;ucell longitude&quot; ;<br>
&gt;<br>
&gt;                       xu_ocean:units = &quot;degrees_E&quot; ;<br>
&gt;<br>
&gt;                       xu_ocean:cartesian_axis = &quot;X&quot; ;<br>
&gt;<br>
&gt;               double yu_ocean(yu_ocean) ;<br>
&gt;<br>
&gt;                       yu_ocean:long_name = &quot;ucell latitude&quot; ;<br>
&gt;<br>
&gt;                       yu_ocean:units = &quot;degrees_N&quot; ;<br>
&gt;<br>
&gt;                       yu_ocean:cartesian_axis = &quot;Y&quot; ;<br>
&gt;<br>
&gt;               float ht(yt_ocean, xt_ocean) ;<br>
&gt;<br>
&gt;                       ht:long_name = &quot;ocean depth on t-cells&quot; ;<br>
&gt;<br>
&gt;                       ht:units = &quot;m&quot; ;<br>
&gt;<br>
&gt;                       ht:valid_range = -1.e+09f, 1.e+09f ;<br>
&gt;<br>
&gt;                       ht:missing_value = -1.e+20f ;<br>
&gt;<br>
&gt;                       ht:_FillValue = -1.e+20f ;<br>
&gt;<br>
&gt;                       ht:cell_methods = &quot;time: point&quot; ;<br>
&gt;<br>
&gt;                       ht:coordinates = &quot;geolon_t geolat_t&quot; ;<br>
&gt;<br>
&gt;                       ht:standard_name = &quot;sea_floor_depth_below_geoid&quot; ;<br>
&gt;<br>
&gt;               float hu(yu_ocean, xu_ocean) ;<br>
&gt;<br>
&gt;                       hu:long_name = &quot;ocean depth on u-cells&quot; ;<br>
&gt;<br>
&gt;                       hu:units = &quot;m&quot; ;<br>
&gt;<br>
&gt;                       hu:valid_range = -1.e+09f, 1.e+09f ;<br>
&gt;<br>
&gt;                       hu:missing_value = -1.e+20f ;<br>
&gt;<br>
&gt;                       hu:_FillValue = -1.e+20f ;<br>
&gt;<br>
&gt;                       hu:cell_methods = &quot;time: point&quot; ;<br>
&gt;<br>
&gt;                       hu:coordinates = &quot;geolon_c geolat_c&quot; ;<br>
&gt;<br>
&gt;               float geolat_c(yu_ocean, xu_ocean) ;<br>
&gt;<br>
&gt;                       geolat_c:long_name = &quot;uv latitude&quot; ;<br>
&gt;<br>
&gt;                       geolat_c:units = &quot;degrees_N&quot; ;<br>
&gt;<br>
&gt;                       geolat_c:valid_range = -91.f, 91.f ;<br>
&gt;<br>
&gt;                       geolat_c:missing_value = 1.e+20f ;<br>
&gt;<br>
&gt;                       geolat_c:_FillValue = 1.e+20f ;<br>
&gt;<br>
&gt;                       geolat_c:cell_methods = &quot;time: point&quot; ;<br>
&gt;<br>
&gt;                       geolat_c:coordinates = &quot;geolon_c geolat_c&quot; ;<br>
&gt;<br>
&gt;               float geolon_c(yu_ocean, xu_ocean) ;<br>
&gt;<br>
&gt;                       geolon_c:long_name = &quot;uv longitude&quot; ;<br>
&gt;<br>
&gt;                       geolon_c:units = &quot;degrees_E&quot; ;<br>
&gt;<br>
&gt;                       geolon_c:valid_range = -281.f, 361.f ;<br>
&gt;<br>
&gt;                       geolon_c:missing_value = 1.e+20f ;<br>
&gt;<br>
&gt;                       geolon_c:_FillValue = 1.e+20f ;<br>
&gt;<br>
&gt;                       geolon_c:cell_methods = &quot;time: point&quot; ;<br>
&gt;<br>
&gt;                       geolon_c:coordinates = &quot;geolon_c geolat_c&quot; ;<br>
&gt;<br>
&gt;               float geolat_t(yt_ocean, xt_ocean) ;<br>
&gt;<br>
&gt;                       geolat_t:long_name = &quot;tracer latitude&quot; ;<br>
&gt;<br>
&gt;                       geolat_t:units = &quot;degrees_N&quot; ;<br>
&gt;<br>
&gt;                       geolat_t:valid_range = -91.f, 91.f ;<br>
&gt;<br>
&gt;                       geolat_t:missing_value = 1.e+20f ;<br>
&gt;<br>
&gt;                       geolat_t:_FillValue = 1.e+20f ;<br>
&gt;<br>
&gt;                       geolat_t:cell_methods = &quot;time: point&quot; ;<br>
&gt;<br>
&gt;                       geolat_t:coordinates = &quot;geolon_t geolat_t&quot; ;<br>
&gt;<br>
&gt;               float geolon_t(yt_ocean, xt_ocean) ;<br>
&gt;<br>
&gt;                       geolon_t:long_name = &quot;tracer longitude&quot; ;<br>
&gt;<br>
&gt;                       geolon_t:units = &quot;degrees_E&quot; ;<br>
&gt;<br>
&gt;                       geolon_t:valid_range = -281.f, 361.f ;<br>
&gt;<br>
&gt;                       geolon_t:missing_value = 1.e+20f ;<br>
&gt;<br>
&gt;                       geolon_t:_FillValue = 1.e+20f ;<br>
&gt;<br>
&gt;                       geolon_t:cell_methods = &quot;time: point&quot; ;<br>
&gt;<br>
&gt;                       geolon_t:coordinates = &quot;geolon_t geolat_t&quot; ;<br>
&gt;<br>
&gt;               float area_t(yt_ocean, xt_ocean) ;<br>
&gt;<br>
&gt;                       area_t:long_name = &quot;tracer cell area&quot; ;<br>
&gt;<br>
&gt;                       area_t:units = &quot;m^2&quot; ;<br>
&gt;<br>
&gt;                       area_t:valid_range = 0.f, 1.e+15f ;<br>
&gt;<br>
&gt;                       area_t:missing_value = 1.e+20f ;<br>
&gt;<br>
&gt;                       area_t:_FillValue = 1.e+20f ;<br>
&gt;<br>
&gt;                       area_t:cell_methods = &quot;time: point&quot; ;<br>
&gt;<br>
&gt;                       area_t:coordinates = &quot;geolon_t geolat_t&quot; ;<br>
&gt;<br>
&gt;               float area_u(yu_ocean, xu_ocean) ;<br>
&gt;<br>
&gt;                       area_u:long_name = &quot;velocity cell area&quot; ;<br>
&gt;<br>
&gt;                       area_u:units = &quot;m^2&quot; ;<br>
&gt;<br>
&gt;                       area_u:valid_range = 0.f, 1.e+15f ;<br>
&gt;<br>
&gt;                       area_u:missing_value = 1.e+20f ;<br>
&gt;<br>
&gt;                       area_u:_FillValue = 1.e+20f ;<br>
&gt;<br>
&gt;                       area_u:cell_methods = &quot;time: point&quot; ;<br>
&gt;<br>
&gt;                       area_u:coordinates = &quot;geolon_c geolat_c&quot; ;<br>
&gt;<br>
&gt;               float kmt(yt_ocean, xt_ocean) ;<br>
&gt;<br>
&gt;                       kmt:long_name = &quot;number of depth levels on t-grid&quot; ;<br>
&gt;<br>
&gt;                       kmt:units = &quot;dimensionless&quot; ;<br>
&gt;<br>
&gt;                       kmt:valid_range = -10.f, 1.e+09f ;<br>
&gt;<br>
&gt;                       kmt:missing_value = -1.e+20f ;<br>
&gt;<br>
&gt;                       kmt:_FillValue = -1.e+20f ;<br>
&gt;<br>
&gt;                       kmt:cell_methods = &quot;time: point&quot; ;<br>
&gt;<br>
&gt;                       kmt:coordinates = &quot;geolon_t geolat_t&quot; ;<br>
&gt;<br>
&gt;               float kmu(yu_ocean, xu_ocean) ;<br>
&gt;<br>
&gt;                       kmu:long_name = &quot;number of depth levels on u-grid&quot; ;<br>
&gt;<br>
&gt;                       kmu:units = &quot;dimensionless&quot; ;<br>
&gt;<br>
&gt;                       kmu:valid_range = -10.f, 1.e+09f ;<br>
&gt;<br>
&gt;                       kmu:missing_value = -1.e+20f ;<br>
&gt;<br>
&gt;                       kmu:_FillValue = -1.e+20f ;<br>
&gt;<br>
&gt;                       kmu:cell_methods = &quot;time: point&quot; ;<br>
&gt;<br>
&gt;                       kmu:coordinates = &quot;geolon_c geolat_c&quot; ;<br>
&gt;<br>
&gt;               float geo_heat(yt_ocean, xt_ocean) ;<br>
&gt;<br>
&gt;                       geo_heat:long_name = &quot;Geothermal heating&quot; ;<br>
&gt;<br>
&gt;                       geo_heat:units = &quot;W/m^2&quot; ;<br>
&gt;<br>
&gt;                       geo_heat:valid_range = -10.f, 1000000.f ;<br>
&gt;<br>
&gt;                       geo_heat:missing_value = -1.e+20f ;<br>
&gt;<br>
&gt;                       geo_heat:_FillValue = -1.e+20f ;<br>
&gt;<br>
&gt;                       geo_heat:cell_methods = &quot;time: point&quot; ;<br>
&gt;<br>
&gt;                       geo_heat:coordinates = &quot;geolon_t geolat_t&quot; ;<br>
&gt;<br>
&gt;                       geo_heat:standard_name =<br>
&gt;     &quot;upward_geothermal_heat_flux_at_sea_floor&quot; ;<br>
&gt;<br>
&gt;               float pbot0(yt_ocean, xt_ocean) ;<br>
&gt;<br>
&gt;                       pbot0:long_name = &quot;reference bottom pressure<br>
&gt;     t-cells&quot; ;<br>
&gt;<br>
&gt;                       pbot0:units = &quot;dbar&quot; ;<br>
&gt;<br>
&gt;                       pbot0:valid_range = -1.e+08f, 1.e+08f ;<br>
&gt;                       pbot0:missing_value = -1.e+20f ;<br>
&gt;                       pbot0:_FillValue = -1.e+20f ;<br>
&gt;                       pbot0:cell_methods = &quot;time: point&quot; ;<br>
&gt;                       pbot0:coordinates = &quot;geolon_t geolat_t&quot; ;<br>
&gt;<br>
&gt;     // global attributes:<br>
&gt;                       :NCX = &quot;1.1.1&quot; ;<br>
&gt;     }<br>
&gt;<br>
&gt;     I hope this is helpful.  I don&#39;t think I can be of much more help, but<br>
&gt;     hope this points you in the right direction.<br>
&gt;<br>
&gt;     Chuck<br>
&gt;<br>
&gt;     On 07/25/2014 12:26 AM, David Stanfel wrote:<br>
&gt;      &gt; Hello,<br>
&gt;      &gt;<br>
&gt;      &gt; I decided to try a different file. Here is the ncdump -c of said<br>
&gt;     file:<br>
&gt;      &gt;<br>
&gt;      &gt; dstanfel:bin ninjawarrior343$ ncdump -c <a href="http://19800101.ocean_neutral.nc" target="_blank">19800101.ocean_neutral.nc</a><br>
&gt;     &lt;<a href="http://19800101.ocean_neutral.nc" target="_blank">http://19800101.ocean_neutral.nc</a>&gt;<br>
&gt;      &gt; &lt;<a href="http://19800101.ocean_neutral.nc" target="_blank">http://19800101.ocean_neutral.nc</a>&gt;<br>
&gt;      &gt;<br>
&gt;      &gt; netcdf \19800101.ocean_neutral {<br>
&gt;      &gt;<br>
&gt;      &gt; dimensions:<br>
&gt;      &gt;<br>
&gt;      &gt; xt_ocean = 24 ;<br>
&gt;      &gt;<br>
&gt;      &gt; yt_ocean = 35 ;<br>
&gt;      &gt;<br>
&gt;      &gt; Time = UNLIMITED ; // (1 currently)<br>
&gt;      &gt;<br>
&gt;      &gt; nv = 2 ;<br>
&gt;      &gt;<br>
&gt;      &gt; st_ocean = 18 ;<br>
&gt;      &gt;<br>
&gt;      &gt; st_edges_ocean = 19 ;<br>
&gt;      &gt;<br>
&gt;      &gt; variables:<br>
&gt;      &gt;<br>
&gt;      &gt; double xt_ocean(xt_ocean) ;<br>
&gt;      &gt;<br>
&gt;      &gt; xt_ocean:long_name = &quot;tcell longitude&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt; xt_ocean:units = &quot;degrees_E&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt; xt_ocean:cartesian_axis = &quot;X&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt; double yt_ocean(yt_ocean) ;<br>
&gt;      &gt;<br>
&gt;      &gt; yt_ocean:long_name = &quot;tcell latitude&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt; yt_ocean:units = &quot;degrees_N&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt; yt_ocean:cartesian_axis = &quot;Y&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt; double Time(Time) ;<br>
&gt;      &gt;<br>
&gt;      &gt; Time:long_name = &quot;Time&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt; Time:units = &quot;days since 1980-01-01 00:00:00&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt; Time:cartesian_axis = &quot;T&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt; Time:calendar_type = &quot;NOLEAP&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt; Time:calendar = &quot;NOLEAP&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt; Time:bounds = &quot;Time_bounds&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt; double nv(nv) ;<br>
&gt;      &gt;<br>
&gt;      &gt; nv:long_name = &quot;vertex number&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt; nv:units = &quot;none&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt; nv:cartesian_axis = &quot;N&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt; double st_ocean(st_ocean) ;<br>
&gt;      &gt;<br>
&gt;      &gt; st_ocean:long_name = &quot;tcell pstar&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt; st_ocean:units = &quot;dbars&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt; st_ocean:cartesian_axis = &quot;Z&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt; st_ocean:positive = &quot;down&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt; st_ocean:edges = &quot;st_edges_ocean&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt; double st_edges_ocean(st_edges_ocean) ;<br>
&gt;      &gt;<br>
&gt;      &gt; st_edges_ocean:long_name = &quot;tcell pstar edges&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt; st_edges_ocean:units = &quot;dbars&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt; st_edges_ocean:cartesian_axis = &quot;Z&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt; st_edges_ocean:positive = &quot;down&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt; float agm(Time, yt_ocean, xt_ocean) ;<br>
&gt;      &gt;<br>
&gt;      &gt; agm:long_name = &quot;GM diffusivity at surface&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt; agm:units = &quot;m^2/sec&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt; agm:valid_range = -10.f, 1.e+10f ;<br>
&gt;      &gt;<br>
&gt;      &gt; agm:missing_value = -1.e+20f ;<br>
&gt;      &gt;<br>
&gt;      &gt; agm:_FillValue = -1.e+20f ;<br>
&gt;      &gt;<br>
&gt;      &gt; agm:cell_methods = &quot;time: mean&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt; agm:time_avg_info = &quot;average_T1,average_T2,average_DT&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt; agm:coordinates = &quot;geolon_t geolat_t&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt; agm:standard_name = &quot;ocean_tracer_bolus_laplacian_diffusivity&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt; float rossby(Time, yt_ocean, xt_ocean) ;<br>
&gt;      &gt;<br>
&gt;      &gt; rossby:long_name = &quot;Rossby radius used in neutral physics&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt; rossby:units = &quot;m&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt; rossby:valid_range = -10.f, 1.e+10f ;<br>
&gt;      &gt;<br>
&gt;      &gt; rossby:missing_value = -1.e+20f ;<br>
&gt;      &gt;<br>
&gt;      &gt; rossby:_FillValue = -1.e+20f ;<br>
&gt;      &gt;<br>
&gt;      &gt; rossby:cell_methods = &quot;time: mean&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt; rossby:time_avg_info = &quot;average_T1,average_T2,average_DT&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt; rossby:coordinates = &quot;geolon_t geolat_t&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt; float slope31(Time, st_ocean, yt_ocean, xt_ocean) ;<br>
&gt;      &gt;<br>
&gt;      &gt; slope31:long_name = &quot;neutral slope -(rho_x/rho_z)&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt; slope31:units = &quot;dimensionless&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt; slope31:valid_range = -1.e+10f, 1.e+10f ;<br>
&gt;      &gt;<br>
&gt;      &gt; slope31:missing_value = -1.e+20f ;<br>
&gt;      &gt;<br>
&gt;      &gt; slope31:_FillValue = -1.e+20f ;<br>
&gt;      &gt;<br>
&gt;      &gt; slope31:cell_methods = &quot;time: mean&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt; slope31:time_avg_info = &quot;average_T1,average_T2,average_DT&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt; slope31:coordinates = &quot;geolon_t geolat_t&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt; float slope32(Time, st_ocean, yt_ocean, xt_ocean) ;<br>
&gt;      &gt;<br>
&gt;      &gt; slope32:long_name = &quot;neutral slope -(rho_y/rho_z)&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt; slope32:units = &quot;dimensionless&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt; slope32:valid_range = -1.e+10f, 1.e+10f ;<br>
&gt;      &gt;<br>
&gt;      &gt; slope32:missing_value = -1.e+20f ;<br>
&gt;      &gt;<br>
&gt;      &gt; slope32:_FillValue = -1.e+20f ;<br>
&gt;      &gt;<br>
&gt;      &gt; slope32:cell_methods = &quot;time: mean&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt; slope32:time_avg_info = &quot;average_T1,average_T2,average_DT&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt; slope32:coordinates = &quot;geolon_t geolat_t&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt; double average_T1(Time) ;<br>
&gt;      &gt;<br>
&gt;      &gt; average_T1:long_name = &quot;Start time for average period&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt; average_T1:units = &quot;days since 1980-01-01 00:00:00&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt; average_T1:missing_value = 1.e+20 ;<br>
&gt;      &gt;<br>
&gt;      &gt; average_T1:_FillValue = 1.e+20 ;<br>
&gt;      &gt;<br>
&gt;      &gt; double average_T2(Time) ;<br>
&gt;      &gt;<br>
&gt;      &gt; average_T2:long_name = &quot;End time for average period&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt; average_T2:units = &quot;days since 1980-01-01 00:00:00&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt; average_T2:missing_value = 1.e+20 ;<br>
&gt;      &gt;<br>
&gt;      &gt; average_T2:_FillValue = 1.e+20 ;<br>
&gt;      &gt;<br>
&gt;      &gt; double average_DT(Time) ;<br>
&gt;      &gt;<br>
&gt;      &gt; average_DT:long_name = &quot;Length of average period&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt; average_DT:units = &quot;days&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt; average_DT:missing_value = 1.e+20 ;<br>
&gt;      &gt;<br>
&gt;      &gt; average_DT:_FillValue = 1.e+20 ;<br>
&gt;      &gt;<br>
&gt;      &gt; double Time_bounds(Time, nv) ;<br>
&gt;      &gt;<br>
&gt;      &gt; Time_bounds:long_name = &quot;Time axis boundaries&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt; Time_bounds:units = &quot;days&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt; Time_bounds:missing_value = 1.e+20 ;<br>
&gt;      &gt;<br>
&gt;      &gt; Time_bounds:_FillValue = 1.e+20 ;<br>
&gt;      &gt;<br>
&gt;      &gt;<br>
&gt;      &gt; // global attributes:<br>
&gt;      &gt;<br>
&gt;      &gt; :filename = &quot;<a href="http://ocean_neutral.nc" target="_blank">ocean_neutral.nc</a> &lt;<a href="http://ocean_neutral.nc" target="_blank">http://ocean_neutral.nc</a>&gt;<br>
&gt;     &lt;<a href="http://ocean_neutral.nc" target="_blank">http://ocean_neutral.nc</a>&gt;&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt; :title = &quot;MOM4 box&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt; :grid_type = &quot;regular&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt; :grid_tile = &quot;N/A&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt; data:<br>
&gt;      &gt;<br>
&gt;      &gt;<br>
&gt;      &gt;   xt_ocean = 1.25, 3.75, 6.25, 8.75, 11.25, 13.75, 16.25, 18.75,<br>
&gt;     21.25,<br>
&gt;      &gt; 23.75,<br>
&gt;      &gt;<br>
&gt;      &gt;      26.25, 28.75, 31.25, 33.75, 36.25, 38.75, 41.25, 43.75,<br>
&gt;     46.25, 48.75,<br>
&gt;      &gt;<br>
&gt;      &gt;      51.25, 53.75, 56.25, 58.75 ;<br>
&gt;      &gt;<br>
&gt;      &gt;<br>
&gt;      &gt;   yt_ocean = 1.25, 3.7488100528717, 6.2393012046814,<br>
&gt;     8.71919345855713,<br>
&gt;      &gt;<br>
&gt;      &gt;      11.1815176010132, 13.6242818832397, 16.0408229827881,<br>
&gt;      &gt; 18.4296073913574,<br>
&gt;      &gt;<br>
&gt;      &gt;      20.784387588501, 23.1042194366455, 25.3833503723145,<br>
&gt;     27.6215190887451,<br>
&gt;      &gt;<br>
&gt;      &gt;      29.813497543335, 31.9597606658936, 34.0555992126465,<br>
&gt;     36.1022338867188,<br>
&gt;      &gt;<br>
&gt;      &gt;      38.0954322814941, 40.0371551513672, 41.9235687255859,<br>
&gt;      &gt; 43.7573394775391,<br>
&gt;      &gt;<br>
&gt;      &gt;      45.5349464416504, 47.259708404541, 48.9283294677734,<br>
&gt;     50.5447235107422,<br>
&gt;      &gt;<br>
&gt;      &gt;      52.1057090759277, 53.6157531738281, 55.0716934204102,<br>
&gt;      &gt; 56.4785079956055,<br>
&gt;      &gt;<br>
&gt;      &gt;      57.8329429626465, 59.1404571533203, 60.3976211547852,<br>
&gt;      &gt; 61.6103477478027,<br>
&gt;      &gt;<br>
&gt;      &gt;      62.7749481201172, 63.8977813720703, 64.9748153686523 ;<br>
&gt;      &gt;<br>
&gt;      &gt;<br>
&gt;      &gt;   Time = 2 ;<br>
&gt;      &gt;<br>
&gt;      &gt;<br>
&gt;      &gt;   nv = 1, 2 ;<br>
&gt;      &gt;<br>
&gt;      &gt;<br>
&gt;      &gt;   st_ocean = 17.75025, 53.25075, 88.75125, 124.25175, 159.75225,<br>
&gt;     195.25275,<br>
&gt;      &gt;<br>
&gt;      &gt;      230.75325, 266.25375, 314.34159289856, 399.171716592407,<br>
&gt;      &gt;<br>
&gt;      &gt;      541.922400494385, 759.079548358154, 1061.10117883301,<br>
&gt;      &gt; 1451.56990550537,<br>
&gt;      &gt;<br>
&gt;      &gt;      1926.90292906494, 2476.64250944824, 3084.30280480957,<br>
&gt;      &gt; 3728.70522619629 ;<br>
&gt;      &gt;<br>
&gt;      &gt;<br>
&gt;      &gt;   st_edges_ocean = 0, 35.5005, 71.001, 106.5015, 142.002, 177.5025,<br>
&gt;      &gt; 213.003,<br>
&gt;      &gt;<br>
&gt;      &gt;      248.5035, 290.29765597229, 356.756654745483, 470.547043066406,<br>
&gt;      &gt;<br>
&gt;      &gt;      650.50097442627, 910.090332641602, 1256.33548026123,<br>
&gt;     1689.23641728516,<br>
&gt;      &gt;<br>
&gt;      &gt;      2201.77290498047, 2780.47265712891, 3406.50401550293,<br>
&gt;     4057.20024763184<br>
&gt;      &gt;<br>
&gt;      &gt;<br>
&gt;      &gt; I wrote a .ctl file:<br>
&gt;      &gt;<br>
&gt;      &gt;<br>
&gt;      &gt; DSET ^<a href="http://19800101.ocean_neutral.nc" target="_blank">19800101.ocean_neutral.nc</a><br>
</div></div>&gt;     &lt;<a href="http://19800101.ocean_neutral.nc" target="_blank">http://19800101.ocean_neutral.nc</a>&gt; &lt;<a href="http://19800101.ocean_neutral.nc" target="_blank">http://19800101.ocean_neutral.nc</a>&gt;<br>

<div><div class="h5">&gt;      &gt;<br>
&gt;      &gt; UNDEF -9.99e8 missing_value<br>
&gt;      &gt;<br>
&gt;      &gt; DTYPE netcdf<br>
&gt;      &gt;<br>
&gt;      &gt; CACHESIZE 100000<br>
&gt;      &gt;<br>
&gt;      &gt; TITLE Sample Model Output<br>
&gt;      &gt;<br>
&gt;      &gt; XDEF 24 LINEAR  1.25   2.5<br>
&gt;      &gt;<br>
&gt;      &gt; YDEF 35 levels  1.25, 3.7488100528717, 6.2393012046814,<br>
&gt;     8.71919345855713,<br>
&gt;      &gt;<br>
&gt;      &gt;      11.1815176010132, 13.6242818832397, 16.0408229827881,<br>
&gt;      &gt; 18.4296073913574,<br>
&gt;      &gt;<br>
&gt;      &gt;      20.784387588501, 23.1042194366455, 25.3833503723145,<br>
&gt;     27.6215190887451,<br>
&gt;      &gt;<br>
&gt;      &gt;      29.813497543335, 31.9597606658936, 34.0555992126465,<br>
&gt;     36.1022338867188,<br>
&gt;      &gt;<br>
&gt;      &gt;      38.0954322814941, 40.0371551513672, 41.9235687255859,<br>
&gt;      &gt; 43.7573394775391,<br>
&gt;      &gt;<br>
&gt;      &gt;      45.5349464416504, 47.259708404541, 48.9283294677734,<br>
&gt;     50.5447235107422,<br>
&gt;      &gt;<br>
&gt;      &gt;      52.1057090759277, 53.6157531738281, 55.0716934204102,<br>
&gt;      &gt; 56.4785079956055,<br>
&gt;      &gt;<br>
&gt;      &gt;      57.8329429626465, 59.1404571533203, 60.3976211547852,<br>
&gt;      &gt; 61.6103477478027,<br>
&gt;      &gt;<br>
&gt;      &gt;      62.7749481201172, 63.8977813720703, 64.9748153686523<br>
&gt;      &gt;<br>
&gt;      &gt; ZDEF 18 levels  17.75025, 53.25075, 88.75125, 124.25175, 159.75225,<br>
&gt;      &gt; 195.25275, 230.75327, 266.25375, 314.34159289856, 399.171716592407,<br>
&gt;      &gt; 541.922400494385, 759.079548358154, 1061.10117883301,<br>
&gt;     1451.56990550537,<br>
&gt;      &gt; 1926.90292906494, 2476.64250944824, 3084.30280480957,<br>
&gt;     3728.70522619629<br>
&gt;      &gt;<br>
&gt;      &gt; TDEF 1 LINEAR 01JAN0001 1mo<br>
&gt;      &gt;<br>
&gt;      &gt; VARS 4<br>
&gt;      &gt;<br>
&gt;      &gt; agm    0  t,y,x  GM diffusivity at surface (m^2/sec)<br>
&gt;      &gt;<br>
&gt;      &gt; rossby 0  t,y,x  Rossby radius used in neutral physics (m)<br>
&gt;      &gt;<br>
&gt;      &gt; slope31 18 t,z,y,x Neutral slope -(rho_x/rho_z)<br>
&gt;      &gt;<br>
&gt;      &gt; slope32 18 t,z,y,x Neutral slope -(rho_y/rho_z)<br>
&gt;      &gt;<br>
&gt;      &gt; ENDVARS<br>
&gt;      &gt;<br>
&gt;      &gt;<br>
&gt;      &gt; When trying to display the variables, I get this error:<br>
&gt;      &gt;<br>
&gt;      &gt;<br>
&gt;      &gt; ga-&gt; d agm<br>
&gt;      &gt;<br>
&gt;      &gt; Cannot contour grid - all undefined values<br>
&gt;      &gt;<br>
&gt;      &gt; ga-&gt; d rossby<br>
&gt;      &gt;<br>
&gt;      &gt; Cannot contour grid - all undefined values<br>
&gt;      &gt;<br>
&gt;      &gt; ga-&gt; d slope31<br>
&gt;      &gt;<br>
&gt;      &gt; Cannot contour grid - all undefined values<br>
&gt;      &gt;<br>
&gt;      &gt; ga-&gt;<br>
&gt;      &gt;<br>
&gt;      &gt;<br>
&gt;      &gt; What could be going wrong here?<br>
&gt;      &gt;<br>
&gt;      &gt; Thanks,<br>
&gt;      &gt;<br>
&gt;      &gt; David<br>
&gt;      &gt;<br>
&gt;      &gt;<br>
&gt;      &gt; On Wed, Jul 23, 2014 at 11:23 AM, David Stanfel<br>
&gt;     &lt;<a href="mailto:dstanfel@ualberta.ca">dstanfel@ualberta.ca</a> &lt;mailto:<a href="mailto:dstanfel@ualberta.ca">dstanfel@ualberta.ca</a>&gt;<br>
</div></div><div class="">&gt;      &gt; &lt;mailto:<a href="mailto:dstanfel@ualberta.ca">dstanfel@ualberta.ca</a> &lt;mailto:<a href="mailto:dstanfel@ualberta.ca">dstanfel@ualberta.ca</a>&gt;&gt;&gt; wrote:<br>
&gt;      &gt;<br>
&gt;      &gt;     Also, how do you determine the CACHESIZE?<br>
&gt;      &gt;<br>
&gt;      &gt;     Thanks,<br>
&gt;      &gt;<br>
&gt;      &gt;     David<br>
&gt;      &gt;<br>
&gt;      &gt;<br>
&gt;      &gt;     On Tue, Jul 22, 2014 at 11:05 PM, David Stanfel<br>
&gt;      &gt;     &lt;<a href="mailto:dstanfel@ualberta.ca">dstanfel@ualberta.ca</a> &lt;mailto:<a href="mailto:dstanfel@ualberta.ca">dstanfel@ualberta.ca</a>&gt;<br>
</div><div><div class="h5">&gt;     &lt;mailto:<a href="mailto:dstanfel@ualberta.ca">dstanfel@ualberta.ca</a> &lt;mailto:<a href="mailto:dstanfel@ualberta.ca">dstanfel@ualberta.ca</a>&gt;&gt;&gt; wrote:<br>
&gt;      &gt;<br>
&gt;      &gt;         Hello,<br>
&gt;      &gt;<br>
&gt;      &gt;         I just ran a climate model which produced some .nc files.<br>
&gt;     I am<br>
&gt;      &gt;         having trouble writing the .ctl file to read said .nc files.<br>
&gt;      &gt;         Here is the ncdump -c for one of the files I am trying to<br>
&gt;     read:<br>
&gt;      &gt;<br>
&gt;      &gt;         dstanfel:bin ninjawarrior343$ ncdump -c<br>
&gt;     <a href="http://19800101.ocean_layer.nc" target="_blank">19800101.ocean_layer.nc</a> &lt;<a href="http://19800101.ocean_layer.nc" target="_blank">http://19800101.ocean_layer.nc</a>&gt;<br>
&gt;      &gt;         &lt;<a href="http://19800101.ocean_layer.nc" target="_blank">http://19800101.ocean_layer.nc</a>&gt;<br>
&gt;      &gt;<br>
&gt;      &gt;         netcdf \19800101.ocean_layer {<br>
&gt;      &gt;<br>
&gt;      &gt;         dimensions:<br>
&gt;      &gt;<br>
&gt;      &gt;         grid_xt_ocean = 24 ;<br>
&gt;      &gt;<br>
&gt;      &gt;         grid_yt_ocean = 35 ;<br>
&gt;      &gt;<br>
&gt;      &gt;         theta = 100 ;<br>
&gt;      &gt;<br>
&gt;      &gt;         theta_edges = 101 ;<br>
&gt;      &gt;<br>
&gt;      &gt;         Time = UNLIMITED ; // (1 currently)<br>
&gt;      &gt;<br>
&gt;      &gt;         nv = 2 ;<br>
&gt;      &gt;<br>
&gt;      &gt;         xt_ocean = 24 ;<br>
&gt;      &gt;<br>
&gt;      &gt;         yt_ocean = 35 ;<br>
&gt;      &gt;<br>
&gt;      &gt;         st_ocean = 18 ;<br>
&gt;      &gt;<br>
&gt;      &gt;         st_edges_ocean = 19 ;<br>
&gt;      &gt;<br>
&gt;      &gt;         potrho = 100 ;<br>
&gt;      &gt;<br>
&gt;      &gt;         potrho_edges = 101 ;<br>
&gt;      &gt;<br>
&gt;      &gt;         variables:<br>
&gt;      &gt;<br>
&gt;      &gt;         double grid_xt_ocean(grid_xt_ocean) ;<br>
&gt;      &gt;<br>
&gt;      &gt;         grid_xt_ocean:long_name = &quot;tcell longitude&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         grid_xt_ocean:units = &quot;degrees_E&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         grid_xt_ocean:cartesian_axis = &quot;X&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         double grid_yt_ocean(grid_yt_ocean) ;<br>
&gt;      &gt;<br>
&gt;      &gt;         grid_yt_ocean:long_name = &quot;tcell latitude&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         grid_yt_ocean:units = &quot;degrees_N&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         grid_yt_ocean:cartesian_axis = &quot;Y&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         double theta(theta) ;<br>
&gt;      &gt;<br>
&gt;      &gt;         theta:long_name = &quot;potential or conservative temperature&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         theta:units = &quot;C&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         theta:cartesian_axis = &quot;Z&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         theta:positive = &quot;up&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         theta:edges = &quot;theta_edges&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         double theta_edges(theta_edges) ;<br>
&gt;      &gt;<br>
&gt;      &gt;         theta_edges:long_name = &quot;potential or conservative<br>
&gt;     temperature<br>
&gt;      &gt;         edges&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         theta_edges:units = &quot;C&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         theta_edges:cartesian_axis = &quot;Z&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         theta_edges:positive = &quot;up&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         double Time(Time) ;<br>
&gt;      &gt;<br>
&gt;      &gt;         Time:long_name = &quot;Time&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         Time:units = &quot;days since 1980-01-01 00:00:00&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         Time:cartesian_axis = &quot;T&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         Time:calendar_type = &quot;NOLEAP&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         Time:calendar = &quot;NOLEAP&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         Time:bounds = &quot;Time_bounds&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         double nv(nv) ;<br>
&gt;      &gt;<br>
&gt;      &gt;         nv:long_name = &quot;vertex number&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         nv:units = &quot;none&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         nv:cartesian_axis = &quot;N&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         double xt_ocean(xt_ocean) ;<br>
&gt;      &gt;<br>
&gt;      &gt;         xt_ocean:long_name = &quot;tcell longitude&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         xt_ocean:units = &quot;degrees_E&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         xt_ocean:cartesian_axis = &quot;X&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         double yt_ocean(yt_ocean) ;<br>
&gt;      &gt;<br>
&gt;      &gt;         yt_ocean:long_name = &quot;tcell latitude&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         yt_ocean:units = &quot;degrees_N&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         yt_ocean:cartesian_axis = &quot;Y&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         double st_ocean(st_ocean) ;<br>
&gt;      &gt;<br>
&gt;      &gt;         st_ocean:long_name = &quot;tcell pstar&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         st_ocean:units = &quot;dbars&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         st_ocean:cartesian_axis = &quot;Z&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         st_ocean:positive = &quot;down&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         st_ocean:edges = &quot;st_edges_ocean&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         double st_edges_ocean(st_edges_ocean) ;<br>
&gt;      &gt;<br>
&gt;      &gt;         st_edges_ocean:long_name = &quot;tcell pstar edges&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         st_edges_ocean:units = &quot;dbars&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         st_edges_ocean:cartesian_axis = &quot;Z&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         st_edges_ocean:positive = &quot;down&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         double potrho(potrho) ;<br>
&gt;      &gt;<br>
&gt;      &gt;         potrho:long_name = &quot;potential density&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         potrho:units = &quot;kg/m^3&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         potrho:cartesian_axis = &quot;Z&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         potrho:positive = &quot;down&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         potrho:edges = &quot;potrho_edges&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         double potrho_edges(potrho_edges) ;<br>
&gt;      &gt;<br>
&gt;      &gt;         potrho_edges:long_name = &quot;potential density edges&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         potrho_edges:units = &quot;kg/m^3&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         potrho_edges:cartesian_axis = &quot;Z&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         potrho_edges:positive = &quot;down&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         float depth_of_theta(Time, theta, grid_yt_ocean,<br>
&gt;     grid_xt_ocean) ;<br>
&gt;      &gt;<br>
&gt;      &gt;         depth_of_theta:long_name = &quot;depth of potential temp<br>
&gt;     surface&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         depth_of_theta:units = &quot;m&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         depth_of_theta:valid_range = 0.f, 1.e+10f ;<br>
&gt;      &gt;<br>
&gt;      &gt;         depth_of_theta:missing_value = -1.e+20f ;<br>
&gt;      &gt;<br>
&gt;      &gt;         depth_of_theta:_FillValue = -1.e+20f ;<br>
&gt;      &gt;<br>
&gt;      &gt;         depth_of_theta:cell_methods = &quot;time: mean&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         depth_of_theta:time_avg_info =<br>
&gt;     &quot;average_T1,average_T2,average_DT&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         float rho(Time, st_ocean, yt_ocean, xt_ocean) ;<br>
&gt;      &gt;<br>
&gt;      &gt;         rho:long_name = &quot;in situ density&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         rho:units = &quot;kg/m^3&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         rho:valid_range = -10.f, 100000.f ;<br>
&gt;      &gt;<br>
&gt;      &gt;         rho:missing_value = -1.e+20f ;<br>
&gt;      &gt;<br>
&gt;      &gt;         rho:_FillValue = -1.e+20f ;<br>
&gt;      &gt;<br>
&gt;      &gt;         rho:cell_methods = &quot;time: mean&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         rho:time_avg_info = &quot;average_T1,average_T2,average_DT&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         rho:coordinates = &quot;geolon_t geolat_t&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         float neutral_rho(Time, st_ocean, yt_ocean, xt_ocean) ;<br>
&gt;      &gt;<br>
&gt;      &gt;         neutral_rho:long_name = &quot;potential density estimate of<br>
&gt;     neutral<br>
&gt;      &gt;         density&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         neutral_rho:units = &quot;kg/m^3&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         neutral_rho:valid_range = -10.f, 100000.f ;<br>
&gt;      &gt;<br>
&gt;      &gt;         neutral_rho:missing_value = -1.e+20f ;<br>
&gt;      &gt;<br>
&gt;      &gt;         neutral_rho:_FillValue = -1.e+20f ;<br>
&gt;      &gt;<br>
&gt;      &gt;         neutral_rho:cell_methods = &quot;time: mean&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         neutral_rho:time_avg_info =<br>
&gt;     &quot;average_T1,average_T2,average_DT&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         neutral_rho:coordinates = &quot;geolon_t geolat_t&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         float pot_rho_0(Time, st_ocean, yt_ocean, xt_ocean) ;<br>
&gt;      &gt;<br>
&gt;      &gt;         pot_rho_0:long_name = &quot;potential density referenced to 0<br>
&gt;     dbar&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         pot_rho_0:units = &quot;kg/m^3&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         pot_rho_0:valid_range = -10.f, 100000.f ;<br>
&gt;      &gt;<br>
&gt;      &gt;         pot_rho_0:missing_value = -1.e+20f ;<br>
&gt;      &gt;<br>
&gt;      &gt;         pot_rho_0:_FillValue = -1.e+20f ;<br>
&gt;      &gt;<br>
&gt;      &gt;         pot_rho_0:cell_methods = &quot;time: mean&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         pot_rho_0:time_avg_info =<br>
&gt;     &quot;average_T1,average_T2,average_DT&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         pot_rho_0:coordinates = &quot;geolon_t geolat_t&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         pot_rho_0:standard_name = &quot;sea_water_potential_density&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         float salt_on_rho(Time, potrho, grid_yt_ocean,<br>
&gt;     grid_xt_ocean) ;<br>
&gt;      &gt;<br>
&gt;      &gt;         salt_on_rho:long_name = &quot;salt on potential density surface&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         salt_on_rho:units = &quot;psu&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         salt_on_rho:valid_range = 0.f, 1000.f ;<br>
&gt;      &gt;<br>
&gt;      &gt;         salt_on_rho:missing_value = -1.e+20f ;<br>
&gt;      &gt;<br>
&gt;      &gt;         salt_on_rho:_FillValue = -1.e+20f ;<br>
&gt;      &gt;<br>
&gt;      &gt;         salt_on_rho:cell_methods = &quot;time: mean&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         salt_on_rho:time_avg_info =<br>
&gt;     &quot;average_T1,average_T2,average_DT&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         float potrho_mix_depth(Time, yt_ocean, xt_ocean) ;<br>
&gt;      &gt;<br>
&gt;      &gt;         potrho_mix_depth:long_name = &quot;Depth of potential density<br>
&gt;     mixed<br>
&gt;      &gt;         layer&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         potrho_mix_depth:units = &quot;m&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         potrho_mix_depth:valid_range = -1000000.f, 1000000.f ;<br>
&gt;      &gt;<br>
&gt;      &gt;         potrho_mix_depth:missing_value = -1.e+20f ;<br>
&gt;      &gt;<br>
&gt;      &gt;         potrho_mix_depth:_FillValue = -1.e+20f ;<br>
&gt;      &gt;<br>
&gt;      &gt;         potrho_mix_depth:cell_methods = &quot;time: mean&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         potrho_mix_depth:time_avg_info =<br>
&gt;      &gt;         &quot;average_T1,average_T2,average_DT&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         potrho_mix_depth:coordinates = &quot;geolon_t geolat_t&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         float potrho_mix_base(Time, yt_ocean, xt_ocean) ;<br>
&gt;      &gt;<br>
&gt;      &gt;         potrho_mix_base:long_name = &quot;Potential density at mixed layer<br>
&gt;      &gt;         base&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         potrho_mix_base:units = &quot;kg/m^3&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         potrho_mix_base:valid_range = -1000000.f, 1000000.f ;<br>
&gt;      &gt;<br>
&gt;      &gt;         potrho_mix_base:missing_value = -1.e+20f ;<br>
&gt;      &gt;<br>
&gt;      &gt;         potrho_mix_base:_FillValue = -1.e+20f ;<br>
&gt;      &gt;<br>
&gt;      &gt;         potrho_mix_base:cell_methods = &quot;time: mean&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         potrho_mix_base:time_avg_info =<br>
&gt;     &quot;average_T1,average_T2,average_DT&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         potrho_mix_base:coordinates = &quot;geolon_t geolat_t&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         float press(Time, st_ocean, yt_ocean, xt_ocean) ;<br>
&gt;      &gt;<br>
&gt;      &gt;         press:long_name = &quot;absolute pressure&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         press:units = &quot;dbar&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         press:valid_range = -10.f, 1000000.f ;<br>
&gt;      &gt;<br>
&gt;      &gt;         press:missing_value = -1.e+20f ;<br>
&gt;      &gt;<br>
&gt;      &gt;         press:_FillValue = -1.e+20f ;<br>
&gt;      &gt;<br>
&gt;      &gt;         press:cell_methods = &quot;time: mean&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         press:time_avg_info = &quot;average_T1,average_T2,average_DT&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         press:coordinates = &quot;geolon_t geolat_t&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         float mld(Time, yt_ocean, xt_ocean) ;<br>
&gt;      &gt;<br>
&gt;      &gt;         mld:long_name = &quot;mixed layer depth determined by density<br>
&gt;     criteria&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         mld:units = &quot;m&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         mld:valid_range = 0.f, 1000000.f ;<br>
&gt;      &gt;<br>
&gt;      &gt;         mld:missing_value = -1.e+20f ;<br>
&gt;      &gt;<br>
&gt;      &gt;         mld:_FillValue = -1.e+20f ;<br>
&gt;      &gt;<br>
&gt;      &gt;         mld:cell_methods = &quot;time: mean&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         mld:time_avg_info = &quot;average_T1,average_T2,average_DT&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         mld:coordinates = &quot;geolon_t geolat_t&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         mld:standard_name =<br>
&gt;      &gt;         &quot;ocean_mixed_layer_thickness_defined_by_sigma_t&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         double average_T1(Time) ;<br>
&gt;      &gt;<br>
&gt;      &gt;         average_T1:long_name = &quot;Start time for average period&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         average_T1:units = &quot;days since 1980-01-01 00:00:00&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         average_T1:missing_value = 1.e+20 ;<br>
&gt;      &gt;<br>
&gt;      &gt;         average_T1:_FillValue = 1.e+20 ;<br>
&gt;      &gt;<br>
&gt;      &gt;         double average_T2(Time) ;<br>
&gt;      &gt;<br>
&gt;      &gt;         average_T2:long_name = &quot;End time for average period&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         average_T2:units = &quot;days since 1980-01-01 00:00:00&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         average_T2:missing_value = 1.e+20 ;<br>
&gt;      &gt;<br>
&gt;      &gt;         average_T2:_FillValue = 1.e+20 ;<br>
&gt;      &gt;<br>
&gt;      &gt;         double average_DT(Time) ;<br>
&gt;      &gt;<br>
&gt;      &gt;         average_DT:long_name = &quot;Length of average period&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         average_DT:units = &quot;days&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         average_DT:missing_value = 1.e+20 ;<br>
&gt;      &gt;<br>
&gt;      &gt;         average_DT:_FillValue = 1.e+20 ;<br>
&gt;      &gt;<br>
&gt;      &gt;         double Time_bounds(Time, nv) ;<br>
&gt;      &gt;<br>
&gt;      &gt;         Time_bounds:long_name = &quot;Time axis boundaries&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         Time_bounds:units = &quot;days&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         Time_bounds:missing_value = 1.e+20 ;<br>
&gt;      &gt;<br>
&gt;      &gt;         Time_bounds:_FillValue = 1.e+20 ;<br>
&gt;      &gt;<br>
&gt;      &gt;<br>
&gt;      &gt;         // global attributes:<br>
&gt;      &gt;<br>
&gt;      &gt;         :filename = &quot;<a href="http://ocean_layer.nc" target="_blank">ocean_layer.nc</a> &lt;<a href="http://ocean_layer.nc" target="_blank">http://ocean_layer.nc</a>&gt;<br>
&gt;     &lt;<a href="http://ocean_layer.nc" target="_blank">http://ocean_layer.nc</a>&gt;&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         :title = &quot;MOM4 box&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         :grid_type = &quot;regular&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         :grid_tile = &quot;N/A&quot; ;<br>
&gt;      &gt;<br>
&gt;      &gt;         data:<br>
&gt;      &gt;<br>
&gt;      &gt;<br>
&gt;      &gt;           grid_xt_ocean = 1.25, 3.75, 6.25, 8.75, 11.25, 13.75,<br>
&gt;     16.25,<br>
&gt;      &gt;         18.75, 21.25,<br>
&gt;      &gt;<br>
&gt;      &gt;              23.75, 26.25, 28.75, 31.25, 33.75, 36.25, 38.75, 41.25,<br>
&gt;      &gt;         43.75, 46.25,<br>
&gt;      &gt;<br>
&gt;      &gt;              48.75, 51.25, 53.75, 56.25, 58.75 ;<br>
&gt;      &gt;<br>
&gt;      &gt;<br>
&gt;      &gt;           grid_yt_ocean = 1.25, 3.7488100528717, 6.2393012046814,<br>
&gt;      &gt;         8.71919345855713,<br>
&gt;      &gt;<br>
&gt;      &gt;              11.1815176010132, 13.6242818832397, 16.0408229827881,<br>
&gt;      &gt;         18.4296073913574,<br>
&gt;      &gt;<br>
&gt;      &gt;              20.784387588501, 23.1042194366455, 25.3833503723145,<br>
&gt;      &gt;         27.6215190887451,<br>
&gt;      &gt;<br>
&gt;      &gt;              29.813497543335, 31.9597606658936, 34.0555992126465,<br>
&gt;      &gt;         36.1022338867188,<br>
&gt;      &gt;<br>
&gt;      &gt;              38.0954322814941, 40.0371551513672, 41.9235687255859,<br>
&gt;      &gt;         43.7573394775391,<br>
&gt;      &gt;<br>
&gt;      &gt;              45.5349464416504, 47.259708404541, 48.9283294677734,<br>
&gt;      &gt;         50.5447235107422,<br>
&gt;      &gt;<br>
&gt;      &gt;              52.1057090759277, 53.6157531738281, 55.0716934204102,<br>
&gt;      &gt;         56.4785079956055,<br>
&gt;      &gt;<br>
&gt;      &gt;              57.8329429626465, 59.1404571533203, 60.3976211547852,<br>
&gt;      &gt;         61.6103477478027,<br>
&gt;      &gt;<br>
&gt;      &gt;              62.7749481201172, 63.8977813720703, 64.9748153686523 ;<br>
&gt;      &gt;<br>
&gt;      &gt;<br>
&gt;      &gt;           theta = -1.84, -1.52, -1.2, -0.88, -0.56, -0.24,<br>
&gt;      &gt;         0.0800000000000003, 0.4,<br>
&gt;      &gt;<br>
&gt;      &gt;              0.72, 1.04, 1.36, 1.68, 2, 2.32, 2.64, 2.96, 3.28, 3.6,<br>
&gt;      &gt;         3.92, 4.24, 4.56,<br>
&gt;      &gt;<br>
&gt;      &gt;              4.88, 5.2, 5.52, 5.84, 6.16, 6.48, 6.8, 7.12, 7.44,<br>
&gt;     7.76,<br>
&gt;      &gt;         8.08, 8.4,<br>
&gt;      &gt;<br>
&gt;      &gt;              8.72, 9.04, 9.36, 9.68000000000001, 10, 10.32, 10.64,<br>
&gt;      &gt;         10.96, 11.28, 11.6,<br>
&gt;      &gt;<br>
&gt;      &gt;              11.92, 12.24, 12.56, 12.88, 13.2, 13.52, 13.84, 14.16,<br>
&gt;      &gt;         14.48, 14.8,<br>
&gt;      &gt;<br>
&gt;      &gt;              15.12, 15.44, 15.76, 16.08, 16.4, 16.72, 17.04, 17.36,<br>
&gt;      &gt;         17.68, 18, 18.32,<br>
&gt;      &gt;<br>
&gt;      &gt;              18.64, 18.96, 19.28, 19.6, 19.92, 20.24, 20.56, 20.88,<br>
&gt;      &gt;         21.2, 21.52,<br>
&gt;      &gt;<br>
&gt;      &gt;              21.84, 22.16, 22.48, 22.8, 23.12, 23.44, 23.76, 24.08,<br>
&gt;      &gt;         24.4, 24.72,<br>
&gt;      &gt;<br>
&gt;      &gt;              25.04, 25.36, 25.68, 26, 26.32, 26.64, 26.96, 27.28,<br>
&gt;     27.6,<br>
&gt;      &gt;         27.92, 28.24,<br>
&gt;      &gt;<br>
&gt;      &gt;              28.56, 28.88, 29.2, 29.52, 29.84 ;<br>
&gt;      &gt;<br>
&gt;      &gt;<br>
&gt;      &gt;           theta_edges = 1028, 1028.1, 1028.2, 1028.3, 1028.4, 1028.5,<br>
&gt;      &gt;         1028.6, 1028.7,<br>
&gt;      &gt;<br>
&gt;      &gt;              1028.8, 1028.9, 1029, 1029.1, 1029.2, 1029.3, 1029.4,<br>
&gt;      &gt;         1029.5, 1029.6,<br>
&gt;      &gt;<br>
&gt;      &gt;              1029.7, 1029.8, 1029.9, 1030, 1030.1, 1030.2, 1030.3,<br>
&gt;      &gt;         1030.4, 1030.5,<br>
&gt;      &gt;<br>
&gt;      &gt;              1030.6, 1030.7, 1030.8, 1030.9, 1031, 1031.1, 1031.2,<br>
&gt;      &gt;         1031.3, 1031.4,<br>
&gt;      &gt;<br>
&gt;      &gt;              1031.5, 1031.6, 1031.7, 1031.8, 1031.9, 1032, 1032.1,<br>
&gt;      &gt;         1032.2, 1032.3,<br>
&gt;      &gt;<br>
&gt;      &gt;              1032.4, 1032.5, 1032.6, 1032.7, 1032.8, 1032.9, 1033,<br>
&gt;      &gt;         1033.1, 1033.2,<br>
&gt;      &gt;<br>
&gt;      &gt;              1033.3, 1033.4, 1033.49999999999, 1033.59999999999,<br>
&gt;      &gt;         1033.69999999999,<br>
&gt;      &gt;<br>
&gt;      &gt;              1033.79999999999, 1033.89999999999, 1033.99999999999,<br>
&gt;      &gt;         1034.09999999999,<br>
&gt;      &gt;<br>
&gt;      &gt;              1034.19999999999, 1034.29999999999, 1034.39999999999,<br>
&gt;      &gt;         1034.49999999999,<br>
&gt;      &gt;<br>
&gt;      &gt;              1034.59999999999, 1034.69999999999, 1034.79999999999,<br>
&gt;      &gt;         1034.89999999999,<br>
&gt;      &gt;<br>
&gt;      &gt;              1034.99999999999, 1035.09999999999, 1035.19999999999,<br>
&gt;      &gt;         1035.29999999999,<br>
&gt;      &gt;<br>
&gt;      &gt;              1035.39999999999, 1035.49999999999, 1035.59999999999,<br>
&gt;      &gt;         1035.69999999999,<br>
&gt;      &gt;<br>
&gt;      &gt;              1035.79999999999, 1035.89999999999, 1035.99999999999,<br>
&gt;      &gt;         1036.09999999999,<br>
&gt;      &gt;<br>
&gt;      &gt;              1036.19999999999, 1036.29999999999, 1036.39999999999,<br>
&gt;      &gt;         1036.49999999999,<br>
&gt;      &gt;<br>
&gt;      &gt;              1036.59999999999, 1036.69999999999, 1036.79999999999,<br>
&gt;      &gt;         1036.89999999999,<br>
&gt;      &gt;<br>
&gt;      &gt;              1036.99999999999, 1037.09999999999, 1037.19999999999,<br>
&gt;      &gt;         1037.29999999999,<br>
&gt;      &gt;<br>
&gt;      &gt;              1037.39999999999, 1037.49999999999, 1037.59999999999,<br>
&gt;      &gt;         1037.69999999999,<br>
&gt;      &gt;<br>
&gt;      &gt;              1037.79999999999, 1037.89999999999, 1037.99999999999 ;<br>
&gt;      &gt;<br>
&gt;      &gt;<br>
&gt;      &gt;           Time = 2 ;<br>
&gt;      &gt;<br>
&gt;      &gt;<br>
&gt;      &gt;           nv = 1, 2 ;<br>
&gt;      &gt;<br>
&gt;      &gt;<br>
&gt;      &gt;           xt_ocean = 1.25, 3.75, 6.25, 8.75, 11.25, 13.75, 16.25,<br>
&gt;     18.75,<br>
&gt;      &gt;         21.25, 23.75,<br>
&gt;      &gt;<br>
&gt;      &gt;              26.25, 28.75, 31.25, 33.75, 36.25, 38.75, 41.25, 43.75,<br>
&gt;      &gt;         46.25, 48.75,<br>
&gt;      &gt;<br>
&gt;      &gt;              51.25, 53.75, 56.25, 58.75 ;<br>
&gt;      &gt;<br>
&gt;      &gt;<br>
&gt;      &gt;           yt_ocean = 1.25, 3.7488100528717, 6.2393012046814,<br>
&gt;      &gt;         8.71919345855713,<br>
&gt;      &gt;<br>
&gt;      &gt;              11.1815176010132, 13.6242818832397, 16.0408229827881,<br>
&gt;      &gt;         18.4296073913574,<br>
&gt;      &gt;<br>
&gt;      &gt;              20.784387588501, 23.1042194366455, 25.3833503723145,<br>
&gt;      &gt;         27.6215190887451,<br>
&gt;      &gt;<br>
&gt;      &gt;              29.813497543335, 31.9597606658936, 34.0555992126465,<br>
&gt;      &gt;         36.1022338867188,<br>
&gt;      &gt;<br>
&gt;      &gt;              38.0954322814941, 40.0371551513672, 41.9235687255859,<br>
&gt;      &gt;         43.7573394775391,<br>
&gt;      &gt;<br>
&gt;      &gt;              45.5349464416504, 47.259708404541, 48.9283294677734,<br>
&gt;      &gt;         50.5447235107422,<br>
&gt;      &gt;<br>
&gt;      &gt;              52.1057090759277, 53.6157531738281, 55.0716934204102,<br>
&gt;      &gt;         56.4785079956055,<br>
&gt;      &gt;<br>
&gt;      &gt;              57.8329429626465, 59.1404571533203, 60.3976211547852,<br>
&gt;      &gt;         61.6103477478027,<br>
&gt;      &gt;<br>
&gt;      &gt;              62.7749481201172, 63.8977813720703, 64.9748153686523 ;<br>
&gt;      &gt;<br>
&gt;      &gt;<br>
&gt;      &gt;           st_ocean = 17.75025, 53.25075, 88.75125, 124.25175,<br>
&gt;     159.75225,<br>
&gt;      &gt;         195.25275,<br>
&gt;      &gt;<br>
&gt;      &gt;              230.75325, 266.25375, 314.34159289856, 399.171716592407,<br>
&gt;      &gt;<br>
&gt;      &gt;              541.922400494385, 759.079548358154, 1061.10117883301,<br>
&gt;      &gt;         1451.56990550537,<br>
&gt;      &gt;<br>
&gt;      &gt;              1926.90292906494, 2476.64250944824, 3084.30280480957,<br>
&gt;      &gt;         3728.70522619629 ;<br>
&gt;      &gt;<br>
&gt;      &gt;<br>
&gt;      &gt;           st_edges_ocean = 0, 35.5005, 71.001, 106.5015, 142.002,<br>
&gt;      &gt;         177.5025, 213.003,<br>
&gt;      &gt;<br>
&gt;      &gt;              248.5035, 290.29765597229, 356.756654745483,<br>
&gt;     470.547043066406,<br>
&gt;      &gt;<br>
&gt;      &gt;              650.50097442627, 910.090332641602, 1256.33548026123,<br>
&gt;      &gt;         1689.23641728516,<br>
&gt;      &gt;<br>
&gt;      &gt;              2201.77290498047, 2780.47265712891, 3406.50401550293,<br>
&gt;      &gt;         4057.20024763184 ;<br>
&gt;      &gt;<br>
&gt;      &gt;<br>
&gt;      &gt;           potrho = 1028.05, 1028.15, 1028.25, 1028.35, 1028.45,<br>
&gt;     1028.55,<br>
&gt;      &gt;         1028.65,<br>
&gt;      &gt;<br>
&gt;      &gt;              1028.75, 1028.85, 1028.95, 1029.05, 1029.15, 1029.25,<br>
&gt;      &gt;         1029.35, 1029.45,<br>
&gt;      &gt;<br>
&gt;      &gt;              1029.55, 1029.65, 1029.75, 1029.85, 1029.95, 1030.05,<br>
&gt;      &gt;         1030.15, 1030.25,<br>
&gt;      &gt;<br>
&gt;      &gt;              1030.35, 1030.45, 1030.55, 1030.65, 1030.75, 1030.85,<br>
&gt;      &gt;         1030.95, 1031.05,<br>
&gt;      &gt;<br>
&gt;      &gt;              1031.15, 1031.25, 1031.35, 1031.45, 1031.55, 1031.65,<br>
&gt;      &gt;         1031.75, 1031.85,<br>
&gt;      &gt;<br>
&gt;      &gt;              1031.95, 1032.05, 1032.15, 1032.25, 1032.35, 1032.45,<br>
&gt;      &gt;         1032.55, 1032.65,<br>
&gt;      &gt;<br>
&gt;      &gt;              1032.75, 1032.85, 1032.95, 1033.05, 1033.15, 1033.25,<br>
&gt;      &gt;         1033.35, 1033.45,<br>
&gt;      &gt;<br>
&gt;      &gt;              1033.54999999999, 1033.64999999999, 1033.74999999999,<br>
&gt;      &gt;         1033.84999999999,<br>
&gt;      &gt;<br>
&gt;      &gt;              1033.94999999999, 1034.04999999999, 1034.14999999999,<br>
&gt;      &gt;         1034.24999999999,<br>
&gt;      &gt;<br>
&gt;      &gt;              1034.34999999999, 1034.44999999999, 1034.54999999999,<br>
&gt;      &gt;         1034.64999999999,<br>
&gt;      &gt;<br>
&gt;      &gt;              1034.74999999999, 1034.84999999999, 1034.94999999999,<br>
&gt;      &gt;         1035.04999999999,<br>
&gt;      &gt;<br>
&gt;      &gt;              1035.14999999999, 1035.24999999999, 1035.34999999999,<br>
&gt;      &gt;         1035.44999999999,<br>
&gt;      &gt;<br>
&gt;      &gt;              1035.54999999999, 1035.64999999999, 1035.74999999999,<br>
&gt;      &gt;         1035.84999999999,<br>
&gt;      &gt;<br>
&gt;      &gt;              1035.94999999999, 1036.04999999999, 1036.14999999999,<br>
&gt;      &gt;         1036.24999999999,<br>
&gt;      &gt;<br>
&gt;      &gt;              1036.34999999999, 1036.44999999999, 1036.54999999999,<br>
&gt;      &gt;         1036.64999999999,<br>
&gt;      &gt;<br>
&gt;      &gt;              1036.74999999999, 1036.84999999999, 1036.94999999999,<br>
&gt;      &gt;         1037.04999999999,<br>
&gt;      &gt;<br>
&gt;      &gt;              1037.14999999999, 1037.24999999999, 1037.34999999999,<br>
&gt;      &gt;         1037.44999999999,<br>
&gt;      &gt;<br>
&gt;      &gt;              1037.54999999999, 1037.64999999999, 1037.74999999999,<br>
&gt;      &gt;         1037.84999999999,<br>
&gt;      &gt;<br>
&gt;      &gt;              1037.94999999999 ;<br>
&gt;      &gt;<br>
&gt;      &gt;<br>
&gt;      &gt;           potrho_edges = 1028, 1028.1, 1028.2, 1028.3, 1028.4,<br>
&gt;     1028.5,<br>
&gt;      &gt;         1028.6, 1028.7,<br>
&gt;      &gt;<br>
&gt;      &gt;              1028.8, 1028.9, 1029, 1029.1, 1029.2, 1029.3, 1029.4,<br>
&gt;      &gt;         1029.5, 1029.6,<br>
&gt;      &gt;<br>
&gt;      &gt;              1029.7, 1029.8, 1029.9, 1030, 1030.1, 1030.2, 1030.3,<br>
&gt;      &gt;         1030.4, 1030.5,<br>
&gt;      &gt;<br>
&gt;      &gt;              1030.6, 1030.7, 1030.8, 1030.9, 1031, 1031.1, 1031.2,<br>
&gt;      &gt;         1031.3, 1031.4,<br>
&gt;      &gt;<br>
&gt;      &gt;              1031.5, 1031.6, 1031.7, 1031.8, 1031.9, 1032, 1032.1,<br>
&gt;      &gt;         1032.2, 1032.3,<br>
&gt;      &gt;<br>
&gt;      &gt;              1032.4, 1032.5, 1032.6, 1032.7, 1032.8, 1032.9, 1033,<br>
&gt;      &gt;         1033.1, 1033.2,<br>
&gt;      &gt;<br>
&gt;      &gt;              1033.3, 1033.4, 1033.49999999999, 1033.59999999999,<br>
&gt;      &gt;         1033.69999999999,<br>
&gt;      &gt;<br>
&gt;      &gt;              1033.79999999999, 1033.89999999999, 1033.99999999999,<br>
&gt;      &gt;         1034.09999999999,<br>
&gt;      &gt;<br>
&gt;      &gt;              1034.19999999999, 1034.29999999999, 1034.39999999999,<br>
&gt;      &gt;         1034.49999999999,<br>
&gt;      &gt;<br>
&gt;      &gt;              1034.59999999999, 1034.69999999999, 1034.79999999999,<br>
&gt;      &gt;         1034.89999999999,<br>
&gt;      &gt;<br>
&gt;      &gt;              1034.99999999999, 1035.09999999999, 1035.19999999999,<br>
&gt;      &gt;         1035.29999999999,<br>
&gt;      &gt;<br>
&gt;      &gt;              1035.39999999999, 1035.49999999999, 1035.59999999999,<br>
&gt;      &gt;         1035.69999999999,<br>
&gt;      &gt;<br>
&gt;      &gt;              1035.79999999999, 1035.89999999999, 1035.99999999999,<br>
&gt;      &gt;         1036.09999999999,<br>
&gt;      &gt;<br>
&gt;      &gt;              1036.19999999999, 1036.29999999999, 1036.39999999999,<br>
&gt;      &gt;         1036.49999999999,<br>
&gt;      &gt;<br>
&gt;      &gt;              1036.59999999999, 1036.69999999999, 1036.79999999999,<br>
&gt;      &gt;         1036.89999999999,<br>
&gt;      &gt;<br>
&gt;      &gt;              1036.99999999999, 1037.09999999999, 1037.19999999999,<br>
&gt;      &gt;         1037.29999999999,<br>
&gt;      &gt;<br>
&gt;      &gt;              1037.39999999999, 1037.49999999999, 1037.59999999999,<br>
&gt;      &gt;         1037.69999999999,<br>
&gt;      &gt;<br>
&gt;      &gt;              1037.79999999999, 1037.89999999999, 1037.99999999999 ;<br>
&gt;      &gt;<br>
&gt;      &gt;         }<br>
&gt;      &gt;<br>
&gt;      &gt;<br>
&gt;      &gt;         I am having trouble deciding which will be the xdef,<br>
&gt;     ydef, zdef,<br>
&gt;      &gt;         and tdef values. Also, how do I find what UNDEF number to<br>
&gt;     use.<br>
&gt;      &gt;         Your help will be greatly appreciated.<br>
&gt;      &gt;<br>
&gt;      &gt;         Thanks,<br>
&gt;      &gt;<br>
&gt;      &gt;         David<br>
&gt;      &gt;<br>
&gt;      &gt;<br>
&gt;      &gt;<br>
&gt;      &gt;<br>
&gt;      &gt;<br>
&gt;      &gt; _______________________________________________<br>
&gt;      &gt; gradsusr mailing list<br>
</div></div>&gt;      &gt; <a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a> &lt;mailto:<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>&gt;<br>
<div class="">&gt;      &gt; <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
&gt;      &gt;<br>
&gt;<br>
&gt;     --<br>
&gt;<br>
&gt;     Please note that <a href="mailto:Charles.Seman@noaa.gov">Charles.Seman@noaa.gov</a><br>
</div>&gt;     &lt;mailto:<a href="mailto:Charles.Seman@noaa.gov">Charles.Seman@noaa.gov</a>&gt; should be considered my NOAA<br>
&gt;     email address, not <a href="mailto:cjs@gfdl.noaa.gov">cjs@gfdl.noaa.gov</a> &lt;mailto:<a href="mailto:cjs@gfdl.noaa.gov">cjs@gfdl.noaa.gov</a>&gt;.<br>
&gt;<br>
&gt;     ********************************************************************<br>
&gt;        Charles Seman <a href="mailto:Charles.Seman@noaa.gov">Charles.Seman@noaa.gov</a> &lt;mailto:<a href="mailto:Charles.Seman@noaa.gov">Charles.Seman@noaa.gov</a>&gt;<br>
<div class="">&gt;        U.S. Department of Commerce / NOAA / OAR<br>
&gt;        Geophysical Fluid Dynamics Laboratory         voice: (609) 452-6547<br>
&gt;        201 Forrestal Road                              fax: (609) 987-5063<br>
&gt;        Princeton, NJ  08540-6649 <a href="http://www.gfdl.noaa.gov/~cjs/" target="_blank">http://www.gfdl.noaa.gov/~cjs/</a><br>
&gt;     ********************************************************************<br>
&gt;<br>
&gt;     &quot;The contents of this message are mine personally and do not reflect any<br>
&gt;     official or unofficial position of the United States Federal Government,<br>
&gt;     the United States Department of Commerce, or NOAA.&quot;<br>
&gt;     _______________________________________________<br>
&gt;     gradsusr mailing list<br>
</div>&gt;     <a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a> &lt;mailto:<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>&gt;<br>
&gt;     <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
<div class="HOEnZb"><div class="h5">&gt;<br>
&gt;<br>
&gt;<br>
&gt;<br>
&gt; _______________________________________________<br>
&gt; gradsusr mailing list<br>
&gt; <a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
&gt; <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
&gt;<br>
<br>
--<br>
<br>
Please note that <a href="mailto:Charles.Seman@noaa.gov">Charles.Seman@noaa.gov</a> should be considered my NOAA<br>
email address, not <a href="mailto:cjs@gfdl.noaa.gov">cjs@gfdl.noaa.gov</a>.<br>
<br>
********************************************************************<br>
  Charles Seman                                <a href="mailto:Charles.Seman@noaa.gov">Charles.Seman@noaa.gov</a><br>
  U.S. Department of Commerce / NOAA / OAR<br>
  Geophysical Fluid Dynamics Laboratory         voice: (609) 452-6547<br>
  201 Forrestal Road                              fax: (609) 987-5063<br>
  Princeton, NJ  08540-6649            <a href="http://www.gfdl.noaa.gov/~cjs/" target="_blank">http://www.gfdl.noaa.gov/~cjs/</a><br>
********************************************************************<br>
<br>
&quot;The contents of this message are mine personally and do not reflect any<br>
official or unofficial position of the United States Federal Government,<br>
the United States Department of Commerce, or NOAA.&quot;<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>
</div></div></blockquote></div><br></div>