教你转Noah-MP陆面模式，获取陆面多种变量

ranmao

andrewsoong 发表于 2020-5-7 09:27
你的GFS数据是2019年6月之前的？

2020年5月5日的，我才下载了FNL（2020-05-06）的数据，运行就没问题，gfs就会报tsk不合理=============================================================
taskid: 0 hostname: pc-157
Ntasks in X            1 , ntasks in Y            1
REAL_EM V4.1.5 PREPROCESSOR
*************************************
Parent domain
ids,ide,jds,jde            1         100           1         124
ims,ime,jms,jme           -4         105          -4         129
ips,ipe,jps,jpe            1         100           1         124
*************************************
DYNAMICS OPTION: Eulerian Mass Coordinate
   alloc_space_field: domain            1 ,             399234548  bytes allocated
Time period #   1 to process = 2020-05-05_00:00:00.
Time period #   2 to process = 2020-05-05_01:00:00.
Time period #   3 to process = 2020-05-05_02:00:00.
Time period #   4 to process = 2020-05-05_03:00:00.
Time period #   5 to process = 2020-05-05_04:00:00.
Time period #   6 to process = 2020-05-05_05:00:00.
Time period #   7 to process = 2020-05-05_06:00:00.
Time period #   8 to process = 2020-05-05_07:00:00.
Time period #   9 to process = 2020-05-05_08:00:00.
Time period #  10 to process = 2020-05-05_09:00:00.
Time period #  11 to process = 2020-05-05_10:00:00.
Time period #  12 to process = 2020-05-05_11:00:00.
Time period #  13 to process = 2020-05-05_12:00:00.
Time period #  14 to process = 2020-05-05_13:00:00.
Time period #  15 to process = 2020-05-05_14:00:00.
Time period #  16 to process = 2020-05-05_15:00:00.
Time period #  17 to process = 2020-05-05_16:00:00.
Time period #  18 to process = 2020-05-05_17:00:00.
Time period #  19 to process = 2020-05-05_18:00:00.
Time period #  20 to process = 2020-05-05_19:00:00.
Time period #  21 to process = 2020-05-05_20:00:00.
Time period #  22 to process = 2020-05-05_21:00:00.
Time period #  23 to process = 2020-05-05_22:00:00.
Time period #  24 to process = 2020-05-05_23:00:00.
Total analysis times to input =   24.

-----------------------------------------------------------------------------

Domain  1: Current date being processed: 2020-05-05_00:00:00.0000, which is loop #   1 out of   24
configflags%julyr, %julday, %gmt:        2020         126   0.00000000   
d01 2020-05-05_00:00:00  Yes, this special data is acceptable to use: OUTPUT FROM METGRID V4.1
d01 2020-05-05_00:00:00  Input data is acceptable to use: met_em.d01.2020-05-05_00:00:00.nc
metgrid input_wrf.F first_date_input = 2020-05-05_00:00:00
metgrid input_wrf.F first_date_nml = 2020-05-05_00:00:00
d01 2020-05-05_00:00:00 Timing for input          0 s.
d01 2020-05-05_00:00:00          flag_soil_layers read from met_em file is  1
d01 2020-05-05_00:00:00 Turning off use of MAX WIND level data in vertical interpolation
d01 2020-05-05_00:00:00 Turning off use of TROPOPAUSE level data in vertical interpolation
Max map factor in domain 1 =  1.00. Scale the dt in the model accordingly.
Using sfcprs3 to compute psfc
d01 2020-05-05_00:00:00 No average surface temperature for use with inland lakes
Assume Noah LSM input
          15          32  water and category            5  both at 50%, landmask is    0.00000000   
          54          67  water and category            7  both at 50%, landmask is    0.00000000   
          59          67  water and category           12  both at 50%, landmask is    0.00000000   
          50          72  water and category           12  both at 50%, landmask is    0.00000000   
          86          73  water and category           12  both at 50%, landmask is    0.00000000   
          90          78  water and category           12  both at 50%, landmask is    0.00000000   
          91          78  water and category           12  both at 50%, landmask is    0.00000000   
          80          86  water and category           12  both at 50%, landmask is    0.00000000   
          52          99  water and category           12  both at 50%, landmask is    0.00000000   
          33         116  water and category            6  both at 50%, landmask is    0.00000000   
LAND  CHANGE =            0
WATER CHANGE =            0
error in the grid%tsk
i,j=          22          30
grid%landmask=   0.00000000   
grid%tsk, grid%sst, grid%tmn=   0.00000000       0.00000000       0.00000000   
-------------- FATAL CALLED ---------------
FATAL CALLED FROM FILE:  <stdin>  LINE:    2977
grid%tsk unreasonable
-------------------------------------------

andrewsoong

ranmao 发表于 2020-5-7 11:50
2020年5月5日的，我才下载了FNL（2020-05-06）的数据，运行就没问题，gfs就会报tsk不合理=============== ...

你这里面包含了湖泊，可能会报错。你把namelist.input中那个下垫面数据改成1,默认是3？记不清了，或者记反了

ranmao

andrewsoong 发表于 2020-5-7 12:40
你这里面包含了湖泊，可能会报错。你把namelist.input中那个下垫面数据改成1,默认是3？记不清了，或者记 ...

&physics                 
mp_physics               = 8,        8,
ra_lw_physics            = 1,        1,
ra_sw_physics            = 1,        1,
radt                     = 10,       10,
sf_sfclay_physics        = 1,        1,
sf_surface_physics       = 2,        2,
bl_pbl_physics           = 1,        1,
bldt                     = 0,        0,
cu_physics               = 0,        0,
cudt                     = 5,        5,
isfflx                   = 1,
ifsnow                   = 0,
icloud                   = 1,
surface_input_source     = 1,
num_soil_layers          = 4,
sf_urban_physics         = 0,        0,
maxiens                  = 1,
maxens                   = 3,
maxens2                  = 3,
maxens3                  = 16,
ensdim                   = 144,
/

sf_surface_physics？我主要是 后续要用CMAQ模型，用Noha LSM比较契合

andrewsoong

ranmao 发表于 2020-5-7 12:51
&physics                 
mp_physics               = 8,        8,
ra_lw_physics            = 1 ...

surface_input_source     = 3 试试

ranmao

andrewsoong 发表于 2020-5-7 16:04
surface_input_source     = 3 试试

还是不得行，不可能替换wrf自带的下垫面数据吧？

柯一

谢谢楼主分享

andrewsoong

ranmao 发表于 2020-5-7 16:55
还是不得行，不可能替换wrf自带的下垫面数据吧？

你换了下垫面？你换成不带湖泊的试试

ranmao

andrewsoong 发表于 2020-5-7 19:36
你换了下垫面？你换成不带湖泊的试试

下垫面替换我还弄不来，初级人员

andrewsoong

ranmao 发表于 2020-5-8 08:55
下垫面替换我还弄不来，初级人员

就是WPS那步，你选不要湖泊的陆地覆盖

ranmao

andrewsoong 发表于 2020-5-8 09:02
就是WPS那步，你选不要湖泊的陆地覆盖

大神，wps里面好像没看到这个选项呢？



=======================================
&share
wrf_core = 'ARW',
max_dom = 2,
start_date = '2020-05-07_00:00:00', '2020-05-07_00:00:00',
end_date   = '2020-05-08_01:00:00', '2020-05-08_01:00:00',
interval_seconds = 3600,
io_form_geogrid = 2,
opt_output_from_geogrid_path = '/sqwork/project/wrf_domains/wrf_test02',
debug_level = 0,
/

&geogrid
parent_id         = 1,1,
parent_grid_ratio = 1,3,
i_parent_start    = 1,28,
j_parent_start    = 1,32,
e_we          = 100,85,
e_sn          = 104,76,
geog_data_res = '30s','2m',
dx = 5800,
dy = 5800,
map_proj =  'lambert',
ref_lat   = 30.184,
ref_lon   = 106.282,
truelat1  = 30.184,
truelat2  = 30.184,
stand_lon = 106.282,
geog_data_path = '/sqwork/data/wps_geo/',
opt_geogrid_tbl_path = '/sqwork/project/wrf_domains/wrf_test02',
ref_x = 50.0,
ref_y = 52.0,
/

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

&metgrid
fg_name = 'FILE',
io_form_metgrid = 2,
opt_output_from_metgrid_path = '/sqwork/project/wrf_domains/wrf_test02',
opt_metgrid_tbl_path = '/sqwork/project/wrf_domains/wrf_test02',
/

&mod_levs
press_pa = 201300 , 200100 , 100000 ,
             95000 ,  90000 ,
             85000 ,  80000 ,
             75000 ,  70000 ,
             65000 ,  60000 ,
             55000 ,  50000 ,
             45000 ,  40000 ,
             35000 ,  30000 ,
             25000 ,  20000 ,
             15000 ,  10000 ,
              5000 ,   1000
/


&domain_wizard
grib_data_path = '/sqwork/data/gfs',
grib_vtable = 'null',
dwiz_name    =wrf-test02
dwiz_desc    =
dwiz_user_rect_x1 =6375
dwiz_user_rect_y1 =1286
dwiz_user_rect_x2 =6509
dwiz_user_rect_y2 =1407
dwiz_show_political =true
dwiz_center_over_gmt =true
dwiz_latlon_space_in_deg =10
dwiz_latlon_linecolor =-8355712
dwiz_map_scale_pct =50.0
dwiz_map_vert_scrollbar_pos =1105
dwiz_map_horiz_scrollbar_pos =6059
dwiz_gridpt_dist_km =5.8
dwiz_mpi_command =null
dwiz_tcvitals =null
dwiz_bigmap =Y
/