[gradsusr] Problem projecting CMIP5 SST
Roberto Mera
RMera at ucsusa.org
Fri Feb 14 16:29:03 EST 2014
Thanks!
________________________________
From: gradsusr-bounces at gradsusr.org [gradsusr-bounces at gradsusr.org] on behalf of Jennifer Adams [jma at cola.iges.org]
Sent: Friday, February 14, 2014 2:41 PM
To: GrADS Users Forum
Subject: Re: [gradsusr] Problem projecting CMIP5 SST
The critical three words here are "rotated pole grid". You're going to need a PDEF file for this. I would write to the contact person (Dr. Tongwen Wu) and ask if they are using GrADS to display this data -- they may already have a PDEF file for you. If not, then you can use the ESMF Tool RegridWeightGen to generate the interpolation weights. Please see ftp://cola.gmu.edu/pub/jma/cmip5/README.
--Jennifer
On Feb 12, 2014, at 3:09 PM, Roberto Mera wrote:
Jennifer,
Basically, what I meant about “falls apart” is that nothing north of -29.5 displays. That’s where latitudes cease being linear.
I think this is what is relevant from the ncdump.
netcdf tos_Omon_bcc-csm1-1_historicalGHG_r1i1p1_185001-201212 {
dimensions:
time = UNLIMITED ; // (1956 currently)
rlat = 232 ;
rlon = 360 ;
bnds = 2 ;
variables:
double time(time) ;
time:bounds = "time_bnds" ;
time:units = "days since 1850-01-01" ;
time:calendar = "noleap" ;
time:axis = "T" ;
time:long_name = "time" ;
time:standard_name = "time" ;
double time_bnds(time, bnds) ;
double rlat(rlat) ;
rlat:bounds = "rlat_bnds" ;
rlat:units = "degrees_north" ;
rlat:axis = "Y" ;
rlat:long_name = "latitude in rotated pole grid" ;
rlat:standard_name = "grid_latitude" ;
double rlat_bnds(rlat, bnds) ;
double rlon(rlon) ;
rlon:bounds = "rlon_bnds" ;
rlon:units = "degrees_east" ;
rlon:axis = "X" ;
rlon:long_name = "longitude in rotated pole grid" ;
rlon:standard_name = "grid_longitude" ;
double rlon_bnds(rlon, bnds) ;
float lat(rlat, rlon) ;
lat:standard_name = "latitude" ;
lat:long_name = "latitude coordinate" ;
lat:units = "degrees_north" ;
float lon(rlat, rlon) ;
lon:standard_name = "longitude" ;
lon:long_name = "longitude coordinate" ;
lon:units = "degrees_east" ;
float tos(time, rlat, rlon) ;
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:original_name = "SST" ;
tos:cell_methods = "time: mean" ;
tos:cell_measures = "area: areacello" ;
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_bcc-csm1-1_historicalGHG_r0i0p0.nc areacello: areacello_fx_bcc-csm1-1_historicalGHG_r0i0p0.nc" ;
tos:coordinates = "lat lon" ;
// global attributes:
:institution = "Beijing Climate Center(BCC),China Meteorological Administration,China" ;
:institute_id = "BCC" ;
:experiment_id = "historicalGHG" ;
:source = "bcc-csm1-1:atmosphere: BCC_AGCM2.1 (T42L26); land: BCC_AVIM1.0;ocean: MOM4_L40 (tripolar, 1 lon x (1-1/3) lat, L40);sea ice: SIS (tripolar,1 lon x (1-1/3) lat)" ;
:model_id = "bcc-csm1-1" ;
:forcing = "GHG" ;
:parent_experiment_id = "piControl" ;
:parent_experiment_rip = "r1i1p1" ;
:branch_time = 450. ;
:contact = "Dr. Tongwen Wu (twwu at cma.gov.cn<mailto:twwu at cma.gov.cn>)" ;
:history = "Output from monthly mean data 2011-10-20T06:00:30Z CMOR rewrote data to comply with CF standards and CMIP5 requirements." ;
:comment = "Historical simulation but with greenhouse gas forcing only, The model initial state is from the 450-year of pre-industrial control simulation. The directions of variables uo and vo are transfered from the native tripolar grid directions to geographical eastward and northward, respectively." ;
From: gradsusr-bounces at gradsusr.org<mailto:gradsusr-bounces at gradsusr.org> [mailto:gradsusr-bounces at gradsusr.org] On Behalf Of Jennifer Adams
Sent: Tuesday, February 11, 2014 7:16 PM
To: GrADS Users Forum
Subject: Re: [gradsusr] Problem projecting CMIP5 SST
What exactly do you mean by "the display falls apart"? What is the ncdump output? What does the rest of your descriptor file look like? What version of GrADS are you using? If the grid coordinates (lat and lon) are 1-dimensional, then you should be able to open this with GrADS without having to create a PDEF file.
---Jennifer
On Feb 11, 2014, at 2:37 PM, Andrew Friedman wrote:
Hi Roberto,
Jennifer Adams has written some scripts which I believe can process such CMIP5 ocean data. See this posting from the listserve:
http://gradsusr.org/pipermail/gradsusr/2013-January/017089.html
Andrew
On Feb 11, 2014, at 7:54 AM, Roberto Mera <RMera at ucsusa.org<mailto:RMera at ucsusa.org>> wrote:
All,
The CMIP5 natural GHG run for BCC model (tos_Omon_bcc-csm1-1_historicalGHG_r1i1p1_185001-201212.nc) has very strange coordinates with 232 ydef (latitude) but a linear 1deg longitude. I have tried defining ydef with levels but the display falls apart after the interval in the latitude stops being every 1 degree. Here’s an ncdump of the latitude data:
ydef 232 levels -81.5 -80.5 -79.5 -78.5 -77.5 -76.5 -75.5 -74.5 -73.5 -72.5 -71.5 -70.5 -69.5 -68.5 -67.5 -66.5 -65.5 -64.5 -63.5 -62.5 -61.5 -60.5 -59.5 -58.5 -57.5 -56.5 -55.5 -54.5 -53.5 -52.5 -51.5 -50.5 -49.5 -48.5 -47.5 -46.5 -45.5 -44.5 -43.5 -42.5 -41.5 -40.5 -39.5 -38.5 -37.5 -36.5 -35.5 -34.5 -33.5 -32.5 -31.5 -30.5 -29.5 -28.500913619995 -27.50820732117 -26.523630142212 -25.552488327026 -24.596292495726 -23.660039901733 -22.744861602783 -21.855302810669 -20.991981506348 -20.158876419067 -19.355987548828 -18.586631774902 -17.850109100342 -17.149007797241 -16.481884002686 -15.850565910339 -15.252840995789 -14.689778327942 -14.158419609070 -13.659105300903 -13.188179016113 -12.745319366455 -12.326250076294 -11.930083274841 -11.552030563355 -11.190753936768 -10.841084480286 -10.50137138367 -10.166210174561 -9.8337898254395 -9.4995431900024 -9.1671237945557 -8.8328762054443 -8.5004568099976 -8.1662101745606 -7.8337903022766 -7.4995431900024 -7.1671237945557 -6.8328766822815 -6.5004572868347 -6.1662101745606 -5.8337907791138 -5.499543666840 -5.1671242713928 -4.8328766822815 -4.5004572868347 -4.1662101745606 -3.8337910175324 -3.499543666840 -3.1671242713928 -2.8328771591187 -2.5004577636719 -2.1662104129791 -1.8337911367416 -1.4995439052582 -1.167124509811 -0.83287727832794 -0.50045788288117 -0.16621069610120 0.16620869934560 0.50045591592789 0.83287531137466 1.1671224832535 1.4995418787003 1.8337891101837 2.1662085056305 2.5004556179047 2.8328750133514 3.1671223640442 3.4995417594910 3.83378887176514 4.16620826721191 4.50045537948608 4.83287477493286 5.16712188720703 5.49954128265381 5.83378887176514 6.16620826721191 6.50045537948608 6.83287477493286 7.16712188720703 7.49954128265381 7.83378839492798 8.16620826721191 8.50045490264893 8.8328742980957 9.16712188720703 9.49954128265381 9.83378791809082 10.1662120819092 10.501368522644 10.8410863876343 11.1907510757446 11.5520324707031 11.9300804138184 12.3262519836426 12.7453165054321 13.1881809234619 13.6591024398804 14.1584215164185 14.6897754669189 15.2528429031372 15.8505630493164 16.4818859100342 17.1490058898926 17.8501110076904 18.5866279602051 19.3559894561768 20.1588745117188 20.9919834136963 21.8552989959717 22.7448635101318 23.6600379943848 24.5962944030762 25.5524845123291 26.5236320495605 27.5082054138184 28.5009155273438 29.4999961853027 30.5000038146973 31.4999961853027 32.5000038146973 33.4999961853027 34.5000038146973 35.4999961853027 36.5000038146973 37.4999961853027 38.5000038146973 39.4999961853027 40.5000038146973 41.4999961853027 42.5000038146973 43.4999961853027 44.5000038146973 45.4999961853027 46.5000038146973 47.4999961853027 48.5000038146973 49.4999961853027 50.5000038146973 51.4999961853027 52.5000038146973 53.4999961853027 54.5000038146973 55.4999961853027 56.500003814697 57.499996185303 58.500003814697 59.499996185303 60.500003814697 61.4999961853027 62.500003814697 63.499996185303 64.5 65.5 66.5 67.5 68.5 69.5 70.5 71.5 72.5 73.5 74.5 75.5 76.5 77.5 78.5 79.5 80.5 81.5 82.5 83.5 84.5 85.5 86.5 87.5 88.5 89.5
xdef is 360 linear -279.5 1
Robert
_______________________________________________
gradsusr mailing list
gradsusr at gradsusr.org<mailto:gradsusr at gradsusr.org>
http://gradsusr.org/mailman/listinfo/gradsusr
_______________________________________________
gradsusr mailing list
gradsusr at gradsusr.org<mailto: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
_______________________________________________
gradsusr mailing list
gradsusr at gradsusr.org<mailto: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/20140214/476522bc/attachment-0001.html
More information about the gradsusr
mailing list