[gradsusr] MOM4pi

Corey Gabriel cjgabriel7 at gmail.com
Thu Mar 20 10:27:02 EDT 2014


Dear Grads users:

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't open the file using sdfopen, so I
tried xdfopen.

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'll see I and j where lon and lat would more typically
appear.

How can I fix my .ctl file to correctly read the data?

Thanks,
Corey

DSET /stu_home/corey/tos_Omon_BNU-ESM_historical_r1i1p1_185001-200512.nc
DTYPE netcdf
TITLE SST (C) Nino3.4
UNDEF 1.e20

*XDEF i 360 linear 60 1YDEF j 180 linear -90 1*
TDEF time 1872 linear 01jan1850 1mo
VARS 1
tos=>sst 0 t,y,x Sea Surface Temperature
ENDVARS


ncdump -c tos_Omon_BNU-ESM_historical_r1i1p1_185001-200512
.nc
netcdf tos_Omon_BNU-ESM_historical_r1i1p1_185001-200512 {
dimensions:
        time = UNLIMITED ; // (1872 currently)
        j = 200 ;
        i = 360 ;
        bnds = 2 ;
        vertices = 4 ;
variables:
        double time(time) ;
                time:bounds = "time_bnds" ;
                time:units = "days since 1850-01-01 00:00:00" ;
                time:calendar = "noleap" ;
                time:axis = "T" ;
                time:long_name = "time" ;
                time:standard_name = "time" ;
        double time_bnds(time, bnds) ;
        int j(j) ;
                j:units = "1" ;
                j:long_name = "cell index along second dimension" ;
        int i(i) ;
                i:units = "1" ;
                i:long_name = "cell index along first dimension" ;
        float lat(j, i) ;
                lat:standard_name = "latitude" ;
                lat:long_name = "latitude coordinate" ;
                lat:units = "degrees_north" ;
                lat:bounds = "lat_vertices" ;
        float lon(j, i) ;
                lon:standard_name = "longitude" ;
                lon:long_name = "longitude coordinate" ;
                lon:units = "degrees_east" ;
                lon:bounds = "lon_vertices" ;
        float lat_vertices(j, i, vertices) ;
                lat_vertices:units = "degrees_north" ;
        float lon_vertices(j, i, vertices) ;
                lon_vertices:units = "degrees_east" ;
        float tos(time, j, i) ;
                tos:standard_name = "sea_surface_temperature" ;
                tos:long_name = "Sea Surface Temperature" ;
                tos:comment = "\"this may differ from \"\"surface
temperature\"\" in regions of sea ice.\"" ;
                tos:units = "K" ;
                tos:cell_methods = "time: mean" ;
                tos:cell_measures = "area: areacello" ;
                tos:history = "2013-03-12T09:10:58Z altered by CMOR:
replaced missing value flag (1e+30) with standard missing value (1e+20)." ;
                tos:missing_value = 1.e+20f ;
                tos:_FillValue = 1.e+20f ;
                tos:associated_files = "baseURL:
http://cmip-pcmdi.llnl.gov/CMIP5/dataLocation gridspecFile:
gridspec_ocean_fx_BNU-ESM_historical_r0i0p0.nc areacello:
areacello_fx_BNU-ESM_historical_r0i0p0.nc" ;
                tos:coordinates = "lat lon" ;
// global attributes:
                :institution = "GCESS,BNU,Beijing,China" ;
                :institute_id = "BNU" ;
                :experiment_id = "historical" ;
                :source = "BNU-ESM" ;
                :model_id = "BNU-ESM" ;
                :forcing = "Nat,Ant" ;
                :parent_experiment_id = "piControl" ;
                :parent_experiment_rip = "r1i1p1" ;
                :branch_time = 0. ;
                :contact = "Ji Duoying (duoyingji at bnu.edu.cn)" ;
                :initialization_method = 1 ;
                :physics_version = 1 ;
                :tracking_id = "748b103b-2782-4c80-bdb6-e7c6d59449d9" ;
                :product = "output" ;
                :experiment = "historical" ;
                :frequency = "mon" ;
                :creation_date = "2013-03-12T09:10:58Z" ;
                :history = "2013-03-12T09:10:58Z CMOR rewrote data to
comply with CF standards and CMIP5 requirements." ;
                :Conventions = "CF-1.4" ;
                :project_id = "CMIP5" ;
                :table_id = "Table Omon (26 July 2011)
7a5904e887d79e2db795ad92d85c7968" ;
                :title = "BNU-ESM model output prepared for CMIP5
historical" ;
                :parent_experiment = "pre-industrial control" ;
                :modeling_realm = "ocean" ;
                :realization = 1 ;
                :cmor_version = "2.7.1" ;
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20140320/7d4d49d0/attachment.html 


More information about the gradsusr mailing list