Advertisement

【WRF-Urban】考虑空间分布城市参数的初始WRF-Urban版本

阅读量:

目录

  • WRF模型运行

    • WRF编译
    • WPS预处理
    • WRF运行
  • 参考

本博客介绍WRF-Urban考虑空间分布城市参数的初始版本,出自GitHub-WRF-distributed-urban
在这里插入图片描述
这是东京工业大学 GUC 实验室使用的 WRF v3.3.1 的修改版本。在此版本中,城市土地利用使用城市形态参数的 2D 地图进行建模。每小时变化的人为热图也可以输入到模型中。
在这里插入图片描述

WRF模型运行

WRF编译

参考WRF原始教程,编译此WRF版本。详细的WRF编译过程可参见另一博客-【WRF教程第二期】WRF编译全过程:以4.5版本为例

WPS预处理

1、下载城市参数(urban parameters)的geogrid文件
下载地址-Present and future 1 km resolution global population density and urban morphological parameters
在这里插入图片描述
2、(可选)下载人为热排放数据
下载地址-Index of /database/AHE/AH4GUC
在这里插入图片描述
3、编译modify_geo_em.f90文件。
需要NetCDF Fortran库。此脚本用于修改geo_em文件,以便可以模拟最主要土地类别和第二主要土地类别。

复制代码
    gfortran modify_geo_em.f90 -o modify_geo_em.exe $(nf-config --fflags) $(nf-config --flibs)
    
    
      
    

4、按照如下方式编辑 WPS 的 GEOGRID.TBL。根据您的目录结构更改 rel_path。
请注意,AHE24H 中的 rel_path 需要根据月份进行更改。

复制代码
    ===============================
    name = DISP # Displacement height (m)
        priority = 1
        dest_type = continuous
        fill_missing = 0
        masked = water
        interp_option =     30s:average_gcell(4.0)+four_pt+average_4pt
        interp_option = default:average_gcell(4.0)+four_pt+average_4pt
        rel_path = default:urban_params/2010/d
    ===============================
    name = Z0_GRD # Roughness length for momentum (m)
        priority = 1
        dest_type = continuous
        fill_missing = 0
        masked = water
        interp_option =     30s:average_gcell(4.0)+four_pt+average_4pt
        interp_option = default:average_gcell(4.0)+four_pt+average_4pt
        rel_path = default:urban_params/2010/z_0
    ===============================
    name = AVE_BH_GRD # Average building height (m)
        priority = 1
        dest_type = continuous
        fill_missing = 0
        masked = water
        interp_option =     30s:average_gcell(4.0)+four_pt+average_4pt
        interp_option = default:average_gcell(4.0)+four_pt+average_4pt
        rel_path = default:urban_params/2010/H_avg
    ===============================
    name = FAI_GRD # Frontal area index (dimensionless)
        priority = 1
        dest_type = continuous
        fill_missing = -1
        masked = water
        interp_option =     30s:average_gcell(4.0)+four_pt+average_4pt
        interp_option = default:average_gcell(4.0)+four_pt+average_4pt
        rel_path = default:urban_params/2010/lambda_f
    ===============================
    name = PAI_GRD # Plan area index (dimensionless)
        priority = 1
        dest_type = continuous
        masked=water
        fill_missing = -1
        interp_option =     30s:average_gcell(4.0)+four_pt+average_4pt
        interp_option = default:average_gcell(4.0)+four_pt+average_4pt
        rel_path = default:urban_params/2010/lambda_p
    ===============================
    name = AHE24H # Anthropogenic heat (W/m^2)
        priority = 1
        dest_type = continuous
        fill_missing = 0
        interp_option =     30s:average_gcell(4.0)+four_pt+average_4pt
        interp_option = default:average_gcell(4.0)+four_pt+average_4pt
        rel_path = default:AHE/01_24 # 01 means January, 24 means 24 hours
        z_dim_name = utc_hour
    ==============================
    
    
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
      
    

5、运行 geogrid.exe 后,在包含 geo_em.d??.nc 文件的目录中,执行

复制代码
    ./modify_geo_em.exe
    
    
      
    

WRF运行

1、修改namelist.input
在 namelist.input 中,修改 domains 部分如下。

复制代码
    &domains
    # Your options (e.g., time_step, e_we, e_sn)
    offset_zd = 1 # Add zero-plane displacement height to topography.
    /
    
    
      
      
      
      
    

2、照常运行 WRF。

参考

全部评论 (0)

还没有任何评论哟~