I currently need to run ICON with a specific data source. I have acquired and remapped all the fields needed, but they are recorded in pressure levels.
I would like to know, is it possible to feed ICON Initial/Boundary condition in pressure coordination? If we can’t use pressure coordination data source, are there any tools to interpolate the data to hybrid coordination? If we can, how should I adjust the namelist?
Thanks for any help.
Did you look into the example of using e.g. IFS / ERA5 data as initial and boundary conditions? Is that helping you? Or do you need something even more specific?
Thank you for your reply. The example of using IFS/ERA5 data taught me to remap the rotated grid to ICON’s triangular grid. But I think IFS/ERA5 data used in test examples are already in hybrid vertical coordinations. What I need to use is an Australian reanalysis dataset, and that dataset is only in the pressure-level.
The horizontal remapping is an important step before the simulation - the vertical is done during. They have hybrid, but hybrid pressure as far as I remember. Did you have a look here? Initial and Boundary Data from IFS — ICON documentation Also be aware, that you have to choose the correct init_mode in the initicon_nml (maybe try 2). If your simulation is crashing, maybe you could provide some information on the error?
Yes, I have looked the ICON tutorial and set the init_mode in initicon_nml to 2. The error was
“mo_initicon: MODE_IFS: perform initialization with IFS analysis
FINISH PE: 0 mo_initicon_io:read_extana_atm netCDF error: NetCDF: Invalid dimension ID or name“
It didn’t state clearly what was the exact error. My guess is the vertical coordination, because we have all required fields.
The ICON documentation have a side note about The IFS Vertical Coordinate. I need to say I don’t fully understand it. Does that mean if we provide the pressure coordination IFS data + LNPS/PS + GEOP_SFC, then ICON model can automatically do the vertical interpolation? Or we need use some other tools? What is the shortname of the pressure (pres or pressure or Pressure or …)? I’m confused because the sample command to get IFS data was “mars4icon_smi -a 70.0/-10.0/40.0/20.0 -r 1279 -l 1/to/137 -g 0.1/0.1 -d 2025100100 -O -L 1 -o ifs_2025100100.grb -p 5 -A”, which selects level 1 to 137. So, it was not dealing with pressure coordination IFS data, but hybrid model levels.
My guess would rather be, that there is some issue with the dimension / name of a variable - and that maybe the horizontal interpolation did not work. Did you run ICON already? Maybe with other forcing where you could compare the dimensions and definitions?
The pressure is usually calculated based on surface pressure and some coefficients, that should work. We could have a closer look, but for that it would be helpful for me to have a look at your data and / or log file. Unfortunately your chat is disabled, maybe you could allow that and we could try to find a solution?
Thank you for your reply. I’m quite confident that the horizontal remapping is correct. I will check other fields once again, before transforming my data to hybrid levels.
How can I enable the “chat”? I just register on this forum for one day, so I’m not familiar with that. Should I leave my email in this post?
Thank you again for your help.