WRFChem restart报错

最新推荐文章于 2024-07-14 11:42:14 发布

肉嘟嘟的科研笔记

最新推荐文章于 2024-07-14 11:42:14 发布

阅读量183

点赞数 2

分类专栏： WRFChem学习文章标签：笔记

本文链接：https://blog.csdn.net/weixin_64799803/article/details/139467654

版权

WRFChem学习专栏收录该内容

6 篇文章 0 订阅

订阅专栏

报错信息

DYNAMICS OPTION: Eulerian Mass Coordinate
   alloc_space_field: domain            1 ,              390974700  bytes allocated
 RESTART run: opening wrfrst_d01_2017-01-14_00:00:00 for reading
   MMINLU error on input
Timing for processing restart file for domain        1:   19.52120 elapsed seconds
Max map factor in domain 1 =  1.10. Scale the dt in the model accordingly.
-------------- FATAL CALLED ---------------
FATAL CALLED FROM FILE:  <stdin>  LINE:    9252
STOP in Noah-MP read_mp_veg_parameters

检查 wrfrst_d01_2017-01-14_00:00:00 发现 MMINLU = "" ; 而wrfinput 的 MMINLU = "MODIFIED_IGBP_MODIS_NOAH"

使用ncl脚本修改 wrfrst_d01_2017-01-14_00:00:00

load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"

begin
  ; 打开现有的文件（假设文件已经存在）
  f = addfile("wrfrst_d01_2017-01-14_00:00:00", "w")

  if (.not.ismissing(f)) then
    ; 定义和设置全局属性
    globalAtt = True
    globalAtt@MMINLU = "MODIFIED_IGBP_MODIS_NOAH" ;;;改成wrfinput中MMINLU的属性

    ; 使用 fileattdef 函数修改文件属性
    fileattdef(f, globalAtt)
  else
    print("Error: The file 'wrfrst_d01_2017-01-14_00:00:00' cannot be opened.")
  end if
end

之后正常运行wrf.exe即可

namelist.input中修改

start_year               = 2017,
start_month              = 01,
start_day                = 14,
start_hour               = 00,
start_minute             = 00,
start_second             = 00,
...............................
restart                  = .true.,