3D(三位变分同化)模式的安装和运行（三）

最新推荐文章于 2023-08-08 11:39:36 发布
最新推荐文章于 2023-08-08 11:39:36 发布
阅读量1w 收藏 19
点赞数 3
CC 4.0 BY-SA版权
文章标签： output input io nested icloud file
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/bernardbai/article/details/7332959
本文详细介绍了WRF模型的三个主要配置文件：namelist.wps、namelist.input及namelist.obsproc的具体参数设置。从地理网格配置到气象物理过程方案选择，再到观测数据预处理，覆盖了从预处理到后处理的整个流程。

摘要生成于 C知道 ，由 DeepSeek-R1 满血版支持， 前往体验 >

三、各个步骤的namelist的配置：

1、WPS/namelist.wps配置：

&share
wrf_core = 'ARW',                                                                                                           #  输出格式定义：WRF还是NMM，'ARW‘表示是WRF。
max_dom = 3,                                                                                                                 #  域/嵌套数
start_date = '2012-03-01_00:00:00','2012-03-01_00:00:00','2012-03-01_00:00:00',     #  时间窗口，开始时间
end_date   = '2012-03-02_00:00:00','2012-03-02_00:00:00','2012-03-02_00:00:00',    #  时间窗口，结束时间

interval_seconds = 10800                                                                                               # 时间间隔（例如美国GFS的3小时间隔的，3*60*60秒）
io_form_geogrid = 2,                                                                                                       # WRF I/O API格式。。。。
/

&geogrid                                                                                                                         # 定义所有模型的大小和位置。
parent_id         =   1,   1,    2,
parent_grid_ratio =   1,   3,    3,
i_parent_start    =   1,  2,    2,
j_parent_start    =   1,  2,    2,

# 北京地区
e_we              =  10,  16,    19,                                                                                        # 开始、结束格点
e_sn              =  10,  16,    19,         
geog_data_res     = '10m','2m',  '30s',
dx = 40000,                                                                                                                    # GFS 0.5 deg 的网格4万米
dy = 40000,
map_proj = 'lambert',                                                                                                     # GFS的映射
ref_lat   =  39.92,                                                                                                           # 中心纬度
ref_lon   = 116.39,                                                                                                         #中心经度
truelat1  =  30.0, 
truelat2  =  60.0,
stand_lon = 116.39,                                                                                                       # 与ref_lon取值相同。

#
geog_data_path = '/software/geog'                                                                               # geog地理信息存放目录
/

&ungrib
out_format = 'WPS',                                                                                                      # ungrib输出格式
prefix = 'FILE',                                                                                                               # 输出文件的前缀
/

&metgrid                                                                                                                       #输出文件的前缀
fg_name = 'FILE'
io_form_metgrid = 2,                                                                                                    # 输出 netcdf格式文件。
/


2、WRFV3/test/em_real/namelist.input配置

&time_control
run_days                            = 0,
run_hours = 24,
run_minutes                         = 0,
run_seconds                         = 0,
start_year                          = 2012,2012,2012,
start_month                         = 03,03,03,
start_day                           = 01,01,01,
start_hour                          = 00,00,00,
start_minute                        = 00,00,00,
start_second                        = 00,00,00,
end_year                            = 2012,2012,2012,
end_month                           = 01,01,01,
end_day                             = 06,06,06,
end_hour                            = 00,00,00,
end_minute                          = 00,00,00,
end_second                          = 00,00,00,
interval_seconds                    = 10800
input_from_file                     = .true.,.true.,.true.,
history_interval                    = 180, 15, 15,  
frames_per_outfile                  = 1000, 1000, 1000,
restart                             = .false.,
restart_interval                    = 5000,
io_form_history                     = 2
io_form_restart                     = 2
io_form_input                       = 2
io_form_boundary                    = 2
debug_level                         = 0
/

&domains
time_step                           = 180,
time_step_fract_num                 = 0,
time_step_fract_den                 = 1,
max_dom                             = 3,
e_we                                =  10,  16,    19,
e_sn                                =  10,  16,    19,
e_vert                              = 28,    28,    28,      #分层层数，0层-27层。
p_top_requested                     = 5000,
num_metgrid_levels                  = 27,
num_metgrid_soil_levels             = 4,
dx                                  = 40000, 13333.3,  4444.4,    #降尺度 除以三。
dy                                  = 40000, 13333.3,  4444.4,
grid_id                             = 1,     2,     3,
parent_id                           = 0,     1,     2,
i_parent_start                      = 1,  2,    2,
j_parent_start                      = 1,  2,    2,
parent_grid_ratio                   = 1,     3,     3,
parent_time_step_ratio              = 1,     3,     3,
feedback                            = 1,
smooth_option                       = 0
/

&physics
mp_physics                          = 3,     3,     3, 
ra_lw_physics                       = 1,     1,     1, 
ra_sw_physics                       = 1,     1,     1,  
radt                                = 10,    10,    10,
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,     1,  
cudt                                = 5,     5,     5,  
isfflx                              = 1,
ifsnow                              = 0,
icloud                              = 1,
surface_input_source                = 1,
num_soil_layers                     = 4,
sf_urban_physics                    = 0,     0,     0,
maxiens                             = 1,
maxens                              = 3,
maxens2                             = 3,
maxens3                             = 16,
ensdim                              = 144,
/

&fdda
/

&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,
/

3、/WRFDA/var/obsproc/namelist.obsproc配置：

&record1
 obs_gts_filename = '/WRF/WRF_3D/GFS_WRF/little_r/obs.data',
 obs_err_filename = 'obserr.txt',
/

&record2
 time_window_min  = '2012-03-01_00:00:00',
 time_analysis    = '2012-03-01_01:30:00',
 time_window_max  = '2012-03-01_03:00:00',
/

&record3
 max_number_of_obs        = 400000,
 fatal_if_exceed_max_obs  = .TRUE.,
/

&record4
 qc_test_vert_consistency = .TRUE.,
 qc_test_convective_adj   = .TRUE.,
 qc_test_above_lid        = .TRUE.,
 remove_above_lid         = .TRUE.,
 domain_check_h           = .true.,
 Thining_SATOB            = false,
 Thining_SSMI             = false,
 Thining_QSCAT            = false,
/

&record5
 print_gts_read           = .TRUE.,
 print_gpspw_read         = .TRUE.,
 print_recoverp           = .TRUE.,
 print_duplicate_loc      = .TRUE.,
 print_duplicate_time     = .TRUE.,
 print_recoverh           = .TRUE.,
 print_qc_vert            = .TRUE.,
 print_qc_conv            = .TRUE.,
 print_qc_lid             = .TRUE.,
 print_uncomplete         = .TRUE.,
/

&record6
 ptop =  1000.0,
 base_pres       = 100000.0,
 base_temp       = 290.0,
 base_lapse      = 50.0,
 base_strat_temp = 215.0,
 base_tropo_pres = 20000.0
/

&record7
 IPROJ = 1,
 PHIC  = 39.92,
 XLONC = 116.39,
 TRUELAT1= 30.0,
 TRUELAT2= 60.0,
 MOAD_CEN_LAT = 39.92,
 STANDARD_LON = 116.39,
/

&record8
 IDD    =   1,
 MAXNES =   1,
 NESTIX =       100,
 NESTJX =       100,
 DIS    =   4.4444,
 NUMC   =   2,
 NESTI  =       6,
 NESTJ  =       6,
 /

&record9
 PREPBUFR_OUTPUT_FILENAME = 'prepbufr_output_filename',
 PREPBUFR_TABLE_FILENAME = 'prepbufr_table_filename',
 OUTPUT_OB_FORMAT = 2
 use_for          = '3DVAR',
 num_slots_past   = 1,
 num_slots_ahead  = 1,
 write_synop = .true.,
 write_ship  = .true.,
 write_metar = .true.,
 write_buoy  = .true.,
 write_pilot = .true.,
 write_sound = .true.,
 write_amdar = .true.,
 write_satem = .true.,
 write_satob = .true.,
 write_airep = .true.,
 write_gpspw = .true.,
 write_gpsztd= .true.,
 write_gpsref= .true.,
 write_gpseph= .true.,
 write_ssmt1 = .true.,
 write_ssmt2 = .true.,
 write_ssmi  = .true.,
 write_tovs  = .true.,
 write_qscat = .true.,
 write_profl = .true.,
 write_bogus = .true.,
 write_airs  = .true.,
 /

 

                                                                                                                                                                                               未完待续

 

 

 

 

 

bernard
关注
【WRF实操】三层/四层嵌套网格设计-以重庆市主城区为例
WW、forever的博客
05-28 86
【WRF实操】三层/四层嵌套网格设计-以重庆市主城区为例
探寻数理物理，体验数学与物理的融合
最新发布
AI天才研究院
07-08 255
数理物理就像一首宏伟的交响曲，数学提供了乐谱和乐器，物理则赋予其旋律和情感。本文将带您踏上一场探索数理物理奇妙世界的旅程，揭示数学与物理如何相互滋养、相互成就。我们将从历史渊源出发，解析核心概念，深入探讨数学工具如何帮助我们理解物理宇宙的奥秘，并通过具体案例展示这种融合如何推动了近代科学的革命性突破。无论您是学生、研究者还是科学爱好者，这篇文章都将为您打开一扇通往数理物理世界的大门，展示人类理性思维的最伟大成就之一。
netcdf2littler:将NetCDF文件转换为Little-R格式的应用程序。 Little-R格式是WRF数据同化系统（WRFDA）预处理器（obsproc）可接受的输入格式。 当前，此应用程序仅支持天气气象站的转换
05-13
NetCDF2LittleR 描述 一个FORTRAN应用程序，用于将NetCDF文件转换为Little-R格式。 Little-R格式是WRF数据同化系统（WRFDA）预处理器（obsproc）可接受的输入格式。 当前，此应用程序仅支持对天气气象站的转换。 安装 NetCDF2LittleR使用cmake构建系统。 NetCDF2LittleR需要构建外部源，例如，通过创建构建目录。 在构建过程中，将创建一个名为convert-littler的二进制文件。 用法 二进制convert-littler需要一个Fortran名称列表作为参数，例如： netcdf2littler --namelist netcdf2littler.namelist 如果未提供任何参数， netcdf2littler尝试使用当前工作目录中的文件netcdf2littler.namelist 。 输入名称列
气象模式 | WRFDA资料同化系统的发展
Earth_Sci_Club的博客
01-08 1279
WRFDA的应用范围广泛非常广泛，从公里尺度的区域和中尺度建模到用于数千公里尺度的全球建模，可用于为WRF模式预报提供初始条件，验证观测与模式之间的差异，进行观测系统的设计、监测和评估等。
3D(三位变分同化)模式的安装和运行（二）
bernardbai的中尺度WRF模型专栏
03-01 7130
二、3D运行:         众所周知,风能是一种清洁的能源.但是风具有间歇性和不连续性的特点,在单位的时间内很难准确的预测风速和风向,更难准确地预测出风能的功率。以下将用 WRF-3D(三维变分同化)模式讲述其模拟的运行过程。         采用WRF模式做模拟的结果往往更贴近数值天气预报，可能是它只对数值天气预报做了一些降低尺度的处理的缘故，WRF可以将时间尺度降到10分钟,与测风设备
3D(三维变分同化)模式的安装和运行(一)
bernardbai的中尺度WRF模型专栏
03-01 5588
一、3D运行环境搭建: 1.首先安装WRF模式的WRFV3.2和WPS3.2(参见《WRF安装过程》)。 2.安装WRFDAV3.2.tar.gz. 环境变量设置：    增加环境变量   export BUFR=1 # tar -zxvf ./WRFDAV3.2.tar.gz #cd WRFDA #./configure wrfda 出现提示如下： -------------
运行WRFDA-3DVAR and 4DVAR
Sun_wenxuan的博客
08-08 972
把wrfbdy_d01文件复制到WRFDA工作文件夹下，因为 da_update_bc.exe.会把该文件覆盖，因此不能简单的链接过来！同样把WRFDA/var/test/update_bc/parame.in 复制到工作文件夹下。运行WRFDA/var/build/da_update_bc.exe。更新 WRF 横向边界条件。更改parame.in。
WRFDA-OBSPROC
Sun_wenxuan的博客
08-04 832
参考WRFDA官网2017年夏季指导课程第一步。
[更新中]文献索引
LeoGIS1987的博客
12-03 9607
2010 2010-01 [1]多个粗差定位的抗掩盖型Bayes方法.pdf [6]模糊神经网络在变形分析与预报中的应用研究.pdf [9]GPS_GLONASS组合精密单点定位研究.pdf [13]基于星间距离变化的动态双向时间同步算法.pdf [17]QZSS对GPS区域性能增强分析.pdf [21]GPS模糊度降相关LLL算法的一种改进.pdf [25]2008年汶川地震前的形变异常及机理解...
WRFDA环境搭建速学：新手的全面教程
本文全面介绍了WRFDA（Weather Research and Forecasting Data Assimilation）环境的搭建和应用。首先，概述了WRFDA环境搭建的必要性，并详细讲解了搭建过程中的软件需求、系统环境准备以及基础环境配置。接着，通过...
WRFDA教程：3DVar同化卫星辐射率资料
islandowner2017的博客
10-15 4399
由于卫星观测的优势，目前各国的资料同化业务化系统中，卫星资料的数量基本是占据第一位。本篇参考WRFDA官网和用户手册，梳理了WRFDA三维变分同化的卫星辐射率资料的运行流程，供初学者参考。
ndvar：GAME的最佳插值（OI），3D和4D变异数据同化
02-16
最佳插值（OI），3D和4D变异数据同化。 文件资料 科学推导 同化技术的科学派生可以在我的理论气象学教科书中找到（德语）： 。 文献资料 该代码的文档可在子目录doc中找到。 建筑 依存关系 netcdf库（Ubuntu：sudo apt-get libnetcdf-dev） CMake（Ubuntu：sudo apt-get install cmake） MPICH（Ubuntu：sudo apt-get install mpich） bzip2（Ubuntu：sudo apt-get install bzip2） python3（Ubuntu：sudo apt-get安装python3） 为了发展 Valgrind（Ubuntu：sudo apt-get install valgrind，用于进行检查） 下载 git clone https://github.com/AUN
wrfda安装和问题分析
09-21
要做资料同化必须安装的wrfda，刚从网上找的好东西，供大家使用！
WRF-DA代码编译与安装(一）——依赖库的编译与安装
胖胖雕的GIS博客
04-02 5234
由于笔者的研究需要用到数据同化技术，所以开始学习WRF相关内容（主要是WRF-DA模块）。本篇内容WRF依赖库的编译与安装
WRFDA资料同化详解
CCfz566的博客
03-08 913
资料同化新方法的快速发展，气象常规资料、卫星遥感观测和大气环境等多种资料日益增加，为资料同化的有效应用奠定了坚实的科学基础，也导致许多新的复杂科学问题，增加了实际应用的难度。为有效提升广大科研、业务人员的资料同化理论基础和应用能力以及研究水平，提高我国气象、海洋和大气环境等领域的数值预报业务的质量和实际预报水平。马老师（副教授）来自高等院校，科研机构的数值预报资料同化高级专家，拥有丰富的科研和实际应用经验，长期从事数值预报、资料同化等国家级科研项目研究与教学工作。1、WRF DA同化系统的配置。
【案例教程】气象数值预报WRF-DA资料同化系统理论、运行与变分、混合同化新方法技术
weixin_49857990的博客
12-12 998
掌握和理解资料同化的理论； 熟悉Linux环境下的配置和应用； 掌握背景误差协方差、物理约束关系的调优和更新； 掌握新观测资料同化的关键技术，如观测算子的构造、伴随程序的发展等； 熟悉同化结果的分析方法； 掌握Bash、GRADS、NCL等后处理工具进行结果处理及分析; 掌握在气象、环境、生态、水文等领域应用技术的拓展和方法; 熟练掌握WRF ETKF-3DVAR DA混合同化系统的应用;
WRF/WRFDA中尺度气象模式和资料同化的理论与应用
walrus_course的博客
07-31 959
本课程围绕当前应用最广泛的WRF/WRFDA，讲解相关的理论基础、运行流程与相关的资料处理步骤，并结合实际应用，介绍数值模式和资料同化应用中的重点与难点。（5）常规气象要素的误差诊断（U10\V10\ T2，Python）（2）卫星资料同化（AMUSA、HIRS、MHS、葵花8等）（5）同化微波辐射资料的注意事项（通道选择、偏差订正等）（4）站点的降水评分（TS\ETS评分，Python）（1）3DVar个例（常规观测、卫星微波辐射等）（1）绘制二维填色图（风温压等，Python）
WRF-DA资料同化系统理论、运行与变分、混合同化新方法
weixin_46747075的博客
04-26 1033
以WRF DA变分资料同化系统为核心，针对实际应用中的重点和难点问题，开展资料同化理论和方法的讲解，并结合实际天气个例，致力于培养和提高学习人员的理论水平和实际应用能力。
三维变分同化
05-30
### 三维变分同化（3D-Var）的概念与实现方法 三维变分同化（3D-Var）是一种用于优化数值天气预报模型初始条件的技术。其核心思想是通过最小化一个成本函数来获得最优的模式初值[^1]。成本函数通常由两部分组成：一是模型预报数据与观测数据之间的差异，二是模型预报数据与先验背景场之间的差异。目标是找到一个初始条件，使得成本函数达到最小值。 在实际应用中，3D-Var方法擅长同化卫星和雷达等非常规资料，并通过构造标量目标函数来实现分析场的优化[^3]。以下是3D-Var方法的具体实现步骤： #### 1. 构造目标函数 目标函数 \( J(x) \) 可以表示为： \[ J(x) = J_b(x) + J_o(x) \] 其中： - \( J_b(x) \) 表示背景场误差项，衡量初始场与背景场之间的差异。 - \( J_o(x) \) 表示观测误差项，衡量初始场与观测数据之间的差异。 #### 2. 最小化目标函数 通过对目标函数 \( J(x) \) 进行极小化，可以得到最优的模式初值。这一过程通常使用梯度下降法或其他优化算法来实现[^1]。 #### 3. 数据处理与生成 在实际运行过程中，需要对观测数据进行预处理。例如，在/WRFDA/var/obsproc目录下运行 obsproc 程序，依据定义处理观测数据文件，生成如 `obs_gts_2012-03-01:00:00:00.3DVAR` 文件[^2]。 #### 4. 配置与运行 在配置阶段，需要设置网格参数、地理数据分辨率以及其他相关参数。例如，对于北京地区的配置可能包括以下内容[^4]： ```python e_we = [10, 16, 19] # 开始、结束格点 e_sn = [10, 16, 19] geog_data_res = ['10m', '2m', '30s'] # 地理数据分辨率 dx = 40000 # 网格间距（米） dy = 40000 map_proj = 'lambert' # 映射投影类型 ref_lat = 39.92 # 中心纬度 ref_lon = 116.39 # 中心经度 truelat1 = 30.0 # 标准纬线1 truelat2 = 60.0 # 标准纬线2 stand_lon = 116.39 # 标准经线 ``` 此外，还需要链接必要的表格文件和可执行程序，例如 LANDUSE.TBL 和 da_wrfvar.exe[^5]。 #### 5. 应用场景 3D-Var 方法广泛应用于气象领域，特别是在数值天气预报中。它能够有效结合背景场和观测数据，生成更精确的初始场，从而提高预报精度。尽管 3D-Var 在时间维度上仅考虑单一时刻的数据，但它为更复杂的 4D-Var 方法奠定了基础[^3]。 ---
评论 4
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值