<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"><<a href="mailto:charles.seman@noaa.gov" target="_blank">charles.seman@noaa.gov</a>></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 "frerun" with an xml file? If so, there is a companion<br>
"frepp" 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>
> It is GFDL Model data using the Molecular Ocean Model MOM_solo type.<br>
> This is for the box1 test case. I did try sdfopen and it came out with<br>
> the following error message:<br>
><br>
> Scanning self-describing file: <a href="http://19800101.ocean_neutral.nc" target="_blank">19800101.ocean_neutral.nc</a><br>
</div>> <<a href="http://19800101.ocean_neutral.nc" target="_blank">http://19800101.ocean_neutral.nc</a>><br>
<div class="">><br>
> SDF Error: 365 day calendars are no longer supported by sdfopen.<br>
><br>
> To open this file with GrADS, use a descriptor file with<br>
><br>
> a complete TDEF entry and OPTIONS 365_day_calendar.<br>
><br>
><br>
> I do not have any pp data and I don't know how to go about creating the<br>
> grid data. How could I do that, putting it into lat, long coordinates?<br>
><br>
> Thanks,<br>
><br>
> David<br>
><br>
><br>
> On Fri, Jul 25, 2014 at 2:39 PM, Charles Seman - NOAA Federal<br>
</div><div class="">> <<a href="mailto:charles.seman@noaa.gov">charles.seman@noaa.gov</a> <mailto:<a href="mailto:charles.seman@noaa.gov">charles.seman@noaa.gov</a>>> wrote:<br>
><br>
> David,<br>
><br>
> This file looks like it's from a GFDL model? Based on the names of the<br>
> coordinates, and the global attributes "grid_type" and "grid_tile", it<br>
> looks like it is a "history" file which is on the model grid and has yet<br>
> to be regridded to a lat/lon grid. I'm not a user of the ocean data,<br>
> but am working with a GFDL history data file from the atmospheric model<br>
> on the cubed-sphere model grid tile4. Be very careful! The grid<br>
> coordinates xt_ocean and yt_ocean may need to be mapped into actual<br>
> global coordinates. To get a view of the data, have you tried opening<br>
> the file in GrADS without using a GrADS ctl file? GrADS is able to open<br>
> and view the tile4 file I'm working with, but one must be aware that the<br>
> map background that GrADS assigns is not physically correct. Can you<br>
> use post-processed "pp" data (if you have this, it should be in a "pp"<br>
> subdirectory)? These "pp" data have been put on a physical grid.<br>
> I talked with someone here at GFDL who works with ocean data, and he<br>
> says Ferret is able to open the datasets he works with...<br>
><br>
> Here's an example of what's in a pp/ocean/<a href="http://ocean.static.nc" target="_blank">ocean.static.nc</a><br>
</div>> <<a href="http://ocean.static.nc" target="_blank">http://ocean.static.nc</a>> file. Note<br>
<div><div class="h5">> how the tracer latitude and longitudes are functions of<br>
> yt_ocean,xt_ocean: geolat_t(yt_ocean,xt_ocean) and<br>
> geolon_t(yt_ocean,xt_ocean). So the lat and lon are themselves<br>
> functions of grid coordinate variables. Here's an example of<br>
> two-dimensional lat, lon coordinates:<br>
> <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>
><br>
> netcdf ocean.static {<br>
><br>
> dimensions:<br>
><br>
> xu_ocean = 360 ;<br>
><br>
> yu_ocean = 200 ;<br>
><br>
> xt_ocean = 360 ;<br>
><br>
> yt_ocean = 200 ;<br>
><br>
> variables:<br>
><br>
> double xt_ocean(xt_ocean) ;<br>
><br>
> xt_ocean:long_name = "tcell longitude" ;<br>
><br>
> xt_ocean:units = "degrees_E" ;<br>
><br>
> xt_ocean:cartesian_axis = "X" ;<br>
><br>
> double yt_ocean(yt_ocean) ;<br>
><br>
> yt_ocean:long_name = "tcell latitude" ;<br>
><br>
> yt_ocean:units = "degrees_N" ;<br>
><br>
> yt_ocean:cartesian_axis = "Y" ;<br>
><br>
> double xu_ocean(xu_ocean) ;<br>
><br>
> xu_ocean:long_name = "ucell longitude" ;<br>
><br>
> xu_ocean:units = "degrees_E" ;<br>
><br>
> xu_ocean:cartesian_axis = "X" ;<br>
><br>
> double yu_ocean(yu_ocean) ;<br>
><br>
> yu_ocean:long_name = "ucell latitude" ;<br>
><br>
> yu_ocean:units = "degrees_N" ;<br>
><br>
> yu_ocean:cartesian_axis = "Y" ;<br>
><br>
> float ht(yt_ocean, xt_ocean) ;<br>
><br>
> ht:long_name = "ocean depth on t-cells" ;<br>
><br>
> ht:units = "m" ;<br>
><br>
> ht:valid_range = -1.e+09f, 1.e+09f ;<br>
><br>
> ht:missing_value = -1.e+20f ;<br>
><br>
> ht:_FillValue = -1.e+20f ;<br>
><br>
> ht:cell_methods = "time: point" ;<br>
><br>
> ht:coordinates = "geolon_t geolat_t" ;<br>
><br>
> ht:standard_name = "sea_floor_depth_below_geoid" ;<br>
><br>
> float hu(yu_ocean, xu_ocean) ;<br>
><br>
> hu:long_name = "ocean depth on u-cells" ;<br>
><br>
> hu:units = "m" ;<br>
><br>
> hu:valid_range = -1.e+09f, 1.e+09f ;<br>
><br>
> hu:missing_value = -1.e+20f ;<br>
><br>
> hu:_FillValue = -1.e+20f ;<br>
><br>
> hu:cell_methods = "time: point" ;<br>
><br>
> hu:coordinates = "geolon_c geolat_c" ;<br>
><br>
> float geolat_c(yu_ocean, xu_ocean) ;<br>
><br>
> geolat_c:long_name = "uv latitude" ;<br>
><br>
> geolat_c:units = "degrees_N" ;<br>
><br>
> geolat_c:valid_range = -91.f, 91.f ;<br>
><br>
> geolat_c:missing_value = 1.e+20f ;<br>
><br>
> geolat_c:_FillValue = 1.e+20f ;<br>
><br>
> geolat_c:cell_methods = "time: point" ;<br>
><br>
> geolat_c:coordinates = "geolon_c geolat_c" ;<br>
><br>
> float geolon_c(yu_ocean, xu_ocean) ;<br>
><br>
> geolon_c:long_name = "uv longitude" ;<br>
><br>
> geolon_c:units = "degrees_E" ;<br>
><br>
> geolon_c:valid_range = -281.f, 361.f ;<br>
><br>
> geolon_c:missing_value = 1.e+20f ;<br>
><br>
> geolon_c:_FillValue = 1.e+20f ;<br>
><br>
> geolon_c:cell_methods = "time: point" ;<br>
><br>
> geolon_c:coordinates = "geolon_c geolat_c" ;<br>
><br>
> float geolat_t(yt_ocean, xt_ocean) ;<br>
><br>
> geolat_t:long_name = "tracer latitude" ;<br>
><br>
> geolat_t:units = "degrees_N" ;<br>
><br>
> geolat_t:valid_range = -91.f, 91.f ;<br>
><br>
> geolat_t:missing_value = 1.e+20f ;<br>
><br>
> geolat_t:_FillValue = 1.e+20f ;<br>
><br>
> geolat_t:cell_methods = "time: point" ;<br>
><br>
> geolat_t:coordinates = "geolon_t geolat_t" ;<br>
><br>
> float geolon_t(yt_ocean, xt_ocean) ;<br>
><br>
> geolon_t:long_name = "tracer longitude" ;<br>
><br>
> geolon_t:units = "degrees_E" ;<br>
><br>
> geolon_t:valid_range = -281.f, 361.f ;<br>
><br>
> geolon_t:missing_value = 1.e+20f ;<br>
><br>
> geolon_t:_FillValue = 1.e+20f ;<br>
><br>
> geolon_t:cell_methods = "time: point" ;<br>
><br>
> geolon_t:coordinates = "geolon_t geolat_t" ;<br>
><br>
> float area_t(yt_ocean, xt_ocean) ;<br>
><br>
> area_t:long_name = "tracer cell area" ;<br>
><br>
> area_t:units = "m^2" ;<br>
><br>
> area_t:valid_range = 0.f, 1.e+15f ;<br>
><br>
> area_t:missing_value = 1.e+20f ;<br>
><br>
> area_t:_FillValue = 1.e+20f ;<br>
><br>
> area_t:cell_methods = "time: point" ;<br>
><br>
> area_t:coordinates = "geolon_t geolat_t" ;<br>
><br>
> float area_u(yu_ocean, xu_ocean) ;<br>
><br>
> area_u:long_name = "velocity cell area" ;<br>
><br>
> area_u:units = "m^2" ;<br>
><br>
> area_u:valid_range = 0.f, 1.e+15f ;<br>
><br>
> area_u:missing_value = 1.e+20f ;<br>
><br>
> area_u:_FillValue = 1.e+20f ;<br>
><br>
> area_u:cell_methods = "time: point" ;<br>
><br>
> area_u:coordinates = "geolon_c geolat_c" ;<br>
><br>
> float kmt(yt_ocean, xt_ocean) ;<br>
><br>
> kmt:long_name = "number of depth levels on t-grid" ;<br>
><br>
> kmt:units = "dimensionless" ;<br>
><br>
> kmt:valid_range = -10.f, 1.e+09f ;<br>
><br>
> kmt:missing_value = -1.e+20f ;<br>
><br>
> kmt:_FillValue = -1.e+20f ;<br>
><br>
> kmt:cell_methods = "time: point" ;<br>
><br>
> kmt:coordinates = "geolon_t geolat_t" ;<br>
><br>
> float kmu(yu_ocean, xu_ocean) ;<br>
><br>
> kmu:long_name = "number of depth levels on u-grid" ;<br>
><br>
> kmu:units = "dimensionless" ;<br>
><br>
> kmu:valid_range = -10.f, 1.e+09f ;<br>
><br>
> kmu:missing_value = -1.e+20f ;<br>
><br>
> kmu:_FillValue = -1.e+20f ;<br>
><br>
> kmu:cell_methods = "time: point" ;<br>
><br>
> kmu:coordinates = "geolon_c geolat_c" ;<br>
><br>
> float geo_heat(yt_ocean, xt_ocean) ;<br>
><br>
> geo_heat:long_name = "Geothermal heating" ;<br>
><br>
> geo_heat:units = "W/m^2" ;<br>
><br>
> geo_heat:valid_range = -10.f, 1000000.f ;<br>
><br>
> geo_heat:missing_value = -1.e+20f ;<br>
><br>
> geo_heat:_FillValue = -1.e+20f ;<br>
><br>
> geo_heat:cell_methods = "time: point" ;<br>
><br>
> geo_heat:coordinates = "geolon_t geolat_t" ;<br>
><br>
> geo_heat:standard_name =<br>
> "upward_geothermal_heat_flux_at_sea_floor" ;<br>
><br>
> float pbot0(yt_ocean, xt_ocean) ;<br>
><br>
> pbot0:long_name = "reference bottom pressure<br>
> t-cells" ;<br>
><br>
> pbot0:units = "dbar" ;<br>
><br>
> pbot0:valid_range = -1.e+08f, 1.e+08f ;<br>
> pbot0:missing_value = -1.e+20f ;<br>
> pbot0:_FillValue = -1.e+20f ;<br>
> pbot0:cell_methods = "time: point" ;<br>
> pbot0:coordinates = "geolon_t geolat_t" ;<br>
><br>
> // global attributes:<br>
> :NCX = "1.1.1" ;<br>
> }<br>
><br>
> I hope this is helpful. I don't think I can be of much more help, but<br>
> hope this points you in the right direction.<br>
><br>
> Chuck<br>
><br>
> On 07/25/2014 12:26 AM, David Stanfel wrote:<br>
> > Hello,<br>
> ><br>
> > I decided to try a different file. Here is the ncdump -c of said<br>
> file:<br>
> ><br>
> > dstanfel:bin ninjawarrior343$ ncdump -c <a href="http://19800101.ocean_neutral.nc" target="_blank">19800101.ocean_neutral.nc</a><br>
> <<a href="http://19800101.ocean_neutral.nc" target="_blank">http://19800101.ocean_neutral.nc</a>><br>
> > <<a href="http://19800101.ocean_neutral.nc" target="_blank">http://19800101.ocean_neutral.nc</a>><br>
> ><br>
> > netcdf \19800101.ocean_neutral {<br>
> ><br>
> > dimensions:<br>
> ><br>
> > xt_ocean = 24 ;<br>
> ><br>
> > yt_ocean = 35 ;<br>
> ><br>
> > Time = UNLIMITED ; // (1 currently)<br>
> ><br>
> > nv = 2 ;<br>
> ><br>
> > st_ocean = 18 ;<br>
> ><br>
> > st_edges_ocean = 19 ;<br>
> ><br>
> > variables:<br>
> ><br>
> > double xt_ocean(xt_ocean) ;<br>
> ><br>
> > xt_ocean:long_name = "tcell longitude" ;<br>
> ><br>
> > xt_ocean:units = "degrees_E" ;<br>
> ><br>
> > xt_ocean:cartesian_axis = "X" ;<br>
> ><br>
> > double yt_ocean(yt_ocean) ;<br>
> ><br>
> > yt_ocean:long_name = "tcell latitude" ;<br>
> ><br>
> > yt_ocean:units = "degrees_N" ;<br>
> ><br>
> > yt_ocean:cartesian_axis = "Y" ;<br>
> ><br>
> > double Time(Time) ;<br>
> ><br>
> > Time:long_name = "Time" ;<br>
> ><br>
> > Time:units = "days since 1980-01-01 00:00:00" ;<br>
> ><br>
> > Time:cartesian_axis = "T" ;<br>
> ><br>
> > Time:calendar_type = "NOLEAP" ;<br>
> ><br>
> > Time:calendar = "NOLEAP" ;<br>
> ><br>
> > Time:bounds = "Time_bounds" ;<br>
> ><br>
> > double nv(nv) ;<br>
> ><br>
> > nv:long_name = "vertex number" ;<br>
> ><br>
> > nv:units = "none" ;<br>
> ><br>
> > nv:cartesian_axis = "N" ;<br>
> ><br>
> > double st_ocean(st_ocean) ;<br>
> ><br>
> > st_ocean:long_name = "tcell pstar" ;<br>
> ><br>
> > st_ocean:units = "dbars" ;<br>
> ><br>
> > st_ocean:cartesian_axis = "Z" ;<br>
> ><br>
> > st_ocean:positive = "down" ;<br>
> ><br>
> > st_ocean:edges = "st_edges_ocean" ;<br>
> ><br>
> > double st_edges_ocean(st_edges_ocean) ;<br>
> ><br>
> > st_edges_ocean:long_name = "tcell pstar edges" ;<br>
> ><br>
> > st_edges_ocean:units = "dbars" ;<br>
> ><br>
> > st_edges_ocean:cartesian_axis = "Z" ;<br>
> ><br>
> > st_edges_ocean:positive = "down" ;<br>
> ><br>
> > float agm(Time, yt_ocean, xt_ocean) ;<br>
> ><br>
> > agm:long_name = "GM diffusivity at surface" ;<br>
> ><br>
> > agm:units = "m^2/sec" ;<br>
> ><br>
> > agm:valid_range = -10.f, 1.e+10f ;<br>
> ><br>
> > agm:missing_value = -1.e+20f ;<br>
> ><br>
> > agm:_FillValue = -1.e+20f ;<br>
> ><br>
> > agm:cell_methods = "time: mean" ;<br>
> ><br>
> > agm:time_avg_info = "average_T1,average_T2,average_DT" ;<br>
> ><br>
> > agm:coordinates = "geolon_t geolat_t" ;<br>
> ><br>
> > agm:standard_name = "ocean_tracer_bolus_laplacian_diffusivity" ;<br>
> ><br>
> > float rossby(Time, yt_ocean, xt_ocean) ;<br>
> ><br>
> > rossby:long_name = "Rossby radius used in neutral physics" ;<br>
> ><br>
> > rossby:units = "m" ;<br>
> ><br>
> > rossby:valid_range = -10.f, 1.e+10f ;<br>
> ><br>
> > rossby:missing_value = -1.e+20f ;<br>
> ><br>
> > rossby:_FillValue = -1.e+20f ;<br>
> ><br>
> > rossby:cell_methods = "time: mean" ;<br>
> ><br>
> > rossby:time_avg_info = "average_T1,average_T2,average_DT" ;<br>
> ><br>
> > rossby:coordinates = "geolon_t geolat_t" ;<br>
> ><br>
> > float slope31(Time, st_ocean, yt_ocean, xt_ocean) ;<br>
> ><br>
> > slope31:long_name = "neutral slope -(rho_x/rho_z)" ;<br>
> ><br>
> > slope31:units = "dimensionless" ;<br>
> ><br>
> > slope31:valid_range = -1.e+10f, 1.e+10f ;<br>
> ><br>
> > slope31:missing_value = -1.e+20f ;<br>
> ><br>
> > slope31:_FillValue = -1.e+20f ;<br>
> ><br>
> > slope31:cell_methods = "time: mean" ;<br>
> ><br>
> > slope31:time_avg_info = "average_T1,average_T2,average_DT" ;<br>
> ><br>
> > slope31:coordinates = "geolon_t geolat_t" ;<br>
> ><br>
> > float slope32(Time, st_ocean, yt_ocean, xt_ocean) ;<br>
> ><br>
> > slope32:long_name = "neutral slope -(rho_y/rho_z)" ;<br>
> ><br>
> > slope32:units = "dimensionless" ;<br>
> ><br>
> > slope32:valid_range = -1.e+10f, 1.e+10f ;<br>
> ><br>
> > slope32:missing_value = -1.e+20f ;<br>
> ><br>
> > slope32:_FillValue = -1.e+20f ;<br>
> ><br>
> > slope32:cell_methods = "time: mean" ;<br>
> ><br>
> > slope32:time_avg_info = "average_T1,average_T2,average_DT" ;<br>
> ><br>
> > slope32:coordinates = "geolon_t geolat_t" ;<br>
> ><br>
> > double average_T1(Time) ;<br>
> ><br>
> > average_T1:long_name = "Start time for average period" ;<br>
> ><br>
> > average_T1:units = "days since 1980-01-01 00:00:00" ;<br>
> ><br>
> > average_T1:missing_value = 1.e+20 ;<br>
> ><br>
> > average_T1:_FillValue = 1.e+20 ;<br>
> ><br>
> > double average_T2(Time) ;<br>
> ><br>
> > average_T2:long_name = "End time for average period" ;<br>
> ><br>
> > average_T2:units = "days since 1980-01-01 00:00:00" ;<br>
> ><br>
> > average_T2:missing_value = 1.e+20 ;<br>
> ><br>
> > average_T2:_FillValue = 1.e+20 ;<br>
> ><br>
> > double average_DT(Time) ;<br>
> ><br>
> > average_DT:long_name = "Length of average period" ;<br>
> ><br>
> > average_DT:units = "days" ;<br>
> ><br>
> > average_DT:missing_value = 1.e+20 ;<br>
> ><br>
> > average_DT:_FillValue = 1.e+20 ;<br>
> ><br>
> > double Time_bounds(Time, nv) ;<br>
> ><br>
> > Time_bounds:long_name = "Time axis boundaries" ;<br>
> ><br>
> > Time_bounds:units = "days" ;<br>
> ><br>
> > Time_bounds:missing_value = 1.e+20 ;<br>
> ><br>
> > Time_bounds:_FillValue = 1.e+20 ;<br>
> ><br>
> ><br>
> > // global attributes:<br>
> ><br>
> > :filename = "<a href="http://ocean_neutral.nc" target="_blank">ocean_neutral.nc</a> <<a href="http://ocean_neutral.nc" target="_blank">http://ocean_neutral.nc</a>><br>
> <<a href="http://ocean_neutral.nc" target="_blank">http://ocean_neutral.nc</a>>" ;<br>
> ><br>
> > :title = "MOM4 box" ;<br>
> ><br>
> > :grid_type = "regular" ;<br>
> ><br>
> > :grid_tile = "N/A" ;<br>
> ><br>
> > data:<br>
> ><br>
> ><br>
> > xt_ocean = 1.25, 3.75, 6.25, 8.75, 11.25, 13.75, 16.25, 18.75,<br>
> 21.25,<br>
> > 23.75,<br>
> ><br>
> > 26.25, 28.75, 31.25, 33.75, 36.25, 38.75, 41.25, 43.75,<br>
> 46.25, 48.75,<br>
> ><br>
> > 51.25, 53.75, 56.25, 58.75 ;<br>
> ><br>
> ><br>
> > yt_ocean = 1.25, 3.7488100528717, 6.2393012046814,<br>
> 8.71919345855713,<br>
> ><br>
> > 11.1815176010132, 13.6242818832397, 16.0408229827881,<br>
> > 18.4296073913574,<br>
> ><br>
> > 20.784387588501, 23.1042194366455, 25.3833503723145,<br>
> 27.6215190887451,<br>
> ><br>
> > 29.813497543335, 31.9597606658936, 34.0555992126465,<br>
> 36.1022338867188,<br>
> ><br>
> > 38.0954322814941, 40.0371551513672, 41.9235687255859,<br>
> > 43.7573394775391,<br>
> ><br>
> > 45.5349464416504, 47.259708404541, 48.9283294677734,<br>
> 50.5447235107422,<br>
> ><br>
> > 52.1057090759277, 53.6157531738281, 55.0716934204102,<br>
> > 56.4785079956055,<br>
> ><br>
> > 57.8329429626465, 59.1404571533203, 60.3976211547852,<br>
> > 61.6103477478027,<br>
> ><br>
> > 62.7749481201172, 63.8977813720703, 64.9748153686523 ;<br>
> ><br>
> ><br>
> > Time = 2 ;<br>
> ><br>
> ><br>
> > nv = 1, 2 ;<br>
> ><br>
> ><br>
> > st_ocean = 17.75025, 53.25075, 88.75125, 124.25175, 159.75225,<br>
> 195.25275,<br>
> ><br>
> > 230.75325, 266.25375, 314.34159289856, 399.171716592407,<br>
> ><br>
> > 541.922400494385, 759.079548358154, 1061.10117883301,<br>
> > 1451.56990550537,<br>
> ><br>
> > 1926.90292906494, 2476.64250944824, 3084.30280480957,<br>
> > 3728.70522619629 ;<br>
> ><br>
> ><br>
> > st_edges_ocean = 0, 35.5005, 71.001, 106.5015, 142.002, 177.5025,<br>
> > 213.003,<br>
> ><br>
> > 248.5035, 290.29765597229, 356.756654745483, 470.547043066406,<br>
> ><br>
> > 650.50097442627, 910.090332641602, 1256.33548026123,<br>
> 1689.23641728516,<br>
> ><br>
> > 2201.77290498047, 2780.47265712891, 3406.50401550293,<br>
> 4057.20024763184<br>
> ><br>
> ><br>
> > I wrote a .ctl file:<br>
> ><br>
> ><br>
> > DSET ^<a href="http://19800101.ocean_neutral.nc" target="_blank">19800101.ocean_neutral.nc</a><br>
</div></div>> <<a href="http://19800101.ocean_neutral.nc" target="_blank">http://19800101.ocean_neutral.nc</a>> <<a href="http://19800101.ocean_neutral.nc" target="_blank">http://19800101.ocean_neutral.nc</a>><br>
<div><div class="h5">> ><br>
> > UNDEF -9.99e8 missing_value<br>
> ><br>
> > DTYPE netcdf<br>
> ><br>
> > CACHESIZE 100000<br>
> ><br>
> > TITLE Sample Model Output<br>
> ><br>
> > XDEF 24 LINEAR 1.25 2.5<br>
> ><br>
> > YDEF 35 levels 1.25, 3.7488100528717, 6.2393012046814,<br>
> 8.71919345855713,<br>
> ><br>
> > 11.1815176010132, 13.6242818832397, 16.0408229827881,<br>
> > 18.4296073913574,<br>
> ><br>
> > 20.784387588501, 23.1042194366455, 25.3833503723145,<br>
> 27.6215190887451,<br>
> ><br>
> > 29.813497543335, 31.9597606658936, 34.0555992126465,<br>
> 36.1022338867188,<br>
> ><br>
> > 38.0954322814941, 40.0371551513672, 41.9235687255859,<br>
> > 43.7573394775391,<br>
> ><br>
> > 45.5349464416504, 47.259708404541, 48.9283294677734,<br>
> 50.5447235107422,<br>
> ><br>
> > 52.1057090759277, 53.6157531738281, 55.0716934204102,<br>
> > 56.4785079956055,<br>
> ><br>
> > 57.8329429626465, 59.1404571533203, 60.3976211547852,<br>
> > 61.6103477478027,<br>
> ><br>
> > 62.7749481201172, 63.8977813720703, 64.9748153686523<br>
> ><br>
> > ZDEF 18 levels 17.75025, 53.25075, 88.75125, 124.25175, 159.75225,<br>
> > 195.25275, 230.75327, 266.25375, 314.34159289856, 399.171716592407,<br>
> > 541.922400494385, 759.079548358154, 1061.10117883301,<br>
> 1451.56990550537,<br>
> > 1926.90292906494, 2476.64250944824, 3084.30280480957,<br>
> 3728.70522619629<br>
> ><br>
> > TDEF 1 LINEAR 01JAN0001 1mo<br>
> ><br>
> > VARS 4<br>
> ><br>
> > agm 0 t,y,x GM diffusivity at surface (m^2/sec)<br>
> ><br>
> > rossby 0 t,y,x Rossby radius used in neutral physics (m)<br>
> ><br>
> > slope31 18 t,z,y,x Neutral slope -(rho_x/rho_z)<br>
> ><br>
> > slope32 18 t,z,y,x Neutral slope -(rho_y/rho_z)<br>
> ><br>
> > ENDVARS<br>
> ><br>
> ><br>
> > When trying to display the variables, I get this error:<br>
> ><br>
> ><br>
> > ga-> d agm<br>
> ><br>
> > Cannot contour grid - all undefined values<br>
> ><br>
> > ga-> d rossby<br>
> ><br>
> > Cannot contour grid - all undefined values<br>
> ><br>
> > ga-> d slope31<br>
> ><br>
> > Cannot contour grid - all undefined values<br>
> ><br>
> > ga-><br>
> ><br>
> ><br>
> > What could be going wrong here?<br>
> ><br>
> > Thanks,<br>
> ><br>
> > David<br>
> ><br>
> ><br>
> > On Wed, Jul 23, 2014 at 11:23 AM, David Stanfel<br>
> <<a href="mailto:dstanfel@ualberta.ca">dstanfel@ualberta.ca</a> <mailto:<a href="mailto:dstanfel@ualberta.ca">dstanfel@ualberta.ca</a>><br>
</div></div><div class="">> > <mailto:<a href="mailto:dstanfel@ualberta.ca">dstanfel@ualberta.ca</a> <mailto:<a href="mailto:dstanfel@ualberta.ca">dstanfel@ualberta.ca</a>>>> wrote:<br>
> ><br>
> > Also, how do you determine the CACHESIZE?<br>
> ><br>
> > Thanks,<br>
> ><br>
> > David<br>
> ><br>
> ><br>
> > On Tue, Jul 22, 2014 at 11:05 PM, David Stanfel<br>
> > <<a href="mailto:dstanfel@ualberta.ca">dstanfel@ualberta.ca</a> <mailto:<a href="mailto:dstanfel@ualberta.ca">dstanfel@ualberta.ca</a>><br>
</div><div><div class="h5">> <mailto:<a href="mailto:dstanfel@ualberta.ca">dstanfel@ualberta.ca</a> <mailto:<a href="mailto:dstanfel@ualberta.ca">dstanfel@ualberta.ca</a>>>> wrote:<br>
> ><br>
> > Hello,<br>
> ><br>
> > I just ran a climate model which produced some .nc files.<br>
> I am<br>
> > having trouble writing the .ctl file to read said .nc files.<br>
> > Here is the ncdump -c for one of the files I am trying to<br>
> read:<br>
> ><br>
> > dstanfel:bin ninjawarrior343$ ncdump -c<br>
> <a href="http://19800101.ocean_layer.nc" target="_blank">19800101.ocean_layer.nc</a> <<a href="http://19800101.ocean_layer.nc" target="_blank">http://19800101.ocean_layer.nc</a>><br>
> > <<a href="http://19800101.ocean_layer.nc" target="_blank">http://19800101.ocean_layer.nc</a>><br>
> ><br>
> > netcdf \19800101.ocean_layer {<br>
> ><br>
> > dimensions:<br>
> ><br>
> > grid_xt_ocean = 24 ;<br>
> ><br>
> > grid_yt_ocean = 35 ;<br>
> ><br>
> > theta = 100 ;<br>
> ><br>
> > theta_edges = 101 ;<br>
> ><br>
> > Time = UNLIMITED ; // (1 currently)<br>
> ><br>
> > nv = 2 ;<br>
> ><br>
> > xt_ocean = 24 ;<br>
> ><br>
> > yt_ocean = 35 ;<br>
> ><br>
> > st_ocean = 18 ;<br>
> ><br>
> > st_edges_ocean = 19 ;<br>
> ><br>
> > potrho = 100 ;<br>
> ><br>
> > potrho_edges = 101 ;<br>
> ><br>
> > variables:<br>
> ><br>
> > double grid_xt_ocean(grid_xt_ocean) ;<br>
> ><br>
> > grid_xt_ocean:long_name = "tcell longitude" ;<br>
> ><br>
> > grid_xt_ocean:units = "degrees_E" ;<br>
> ><br>
> > grid_xt_ocean:cartesian_axis = "X" ;<br>
> ><br>
> > double grid_yt_ocean(grid_yt_ocean) ;<br>
> ><br>
> > grid_yt_ocean:long_name = "tcell latitude" ;<br>
> ><br>
> > grid_yt_ocean:units = "degrees_N" ;<br>
> ><br>
> > grid_yt_ocean:cartesian_axis = "Y" ;<br>
> ><br>
> > double theta(theta) ;<br>
> ><br>
> > theta:long_name = "potential or conservative temperature" ;<br>
> ><br>
> > theta:units = "C" ;<br>
> ><br>
> > theta:cartesian_axis = "Z" ;<br>
> ><br>
> > theta:positive = "up" ;<br>
> ><br>
> > theta:edges = "theta_edges" ;<br>
> ><br>
> > double theta_edges(theta_edges) ;<br>
> ><br>
> > theta_edges:long_name = "potential or conservative<br>
> temperature<br>
> > edges" ;<br>
> ><br>
> > theta_edges:units = "C" ;<br>
> ><br>
> > theta_edges:cartesian_axis = "Z" ;<br>
> ><br>
> > theta_edges:positive = "up" ;<br>
> ><br>
> > double Time(Time) ;<br>
> ><br>
> > Time:long_name = "Time" ;<br>
> ><br>
> > Time:units = "days since 1980-01-01 00:00:00" ;<br>
> ><br>
> > Time:cartesian_axis = "T" ;<br>
> ><br>
> > Time:calendar_type = "NOLEAP" ;<br>
> ><br>
> > Time:calendar = "NOLEAP" ;<br>
> ><br>
> > Time:bounds = "Time_bounds" ;<br>
> ><br>
> > double nv(nv) ;<br>
> ><br>
> > nv:long_name = "vertex number" ;<br>
> ><br>
> > nv:units = "none" ;<br>
> ><br>
> > nv:cartesian_axis = "N" ;<br>
> ><br>
> > double xt_ocean(xt_ocean) ;<br>
> ><br>
> > xt_ocean:long_name = "tcell longitude" ;<br>
> ><br>
> > xt_ocean:units = "degrees_E" ;<br>
> ><br>
> > xt_ocean:cartesian_axis = "X" ;<br>
> ><br>
> > double yt_ocean(yt_ocean) ;<br>
> ><br>
> > yt_ocean:long_name = "tcell latitude" ;<br>
> ><br>
> > yt_ocean:units = "degrees_N" ;<br>
> ><br>
> > yt_ocean:cartesian_axis = "Y" ;<br>
> ><br>
> > double st_ocean(st_ocean) ;<br>
> ><br>
> > st_ocean:long_name = "tcell pstar" ;<br>
> ><br>
> > st_ocean:units = "dbars" ;<br>
> ><br>
> > st_ocean:cartesian_axis = "Z" ;<br>
> ><br>
> > st_ocean:positive = "down" ;<br>
> ><br>
> > st_ocean:edges = "st_edges_ocean" ;<br>
> ><br>
> > double st_edges_ocean(st_edges_ocean) ;<br>
> ><br>
> > st_edges_ocean:long_name = "tcell pstar edges" ;<br>
> ><br>
> > st_edges_ocean:units = "dbars" ;<br>
> ><br>
> > st_edges_ocean:cartesian_axis = "Z" ;<br>
> ><br>
> > st_edges_ocean:positive = "down" ;<br>
> ><br>
> > double potrho(potrho) ;<br>
> ><br>
> > potrho:long_name = "potential density" ;<br>
> ><br>
> > potrho:units = "kg/m^3" ;<br>
> ><br>
> > potrho:cartesian_axis = "Z" ;<br>
> ><br>
> > potrho:positive = "down" ;<br>
> ><br>
> > potrho:edges = "potrho_edges" ;<br>
> ><br>
> > double potrho_edges(potrho_edges) ;<br>
> ><br>
> > potrho_edges:long_name = "potential density edges" ;<br>
> ><br>
> > potrho_edges:units = "kg/m^3" ;<br>
> ><br>
> > potrho_edges:cartesian_axis = "Z" ;<br>
> ><br>
> > potrho_edges:positive = "down" ;<br>
> ><br>
> > float depth_of_theta(Time, theta, grid_yt_ocean,<br>
> grid_xt_ocean) ;<br>
> ><br>
> > depth_of_theta:long_name = "depth of potential temp<br>
> surface" ;<br>
> ><br>
> > depth_of_theta:units = "m" ;<br>
> ><br>
> > depth_of_theta:valid_range = 0.f, 1.e+10f ;<br>
> ><br>
> > depth_of_theta:missing_value = -1.e+20f ;<br>
> ><br>
> > depth_of_theta:_FillValue = -1.e+20f ;<br>
> ><br>
> > depth_of_theta:cell_methods = "time: mean" ;<br>
> ><br>
> > depth_of_theta:time_avg_info =<br>
> "average_T1,average_T2,average_DT" ;<br>
> ><br>
> > float rho(Time, st_ocean, yt_ocean, xt_ocean) ;<br>
> ><br>
> > rho:long_name = "in situ density" ;<br>
> ><br>
> > rho:units = "kg/m^3" ;<br>
> ><br>
> > rho:valid_range = -10.f, 100000.f ;<br>
> ><br>
> > rho:missing_value = -1.e+20f ;<br>
> ><br>
> > rho:_FillValue = -1.e+20f ;<br>
> ><br>
> > rho:cell_methods = "time: mean" ;<br>
> ><br>
> > rho:time_avg_info = "average_T1,average_T2,average_DT" ;<br>
> ><br>
> > rho:coordinates = "geolon_t geolat_t" ;<br>
> ><br>
> > float neutral_rho(Time, st_ocean, yt_ocean, xt_ocean) ;<br>
> ><br>
> > neutral_rho:long_name = "potential density estimate of<br>
> neutral<br>
> > density" ;<br>
> ><br>
> > neutral_rho:units = "kg/m^3" ;<br>
> ><br>
> > neutral_rho:valid_range = -10.f, 100000.f ;<br>
> ><br>
> > neutral_rho:missing_value = -1.e+20f ;<br>
> ><br>
> > neutral_rho:_FillValue = -1.e+20f ;<br>
> ><br>
> > neutral_rho:cell_methods = "time: mean" ;<br>
> ><br>
> > neutral_rho:time_avg_info =<br>
> "average_T1,average_T2,average_DT" ;<br>
> ><br>
> > neutral_rho:coordinates = "geolon_t geolat_t" ;<br>
> ><br>
> > float pot_rho_0(Time, st_ocean, yt_ocean, xt_ocean) ;<br>
> ><br>
> > pot_rho_0:long_name = "potential density referenced to 0<br>
> dbar" ;<br>
> ><br>
> > pot_rho_0:units = "kg/m^3" ;<br>
> ><br>
> > pot_rho_0:valid_range = -10.f, 100000.f ;<br>
> ><br>
> > pot_rho_0:missing_value = -1.e+20f ;<br>
> ><br>
> > pot_rho_0:_FillValue = -1.e+20f ;<br>
> ><br>
> > pot_rho_0:cell_methods = "time: mean" ;<br>
> ><br>
> > pot_rho_0:time_avg_info =<br>
> "average_T1,average_T2,average_DT" ;<br>
> ><br>
> > pot_rho_0:coordinates = "geolon_t geolat_t" ;<br>
> ><br>
> > pot_rho_0:standard_name = "sea_water_potential_density" ;<br>
> ><br>
> > float salt_on_rho(Time, potrho, grid_yt_ocean,<br>
> grid_xt_ocean) ;<br>
> ><br>
> > salt_on_rho:long_name = "salt on potential density surface" ;<br>
> ><br>
> > salt_on_rho:units = "psu" ;<br>
> ><br>
> > salt_on_rho:valid_range = 0.f, 1000.f ;<br>
> ><br>
> > salt_on_rho:missing_value = -1.e+20f ;<br>
> ><br>
> > salt_on_rho:_FillValue = -1.e+20f ;<br>
> ><br>
> > salt_on_rho:cell_methods = "time: mean" ;<br>
> ><br>
> > salt_on_rho:time_avg_info =<br>
> "average_T1,average_T2,average_DT" ;<br>
> ><br>
> > float potrho_mix_depth(Time, yt_ocean, xt_ocean) ;<br>
> ><br>
> > potrho_mix_depth:long_name = "Depth of potential density<br>
> mixed<br>
> > layer" ;<br>
> ><br>
> > potrho_mix_depth:units = "m" ;<br>
> ><br>
> > potrho_mix_depth:valid_range = -1000000.f, 1000000.f ;<br>
> ><br>
> > potrho_mix_depth:missing_value = -1.e+20f ;<br>
> ><br>
> > potrho_mix_depth:_FillValue = -1.e+20f ;<br>
> ><br>
> > potrho_mix_depth:cell_methods = "time: mean" ;<br>
> ><br>
> > potrho_mix_depth:time_avg_info =<br>
> > "average_T1,average_T2,average_DT" ;<br>
> ><br>
> > potrho_mix_depth:coordinates = "geolon_t geolat_t" ;<br>
> ><br>
> > float potrho_mix_base(Time, yt_ocean, xt_ocean) ;<br>
> ><br>
> > potrho_mix_base:long_name = "Potential density at mixed layer<br>
> > base" ;<br>
> ><br>
> > potrho_mix_base:units = "kg/m^3" ;<br>
> ><br>
> > potrho_mix_base:valid_range = -1000000.f, 1000000.f ;<br>
> ><br>
> > potrho_mix_base:missing_value = -1.e+20f ;<br>
> ><br>
> > potrho_mix_base:_FillValue = -1.e+20f ;<br>
> ><br>
> > potrho_mix_base:cell_methods = "time: mean" ;<br>
> ><br>
> > potrho_mix_base:time_avg_info =<br>
> "average_T1,average_T2,average_DT" ;<br>
> ><br>
> > potrho_mix_base:coordinates = "geolon_t geolat_t" ;<br>
> ><br>
> > float press(Time, st_ocean, yt_ocean, xt_ocean) ;<br>
> ><br>
> > press:long_name = "absolute pressure" ;<br>
> ><br>
> > press:units = "dbar" ;<br>
> ><br>
> > press:valid_range = -10.f, 1000000.f ;<br>
> ><br>
> > press:missing_value = -1.e+20f ;<br>
> ><br>
> > press:_FillValue = -1.e+20f ;<br>
> ><br>
> > press:cell_methods = "time: mean" ;<br>
> ><br>
> > press:time_avg_info = "average_T1,average_T2,average_DT" ;<br>
> ><br>
> > press:coordinates = "geolon_t geolat_t" ;<br>
> ><br>
> > float mld(Time, yt_ocean, xt_ocean) ;<br>
> ><br>
> > mld:long_name = "mixed layer depth determined by density<br>
> criteria" ;<br>
> ><br>
> > mld:units = "m" ;<br>
> ><br>
> > mld:valid_range = 0.f, 1000000.f ;<br>
> ><br>
> > mld:missing_value = -1.e+20f ;<br>
> ><br>
> > mld:_FillValue = -1.e+20f ;<br>
> ><br>
> > mld:cell_methods = "time: mean" ;<br>
> ><br>
> > mld:time_avg_info = "average_T1,average_T2,average_DT" ;<br>
> ><br>
> > mld:coordinates = "geolon_t geolat_t" ;<br>
> ><br>
> > mld:standard_name =<br>
> > "ocean_mixed_layer_thickness_defined_by_sigma_t" ;<br>
> ><br>
> > double average_T1(Time) ;<br>
> ><br>
> > average_T1:long_name = "Start time for average period" ;<br>
> ><br>
> > average_T1:units = "days since 1980-01-01 00:00:00" ;<br>
> ><br>
> > average_T1:missing_value = 1.e+20 ;<br>
> ><br>
> > average_T1:_FillValue = 1.e+20 ;<br>
> ><br>
> > double average_T2(Time) ;<br>
> ><br>
> > average_T2:long_name = "End time for average period" ;<br>
> ><br>
> > average_T2:units = "days since 1980-01-01 00:00:00" ;<br>
> ><br>
> > average_T2:missing_value = 1.e+20 ;<br>
> ><br>
> > average_T2:_FillValue = 1.e+20 ;<br>
> ><br>
> > double average_DT(Time) ;<br>
> ><br>
> > average_DT:long_name = "Length of average period" ;<br>
> ><br>
> > average_DT:units = "days" ;<br>
> ><br>
> > average_DT:missing_value = 1.e+20 ;<br>
> ><br>
> > average_DT:_FillValue = 1.e+20 ;<br>
> ><br>
> > double Time_bounds(Time, nv) ;<br>
> ><br>
> > Time_bounds:long_name = "Time axis boundaries" ;<br>
> ><br>
> > Time_bounds:units = "days" ;<br>
> ><br>
> > Time_bounds:missing_value = 1.e+20 ;<br>
> ><br>
> > Time_bounds:_FillValue = 1.e+20 ;<br>
> ><br>
> ><br>
> > // global attributes:<br>
> ><br>
> > :filename = "<a href="http://ocean_layer.nc" target="_blank">ocean_layer.nc</a> <<a href="http://ocean_layer.nc" target="_blank">http://ocean_layer.nc</a>><br>
> <<a href="http://ocean_layer.nc" target="_blank">http://ocean_layer.nc</a>>" ;<br>
> ><br>
> > :title = "MOM4 box" ;<br>
> ><br>
> > :grid_type = "regular" ;<br>
> ><br>
> > :grid_tile = "N/A" ;<br>
> ><br>
> > data:<br>
> ><br>
> ><br>
> > grid_xt_ocean = 1.25, 3.75, 6.25, 8.75, 11.25, 13.75,<br>
> 16.25,<br>
> > 18.75, 21.25,<br>
> ><br>
> > 23.75, 26.25, 28.75, 31.25, 33.75, 36.25, 38.75, 41.25,<br>
> > 43.75, 46.25,<br>
> ><br>
> > 48.75, 51.25, 53.75, 56.25, 58.75 ;<br>
> ><br>
> ><br>
> > grid_yt_ocean = 1.25, 3.7488100528717, 6.2393012046814,<br>
> > 8.71919345855713,<br>
> ><br>
> > 11.1815176010132, 13.6242818832397, 16.0408229827881,<br>
> > 18.4296073913574,<br>
> ><br>
> > 20.784387588501, 23.1042194366455, 25.3833503723145,<br>
> > 27.6215190887451,<br>
> ><br>
> > 29.813497543335, 31.9597606658936, 34.0555992126465,<br>
> > 36.1022338867188,<br>
> ><br>
> > 38.0954322814941, 40.0371551513672, 41.9235687255859,<br>
> > 43.7573394775391,<br>
> ><br>
> > 45.5349464416504, 47.259708404541, 48.9283294677734,<br>
> > 50.5447235107422,<br>
> ><br>
> > 52.1057090759277, 53.6157531738281, 55.0716934204102,<br>
> > 56.4785079956055,<br>
> ><br>
> > 57.8329429626465, 59.1404571533203, 60.3976211547852,<br>
> > 61.6103477478027,<br>
> ><br>
> > 62.7749481201172, 63.8977813720703, 64.9748153686523 ;<br>
> ><br>
> ><br>
> > theta = -1.84, -1.52, -1.2, -0.88, -0.56, -0.24,<br>
> > 0.0800000000000003, 0.4,<br>
> ><br>
> > 0.72, 1.04, 1.36, 1.68, 2, 2.32, 2.64, 2.96, 3.28, 3.6,<br>
> > 3.92, 4.24, 4.56,<br>
> ><br>
> > 4.88, 5.2, 5.52, 5.84, 6.16, 6.48, 6.8, 7.12, 7.44,<br>
> 7.76,<br>
> > 8.08, 8.4,<br>
> ><br>
> > 8.72, 9.04, 9.36, 9.68000000000001, 10, 10.32, 10.64,<br>
> > 10.96, 11.28, 11.6,<br>
> ><br>
> > 11.92, 12.24, 12.56, 12.88, 13.2, 13.52, 13.84, 14.16,<br>
> > 14.48, 14.8,<br>
> ><br>
> > 15.12, 15.44, 15.76, 16.08, 16.4, 16.72, 17.04, 17.36,<br>
> > 17.68, 18, 18.32,<br>
> ><br>
> > 18.64, 18.96, 19.28, 19.6, 19.92, 20.24, 20.56, 20.88,<br>
> > 21.2, 21.52,<br>
> ><br>
> > 21.84, 22.16, 22.48, 22.8, 23.12, 23.44, 23.76, 24.08,<br>
> > 24.4, 24.72,<br>
> ><br>
> > 25.04, 25.36, 25.68, 26, 26.32, 26.64, 26.96, 27.28,<br>
> 27.6,<br>
> > 27.92, 28.24,<br>
> ><br>
> > 28.56, 28.88, 29.2, 29.52, 29.84 ;<br>
> ><br>
> ><br>
> > theta_edges = 1028, 1028.1, 1028.2, 1028.3, 1028.4, 1028.5,<br>
> > 1028.6, 1028.7,<br>
> ><br>
> > 1028.8, 1028.9, 1029, 1029.1, 1029.2, 1029.3, 1029.4,<br>
> > 1029.5, 1029.6,<br>
> ><br>
> > 1029.7, 1029.8, 1029.9, 1030, 1030.1, 1030.2, 1030.3,<br>
> > 1030.4, 1030.5,<br>
> ><br>
> > 1030.6, 1030.7, 1030.8, 1030.9, 1031, 1031.1, 1031.2,<br>
> > 1031.3, 1031.4,<br>
> ><br>
> > 1031.5, 1031.6, 1031.7, 1031.8, 1031.9, 1032, 1032.1,<br>
> > 1032.2, 1032.3,<br>
> ><br>
> > 1032.4, 1032.5, 1032.6, 1032.7, 1032.8, 1032.9, 1033,<br>
> > 1033.1, 1033.2,<br>
> ><br>
> > 1033.3, 1033.4, 1033.49999999999, 1033.59999999999,<br>
> > 1033.69999999999,<br>
> ><br>
> > 1033.79999999999, 1033.89999999999, 1033.99999999999,<br>
> > 1034.09999999999,<br>
> ><br>
> > 1034.19999999999, 1034.29999999999, 1034.39999999999,<br>
> > 1034.49999999999,<br>
> ><br>
> > 1034.59999999999, 1034.69999999999, 1034.79999999999,<br>
> > 1034.89999999999,<br>
> ><br>
> > 1034.99999999999, 1035.09999999999, 1035.19999999999,<br>
> > 1035.29999999999,<br>
> ><br>
> > 1035.39999999999, 1035.49999999999, 1035.59999999999,<br>
> > 1035.69999999999,<br>
> ><br>
> > 1035.79999999999, 1035.89999999999, 1035.99999999999,<br>
> > 1036.09999999999,<br>
> ><br>
> > 1036.19999999999, 1036.29999999999, 1036.39999999999,<br>
> > 1036.49999999999,<br>
> ><br>
> > 1036.59999999999, 1036.69999999999, 1036.79999999999,<br>
> > 1036.89999999999,<br>
> ><br>
> > 1036.99999999999, 1037.09999999999, 1037.19999999999,<br>
> > 1037.29999999999,<br>
> ><br>
> > 1037.39999999999, 1037.49999999999, 1037.59999999999,<br>
> > 1037.69999999999,<br>
> ><br>
> > 1037.79999999999, 1037.89999999999, 1037.99999999999 ;<br>
> ><br>
> ><br>
> > Time = 2 ;<br>
> ><br>
> ><br>
> > nv = 1, 2 ;<br>
> ><br>
> ><br>
> > xt_ocean = 1.25, 3.75, 6.25, 8.75, 11.25, 13.75, 16.25,<br>
> 18.75,<br>
> > 21.25, 23.75,<br>
> ><br>
> > 26.25, 28.75, 31.25, 33.75, 36.25, 38.75, 41.25, 43.75,<br>
> > 46.25, 48.75,<br>
> ><br>
> > 51.25, 53.75, 56.25, 58.75 ;<br>
> ><br>
> ><br>
> > yt_ocean = 1.25, 3.7488100528717, 6.2393012046814,<br>
> > 8.71919345855713,<br>
> ><br>
> > 11.1815176010132, 13.6242818832397, 16.0408229827881,<br>
> > 18.4296073913574,<br>
> ><br>
> > 20.784387588501, 23.1042194366455, 25.3833503723145,<br>
> > 27.6215190887451,<br>
> ><br>
> > 29.813497543335, 31.9597606658936, 34.0555992126465,<br>
> > 36.1022338867188,<br>
> ><br>
> > 38.0954322814941, 40.0371551513672, 41.9235687255859,<br>
> > 43.7573394775391,<br>
> ><br>
> > 45.5349464416504, 47.259708404541, 48.9283294677734,<br>
> > 50.5447235107422,<br>
> ><br>
> > 52.1057090759277, 53.6157531738281, 55.0716934204102,<br>
> > 56.4785079956055,<br>
> ><br>
> > 57.8329429626465, 59.1404571533203, 60.3976211547852,<br>
> > 61.6103477478027,<br>
> ><br>
> > 62.7749481201172, 63.8977813720703, 64.9748153686523 ;<br>
> ><br>
> ><br>
> > st_ocean = 17.75025, 53.25075, 88.75125, 124.25175,<br>
> 159.75225,<br>
> > 195.25275,<br>
> ><br>
> > 230.75325, 266.25375, 314.34159289856, 399.171716592407,<br>
> ><br>
> > 541.922400494385, 759.079548358154, 1061.10117883301,<br>
> > 1451.56990550537,<br>
> ><br>
> > 1926.90292906494, 2476.64250944824, 3084.30280480957,<br>
> > 3728.70522619629 ;<br>
> ><br>
> ><br>
> > st_edges_ocean = 0, 35.5005, 71.001, 106.5015, 142.002,<br>
> > 177.5025, 213.003,<br>
> ><br>
> > 248.5035, 290.29765597229, 356.756654745483,<br>
> 470.547043066406,<br>
> ><br>
> > 650.50097442627, 910.090332641602, 1256.33548026123,<br>
> > 1689.23641728516,<br>
> ><br>
> > 2201.77290498047, 2780.47265712891, 3406.50401550293,<br>
> > 4057.20024763184 ;<br>
> ><br>
> ><br>
> > potrho = 1028.05, 1028.15, 1028.25, 1028.35, 1028.45,<br>
> 1028.55,<br>
> > 1028.65,<br>
> ><br>
> > 1028.75, 1028.85, 1028.95, 1029.05, 1029.15, 1029.25,<br>
> > 1029.35, 1029.45,<br>
> ><br>
> > 1029.55, 1029.65, 1029.75, 1029.85, 1029.95, 1030.05,<br>
> > 1030.15, 1030.25,<br>
> ><br>
> > 1030.35, 1030.45, 1030.55, 1030.65, 1030.75, 1030.85,<br>
> > 1030.95, 1031.05,<br>
> ><br>
> > 1031.15, 1031.25, 1031.35, 1031.45, 1031.55, 1031.65,<br>
> > 1031.75, 1031.85,<br>
> ><br>
> > 1031.95, 1032.05, 1032.15, 1032.25, 1032.35, 1032.45,<br>
> > 1032.55, 1032.65,<br>
> ><br>
> > 1032.75, 1032.85, 1032.95, 1033.05, 1033.15, 1033.25,<br>
> > 1033.35, 1033.45,<br>
> ><br>
> > 1033.54999999999, 1033.64999999999, 1033.74999999999,<br>
> > 1033.84999999999,<br>
> ><br>
> > 1033.94999999999, 1034.04999999999, 1034.14999999999,<br>
> > 1034.24999999999,<br>
> ><br>
> > 1034.34999999999, 1034.44999999999, 1034.54999999999,<br>
> > 1034.64999999999,<br>
> ><br>
> > 1034.74999999999, 1034.84999999999, 1034.94999999999,<br>
> > 1035.04999999999,<br>
> ><br>
> > 1035.14999999999, 1035.24999999999, 1035.34999999999,<br>
> > 1035.44999999999,<br>
> ><br>
> > 1035.54999999999, 1035.64999999999, 1035.74999999999,<br>
> > 1035.84999999999,<br>
> ><br>
> > 1035.94999999999, 1036.04999999999, 1036.14999999999,<br>
> > 1036.24999999999,<br>
> ><br>
> > 1036.34999999999, 1036.44999999999, 1036.54999999999,<br>
> > 1036.64999999999,<br>
> ><br>
> > 1036.74999999999, 1036.84999999999, 1036.94999999999,<br>
> > 1037.04999999999,<br>
> ><br>
> > 1037.14999999999, 1037.24999999999, 1037.34999999999,<br>
> > 1037.44999999999,<br>
> ><br>
> > 1037.54999999999, 1037.64999999999, 1037.74999999999,<br>
> > 1037.84999999999,<br>
> ><br>
> > 1037.94999999999 ;<br>
> ><br>
> ><br>
> > potrho_edges = 1028, 1028.1, 1028.2, 1028.3, 1028.4,<br>
> 1028.5,<br>
> > 1028.6, 1028.7,<br>
> ><br>
> > 1028.8, 1028.9, 1029, 1029.1, 1029.2, 1029.3, 1029.4,<br>
> > 1029.5, 1029.6,<br>
> ><br>
> > 1029.7, 1029.8, 1029.9, 1030, 1030.1, 1030.2, 1030.3,<br>
> > 1030.4, 1030.5,<br>
> ><br>
> > 1030.6, 1030.7, 1030.8, 1030.9, 1031, 1031.1, 1031.2,<br>
> > 1031.3, 1031.4,<br>
> ><br>
> > 1031.5, 1031.6, 1031.7, 1031.8, 1031.9, 1032, 1032.1,<br>
> > 1032.2, 1032.3,<br>
> ><br>
> > 1032.4, 1032.5, 1032.6, 1032.7, 1032.8, 1032.9, 1033,<br>
> > 1033.1, 1033.2,<br>
> ><br>
> > 1033.3, 1033.4, 1033.49999999999, 1033.59999999999,<br>
> > 1033.69999999999,<br>
> ><br>
> > 1033.79999999999, 1033.89999999999, 1033.99999999999,<br>
> > 1034.09999999999,<br>
> ><br>
> > 1034.19999999999, 1034.29999999999, 1034.39999999999,<br>
> > 1034.49999999999,<br>
> ><br>
> > 1034.59999999999, 1034.69999999999, 1034.79999999999,<br>
> > 1034.89999999999,<br>
> ><br>
> > 1034.99999999999, 1035.09999999999, 1035.19999999999,<br>
> > 1035.29999999999,<br>
> ><br>
> > 1035.39999999999, 1035.49999999999, 1035.59999999999,<br>
> > 1035.69999999999,<br>
> ><br>
> > 1035.79999999999, 1035.89999999999, 1035.99999999999,<br>
> > 1036.09999999999,<br>
> ><br>
> > 1036.19999999999, 1036.29999999999, 1036.39999999999,<br>
> > 1036.49999999999,<br>
> ><br>
> > 1036.59999999999, 1036.69999999999, 1036.79999999999,<br>
> > 1036.89999999999,<br>
> ><br>
> > 1036.99999999999, 1037.09999999999, 1037.19999999999,<br>
> > 1037.29999999999,<br>
> ><br>
> > 1037.39999999999, 1037.49999999999, 1037.59999999999,<br>
> > 1037.69999999999,<br>
> ><br>
> > 1037.79999999999, 1037.89999999999, 1037.99999999999 ;<br>
> ><br>
> > }<br>
> ><br>
> ><br>
> > I am having trouble deciding which will be the xdef,<br>
> ydef, zdef,<br>
> > and tdef values. Also, how do I find what UNDEF number to<br>
> use.<br>
> > Your help will be greatly appreciated.<br>
> ><br>
> > Thanks,<br>
> ><br>
> > David<br>
> ><br>
> ><br>
> ><br>
> ><br>
> ><br>
> > _______________________________________________<br>
> > gradsusr mailing list<br>
</div></div>> > <a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a> <mailto:<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>><br>
<div class="">> > <a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
> ><br>
><br>
> --<br>
><br>
> Please note that <a href="mailto:Charles.Seman@noaa.gov">Charles.Seman@noaa.gov</a><br>
</div>> <mailto:<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> <mailto:<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> <mailto:<a href="mailto:Charles.Seman@noaa.gov">Charles.Seman@noaa.gov</a>><br>
<div class="">> 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>
> "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."<br>
> _______________________________________________<br>
> gradsusr mailing list<br>
</div>> <a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a> <mailto:<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>
<div class="HOEnZb"><div class="h5">><br>
><br>
><br>
><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>
<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>
"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."<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>