[Wrf-users] Using HadGEM2-ES cmip5 output in WRF

[Dominikus Heinzeller]

Hi,

On the website you referred to http://www.meteo.unican.es/wiki/cordexwrf/SoftwareTools/CmorPreprocessor see section 2. Transforming hybrid levels to pressure levels. This method converts the data before running it through WPS.

An alternative solution is to use the hybrid level data and let WPS do the transformation. For this you need to obtain the coefficients a and b (see the CmorPreprocessor website) and follow the guidelines here: http://www2.mmm.ucar.edu/wrf/users/docs/user_guide/users_guide_chap3.htm#_WPS_Utility_Programs C. calc_ecmwf_p.exe.

Hope this helps. Regards,

Dom

On 23/05/2014, at 2:25 pm, jagan TNAU <jagan at tnau.ac.in> wrote:

> Dear All,
> 
> I have cmip5 HadGEM2-ES netcdf file which has z axis as hybrid vertical height coordinate (meters). I want this to be interpolated to required pressure level so that I can ingest in to the WRF. The CMOR pre-processor available at http://www.meteo.unican.es/wiki/cordexwrf/SoftwareTools/CmorPreprocessor does not provide procedure for this kind of files. Is there any procedure available for using this data in WRF.
> 
> Please find below the ncdump output for netcdf file to be interpolated.
> 
> Can I get some help on this please.
> 
> --------------------------------------------------
> netcdf va_6hrLev_HadGEM2-ES_rcp45_r1i1p1_200512010600-200603010000 {
> dimensions:
>     time = UNLIMITED ; // (360 currently)
>     lev = 38 ;
>     lat = 144 ;
>     lon = 192 ;
>     bnds = 2 ;
> variables:
>     double time(time) ;
>         time:units = "days since 1859-12-01" ;
>         time:calendar = "360_day" ;
>         time:axis = "T" ;
>         time:long_name = "time" ;
>         time:standard_name = "time" ;
>     double lev(lev) ;
>         lev:bounds = "lev_bnds" ;
>         lev:units = "m" ;
>         lev:axis = "Z" ;
>         lev:positive = "up" ;
>         lev:long_name = "hybrid height coordinate" ;
>         lev:standard_name = "atmosphere_hybrid_height_coordinate" ;
>         lev:formula = "z(k,j,i) = a(k) + b(k)*orog(j,i)" ;
>         lev:formula_terms = "a: lev b: b orog: orog" ;
>     double lev_bnds(lev, bnds) ;
>         lev_bnds:formula = "z(k,j,i) = a(k) + b(k)*orog(j,i)" ;
>         lev_bnds:standard_name = "atmosphere_hybrid_height_coordinate" ;
>         lev_bnds:units = "m" ;
>         lev_bnds:formula_terms = "a: lev_bnds b: b_bnds orog: orog" ;
>     double b(lev) ;
>         b:long_name = "vertical coordinate formula term: b(k)" ;
>     float orog(lat, lon) ;
>         orog:standard_name = "surface_altitude" ;
>         orog:long_name = "Surface Altitude" ;
>         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" ;
>         orog:units = "m" ;
>     double b_bnds(lev, bnds) ;
>         b_bnds:long_name = "vertical coordinate formula term: b(k+1/2)" ;
>     double lat(lat) ;
>         lat:bounds = "lat_bnds" ;
>         lat:units = "degrees_north" ;
>         lat:axis = "Y" ;
>         lat:long_name = "latitude" ;
>         lat:standard_name = "latitude" ;
>     double lat_bnds(lat, bnds) ;
>     double lon(lon) ;
>         lon:bounds = "lon_bnds" ;
>         lon:units = "degrees_east" ;
>         lon:axis = "X" ;
>         lon:long_name = "longitude" ;
>         lon:standard_name = "longitude" ;
>     double lon_bnds(lon, bnds) ;
>     float va(time, lev, lat, lon) ;
>         va:standard_name = "northward_wind" ;
>         va:long_name = "Northward Wind" ;
>         va:comment = "on all model levels" ;
>         va:units = "m s-1" ;
>         va:original_name = "mo: m01s00i003" ;
>         va:history = "2010-09-08T14:00:25Z altered by CMOR: replaced missing value flag (-1.07374e+09) with standard missing value (1e+20)." ;
>         va:missing_value = 1.e+20f ;
>         va:_FillValue = 1.e+20f ;
>         va:associated_files = "baseURL: http://cmip-pcmdi.llnl.gov/CMIP5/dataLocation gridspecFile: gridspec_fx_HadGEM2-ES_rcp45_r0i0p0.nc" ;
> 
> 
> -- 
> With regards
> 
> Dr.R.Jagannathan                
> Professor of Agronomy
> Agro Climate Research Centre
> Tamil Nadu Agricultural University 
> Coimbatore - 641 003 India        
> 
> PHONE:  Mob: +91 94438 89891
> 
> DO NOT PRINT THIS E-MAIL UNLESS NECESSARY. THE ENVIRONMENT CONCERNS US ALL.
> _______________________________________________
> Wrf-users mailing list
> Wrf-users at ucar.edu
> http://mailman.ucar.edu/mailman/listinfo/wrf-users

-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/wrf-users/attachments/20140523/70531abf/attachment.html