[gradsusr] MOM4pi

[Jennifer Adams]

You will need a PDEF file to interpolate this data to a regular lat/lon grid. I have  posted my technique and scripts for doing this in ftp://cola.gmu.edu/pub/jma/cmip5. --Jennifer


On Mar 20, 2014, at 10:27 AM, Corey Gabriel wrote:

> 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 1
> YDEF 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" ;
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr

--
Jennifer M. Adams
Center for Ocean-Land-Atmosphere Studies (COLA)
111 Research Hall, Mail Stop 2B3
George Mason University
4400 University Drive
Fairfax, VA 22030 





-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20140325/d23a35ef/attachment.html