Problems creating initial conditions for a global run from ERA5 data

C

cri09

New member
Hi there,
I'm Cristiano D'Amico and I attended the St Andrews MPAS-A tutorial last year.
After successfully installing and compiling MPAS-A v8.3.1 on ECMWF ATOS HPC, I was simply trying to replicate the test case of creating the static and initial conditions for a global 240-km quasi-uniform grid. I succeeded in creating the static.nc file, but when I tried to create the init.nc file starting from the ERA5:2014-09-10_00 intermediate file, the init_atmosphere program started and ended after less than 2 minutes, without reporting a specific error in the log.* files but also failing in producing the init.nc file.

From the log.init_atmosphere.0000.out file, the last message is:
...
Recomputing RH w.r.t. ice below freezing
###################################
Besides that, in the log file from the bash script I run on the HPC I only see this error message:
...
Caught signal 8 (Floating point exception: floating-point overflow)
==== backtrace (tid:2791658) ====
0 0x0000000000012990 __funlockfile() :0
1 0x000000000007317e __iscanonicall() ???:0
2 0x0000000000516236 __init_atm_cases_MOD_init_atm_case_gfs() /etc/ecmwf/nfs/dh2_2_perm_b/ecme1994/MPAS-Model/src/core_init_atmosphere/mpas_init_atm_cases.F:5207
3 0x0000000000579463 __init_atm_cases_MOD_init_atm_setup_case() /etc/ecmwf/nfs/dh2_2_perm_b/ecme1994/MPAS-Model/src/core_init_atmosphere/mpas_init_atm_cases.F:175
4 0x00000000004a1101 __init_atm_core_MOD_init_atm_core_run() /etc/ecmwf/nfs/dh2_2_perm_b/ecme1994/MPAS-Model/src/core_init_atmosphere/mpas_init_atm_core.F:92
5 0x00000000004092a1 __mpas_subdriver_MOD_mpas_run() /etc/ecmwf/nfs/dh2_2_perm_b/ecme1994/MPAS-Model/src/driver/mpas_subdriver.F:416
6 0x000000000040796e MAIN__() /etc/ecmwf/nfs/dh2_2_perm_b/ecme1994/MPAS-Model/src/driver/mpas.F:20
7 0x00000000004079d9 main() /etc/ecmwf/nfs/dh2_2_perm_b/ecme1994/MPAS-Model/src/driver/mpas.F:10
8 0x000000000003a865 __libc_start_main() ???:0
9 0x000000000040784e _start() ???:0
#####################################

I'm trying to attach below the files I believe could be the most useful for you to help me solving this problem, but at the moment of attaching the files to this thread I read the message "the upload failed because the temporary directory was missing. The site administrator will need to resolve this before any files can be uploaded" on every attached file.
Could you please allow me to attach them (otherwise it's difficult for me to clarify the situation originating the error)?
Thank you so much in advance.

Best regards,
Cristiano D'Amico
 
Hi Cristiano,

I believe this is a data issue. I am not sure where you download ERA5 data and how you process this data for MPAS initialization. Indeed we have troubles before when we try to use ERA5 as input for MPAS run.

To address this issue, we have created a python script specifically for ERA5 data processing. Please download the source code from :
github.com

GitHub - NCAR/era5_to_int: A simple Python script for converting ERA5 model-level netCDF files to the WPS intermediate format

A simple Python script for converting ERA5 model-level netCDF files to the WPS intermediate format - NCAR/era5_to_int
github.com github.com

Then follow the instruction to process the data.

You can download ERA5 from NSF NCAR GDEX Dataset d633000.

Please let me know if you have more questions.
 
Hi Ming,

thank you so much for your explanation and for the suggestion of using your era5_to_int script, which works perfectly!
However, now I'm facing another problem regarding ungrib.exe...
I would like to run an MPAS-A global simulation initialized with a ECMWF-IFS global analysis (at a single time instant), which I've downloaded from MARS archive in GRIB1 format (the only one available to me for the case study I'm analyzing). However, when I try to "ungrib" this file, what I obtain from the ungrib.exe program is the following error:

*** Starting program ungrib.exe ***
Start_date = 2014-11-07_00:00:00 , End_date = 2014-11-07_00:00:00
output format is WPS
Path to intermediate files is ./
End-of-record mark (7777) not found 0
Sec0(1) = 8492361 0
*** stopping in findgrib in gribcode ***\n
\tI could not find the GRIB string in the input file
\tafter testing the first 100,000 bytes.
\tThe file may be corrupt or it is not a GRIB file.
\tPerhaps a gzipped GRIB file or netcdf?\n
STOP findgrib
##########################################################

I checked the file and it is actually a GRIB file as when I type "file AN_IFS_20141107_0000_0.grib1" I obtain:
AN_IFS_20141107_0000_0.grib1: Gridded binary (GRIB) version 1.
I checked the fields in the file with "grib_ls" and "cdo sinfov" and they seem to be all okay. I also tried to dowload the IFS global analysis file using different grid types but none of these files is even read by the ungrib.exe program. In addition, I exclude that there could have been troubles in the compilation of ungrib since ungrib.exe worked well when I ran it on a NCEP-GFS FNL analysis file.

Do you have an idea of what could be the problem or have any suggestion to try to solve it? Or maybe did someone experience a similar problem in the past?
Thannk you in advance.
Best regards,
Cristiano D'Amico
 
ECMWF-IFS global analysis has some special issues and we cannot process IFS at present. As far as I know, WPS doesn't work although great efforts have been made years ago. Because of the hard situation NCAR is facing right now and because this is not our top priority, we suggest users to avoid using IFS data at present.
 
Last edited:
You must log in or register to reply here.
Back
Top