<div dir="ltr"><div>Dear Grads users:</div><div><br></div><div>I am attempting to plot a time series of sea surface temperatures using a .nc file.  The ocean domain of the model output is in MOM4pi format with 200lat*360lon data points.  GrADS won&#39;t open the file using sdfopen, so I tried xdfopen.  </div>
<div><br></div><div>I have copied the ncdump as well as my .ctl file for this model.  Basically, what I am asking is how to arrange the XDEF and YDEF in the .ctl file so that I can read the file, and plot the data I need over the appropriate geographical domain.  If you look at the bolded part of the .ctl file, you&#39;ll see I and j where lon and lat would more typically appear.  </div>
<div><br></div><div>How can I fix my .ctl file to correctly read the data?</div><div><br></div><div>Thanks,</div><div>Corey </div><div><br></div><div>DSET /stu_home/corey/tos_Omon_BNU-ESM_historical_r1i1p1_185001-200512.nc<br>
DTYPE netcdf<br>TITLE SST (C) Nino3.4<br>UNDEF 1.e20 <br><strong>XDEF i 360 linear 60 1<br>YDEF j 180 linear -90 1</strong> <br>TDEF time 1872 linear 01jan1850 1mo<br>VARS 1 <br>tos=&gt;sst 0 t,y,x Sea Surface Temperature<br>
ENDVARS</div><div><br></div><div><br></div><div>ncdump -c tos_Omon_BNU-ESM_historical_r1i1p1_185001-200512<br>.nc<br>netcdf tos_Omon_BNU-ESM_historical_r1i1p1_185001-200512 {<br>dimensions:<br>        time = UNLIMITED ; // (1872 currently)<br>
        j = 200 ;<br>        i = 360 ;<br>        bnds = 2 ;<br>        vertices = 4 ;<br>variables:<br>        double time(time) ;<br>                time:bounds = &quot;time_bnds&quot; ;<br>                time:units = &quot;days since 1850-01-01 00:00:00&quot; ;<br>
                time:calendar = &quot;noleap&quot; ;<br>                time:axis = &quot;T&quot; ;<br>                time:long_name = &quot;time&quot; ;<br>                time:standard_name = &quot;time&quot; ;<br>        double time_bnds(time, bnds) ;<br>
        int j(j) ;<br>                j:units = &quot;1&quot; ;<br>                j:long_name = &quot;cell index along second dimension&quot; ;<br>        int i(i) ;<br>                i:units = &quot;1&quot; ;<br>                i:long_name = &quot;cell index along first dimension&quot; ;<br>
        float lat(j, i) ;<br>                lat:standard_name = &quot;latitude&quot; ;<br>                lat:long_name = &quot;latitude coordinate&quot; ;<br>                lat:units = &quot;degrees_north&quot; ;<br>                lat:bounds = &quot;lat_vertices&quot; ;<br>
        float lon(j, i) ;<br>                lon:standard_name = &quot;longitude&quot; ;<br>                lon:long_name = &quot;longitude coordinate&quot; ;<br>                lon:units = &quot;degrees_east&quot; ;<br>                lon:bounds = &quot;lon_vertices&quot; ;<br>
        float lat_vertices(j, i, vertices) ;<br>                lat_vertices:units = &quot;degrees_north&quot; ;<br>        float lon_vertices(j, i, vertices) ;<br>                lon_vertices:units = &quot;degrees_east&quot; ;<br>
        float tos(time, j, i) ;<br>                tos:standard_name = &quot;sea_surface_temperature&quot; ;<br>                tos:long_name = &quot;Sea Surface Temperature&quot; ;<br>                tos:comment = &quot;\&quot;this may differ from \&quot;\&quot;surface temperature\&quot;\&quot; in regions of sea ice.\&quot;&quot; ;<br>
                tos:units = &quot;K&quot; ;<br>                tos:cell_methods = &quot;time: mean&quot; ;<br>                tos:cell_measures = &quot;area: areacello&quot; ;<br>                tos:history = &quot;2013-03-12T09:10:58Z altered by CMOR: replaced missing value flag (1e+30) with standard missing value (1e+20).&quot; ;<br>
                tos:missing_value = 1.e+20f ;<br>                tos:_FillValue = 1.e+20f ;<br>                tos:associated_files = &quot;baseURL: <a href="http://cmip-pcmdi.llnl.gov/CMIP5/dataLocation">http://cmip-pcmdi.llnl.gov/CMIP5/dataLocation</a> gridspecFile: gridspec_ocean_fx_BNU-ESM_historical_r0i0p0.nc areacello: areacello_fx_BNU-ESM_historical_r0i0p0.nc&quot; ;<br>
                tos:coordinates = &quot;lat lon&quot; ;</div><div>// global attributes:<br>                :institution = &quot;GCESS,BNU,Beijing,China&quot; ;<br>                :institute_id = &quot;BNU&quot; ;<br>                :experiment_id = &quot;historical&quot; ;<br>
                :source = &quot;BNU-ESM&quot; ;<br>                :model_id = &quot;BNU-ESM&quot; ;<br>                :forcing = &quot;Nat,Ant&quot; ;<br>                :parent_experiment_id = &quot;piControl&quot; ;<br>
                :parent_experiment_rip = &quot;r1i1p1&quot; ;<br>                :branch_time = 0. ;<br>                :contact = &quot;Ji Duoying (<a href="mailto:duoyingji@bnu.edu.cn">duoyingji@bnu.edu.cn</a>)&quot; ;<br>
                :initialization_method = 1 ;<br>                :physics_version = 1 ;<br>                :tracking_id = &quot;748b103b-2782-4c80-bdb6-e7c6d59449d9&quot; ;<br>                :product = &quot;output&quot; ;<br>
                :experiment = &quot;historical&quot; ;<br>                :frequency = &quot;mon&quot; ;<br>                :creation_date = &quot;2013-03-12T09:10:58Z&quot; ;<br>                :history = &quot;2013-03-12T09:10:58Z CMOR rewrote data to comply with CF standards and CMIP5 requirements.&quot; ;<br>
                :Conventions = &quot;CF-1.4&quot; ;<br>                :project_id = &quot;CMIP5&quot; ;<br>                :table_id = &quot;Table Omon (26 July 2011) 7a5904e887d79e2db795ad92d85c7968&quot; ;<br>                :title = &quot;BNU-ESM model output prepared for CMIP5 historical&quot; ;<br>
                :parent_experiment = &quot;pre-industrial control&quot; ;<br>                :modeling_realm = &quot;ocean&quot; ;<br>                :realization = 1 ;<br>                :cmor_version = &quot;2.7.1&quot; ;<br>
</div></div>