WRF WPS : namelist 学习笔记

news2025/2/24 13:22:31

WPS

& share

采用ARW方式进行模拟,除了ARW还有NMM,不过科研上常用ARW:

wrf_core = ‘ARW’

最大的嵌套层数为3层,初学者一般是从一层开始逐步加多:

max_dom = 3
# max_dom = 2

设置模式开始和结束 的时间,从左到右依次是第一层第二层和第三层;刚安装好WRF时候,里边只有两层,如果您需要转三层的模式,需要自行添加)

start_date = 2012-07-21_00:00:00,2012-07-21_00:00:00,2012-07-21_00:00:00,
start_date = 2012-07-21_00:00:00,2012-07-21_00:00:00

end_date = 2012-07-22_00:00:00,2012-07-22_00:00:00,2012-07-22_00:00:00,
end_date = 2012-07-22_00:00:00,2012-07-22_00:00:00


# 设置输出文件的位置。

opt_output_from_geogrid_path='/WPS/'

前处理的两次分析时间之间的时间间隔(在WRF运行过程中会发现具体的情况,以秒为单位。也就是模式实际输入数据的时间间隔,一般为输入边界条件的文件的时间间隔。没有默认值;输入数据的时间间隔)

interval_seconds = 21600 

这里是一个输出格式的选项,默认为2;当输入1时为binary,后缀是.int文件,当输入2时为NetCDF,后缀为.nc文件,当输入3时为GRIB1,后缀为.gr1文件。

io_form_geogrid = 2, 

& geogrid

嵌套区域的母区域的标号。模式本身没有母区域,因为此项第一列一般设为1,第二列必须等于1。总列数必须等于NUM_DOMAINS

parent_id = 0, 1, 2,
# parent_id = 1 , 1

嵌套时母网格与子网格之间的比例,没有默认值,例如这里就是1:3:9的比例,最外层为27公里,中间层为9公里,最里层为3公里

parent_grid_ratio = 1, 3, 3,
# parent_grid_ratio = 1, 3 ,

水平网格各层之间的开始的隔点数,从网格左下角开始

i_parent_start = 1, 31, 35,

各层网格的结束点,因为起始点从1开始,所以也是总隔点数

e_we = 101, 109, 109, 

e_sn = 91, 109, 109,

(区域所对应选择的地表面静态数据)

geog_data_res = 10m,2m,2m, 

(xy方向最外层网格的格距,与输入数据保持格距一致)

dx = 27000,
dy = 27000,

选择的投影方式lambert, polar, mercator, and lat-lon

map_proj = lambert, 

中心纬度,也就是模式定最外层所确定的中心位置

ref_lat = 39.9, 

中心经度,也就是模式定最外层所确定的中心位置

ref_lon = 116.3,  

投影的标准纬

truelat1 = 30.0, 
truelat2 = 60.0,

标准经度,实际中图形以这条轴展开

stand_lon = 116.3,  

设置地表静态数据的路径

geog_data_path = '/WPS/geog'

文件输出路径:

opt_geogrid_tbl_path = '/domain/output/'

& ungrib

输出文件的形式,有==MM5,WPS和SI==三种形式

out_format = WPS, 

输出文件的前缀,默认为FILE形式

prefix = 'FILE', 

& metgrid

fg_name = 'FILE'
io_form_metgrid = 2,

输出路径:

opt_output_from_metgrid_path = '/metgrid/output',

WRF

根据namelist.wps来设置namelist.input,对于domain,time要和namelist.wps保持一致。

&time_control

时间是累计加起来的。

 运行时间天、时、分、秒
 run_days                            = 2,
 run_hours                           = 48,
 run_minutes                         = 0,
 run_seconds                         = 0,

# 起始/结束 年份、月份、日期、小时、分钟、秒钟 这里的每一列需要根据你选取的嵌套层进行设置,与WPS一致

 start_year                          = 2016,2016,   
 start_month                         = 10,10,
 start_day                           = 06,06,
 start_hour                          = 00,00,
 end_year                            = 2016,2016,
 end_month                           = 10,10,
 end_day                             = 08,08,
 end_hour                            = 00,00,
 #前处理程序的两次分析时间的时间间隔
 interval_seconds                    = 21600,
 #嵌套初始场输入选项。嵌套时,指定嵌套网格是否使用不同的初始场文件
 input_from_file                     = .true.,.true.,
 #指定模式结果输出的时间间隔,以分钟为单位。默认单位是分钟,如果在变量名后面加入'_s'或者 '_h'(例如, history_interval_s),可以将它换成秒或者小时。这个值在不同区域可以不同。
 history_interval                    = 180, 180, 
 #每一个wrfout文件包含多少帧文件,每个wrfout里面记录多少个时次的数据
 frames_per_outfile                  = 1,1,
 #是否重新启动,仅在有restart文件(wrfrst*,之前运行时生成) 时使用。
 restart                             = .false.,
 #重启时间间隔
 restart_interval                    = 1440,
 #指定模式断电重启输出的格式,2为nc文件;1为二进制文件,4为phd5格式,5为grib1格式,
 io_form_history                     = 2
 #restart文件(wrfrst_d0*)将会写入的格式
 io_form_restart                     = 2
 #输入文件 (met_em_d0*) 的格式.
 io_form_input                       = 2
 # 边界文件 (wrfbdy_d01)写入格式.
 io_form_boundary                    = 2
 #输出运行记录的级别,数值越大,输出的信息越详细
 debug_level  =0

/

需要与wps中一致

 start_year                          = 2016,2016,   
 start_month                         = 10,10,
 start_day                           = 06,06,
 start_hour                          = 00,00,
 end_year                            = 2016,2016,
 end_month                           = 10,10,
 end_day                             = 08,08,
 end_hour                            = 00,00,
 interval_seconds                    = 21600,
 input_from_file                     = .true.,.true.,
 history_interval                    = 180, 180, 

&domains

#积分的时间步长,整数型,单位为秒
 time_step                           = 150,
 最大区域层数
 max_dom                             = 2,
 e_we                                = 91,91,
 e_sn                                = 100,100,
 e_vert                              = 45,45,
 s_vert                              = 1,1,
# 模型中要使用的气压上限(单位:Pa)。在传入的WPS数据中必须要有这个层次的数据。默认值和建议值为5000 Pa,不建议气压顶使用低于值(在大气中)。
 p_top_requested                     = 5000,
# 传入(来自WPS输入数据)垂直层次的数量。这由您使用的输入数据决定。
 num_metgrid_levels                  = 32,
 num_metgrid_soil_levels             = 4,
# x和y方向的格距,与WPS不同(在WPS中,嵌套值是基于parent_grid_ratio计算的),必须指定这些值,但必须根据parent_grid_ratio保持准确。
 dx                                  = 27000,
 dy                                  = 27000,
# 计算区域的编号,一般从1开始
 grid_id                             = 1,     1,
# 嵌套网格的上一级网格(母网格)的编号,一般从0开始
 parent_id                           = 1,     3,
 
 i_parent_start                      = 1,         30,
 j_parent_start                      = 1,         30,
# 相对于域父级的嵌套比率。对于最粗糙的域,父网格比率应设置为1。这些值应与您在namelist.wps中使用的值相同。
 parent_grid_ratio                   = 1,     3,  
# 这就是模型为嵌套域定义时间步长的方式。对于最粗糙的域,该值应设置为1。
 parent_time_step_ratio              = 1,     3,    
# 这将确定在使用嵌套时是否使用feedback,并应用于是单向嵌套运行还是双向嵌套运行。
 feedback                            = 1,
# 确定如果启用反馈(feedback=1),是否对嵌套区域中的父域使用平滑。有三种选择:
	0=无平滑
	1=1-2-1平滑
	2=平滑去光滑。通常建议将此选项设置为0。
 smooth_option                       = 0
 /

栅格空间示例:

s_we: x方向(西-东)方向的起始格点值,通常为1
e_we:x方向(西-东)方向的终止格点值,通常为x方向上的格点数
s_sn:y方向的起始点值
e_sn:y方向的终止格点值
s_vert:z方向的起始格点值
e_vert:z方向的终止格点值,即全垂直eta层的总层数,该层数在各嵌套网格中必须一致



e_we = 201, 304,
e_sn = 188, 289,

image.png

e_sn=度数(经纬度)×111(km)➗分辨率(36km)

e_we同理

i_parent_start=1,最外层的度数×111km除以这一层的分辨率(12km)

grid_id:域名。

最粗糙的栅格应设置为1。在下面的示例中,我们展示了一个更复杂的示例,其中我们有4个域(如图所示编号-注意,我们可以对它们进行不同的编号)。对于此示例:

image.png

parent_id:嵌套父级的域号。

最粗糙的域应设置为1。在下面的示例中,我们展示了一个更复杂的示例,其中我们有4个域(如图所示编号-注意,我们可以对它们进行不同的编号)。例如

grid_id = 1, 2, 3, 4
parent_id = 1, 1, 2, 1

可以看到d02和d04将d01作为父级,而d03将d02作为父级:

i_parent_start
j_parent_start:

父域中嵌套左下角的x和y坐标。对于最粗糙的域,应为i和j指定值1。这些值与您在WPS中指定的值相同。

i_parent_start = 1, 31,
j_parent_start = 1, 17,

image.png

需要与wps中一致

 e_we                                = 91,91,
 e_sn                                = 100,100,
 e_vert                              = 45,45,
 s_vert                              = 1,1,
  grid_id                             = 1,     1,
 parent_id                           = 1,     3,
 i_parent_start                      = 1,         30,
 j_parent_start                      = 1,         30,
 parent_grid_ratio                   = 1,     3,    
 parent_time_step_ratio              = 1,     3, 

&physics

 physics_suite                       = 'CONUS'
 mp_physics                          = -1,    -1,    -1,
 cu_physics                          = -1,    -1,     0,
 ra_lw_physics                       = -1,    -1,    -1,
 ra_sw_physics                       = -1,    -1,    -1,
 bl_pbl_physics                      = -1,    -1,    -1,
 sf_sfclay_physics                   = -1,    -1,    -1,
 sf_surface_physics                  = -1,    -1,    -1,
 radt                                = 30,    30,    30,
 bldt                                = 0,     0,     0,
 cudt                                = 5,     5,     5,
 icloud                              = 1,
 num_land_cat                        = 21,
 sf_urban_physics                    = 0,     0,     0,
 /

&fdda

 /

&dynamics

 hybrid_opt                          = 2, 
 w_damping                           = 0,
 diff_opt                            = 1,      1,      1,
 km_opt                              = 4,      4,      4,
 diff_6th_opt                        = 0,      0,      0,
 diff_6th_factor                     = 0.12,   0.12,   0.12,
 base_temp                           = 290.
 damp_opt                            = 3,
 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,     
 gwd_opt                             = 1,
 /

&bdy_control

 spec_bdy_width                      = 5,
 specified                           = .true.
 /

&grib2

&namelist_quilt

nio_tasks_per_group = 0,
nio_groups = 1,
 /

./real.exe 运行成功会出现:

image.png

常用linux 命令

cd 命令

cd directory 改变工作目录
cd .. 退到上一层目录

② ls 命令

-a 显示目录下所有子目录与文件(包括隐藏文件)
-l 显示文件的详细信息

③ cp命令

cp –r default OnlineTut
-r 递归复制源目录下所有的子目录和文件

④ vi 命令

vi file_name 开始编辑或者创建一个文件
编辑命令: <Esc> 模式切换
x 删除光标所在文字
dd 删除光标所在行
a 在光标后新增文字
i 在光标前新增文字
o 在光标下方新增一行
O 在光标上方新增一行
:wq 以原档案名保存并退出
:q! 不保存文档强制退出

⑤ rm命令

rm file_name 删除文件
rm –r deirectory_name 递归删除全部目录和子目录

⑥ mkdir命令

mkdir directory_name 创建新目录

⑦ rmdir 命令

rmdir directory_name 删除空目录

⑧ gunzip命令

gunzip xxxxx.tar.gz 解压缩。

⑨ tar命令

tar –xvf xxxxx.tar 解压文件
x 从档案文件中释放文件
v 详细报告tar处理的文件信息
f 使用档案文件或设备,这个选项通常是必选的
可以和上面的⑧合写成 tar xvfz xxxxx.tar.gz

⑫ pwd 命令

pwd 显示出当前工作目录的绝对路径

⑬ Ctrl + C

强制放弃正在执行的任务

⑭ exit 命令

退出UNIX系统(包括退出SSH)

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.coloradmin.cn/o/1464765.html

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈,一经查实,立即删除!

相关文章

轻松掌握opencv的8种图像变换

文章目录 opencv的8种图像变换1. 图像放大、缩小2. 图像平移3. 图像旋转4. 图像仿射变换5. 图像裁剪6. 图像的位运算&#xff08;AND, OR, XOR&#xff09;7. 图像的分离和融合8. 图像的颜色空间 opencv的8种图像变换 1. 图像放大、缩小 我们先看下原图 import cv2 import ma…

文献速递:GAN医学影像合成--基于生成对抗网络的肺部图像分类的多域医学图像翻译生成

文献速递&#xff1a;GAN医学影像合成–基于生成对抗网络的肺部图像分类的多域医学图像翻译生成 01 文献速递介绍 在2019年底&#xff0c;一种称为2019冠状病毒病&#xff08;COVID-19&#xff09;的新型冠状病毒肺炎出现&#xff0c;迅速成为全球性大流行。感染COVID-19可以…

本机防攻击简介

定义 在网络中&#xff0c;存在着大量针对CPU&#xff08;Central Processing Unit&#xff09;的恶意攻击报文以及需要正常上送CPU的各类报文。针对CPU的恶意攻击报文会导致CPU长时间繁忙的处理攻击报文&#xff0c;从而引发其他业务的中断甚至系统的中断&#xff1b;大量正常…

IOS和Android系统架构

IOS的系统架构 iOS的为Objective-C和Swift&#xff0c;Objective-C的优势是效率高但比较“唯一”。 响应顺序&#xff1a;Touch--Media--Service--Core架构 分为四个层次&#xff1a;核心操作系统层&#xff08;core OS layer&#xff09;、核心服务层&#xff08;Core Serv…

STM32Cubemx TB6612直流电机驱动

一、TB6612FNG TB6612是一个支持双电机的驱动模块&#xff0c;支持PWM调速。PWMA、AIN1、AIN2 为一组控制引脚&#xff0c;PWMA 为 PWM 速度控制引脚&#xff0c;AIN1、AIN2 为方向控制引脚&#xff1b;PWMB、BIN1、BIN2 为一组控制引脚&#xff0c;PWMB 为 PWM 速度控制引脚&…

一种基于动态水位值的Flink调度优化算法(flink1.5以前),等同于实现flink的Credit-based反压原理

优化flink反压 说明1 flink反压介绍1.1 介绍1.2 大数据系统反压现状1.4 flink task与task之间的反压1.5 netty水位机制作用分析 2 反压优化算法3 重点&#xff01; 但是 可但是 flink1.5以后的反压过程。4 flink反压问题的查找瓶颈办法 说明 首先说明&#xff0c;偶然看了个论…

短剧小程序系统,重塑视频观看体验的科技革命

随着科技的飞速发展&#xff0c;人们对于数字化内容的消费需求也在不断增长。在这个大背景下&#xff0c;短剧小程序作为一种新型的视频观看方式&#xff0c;正逐渐受到大众的青睐。本文将探讨短剧小程序的发展背景、特点以及市场前景&#xff0c;分析其在重塑视频观看体验方面…

哪个牌子的护眼台灯比较好用?纯干货护眼台灯品牌推荐

有些家长陪孩子写作业的时候发现他们总是在揉眼睛&#xff0c;学习时间久了还会用力眨眼睛。其实无论是白天还是晚上&#xff0c;孩子在家学习&#xff0c;看书&#xff0c;搭积木等&#xff0c;如果灯光不给力&#xff0c;一定要用台灯来给孩子补光&#xff0c;避免因为光线环…

Linux系列讲解 —— 【Vim编辑器】在Ubuntu18.04中安装新版Vim

平时用的电脑系统是Ubuntu18.04&#xff0c;使用apt安装VIM的默认版本是8.0。如果想要安装新版的Vim编辑器&#xff0c;只能下载Vim源码后进行编译安装。 目录 1. 下载Vim源码2. 编译3. 安装4. 遇到的问题4.1 打开vim后&#xff0c;文本开头有乱码现象。4.2 在Vim编辑器中&…

测试环境搭建整套大数据系统(三:搭建集群zookeeper,hdfs,mapreduce,yarn,hive)

一&#xff1a;搭建zk https://blog.csdn.net/weixin_43446246/article/details/123327143 二&#xff1a;搭建hadoop&#xff0c;yarn&#xff0c;mapreduce。 1. 安装hadoop。 sudo tar -zxvf hadoop-3.2.4.tar.gz -C /opt2. 修改java配置路径。 cd /opt/hadoop-3.2.4/etc…

【Spring】SpringBoot 热部署

目 录 一.添加热部署框架支持二.Settings 开启项目自动编译三.设置运行项目中的热部署( idea 2021.2版本)四.使用 debug 方式运行项目代码示例&#xff1a; 一.添加热部署框架支持 <dependency><groupId>org.springframework.boot</groupId><artifactId&…

适合中国人体质的低成本创业项目,抖音小店抓住小钱到大钱的之路

大家好&#xff0c;我是电商花花。 人啊&#xff0c;这一辈子想要赚钱&#xff0c;想要脱贫致富&#xff0c;小钱靠勤&#xff0c;中财靠运&#xff0c;大富靠命。 我还依稀记得母亲说的一句话&#xff0c;小钱不勤不聚&#xff0c;还谈何赚大钱&#xff0c;所有的大钱也都是…

Unity3d Mesh篇(二)— 创建Unity Logo平面

文章目录 前言一、Mesh组成二、使用步骤GetVertices方法GetNormal方法GetTriangles方法OnDrawGizmos方法 三、效果四、总结 前言 本篇将使用C#脚本实现在Unity中创建平面&#xff0c;并通过调整顶点、UV坐标和三角形来生成Unity Logo 的效果。 一、Mesh组成 顶点&#xff08;…

了解红帽认证,看这篇就够了!

红帽公司成立于1993年&#xff0c;是全球首家收入超10亿美元的开源公司&#xff0c;总部位于美国&#xff0c;分支机构遍布全球。 红帽公司作为全球领先的开源和Linux系统提供商&#xff0c;其产品已被业界广泛认可并使用&#xff0c;尤其是RHEL系统在业内拥有超高的Linux系统…

css3d制作正方体

使用css3d技术 &#xff0c;制作一个可以动态动画的正方体模型 效果图&#xff1a; 代码如下&#xff1a; <!DOCTYPE html> <html> <head><style>/* 设置高度宽度100%并且左右居中、上下居中 */html,body {width: 100%;height: 100%;display: flex…

C#串口 Modbus通讯工具类

一、安装Modbus包 二、创建modbushelper类 1、打开串口 public bool IfCOMOpend; //用于实例内的COM口的状态 public SerialPort OpenedCOM;//用于手动输入的COM转成SERIAL PORT /// <summary> /// 打开串口 /// </summary> /// <param name="COMname&quo…

2-22算法习题总结

贪心问题 凌乱的yyy / 线段覆盖 题目背景 快 noip 了&#xff0c;yyy 很紧张&#xff01; 题目描述 现在各大 oj 上有 n n n 个比赛&#xff0c;每个比赛的开始、结束的时间点是知道的。 yyy 认为&#xff0c;参加越多的比赛&#xff0c;noip 就能考的越好&#xff08;假…

Day04-课后练习以及参考答案(流程控制语句_循环结构)

文章目录 巩固题1、5个一行输出1-100之间的偶数2、趣味折纸3、实现输出如下任一个数字三角形4、计算这一天是这一年的第几天 拔高题1、计算这一天是这一年的第几天2、计算这一天是在打鱼还是晒网3、打印『X』对称图形4、打印空心菱形 巩固题 1、5个一行输出1-100之间的偶数 &…

【更新】ARCGIS之成片区开发方案报备坐标txt格式批量导出工具(定制开发版)

序言 之前开发的成片区开发方案报备格式是按湖北省的标准定制的&#xff0c;目前&#xff0c;自然资源部又有了新的格式要求&#xff0c;现在新增国标版的成片区开发方案报备格式导出。 之前版本软件详见&#xff1a;软件介绍 一、软件简介 本软件是基于arcgis二次开发的工具&…

YOLOv9 | 利用YOLOv9训练自己的数据集 -> 推理、验证(源码解读 + 手撕结构图)

一、本文介绍 本文给大家带来的是全新的SOTA模型YOLOv9的基础使用教程&#xff0c;需要注意的是YOLOv9发布时间为2024年2月21日&#xff0c;截至最近的日期也没有过去几天&#xff0c;从其实验结果上来看&#xff0c;其效果无论是精度和参数量都要大于过去的一些实时检测模型&…