PRECIS support forum
http://forum.precisrcm.com/

plotting converted netcdf file
http://forum.precisrcm.com/viewtopic.php?f=3&t=48
Page 1 of 1

Author:  gkebede12 [ Tue Oct 28, 2014 7:34 pm ]
Post subject:  plotting converted netcdf file

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 .

Author:  sotucker [ Wed Oct 29, 2014 10:26 am ]
Post subject:  Re: plotting converted netcdf file

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.

Page 1 of 1 All times are UTC
Powered by phpBB® Forum Software © phpBB Group
https://www.phpbb.com/