Scheduled Downtime
On Friday 21 April 2023 @ 5pm MT, this website will be down for maintenance and expected to return online the morning of 24 April 2023 at the latest

Error while Running real.exe with ERA5 data in WRF-4.3.3

S

srikar_iist

New member
View attachment NAMELIST-WPS.txtHi,

I have reently installed WRF-4.3.3 along with WPS-4.3.1 in Intel HPC, and to test run the same in have take data from ERA5

Single levels: Product type:Reanalysis Variable:10m u-component of wind, 10m v-component of wind, 2m dewpoint temperature, 2m temperature, Land-sea mask, Mean sea level pressure, Sea surface temperature, Sea-ice cover, Skin temperature, Snow depth, Soil temperature level 1, Soil temperature level 2, Soil temperature level 3, Soil temperature level 4, Volumetric soil water layer 1, Volumetric soil water layer 2, Volumetric soil water layer 3, Volumetric soil water layer 4

Pressure Levels Product type:ReanalysisVariable:Divergence, Geopotential, Relative humidity, Specific humidity, Temperature, U-component of wind, V-component of windPressure level:1 hPa, 2 hPa, 3 hPa, 5 hPa, 7 hPa, 10 hPa, 20 hPa, 30 hPa, 50 hPa, 70 hPa, 100 hPa, 125 hPa, 150 hPa, 175 hPa, 200 hPa, 225 hPa, 250 hPa, 300 hPa, 350 hPa, 400 hPa, 450 hPa, 500 hPa, 550 hPa, 600 hPa, 650 hPa, 700 hPa, 750 hPa, 775 hPa, 800 hPa, 825 hPa, 850 hPa, 875 hPa, 900 hPa, 925 hPa, 950 hPa, 975 hPa, 1000 hPa

Year:2016 Month:August Day:15, 16, 17, 18, 19
Time:00:00, 01:00, 02:00, 03:00, 04:00, 05:00, 06:00, 07:00, 08:00, 09:00, 10:00, 11:00, 12:00, 13:00, 14:00, 15:00, 16:00, 17:00, 18:00, 19:00, 20:00, 21:00, 22:00, 23:00
Sub-region extraction:North 40°, West 50°, South 0°, East 110°
Format:GRIB

Attached is my namelist.wps and namelist.input files.

All metgrid files were generated and Im getting the following errors in execution of real.

module_io_quilt_old.F 2931 T
Ntasks in X 1 , ntasks in Y 1
Domain # 1: dx = 12000.000 m
Domain # 2: dx = 4000.000 m
REAL_EM V4.3.3 PREPROCESSOR
No git found or not a git repository, git commit version not available.

Parent domain
ids,ide,jds,jde 1 247 1 275
ims,ime,jms,jme -4 252 -4 280
ips,ipe,jps,jpe 1 247 1 275

DYNAMICS OPTION: Eulerian Mass Coordinate
alloc_space_field: domain 1 , 3090297656 bytes allocated
d01 2016-08-16_12:00:00 Yes, this special data is acceptable to use: OUTPUT FROM METGRID V4.3.1
d01 2016-08-16_12:00:00 Input data is acceptable to use: met_em.d01.2016-08-16_12:00:00.nc
metgrid input_wrf.F first_date_input = 2016-08-16_12:00:00
metgrid input_wrf.F first_date_nml = 2016-08-16_12:00:00
d01 2016-08-16_12:00:00 Timing for input 0 s.
d01 2016-08-16_12:00:00 flag_soil_layers read from met_em file is 1
Max map factor in domain 1 = 1.03. Scale the dt in the model accordingly.
Using sfcprs to compute psfc
forrtl: severe (174): SIGSEGV, segmentation fault occurred
Image PC Routine Line Source
real.exe 0000000002BA0BBA for__signal_handl Unknown Unknown
libpthread-2.17.s 00007FEE860295D0 Unknown Unknown Unknown
real.exe 0000000000450759 Unknown Unknown Unknown
real.exe 0000000000490D5D Unknown Unknown Unknown
real.exe 000000000049B97D Unknown Unknown Unknown
real.exe 000000000041A56F Unknown Unknown Unknown
real.exe 0000000000418D12 Unknown Unknown Unknown
libc-2.17.so 00007FEE85C6F3D5 __libc_start_main Unknown Unknown
real.exe 0000000000418C29 Unknown

NAMELIST-INPUT.txt
NAMELIST-WPS.txt

(tried for timestep of both 15s and 30s)

Kindly Help.
 

Attachments

  • NAMELIST-INPUT.txt
    3.8 KB · Views: 26
I am suspicious that you don't have the file that contains invariant variables, especially landmask and terrain height, which are two important variables for WRF/WPS.

Please let me know if I am wrong. Otherwise, please download the following two files (from NCAR CISL)

e5.oper.invariant.128_172_lsm.regn320sc.2016010100_2016010100.grb
e5.oper.invariant.128_129_z.regn320sc.2016010100_2016010100.grb
 
Last edited:
Dear Mr. Chen,

I will check. Thank you

Are these files in-built with the installation of WRF and WPS, or have to be installed for simulation using ERA5 data.
If I have to check for them, would you please recommend which folder the files mentioned have to be present in?

Thank you
 
Hi,
These files are not 'built-in' files. You have to download them, then go through the ungrib and metgrid processes.

Note that the two files need to be ungribbed separately (with start_date and end_date = '2016-01-01_00:00:00', and prefix = 'FILE2' in &unribThese files

Then rename FILE2:2016-01-01_00 to be FILE2, and set
&metgrid
fg_name = 'FILE1',
io_form_metgrid = 2,
constants_name = './FILE2'

From here you can continue to run metgrid.exe.
 
Hi ,

I couldn't find the files on the website, I might have overlooked them.

Mind sending the links to download the files.

And are these files required to be 'ungribbed' for simulations of different dates, as you have specifically mentioned that 'the two files need to be ungribbed separately (with start_date and end_date = '2016-01-01_00:00:00', and prefix = 'FILE2'

Thanks,
 
Mr Chen, There is no reply whatsoever from the above email, been a while. Is there any possible work around,
Would you mind trying the simulation at your end, the WPS- version that i have used is WPS-4.3
 
Ming Chen said:
Please look at the website here:

NCAR RDA Dataset ds630.0

NCAR RDA
rda.ucar.edu rda.ucar.edu
Click to expand...
Hi, Seems that they have changed the dataset to

NCAR RDA Dataset ds633.0

NCAR RDA
rda.ucar.edu rda.ucar.edu

I have downloaded the following files:
e5.oper.invariant.128_129_z.ll025sc.1979010100_1979010100.grb
e5.oper.invariant.128_172_lsm.ll025sc.1979010100_1979010100.grb

and my namelist.wps is as follows: (change of date from previous - August 2016 to May 2013)

&share
wrf_core = 'ARW',
max_dom = 2,
start_date = '1979-01-01_00:00:00','1979-01-01_00:00:00',
end_date = '1979-01-01_00:00:00','1979-01-01_00:00:00',
interval_seconds = 10800
io_form_geogrid = 2,
debug_level = 0,
/

&geogrid
parent_id = 1, 1,
parent_grid_ratio = 1, 3,
i_parent_start = 1, 75,
j_parent_start = 1, 113,
e_we = 248, 499,
e_sn = 297, 358,
s_we = 1, 1,
s_sn = 1, 1,
geog_data_res = 'default','default',
dx = 12000,
dy = 12000,
map_proj = 'lambert',
ref_lat = 21.251,
ref_lon = 83.529,
truelat1 = 21.251,
truelat2 = 21.251,
stand_lon = 83.529,
geog_data_path = '/opt/ohpc/pub/nWRF/WPS_GEOG/'
/

&ungrib
out_format = 'WPS',
prefix = 'FILE2',
/

&metgrid
fg_name = 'FILE1'
io_form_metgrid = 2,
constants_name = './FILE2',
/

OUTPUT:

FILE2:1979-01-01_00 - NAME CHANGED TO FILE2
$ mv FILE2\:1979-01-01_00 FILE2

$ ./metgrid.exe

Processing domain 1 of 2
./FILE2
Processing 1979-01-01_00
FILE1
WARNING: Couldn't open file FILE1:1979-01-01_00 for input.
ERROR: The mandatory field TT was not found in any input data.
application called MPI_Abort(MPI_COMM_WORLD, 61092320) - process 0
[unset]: write_line error; fd=-1 buf=:cmd=abort exitcode=61092320
:
system msg for write_line failure : Bad file descriptor

Would you please guide me in detail how to proceed?
Your previous help is deeply appreciated Mr. Chen
 
Last edited:
Please see the attached document that describes how to run WRF/WPS using ERA5 as input data.
 

Attachments

  • ERA5_download_instruction.pdf
    539.2 KB · Views: 71
Ming Chen said:
Please see the attached document that describes how to run WRF/WPS using ERA5 as input data.
Click to expand...
Sorry Mr. Chen, the issue is not resolved.

I am sharing my namelists here, pls try out if possible and revert back. Thanks

&share
wrf_core = 'ARW',
max_dom = 2,
start_date = '2013-05-27_00:00:00','2013-05-27_00:00:00',
end_date = '2013-06-02_00:00:00','2013-06-02_00:00:00',
interval_seconds = 10800,
io_form_geogrid = 2,
debug_level = 0,
/

&geogrid
parent_id = 1, 1, 1, 1, 1,
parent_grid_ratio = 1, 3, 3, 3, 3,
i_parent_start = 1, 33, 70, 55, 185,
j_parent_start = 1, 118, 125, 195, 148,
s_we = 1, 1, 1, 1, 1,
e_we = 248, 499, 103, 112, 103,
s_sn = 1, 1, 1, 1, 1,
e_sn = 297, 358, 103, 112, 103,
geog_data_res = 'default', 'default', 'default', 'default', 'default',
dx = 12000,
dy = 12000,
map_proj = 'lambert',
ref_lat = 21.251,
ref_lon = 83.529,
truelat1 = 21.251,
truelat2 = 21.251,
stand_lon = 83.529,
geog_data_path = '/opt/ohpc/pub/nWRF/WPS_GEOG/'


NAMELIST.INPUT:

&time_control
run_days = 6,
run_hours = 0,
run_minutes = 0,
run_seconds = 0,
start_year = 2013, 2013, 2014, 2014, 2014,
start_month = 05, 05, 03, 03, 03,
start_day = 27, 27, 24, 24, 24,
start_hour = 00, 00, 00, 00, 00,
start_minute = 00, 00, 00, 00, 00,
start_second = 00, 00, 00, 00, 00,
end_year = 2013, 2013, 2014, 2014, 2014,
end_month = 06, 06, 05, 05, 05,
end_day = 02, 02, 31, 31, 31,
end_hour = 00, 00, 00, 00, 00,
end_minute = 00, 00, 00, 00, 00,
end_second = 00, 00, 00, 00, 00,
interval_seconds = 10800,
input_from_file = .true., .true., .true., .true., .true.,
history_interval = 60, 60, 720, 720, 720,
frames_per_outfile = 1, 1, 1, 1, 1,
restart = .false.,
restart_interval = 360,
io_form_history = 2,
io_form_restart = 2,
io_form_input = 2,
io_form_boundary = 2,
debug_level = 0,
/

&domains
time_step = 36,
time_step_fract_num = 0,

time_step_fract_den = 1,
max_dom = 1,
e_we = 248, 499, 103, 112, 103,
e_sn = 297, 358, 103, 112, 103,
e_vert = 35, 35, 35, 35, 35,
num_metgrid_levels = 38,
num_metgrid_soil_levels = 3,

eta_levels = 1.000, 0.993, 0.983, 0.970, 0.954,
0.934, 0.909, 0.880, 0.845, 0.807,
0.765, 0.719, 0.672, 0.622, 0.571,
0.520, 0.468, 0.420, 0.376, 0.335,
0.298, 0.263, 0.231, 0.202, 0.175,
0.150, 0.127, 0.106, 0.088, 0.070,
0.055, 0.040, 0.026, 0.013, 0.000,
p_top_requested = 5000,
dx = 12000, 4000, 5000, 5000, 5000,
dy = 12000, 4000, 5000, 5000, 5000,
grid_id = 1, 2, 3, 4, 5,
parent_id = 1, 1, 1, 1, 1,
parent_grid_ratio = 1, 3, 3, 3, 3,
i_parent_start = 1, 75, 70, 55, 185,
j_parent_start = 1,113, 125, 195, 148,

parent_time_step_ratio = 1, 3, 3, 3, 3,
feedback = 0,
smooth_option = 0


&physics
mp_physics = 7, 7,
gsfcgce_hail = 0,
gsfcgce_2ice = 0,
mp_zero_out = 2,
mp_zero_out_thresh = 1.e-12,
ra_lw_physics = 5, 5,
ra_sw_physics = 5, 5,
radt = 15, 10,
sf_sfclay_physics = 5, 5,
sf_surface_physics = 2, 2,
sf_urban_physics = 0, 0,
bl_pbl_physics = 6, 6,

bldt = 0, 0,
cu_physics = 1, 0,
cudt = 15, 10,
isfflx = 1,
ifsnow = 0,
icloud = 1,
surface_input_source = 1,
num_soil_layers = 4,
maxiens = 1,
maxens = 3,
maxens2 = 3,
maxens3 = 16,
ensdim =144,
ishallow = 0,
cu_rad_feedback = .true.,
num_land_cat = 21,
rdlai2d = .true.,
bl_mynn_edmf = 1
/

&fdda
grid_fdda = 0, 0, 0, 0, 0,
gfdda_inname = "wrffdda_d<domain>",
gfdda_end_h = 48, 48,
gfdda_interval_m = 60, 60,
fgdt = 0, 0,
if_no_pbl_nudging_uv = 0, 0,
if_no_pbl_nudging_t = 0, 0,
if_no_pbl_nudging_q = 0, 0,
if_zfac_uv = 0, 0,
k_zfac_uv = 1, 1,
if_zfac_t = 0, 0,
k_zfac_t = 1, 1,
if_zfac_q = 0, 0,
k_zfac_q = 1, 1,
guv = 0.0003, 0.0003,
gt = 0.000, 0.000,
gq = 0.000, 0.000,
if_ramping = 0,
dtramp_min = 60.0,

io_form_gfdda = 2,

/

&dynamics
w_damping = 1,
diff_opt = 1, 1,
km_opt = 4, 4,
diff_6th_opt = 2, 2,
diff_6th_factor = 0.12, 0.12,
epssm = 0.1, 0.1,
base_temp = 290.
damp_opt = 3,
zdamp = 5000., 5000.,
dampcoef = 0.2, 0.2,
khdif = 0, 0,
kvdif = 0, 0,
non_hydrostatic = .true., .true.,
moist_adv_opt = 1, 1,
scalar_adv_opt = 1, 1,
v_mom_adv_order = 5, 5,
v_sca_adv_order = 5, 5,
h_mom_adv_order = 5, 5,
h_sca_adv_order = 5, 5,
/

&bdy_control
spec_bdy_width = 5,
spec_zone = 1,
relax_zone = 4,
specified = .true., .false.,
nested = .false., .true.,
/

&grib2
/

&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/
 
I am sorry to answer late. This post somehow was skipped. Just wonder whether you have solved your problem?
One obvious mistake is that, in your namelist.wps, you set

start_date = '1979-01-01_00:00:00','1979-01-01_00:00:00',
end_date = '1979-01-01_00:00:00','1979-01-01_00:00:00',

But in your namelist.input, you have a different time settings:

start_year = 2013, 2013, 2014, 2014, 2014,
start_month = 05, 05, 03, 03, 03,
start_day = 27, 27, 24, 24, 24,
start_hour = 00, 00, 00, 00, 00,
start_minute = 00, 00, 00, 00, 00,
start_second = 00, 00, 00, 00, 00,
end_year = 2013, 2013, 2014, 2014, 2014,
end_month = 06, 06, 05, 05, 05,
end_day = 02, 02, 31, 31, 31,
end_hour = 00, 00, 00, 00, 00,
end_minute = 00, 00, 00, 00, 00,
end_second = 00, 00, 00, 00, 00,

Note that the time in WPS and in WRF must match.
 
I have success in using ERA5 as input with WRF 4.3.2. Following the guide in the page below.

Initializing the WRF model with ERA5 (pressure level)

Note This is the tutorial of using WRF with ERA5 pressure level (38 vertical levels) data. If you prefer model level (138 vertical levels), please check another tutorial. I recommend the pressure level data, because the download speed is much faster. If you don’t have cdsapi, please check this...
dreambooker.site dreambooker.site
 
You must log in or register to reply here.
Top