[Wrf-users] ERA-Interim
Atsuyoshi MANDA
am at bio.mie-u.ac.jp
Mon Oct 3 17:25:54 MDT 2016
I post the following email again since I have not received the
confirmation email. I am sorry for that if you receive this email twice.
Dear Juliet,
Thank you very much for your reply. I have checked my met_em file:
$ ncdump -c met_em.d01.2011-01-19_18:00:00.nc
netcdf met_em.d01.2011-01-19_18\:00\:00 {
dimensions:
Time = UNLIMITED ; // (1 currently)
DateStrLen = 19 ;
west_east = 299 ;
south_north = 259 ;
num_metgrid_levels = 33 ;
num_st_layers = 4 ;
num_sm_layers = 4 ;
south_north_stag = 260 ;
west_east_stag = 300 ;
I have omitted the pressure level data higher than 10 hPa. I think
that's why num_metgrid_levels in my met_em file equals 33.
I have set "num_metgrid_levels option" in my namelist.input to 33. Is it
OK?
I totally agree that I should use python script to download the ERA-Interim data, but I fails to install phyton api. My OpenSSL libbary seems to cause this issue. My simulation time is about 5 days and is not so long, thus I tentativly use the web interface (http://apps.ecmwf.int/datasets/data/interim-full-daily/levtype=pl/).
Thank you for your suggestion.
Atsuyoshi
2016-10-03 13:26 GMT+09:00 afwande juliet <afwandej965 at gmail.com>:
> Thanks for this email. This had to do with metgrid levels, One thing I
> realized with your Namelist. Input, you have set metgrid levels =33,as much
> as I know it, erainterim metgrid level should be =38, you need to check your
> metgrid files and confirm this.
> In most cases because of huge files needed from erainterim for wrf model, I
> use a python script to download all the data needed at one go. I can share
> this script if you needed it. But before then check that your metgrid levels
> in metgrid Data is correct
Hello, wrf-users,
I am having a problem to run the WRF3.4.1 using ERA-Interim. Ungrib and
metgrid seem to work, but real.exe fails:
Timing for input 0 s.
flag_soil_layers read from met_em file is 1
Max map factor in domain 1 = 1.06. Scale the dt in the model
accordingly.
ERROR in psfc: flag_psfc = 1, flag_soilhgt = 0, flag_slp = 0,
sfcp_to_sfcp =
F
-------------- FATAL CALLED ---------------
FATAL CALLED FROM FILE: <stdin> LINE: 990
not enough info for a p sfc computation
-------------------------------------------
-------------- FATAL CALLED ---------------
FATAL CALLED FROM FILE: <stdin> LINE: 990
not enough info for a p sfc computation
-------------------------------------------
I have checked met_em files and found it does not have flag_soilhgt
and flag_slp.
I have used the following Vtable files in ungrib/Variable_Tables
directory Vtable.ERA-interim.P (for pressure level data)
Vtable.ERA-interim.S1 (for surface data)
Vtable.ERA-interim.S2 (for time invariant field)
Any suggestions to solve the problem is appreciated.
Thank you,
Atsuyoshi
P.S. I am attaching some log files at the end of this email. If more
information is needed, please let me know.
--
Atsuyoshi Manda, PhD
Mie University, Japan
am at bio.mie-u.ac.jp
[1] Outputs of the g1print.exe:
[Pressure level data]
$ g1print.exe /work4/data/ERA-Interim/ERA-I.P.110116-24.grib |head
-10; echo; echo; echo; g1print.exe
/work4/data/ERA-Interim/ERA-I.P.110116-24.grib| tail -15
----------------------------------------------------
rec GRIB GRIB Lvl Lvl Lvl Time Fcst
Num Code name Code one two hour
----------------------------------------------------
1 129 Z 100 10 0 2011-01-16_00:00 +00
2 130 T 100 10 0 2011-01-16_00:00 +00
3 133 Q 100 10 0 2011-01-16_00:00 +00
4 157 R 100 10 0 2011-01-16_00:00 +00
5 131 U 100 10 0 2011-01-16_00:00 +00
6 132 V 100 10 0 2011-01-16_00:00 +00
6905 131 U 100 975 0 2011-01-24_18:00 +00
6906 132 V 100 975 0 2011-01-24_18:00 +00
6907 129 Z 100 1000 0 2011-01-24_18:00 +00
6908 130 T 100 1000 0 2011-01-24_18:00 +00
6909 133 Q 100 1000 0 2011-01-24_18:00 +00
6910 157 R 100 1000 0 2011-01-24_18:00 +00
6911 131 U 100 1000 0 2011-01-24_18:00 +00
6912 132 V 100 1000 0 2011-01-24_18:00 +00
***** End-Of-File on C unit 3
Copen: File = /work4/data/ERA-Interim/ERA-I.P.110116-24.grib
Fortran Unit = 0
UNIX File descriptor: 3
[Surface data]
$ g1print.exe /work4/data/ERA-Interim/ERA-I.S1.110116-24.grib |head
-10; echo; echo; echo; g1print.exe
/work4/data/ERA-Interim/ERA-I.S1.110116-24.grib| tail -15
----------------------------------------------------
rec GRIB GRIB Lvl Lvl Lvl Time Fcst
Num Code name Code one two hour
----------------------------------------------------
1 31 CI 1 0 0 2011-01-16_00:00 +00
2 32 ASN 1 0 0 2011-01-16_00:00 +00
3 33 RSN 1 0 0 2011-01-16_00:00 +00
4 34 SSTK 1 0 0 2011-01-16_00:00 +00
5 39 SWVL1 112 0 7 2011-01-16_00:00 +00
6 40 SWVL2 112 7 28 2011-01-16_00:00 +00
749 166 10V 1 0 0 2011-01-24_18:00 +00
750 167 2T 1 0 0 2011-01-24_18:00 +00
751 168 2D 1 0 0 2011-01-24_18:00 +00
752 170 STL2 112 7 28 2011-01-24_18:00 +00
753 183 STL3 112 28 100 2011-01-24_18:00 +00
754 235 SKT 1 0 0 2011-01-24_18:00 +00
755 236 STL4 112 100 255 2011-01-24_18:00 +00
756 238 TSN 1 0 0 2011-01-24_18:00 +00
***** End-Of-File on C unit 3
Copen: File = /work4/data/ERA-Interim/ERA-I.S1.110116-24.grib
Fortran Unit = 0
UNIX File descriptor: 3
[Time invariant field]
$ g1print.exe /work4/data/ERA-Interim/ERA-I.TINV.110116-24.grib |head
-10; echo; echo; echo; g1print.exe
/work4/data/ERA-Interim/ERA-I.TINV.110116-24.grib| tail -15
----------------------------------------------------
rec GRIB GRIB Lvl Lvl Lvl Time Fcst
Num Code name Code one two hour
----------------------------------------------------
1 29 TVL 1 0 0 1989-01-01_12:00 +00
2 129 Z 1 0 0 1989-01-01_12:00 +00
3 172 LSM 1 0 0 1989-01-01_12:00 +00
***** End-Of-File on C unit 3
----------------------------------------------------
rec GRIB GRIB Lvl Lvl Lvl Time Fcst
Num Code name Code one two hour
----------------------------------------------------
1 29 TVL 1 0 0 1989-01-01_12:00 +00
2 129 Z 1 0 0 1989-01-01_12:00 +00
3 172 LSM 1 0 0 1989-01-01_12:00 +00
***** End-Of-File on C unit 3
Copen: File = /work4/data/ERA-Interim/ERA-I.TINV.110116-24.grib
Fortran Unit = 0
UNIX File descriptor: 3
[2] namelist.wps
$ cat namelist.wps
!
! PL1101_R01.00
!
! /work1/am/WRF3.4.1.par/WPS
! [~/WRF3.4.1.par/WPS]
&share
wrf_core = 'ARW',
max_dom = 1,
start_date = "2011-01-19_18:00:00", "2011-01-19_18:00:00",
"2011-01-19_18:00:00"
end_date = "2011-01-24_00:00:00", "2016-01-24_00:00:00",
"2011-01-24_00:00:00"
interval_seconds = 21600
io_form_geogrid = 2,
/
&geogrid
parent_id = 1, 1, 2
parent_grid_ratio = 1, 4, 4
i_parent_start = 1, 21, 39
j_parent_start = 1, 8, 58
e_we = 300, 113, 241
e_sn = 260, 141, 269
geog_data_res = '10m','2m','30s'
dx = 10000 !, 10000, 2500
dy = 10000 !, 10000, 2500
map_proj = 'polar',
ref_lat = 72.00,
ref_lon = -5.00,
truelat1 = 72.00,
stand_lon = -5.00,
geog_data_path = "/work4/data/wrf/geog"
/
&ungrib
out_format = "WPS",
prefix = "LAND" !"SURF", !"PRES",
/
&metgrid
fg_name = "PRES", "SURF",
constants_name='LAND:1989-01-01_12'
io_form_metgrid = 2,
/
[3] namelist.input
$ cat namelist.input
!
! PL1101_R01.00
!
!
&time_control
run_days = 4,
run_hours = 6,
run_minutes = 0,
run_seconds = 0,
start_year = 2011,2011,2011
start_month = 01,01,01
start_day = 19,19,19
start_hour = 18,18,18
start_minute = 00,00,00
start_second = 00,00,00
end_year = 2011,2011,2011
end_month = 01,01,01,
end_day = 24,24,24
end_hour = 0, 0,0
end_minute = 00,00,00
end_second = 00,00,00
interval_seconds = 21600
input_from_file = .true.,.true.,.true.,
history_interval = 180,180,60,
frames_per_outfile = 10000, 10000, 10000,
restart = .false.,
restart_interval = 500000, !restart_interval
io_form_history = 2
io_form_restart = 2
io_form_input = 2
io_form_boundary = 2
debug_level = 0
auxinput1_inname = "../../../WPS/met_em.d<domain>.
<date>"
auxinput4_inname = "wrflowinp_d<domain>",
auxinput4_interval = 360,
io_form_auxinput4 = 2,
/
&domains
time_step = 60, !60,
time_step_fract_num = 0,
time_step_fract_den = 1,
max_dom = 3,
e_we = 300, 113, 241
e_sn = 260, 141, 269
p_top_requested = 5000,
num_metgrid_levels = 33 !32, !fnl(2015 or later), !27->fnl,
22->MSM, 18->GSM,
num_metgrid_soil_levels = 4,
dx = 10000 ,10000, 2500
dy = 10000 ,10000, 2500
grid_id = 1, 2, 3
parent_id = 0, 1, 2
i_parent_start = 1, 21, 39
j_parent_start = 1, 8, 58
parent_grid_ratio = 1, 4, 4,
parent_time_step_ratio = 1, 4, 4,
feedback = 0,
smooth_option = 0
/
&physics
mp_physics = 3, 3, 3,
ra_lw_physics = 1, 1, 1,
ra_sw_physics = 1, 1, 1,
radt = 30, 30, 30,
sf_sfclay_physics = 1, 1, 1,
sf_surface_physics = 2, 2, 2,
bl_pbl_physics = 1, 1, 1,
bldt = 0, 0, 0,
cu_physics = 1, 1, 0,
cudt = 5, 5, 5,
isfflx = 1,
ifsnow = 0,
icloud = 1,
surface_input_source = 1,
num_soil_layers = 4,
sf_urban_physics = 0, 0, 0,
sst_update = 1,
/
! Set up the analysis nudging fdda namelist
! http://www2.mmm.ucar.edu/wrf/users/wrfv2/How_to_run_grid_fdda.html
&fdda
grid_fdda = 0,0,0,
gfdda_inname = "wrffdda_d<domain>"
gfdda_end_h = 1000,1000,0,
gfdda_interval_m = 360,360,360,
io_form_gfdda = 2,
fgdt = 0
if_no_pbl_nudging_uv = 0,1,1,
if_no_pbl_nudging_t = 0,1,1,
if_no_pbl_nudging_q = 0,1,1,
if_zfac_uv = 0,0,0,
k_zfac_uv = 0,0,0,
if_zfac_t = 0,0,0,
k_zfac_t = 0,0,0
if_zfac_q = 0,0,0,
k_zfac_q = 10,10,10,
guv = 0.0003,0.0003,0.0003,
gt = 0.0003,0.0003,0.0003,
gq = 0.0003,0.0003,0.0003,
if_ramping = 0
dtramp_min = 60.0,
/
&dynamics
w_damping = 0,
diff_opt = 1,
km_opt = 4,
diff_6th_opt = 0, 0, 0,
diff_6th_factor = 0.12, 0.12, 0.12,
base_temp = 290.
damp_opt = 0,
zdamp = 5000., 5000., 5000.,
dampcoef = 0.2, 0.2, 0.2
khdif = 0, 0, 0,
kvdif = 0, 0, 0,
non_hydrostatic = .true., .true., .true.,
moist_adv_opt = 1, 1, 1,
scalar_adv_opt = 1, 1, 1,
/
&bdy_control
spec_bdy_width = 5,
spec_zone = 1,
relax_zone = 4,
specified = .true., .false.,.false.,
nested = .false., .true., .true.,
/
&grib2
/
&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/
More information about the Wrf-users
mailing list