[RegCNET] convert netcdf files in Hybrid height coordinate to pressure level

[chao.tang]

Dear all,


I am searching for the way to convert netcdf files in hybrid height coordinate(units: m, output of HadGEM2-ES) to pressure level (EIN75).

As examples (Specific Humidity) bellow I attached two files with just one time step.




the first file (HadGEM2 one ) is in hybrid height coordinate, the second (EIN75) is in pressure levels.
Please note that they cover different domains, the EIN75 one was cutted to reduce the size for easier transfer by email.


dumps:
===================================================
netcdf hus_6hrLev_HadGEM2-ES_historical_r1i1p1_198412010600 {
dimensions:
	lon = 192 ;
	nb2 = 2 ;
	lat = 145 ;
	x = 2 ;
	lev = 38 ;
	time = UNLIMITED ; // (1 currently)
variables:
	double lon(lon) ;
		lon:standard_name = "longitude" ;
		lon:long_name = "longitude" ;
		lon:units = "degrees_east" ;
		lon:axis = "X" ;
		lon:bounds = "lon_bnds" ;
	double lon_bnds(lon, nb2) ;
	double lat(lat) ;
		lat:standard_name = "latitude" ;
		lat:long_name = "latitude" ;
		lat:units = "degrees_north" ;
		lat:axis = "Y" ;
		lat:bounds = "lat_bnds" ;
	double lat_bnds(lat, nb2) ;
	double lev(lev) ;
		lev:long_name = "hybrid height coordinate" ;
		lev:units = "m" ;
		lev:positive = "up" ;
		lev:axis = "Z" ;
		lev:bounds = "lev_bnds" ;
	double lev_bnds(lev, nb2) ;
	double time(time) ;
		time:standard_name = "time" ;
		time:long_name = "time" ;
		time:units = "days since 1859-12-01 00:00:00" ;
		time:calendar = "360_day" ;
	double orog(lat, lon) ;
		orog:standard_name = "surface_altitude" ;
		orog:long_name = "Surface Altitude" ;
		orog:units = "m" ;
		orog:comment = "height above the geoid; as defined here, \"\"the geoid\"\" is a surface of constant geopotential that, if the ocean were at rest, would coincide with mean sea level. Under this definition, the geoid changes as the mean volume of the ocean changes (e.g., due to glacial melt, or global warming of the ocean).  Report here the height above the present-day geoid.  Over ocean, report as 0.0" ;
	double b_bnds(lev, x) ;
		b_bnds:long_name = "vertical coordinate formula term: b(k+1/2)" ;
	double hus(time, lev, lat, lon) ;
		hus:standard_name = "specific_humidity" ;
		hus:long_name = "Specific Humidity" ;
		hus:units = "1" ;
		hus:_FillValue = 1.00000002004088e+20 ;
		hus:missing_value = 1.00000002004088e+20 ;
		hus:comment = "on all model levels" ;
		hus:original_name = "mo: m01s00i010" ;
		hus:history = "2010-09-03T11:46:37Z altered by CMOR: replaced missing value flag (-1.07374e+09) with standard missing value (1e+20)." ;
		hus:associated_files = "baseURL: http://cmip-pcmdi.llnl.gov/CMIP5/dataLocation gridspecFile: gridspec_fx_HadGEM2-ES_historical_r0i0p0.nc areacella: areacella_fx_HadGEM2-ES_historical_r0i0p0.nc" ;



===================================================
netcdf q_EIN75.198412010600 {
dimensions:
	longitude = 134 ;
	latitude = 54 ;
	level = 37 ;
	time = UNLIMITED ; // (1 currently)
variables:
	float longitude(longitude) ;
		longitude:standard_name = "longitude" ;
		longitude:long_name = "longitude" ;
		longitude:units = "degrees_east" ;
		longitude:axis = "X" ;
	float latitude(latitude) ;
		latitude:standard_name = "latitude" ;
		latitude:long_name = "latitude" ;
		latitude:units = "degrees_north" ;
		latitude:axis = "Y" ;
	double level(level) ;
		level:standard_name = "air_pressure" ;
		level:long_name = "pressure_level" ;
		level:units = "millibars" ;
		level:positive = "down" ;
		level:axis = "Z" ;
	double time(time) ;
		time:standard_name = "time" ;
		time:long_name = "time" ;
		time:units = "hours since 1900-01-01 00:00:00" ;
		time:calendar = "standard" ;
	double q(time, level, latitude, longitude) ;
		q:standard_name = "specific_humidity" ;
		q:long_name = "Specific humidity" ;
		q:units = "kg kg**-1" ;
		q:_FillValue = -32767. ;
		q:missing_value = -32767. ;
===================================================

Any idea would be really appreciated.

Best regards,
Chao TANG
-------------------------------------------------------------------------------------------------
Laboratoire d'Energétique, d'Electronique et Procédés 
Unité de Recherche EA 4079 
Université de La Réunion, LE2P
15 avenue René Cassin CS 92003
97744 ST DENIS CEDEX 9
Mobile: 0693 91 57 38 <mailto:secretariat.le2p at univ-reunion.fr>
Courriel: chao.tang at univ-reunion.fr  <mailto:secretariat.le2p at univ-reunion.fr>
-------------------------------------------------------------------------------------------------

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.ictp.it/pipermail/regcnet/attachments/20150622/e00431dd/attachment-0003.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: hus_6hrLev_HadGEM2-ES_historical_r1i1p1_198412010600.nc
Type: application/octet-stream
Size: 8738962 bytes
Desc: not available
URL: <https://lists.ictp.it/pipermail/regcnet/attachments/20150622/e00431dd/attachment-0002.obj>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.ictp.it/pipermail/regcnet/attachments/20150622/e00431dd/attachment-0004.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: q_EIN75.198412010600.nc
Type: application/octet-stream
Size: 2144844 bytes
Desc: not available
URL: <https://lists.ictp.it/pipermail/regcnet/attachments/20150622/e00431dd/attachment-0003.obj>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.ictp.it/pipermail/regcnet/attachments/20150622/e00431dd/attachment-0005.html>