[Wrf-users] ERA-Interim

Atsuyoshi MANDA am at bio.mie-u.ac.jp
Sun Oct 2 17:43:23 MDT 2016


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