HadGEM2-ES: hybrid height coordinate to levels pressure
Added by alan llacza about 10 years ago
Hello everybody.
I am using data HadGEM2-ES to simulate with WRF model, however not achieving convert coordinates hydrid to pressure levels needed to WRF. I tried with the "ml2pl" command, but send me the following error:
cdo ml2hl (Warning): No data on hybrid model level found!
Any ideas or steps for this?
NOTE: Attached 2 files, the first with the original data and the second data surface pressure.
Alan Llacza
| ta_orig_20070301.nc (4.16 MB) ta_orig_20070301.nc | Original date | ||
| ta_ps_20070301.nc (4.27 MB) ta_ps_20070301.nc | Temperature and Surface pressure |
Replies (7)
RE: HadGEM2-ES: hybrid height coordinate to levels pressure - Added by Uwe Schulzweida about 10 years ago
The CDO operator mp2pl supports only 3D model data on hybrid sigma pressure coordinates. Your data contains 3D data on hybrid height coordinates. It is not possible to process this data with mp2pl, sorry.
RE: HadGEM2-ES: hybrid height coordinate to levels pressure - Added by Viktoria mohr about 9 years ago
Hi there,
I'd like to pick up the question from Alan Llazca. I would like to work with HadGEM2-ES data, which uses hybrid height in [m] as vertical coordinate, and convert it into pressure levels.
ncdump gives me the lev formula "z = a + b*orog " , so oviously ml2pl won't work due to missing p0 and ps which is given in Datasets which use hybrid sigma pressure coordinate as vertical level...
So is there maybe now (as the question from Alan was posted almost a year ago) any chance to convert hybrid height to pl with CDO? In the CDO users guide I found the operator mh2pl which appeared in one sentence only but is not further described (I also tried to use it but it doesn't work). Any other ideas?
Many Thanks
Viktoria
RE: HadGEM2-ES: hybrid height coordinate to levels pressure - Added by alan llacza about 9 years ago
Hello !
I send references and scripts (examples).
There is 2 steps to interpolate to levels pressure:
a) Find pressure value over height coordinate
REFERENCE: Program_documentation_hadcmip2caf_v13.pf => hadcmip2caf.f90
Script (GrADS): mkNP_th.gs (Only Temperature and Humidity). The script to components of the wind is similar.
NOTE: Change value of "ak_tq" and "bk_tq" to components of the wind
NOTE2: Horizontal Interpolation of components of wind to grid Temperature.
b) Interpolation between pressure levels
REFERENCE: http://moe.met.fsu.edu/~rhart/software/pinterp.gs
SCRIPT: int2_th.gs
prs_th.ctl (File of to help script GrADS).
Alan LLacza
| EXAMPLE.tar.gz (4.62 KB) EXAMPLE.tar.gz |
RE: HadGEM2-ES: hybrid height coordinate to levels pressure - Added by Viktoria mohr about 9 years ago
Hi Alan,
thanks a lot I will try this!
cheers
Viktoria
RE: HadGEM2-ES: hybrid height coordinate to levels pressure - Added by KONSTANTINA POLITI about 9 years ago
Hi guys!
I am trying to work with HadGEM2-ES and I don't know how to make the conversion of coordinates hydrid to pressure levels needed to WRF, as I noticed the problem that Victoria indicated.
I am not familiar with GrADS, is this a complex procedure? I work in ubuntu.
I have to download the software and make some changes in the files of the Alan's Example folder?
I wish you could help me just a little...thank you very much in advance.
Nadia
RE: HadGEM2-ES: hybrid height coordinate to levels pressure - Added by Viktoria mohr about 9 years ago
Hi Nadia,
I am sorry, I still can't tell you any experience as I had to work on another dataset before I can process the HadGEM data...
so as soon as I can tell anything, i will let you know!
hope you won't become desperate...
regards
Viktoria
RE: HadGEM2-ES: hybrid height coordinate to levels pressure - Added by KONSTANTINA POLITI about 9 years ago
Hello,
Thank you Victoria!
I have already installed Grads, but I am working too to another dataset and preprocessing.
When I ll be focus on Hadgem and Grads, if I have some progress I will also inform you!
Cheers,
Nadia