[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