PRECIS support forum

It is currently Tue Mar 19, 2024 2:05 am

All times are UTC




Post new topic Reply to topic  [ 2 posts ] 
Author Message
PostPosted: Tue Oct 28, 2014 7:34 pm 
Offline

Joined: Tue Oct 28, 2014 6:27 pm
Posts: 15
this is short info of netCDF
-1 : Institut Source Ttype Levels Num Gridsize Num Dtype : Parameter name
1 : unknown unknown instant 1 1 4080 1 F32 : precip
Grid coordinates :
1 : lonlat > size : dim = 4080 nx = 68 ny = 60
x : first = 353.559998 last = 368.299988 inc = 0.219999854 degrees_east
y : first = 5.04000044 last = -7.93999958 inc = -0.220000267 degrees_north

and this is the ncdump -h of the netcdf:
dimensions:
x = 84 ;
y = 76 ;
ht = 1 ;
t = UNLIMITED ; // (93 currently)
variables:
float x(x) ;
x:standard_name = "longitude" ;
x:long_name = "longitude" ;
x:units = "degrees_east" ;
x:axis = "X" ;
float y(y) ;
y:standard_name = "latitude" ;
y:long_name = "latitude" ;
y:units = "degrees_north" ;
y:axis = "Y" ;
float ht(ht) ;
ht:units = "m" ;
ht:positive = "up" ;
ht:axis = "Z" ;
double t(t) ;
t:standard_name = "time" ;
t:units = "days since 1859-12-01 00:00:00" ;
t:calendar = "360_day" ;
float temp(t, ht, y, x) ;
temp:long_name = "TEMPERATURE AT 1.5M" ;
temp:units = "K" ;
temp:_FillValue = 2.e+20f ;
temp:missing_value = 2.e+20f ;
temp:source = "Unified Model Output:" ;
temp:name = "temp" ;
temp:title = "TEMPERATURE AT 1.5M" ;
temp:date = "01/12/59" ;
temp:time = "00:00" ;
temp:north_pole = 95.5f, 57.76f ;

// global attributes:
:CDI = "Climate Data Interface version 1.6.2 (http://code.zmaw.de/projects/cdi)" ;
:Conventions = "CF-1.4" ;
:history = "Sun Oct 19 17:19:42 2014: cdo addc,-273.15 jja_temp_75_05.nc jja_tempc_75_05.nc\n",
"Sun Oct 19 16:58:39 2014: cdo selmon,6,7,8, tmp_mon_1975-2005.NC jja_temp_75_05.nc\n",
"Sun Oct 19 16:56:42 2014: cdo selyear,1975/2005 tmonthly tmp_mon_1975-2005.NC\n",
"Sun Sep 7 10:51:41 CDT 2014 - XCONV V1.91 16-February-2006" ;
:CDO = "Climate Data Operators version 1.6.2 (http://code.zmaw.de/projects/cdo)" ;
}
Using matlab the plotting came out inseficient way. Simon and Dave looked in to it advised to post it. basd on that my quastion is, If rotated pole requires any conversion to regular lon lat or if any problem in the information that is provided, any advise is appreciated .


Top
 Profile  
 
PostPosted: Wed Oct 29, 2014 10:26 am 
Offline

Joined: Thu Mar 27, 2014 5:03 pm
Posts: 73
I think that as you suspected, you are having problems since your plotting software is not recognising that this is rotated pole data. The issue is probably due to the NetCDF meta-data (the data that describes what is in the file) is not cf compliant. See http://cfconventions.org/ for more information.

In order to create cf compliant NetCDF files from PP files, I recommend using the cfa tool. For information about it including how to install see http://cfpython.bitbucket.org/.

Unfortunately we have also found that some CDO commands strip meta-data from output files.

_________________
Simon Tucker
Met Office Hadley Centre


Top
 Profile  
 
Display posts from previous:  Sort by  
Post new topic Reply to topic  [ 2 posts ] 

All times are UTC


Who is online

Users browsing this forum: No registered users and 2 guests


You cannot post new topics in this forum
You cannot reply to topics in this forum
You cannot edit your posts in this forum
You cannot delete your posts in this forum
You cannot post attachments in this forum

Search for:
Jump to:  
cron
Powered by phpBB® Forum Software © phpBB Group