相关文章推荐
礼貌的针织衫  ·  Java8函数之旅 (八) - ...·  7 月前    · 
仗义的大蒜  ·  TS 中 never,void,unknow ...·  1 年前    · 
拉风的领结  ·  JDBC02 ...·  1 年前    · 
风流倜傥的充电器  ·  line.strip().split ...·  1 年前    · 
  • 在学院上使用ERA5数据驱动WRF模式,使用一层嵌套网格
  • 数据时间范围为:2004年6月19日–2004年6月30日
  • 经纬度范围:南北纬45°左右,经度全部包括。中心经度为180°,中心纬度为0°
  • 输出的数据时间间隔为:每6小时一次
  • 分辨率为:36000x36000


步骤还是老样子,先进行 WPS 前处理,再运行WRF模式


WPS前处理



1、链接运行数据的类型变量


cd命令 调整到WPS目录下,先链接到ERA5驱动的Vtable,运行完,会发现在当前目录下得到 Vtable文件


ln -sf ungrib/Variable_Tables/Vtable.ECMWF Vtable


121824f256e341b58dcce89e7ee2dabc.png


2、链接下载的数据


链接你下载的ERA5数据,根据你数据放置的文件夹进行链接


./link_grib.csh  ./ERA5_data/*


以上两步完成后,修改namelist.wps中的参数,参数如下设置:


&share
 wrf_core = 'ARW',
 start_date = '2004-06_19_00:00:00','2004-06_24_00:00:00','2004-06_24_00:00:00',  
 end_date = '2004-06_30_18:00:00','2004-06_25_00:00:00','2004-06_25_00:00:00',
 interval_seconds = 21600,   #输出文件的时间间隔
 max_dom = 1,               #嵌套层数
 io_form_geogrid = 2,        #输出文件类型,2表示nc文件
 opt_output_from_geogrid_path = '/Models/WRF4.3.1/WPS/',   #输出文件路径
&geogrid
 parent_id         =   1,      1,      2,  #嵌套区域名称
 parent_grid_ratio =   1,      3,      3,  #嵌套区域比例,主要与分辨率相关
 i_parent_start    =   1,     21,    27,   #起始点
 j_parent_start    =   1,     21,    26,   #起始点
 e_we                     = 1110,      574,    1471,    #根据你的经纬度间隔计算的而来
 e_sn                     = 276,      370,      889,
 geog_data_res     = 'default','default','default',
 dx = 36000, #分辨率
 dy = 36000,
 map_proj  = 'mercator', #投影方式
 ref_lat   =  0.00, #中心纬度
 ref_lon   = 180.00, #中心经度
 truelat1  =  10.0,
 truelat2  =  20.0,
 stand_lon =180.00,  #以哪个经度进行展开对称
 geog_data_path = '/Pakages/WPS_GEOG/', #存放的地理信息数据的路径
&ungrib
 out_format = 'WPS',
 prefix = '/Models/WRF4.3.1/WPS/ERA',   #输出路径以及前缀
&metgrid
 fg_name = 'ERA',
 io_form_metgrid = 2, 
 opt_output_from_metgrid_path = '/Models/WRF4.3.1/WPS/',
/


3、ungrid


修改完成后,输入命令运行 ./ungrid.exe ,如下运行成功,而且发现生成得到了前缀文件


c4eaf54e3b33496fbce4d84aca87d9dd.png


9d7df446a94e4316b9b66e6d16a82218.png


4、geogrid


同理,继续输入命令运行 ./geogrid.exe ,同样得到如下文件:


dfc84bc7662f4f64bd09af8d3da1aa47.png

8306deb8fbf444aeb8c3871b62fd3131.png


5、megrid


同理,继续输入命令运行 ./metgrid.exe ,同样得到如下文件:


ad1993f62524403da133e19e86bfcebd.png

fab699287dd34eb795443698bcb7255d.png


得到生产的met*.nc文件,这里WPS前处理结束。下面进行 WRF 模式处理即可


WRF



一样,先使用cd命令调整目录到 WRF/run/ 下,然后链接刚刚在wps前处理下生成的 met*nc文件

然后,根据你的namelist.wps的时间设置,时间间隔,嵌套层数,等信息,调整你的namelist.input,与其保持一致。这里将需要保持一致的参数大致列一下:


其中,本次实验我修改的 namelist.input 如下所示:


real.exe


上述namelist.input修改好后,运行 ./real.exe ,生成得到以下文件:



可以看一下rsl.out.0000,显示success就成功了


wrf.exe


下面就是提交作业到集群上跑起来了,需要先创建一个 .sh 脚本,下面给出一个作业提交的例子:



通过命令:

qsub job.sh ,进行提交

qstat -f ,查看运行状态

qdel 作业ID ” ,停止作业。

耐心等待即可


使用shell在Linux系统下下载cmip6文件出现报错:No ESG Credentials found in /Users/daniele/.esg/credentials.pem python--海温、OLR数据分布做显著性检验,绘制空间分布并打点 python--使用convolve 对二维数据进行平滑 python--测试使用不同的方式计算位涡平流项的差异 cdo (Climate Data Operators ) 常用命令介绍:(查看文件信息、多文件合并、数据裁剪、数据插值、数据计算。。)
cdo (Climate Data Operators ) 常用命令介绍:(查看文件信息、多文件合并、数据裁剪、数据插值、数据计算。。)
西门子S7-200 SMART如何使用状态图表,如何创建、监视、强制、趋势显示
上篇文章中我们学习了S7-200 SMART系统块的组态,本篇我们来介绍在编程软件STEP7-Micro/WIN SMART中如何使用状态图表,以及如何创建、监视、强制、趋势显示。在STEP7-Micro/WIN SMART与PLC之间成功建立通信,并且将程序下载到PLC后,就可以监控和调试程序了。程序状态监控可以监视程序的运行情况,但是如果需要监控的变量较多,不能在程序编辑器中同时显示的时候就需要使用状态图表监控。接下来我们来介绍在STEP7-Micro/WIN SMART如何使用状态图表监控和调试程序。
西门子S7-200 SMART位逻辑指令概述及应用举例
本篇文章我来带领大家学习西门子S7-200 SMART的位逻辑指令。位逻辑指令是PLC编程中最基本、使用最频繁的指令,按不同的功能和用途具有不同的形式,总的来说可以分为下述几大类:标准位逻辑指令、置位/复位指令、立即位逻辑指令、其他位逻辑指令。
西门子S7-200 SMART PID回路控制,如何创建状态图表,下载并测试程序
上篇文章中我们学习了调用向导生成的子程序编程来实现西门子S7-200 SMART PID回路控制,今天我们来学习如何创建状态图表并利用状态图表来测试程序。首先我们来看一下如何创建状态图表。
使用shell在Linux系统下下载cmip6文件出现报错:No ESG Credentials found in /Users/daniele/.esg/credentials.pem python--海温、OLR数据分布做显著性检验,绘制空间分布并打点 python--使用convolve 对二维数据进行平滑 python--测试使用不同的方式计算位涡平流项的差异