WRF模式上机手册

1、WRF模式上机手册WRF模式上机手册一.安装1登陆系统连接服务器：telnet 00（如果是用客户端软件，则直接用客户端软件进行登陆）输入用户名：*输入密码：*创建自己的用户目录（如huangq）：mkdir huangq进到用户自己的目录（如huangq）：cd huangq2编译安装WRF模式主体1)获取源程序包(获取源程序代码可从WRF的官方网站下载,/doc/4e5e9bc0580216fc710afddf.html/wrf/users/download/get_sources.html2)cp /public1/Mo

2、del/WRFV3.7.TAR.gz ./3)解压源程序压缩包tar xzvf WRFV3.7.TAR4)进入释放后的源程序目录cd WRFV35)设置环境变量NETCDFexport NETCDF= /public/software/mathlib/netcdf/4.3.0/intel/6)配置编译环境./configure出现如下的选择列表：checking for perl5. nochecking for perl. found /usr/bin/perl (perl)Will use NETCDF in dir: /public/software/mathlib/netcdf/4.3

3、.0/intel/PHDF5 not set in environment. Will configure WRF for use without.Will use time to report timing informationIf you REALL Y want Grib2 output from WRF, modify the arch/Config_new.pl script.Right now you are not getting the Jasper lib, from the environment, compiled into WRF.-Please select fro

4、m among the following Linux x86_64 options:1. (serial)2. (smpar)3. (dmpar)4. (dm+sm) PGI (pgf90/gcc)5. (serial)6. (smpar)7. (dmpar)8. (dm+sm) PGI (pgf90/pgcc): SGI MPT9. (serial) 10. (smpar) 11. (dmpar) 12. (dm+sm) PGI (pgf90/gcc): PGI accelerator13. (serial) 14. (smpar) 15. (dmpar) 16. (dm+sm) INTE

5、L (ifort/icc)17. (dm+sm) INTEL (ifort/icc): XeonPhi (MIC architecture)18. (serial) 19. (smpar) 20. (dmpar) 21. (dm+sm) INTEL (ifort/icc): Xeon (SNBwith AVX mods)22. (serial) 23. (smpar) 24. (dmpar) 25. (dm+sm) INTEL (ifort/icc): SGI MPT26. (serial) 27. (smpar) 28. (dmpar) 29. (dm+sm) INTEL (ifort/ic

6、c): IBM POE30. (serial) 31. (dmpar) PA THSCALE (pathf90/pathcc)32. (serial) 33. (smpar) 34. (dmpar) 35. (dm+sm) GNU (gfortran/gcc)36. (serial) 37. (smpar) 38. (dmpar) 39. (dm+sm) IBM (xlf90_r/cc_r)40. (serial) 41. (smpar) 42. (dmpar) 43. (dm+sm) PGI (ftn/gcc): Cray XC CLE44. (serial) 45. (smpar) 46.

7、 (dmpar) 47. (dm+sm) CRAY CCE (ftn/gcc): Cray XEand XC48. (serial) 49. (smpar) 50. (dmpar) 51. (dm+sm) INTEL (ftn/icc): Cray XC52. (serial) 53. (smpar) 54. (dmpar) 55. (dm+sm) PGI (pgf90/pgcc)56. (serial) 57. (smpar) 58. (dmpar) 59. (dm+sm) PGI (pgf90/gcc): -f90=pgf9060. (serial) 61. (smpar) 62. (dm

8、par) 63. (dm+sm) PGI (pgf90/pgcc): -f90=pgf90Enter selection 1-63 : 15 (建议选择16)7)编译模式主体./compile em_real编译成功后，在main目录下有real.exe和wrf.exe。如果没有，则说明编译失败。3编译WPS1)获取源程序包(获取源程序代码可从WRF的官方网站下载)cp /public1/Model/WPSV3.7.TAR.gz ./2)解压源程序压缩包tar xzvf WPSV3.7.TAR.gz3)进入释放后的源程序目录cd WPS4)设置环境变量NETCDFexport NETCDF=

9、/public/software/mathlib/netcdf/4.3.0/intel/5)配置编译环境./configure出现如下的选择列表：-Will use NETCDF in dir: /public/software/mathlib/netcdf/4.3.0/intel/Found Jasper environment variables for GRIB2 support.$JASPERLIB = /public/software/mathlib/jasper/lib$JASPERINC = /public/software/mathlib/jasper/include-Plea

10、se select from among the following supported platforms.1. Linux x86_64, gfortran (serial)2. Linux x86_64, gfortran (serial_NO_GRIB2)3. Linux x86_64, gfortran (dmpar)4. Linux x86_64, gfortran (dmpar_NO_GRIB2)5. Linux x86_64, PGI compiler (serial)6. Linux x86_64, PGI compiler (serial_NO_GRIB2)7. Linux

11、 x86_64, PGI compiler (dmpar)8. Linux x86_64, PGI compiler (dmpar_NO_GRIB2)9. Linux x86_64, PGI compiler, SGI MPT (serial)10. Linux x86_64, PGI compiler, SGI MPT (serial_NO_GRIB2)11. Linux x86_64, PGI compiler, SGI MPT (dmpar)12. Linux x86_64, PGI compiler, SGI MPT (dmpar_NO_GRIB2)13. Linux x86_64,

12、IA64 and Opteron (serial)14. Linux x86_64, IA64 and Opteron (serial_NO_GRIB2)15. Linux x86_64, IA64 and Opteron (dmpar)16. Linux x86_64, IA64 and Opteron (dmpar_NO_GRIB2)17. Linux x86_64, Intel compiler (serial)18. Linux x86_64, Intel compiler (serial_NO_GRIB2)19. Linux x86_64, Intel compiler (dmpar

13、)20. Linux x86_64, Intel compiler (dmpar_NO_GRIB2)21. Linux x86_64, Intel compiler, SGI MPT (serial)22. Linux x86_64, Intel compiler, SGI MPT (serial_NO_GRIB2)23. Linux x86_64, Intel compiler, SGI MPT (dmpar)24. Linux x86_64, Intel compiler, SGI MPT (dmpar_NO_GRIB2)25. Linux x86_64, Intel compiler,

14、IBM POE (serial)26. Linux x86_64, Intel compiler, IBM POE (serial_NO_GRIB2)27. Linux x86_64, Intel compiler, IBM POE (dmpar)28. Linux x86_64, Intel compiler, IBM POE (dmpar_NO_GRIB2)29. Linux x86_64 g95 compiler (serial)30. Linux x86_64 g95 compiler (serial_NO_GRIB2)31. Linux x86_64 g95 compiler (dm

15、par)32. Linux x86_64 g95 compiler (dmpar_NO_GRIB2)33. Cray XE/XC CLE/Linux x86_64, Cray compiler (serial)34. Cray XE/XC CLE/Linux x86_64, Cray compiler (serial_NO_GRIB2)35. Cray XE/XC CLE/Linux x86_64, Cray compiler (dmpar)36. Cray XE/XC CLE/Linux x86_64, Cray compiler (dmpar_NO_GRIB2)37. Cray XC CL

16、E/Linux x86_64, Intel compiler (serial)38. Cray XC CLE/Linux x86_64, Intel compiler (serial_NO_GRIB2)39. Cray XC CLE/Linux x86_64, Intel compiler (dmpar)40. Cray XC CLE/Linux x86_64, Intel compiler (dmpar_NO_GRIB2)6)Enter selection 1-40 : (建议选17)7)编译WPS./compile wps编译成功后，在WPS的主目录下有geogrid.exe、ungrib

17、.exe和metgrid.exe。如果没有，则说明编译失败。8)编译WPS的辅助工具compile util编译成功后，在util目录下有g1print.exe、mod_levs.exe、plotgrids.exe、avg_tsfc.exe、g2print.exe、plotfmt.exe和rd_intermediate.exe。如果没有，则说明编译失败。配置编译环境是选择了不支持GRIB2的选项，所以g2print.exe是没有产生。二.运行WPS，进行数据前处理1运行geogrid.exe，生成模式格点和处理静态数据1)获取静态资料用户可以从WRF的官方网站获取WRF模式运行所需的静态资料。

18、高分辨率的静态数据压缩包有371M，解压后大约有10G，低分辩率的静态数据包有11.2M，解压后大约有264M。在/public1/WRF-DATA/geog目录中有已经解压后的静态数据，我们可以做一个连接后直接使用。cd /public/home/zcdmn/wangln sf /public1/WRF-DATA/geog geog到geog目录下查看一下：cd geogls2)编辑参数列表文件namelist.wps回到WPS的主目录，进行编辑操作。&sharewrf_core = ARW,max_dom = 1,io_form_geogrid = 2,opt_output_from_ge

19、ogrid_path = ./,debug_level = 0/&geogridparent_id = 1,parent_grid_ratio = 1,i_parent_start = 1,j_parent_start = 1,s_we = 1,e_we = 74,s_sn = 1,e_sn = 61,geog_data_res = 10m,2m,dx = 60000,dy = 60000,map_proj = lambert,ref_lat = 40.0ref_lon = 102truelat1 = 30.0,truelat2 = 60.0,stand_lon =102.geog_data_

20、path = /public1/WRF-DATA/geogopt_geogrid_tbl_path = geogrid/提示：这里可以用wps的辅助工具plotgrids程序来调整自己定义的区域参数。具体操作是：a)运行plotgrids.exe程序util/plotgrids.exeb)运行图形查看软件idt/usr/local/bin/ idt gmetac)弹出窗口里，点击“”键。3)创建运行geogrid.exe程序的作业卡vi geogrid.cmd创建后的geogrid.cmd的文件内容格式如下：# run the program geogrid.exe of WPS# job_t

21、ype = parallel# class = normald# comment=WRF# initialdir = /dgpfs/fs2/user/WRF/WPS# executable = geogrid.exe# input = /dev/null# output =geogrid.out# error =geogrid.err# node = 1# tasks_per_node = 8# node_usage = shared# network.MPI = css0,shared,US# notification = never# wall_clock_limit = 00:10:00

22、# queue4)提交作业卡geogrid.cmdllsubmit geogrid.cmd运行成功之后，在当前目录下会产生一个静态数据文件geo_em.d01.nc。2运行ungrib.exe，解码GRIB数据1)获取GRIB码数据文件用户可以从网络中获取全球模式GRIB格式的气象场资料。在IBM机器上，可以到以下目录下进行连接：/dgpfs/fs2/wrf_xp/grib/ln sf /dgpfs/fs2/wrf_xp/grib/ grib2)根据数据来源模式具体情况，选择与此模式相应的变量列表文件。这里以GFS模式的气象要素场为例：ln sf ungrib/Variable_Tables/

23、Vtable.GFS Vtable3)连接GRIB码数据文件成GRIB.？格式l ink_grib.csh grib/gfs4)编辑参数列表文件namelist.wps&sharewrf_core = ARW,max_dom = 1,start_date = 2021-08-16_12:00:00,end_date = 2021-08-16_18:00:00,interval_seconds = 21600io_form_geogrid = 2,opt_output_from_geogrid_path = ./,debug_level = 0/&ungribout_format = WPS/提

24、示：查看GRIB码数据的头信息可以用g1print.exe命令：util/g1print.exe GRIB.AAA5)创建运行ungrib.exe程序的作业卡vi ungrib.cmd创建后的ungrib.cmd的文件内容格式如下：# run the program ungrib.exe of WPS# job_type = serial# class = interactive# initialdir = /dgpfs/fs2/wrftc?/denglt/WPS# executable = ungrib.exe# input = /dev/null# output =ungrib.out#

25、error =ungrib.err# notification = never# wall_clock_limit = 00:10:00# queue6)提交作业卡ungrib.cmdllsubmit ungrib.cmd提示：a)ungrib的结果文件“FILE:”可以用WPS的工具软件rd_intermediate.exe程序进行查看，以验证数据的正确性。查看时，只需执行rd_intermediate.exe文件，并在这个可执行文件后面加上要查看的数据文件名。例如：util/rd_intermediate.exe FILE:2021-07-30_12b)用WPS的工具软件plotfmt.e

26、xe可以绘制中间文件的图形。具体操作是：执行plotfmt.exe程序，并带上数据文件名作为输入参数。例如：util/plotfmt.exe FILE:2021-07-30_12运行图形查看软件idtidt gmeta弹出窗口里，点击“”键。3运行metgrid.exe，对解码后的数据进行水平插值1) 编辑参数列表文件namelist.wps&sharewrf_core = ARW,max_dom = 1,start_date = 2021-08-16_12:00:00,end_date = 2021-08-16_18:00:00,interval_seconds = 21600io_form

27、_geogrid = 2,opt_output_from_geogrid_path = ./,debug_level = 0/&metgridfg_name = ./FILEio_form_metgrid = 2,opt_output_from_metgrid_path = ./,opt_metgrid_tbl_path = metgrid/,/2) 创建运行metgrid.exe程序的作业卡vi metgrid.cmd创建后，metgrid.cmd文件的内容格式为：# run the program metgrid.exe of WPS# job_type = parallel# class

28、 = normald# initialdir = /dgpfs/fs2/wrftc?/denglt/WPS# executable = metgrid.exe# input = /dev/null# output =metgrid.out# error =metgrid.err# node = 1# tasks_per_node = 8# node_usage = shared# network.MPI = css0,shared,US# notification = never# wall_clock_limit = 00:10:00# queue3）提交作业卡metgrid.cmdllsu

29、bmit metgrid.cmd运行成功后，在opt_output_from_metgrid_path设置的路径下会产生met_em.d01.yyyy-mm-dd_hh:00:00.nc格式的文件。三.运行WRF 模式主体1运行初始化程序real1)编辑参数列表文件namelist.inputWRF模式有专门运行real.exe和wrf.exe程序的路径（WRFV3/run或者WRFV3/test/em_real）。我们转到运行路径下cd WRFV3/test/em_real/然后就可以编辑参数列表文件vi namelist.input2)连接real.exe的输入场数据real.exe的输入

30、场数据文件就是metgrid的输出文件，文件名格式为“met_em.d01.*”。在运行real.exe前，必须将这些文件拷贝或者连接到real.exe的运行路径(test/em_real) 下面。ln -sf ././WPS/met_em.d01.* .3)创建运行real.exe程序的作业卡vi real.exe提示：根据编译WRF模式时的选项情况，我们可知real.exe编译后是串行还是并行程序。根据real.exe的编译情况，写出相应的串行或者并行作业卡。4)提交作业./real.exe2运行WRF模式1）编辑参数列表文件namelist.input进入到WRF的运行路径（WRFV2/

31、run或者test/em_real）cd WRFV3/run编辑参数列表文件（一般是在real.exe运行完以后，wrf.exe的参数列表就不做任何变动）vi namelist.input2）创建运行wrf.exe程序的作业卡job-submission.sh#!/bin/bash#PBS -q high#PBS -l nodes=2:ppn=24#PBS -j oecd $PBS_O_WORKDIRcat $PBS_NODEFILE |uniq nodesNP=cat $PBS_NODEFILE | wc -lmpirun -np 48 -machinefile nodes wrf.exe运

32、行成功后，在real.exe的运行路径(test/em_real) 下面会产生wrfinput_d01和wrfbdy_d01格式的文件。提示：根据编译WRF模式时的选项情况，我们可知wrf.exe编译后是串行还是并行程序。根据wrf.exe的编译情况，写出相应的串行或者并行作业卡。提交作业卡wrf.cmdqsub job-submission.sh运行成功后，在wrf.exe的运行路径(test/em_real) 下面会产生wrfout_d01_*格式的文件。四.运行后处理程序1编译安装程序1）获取程序进入到自己的运行目录cd /dgpfs/fs2/wrftc?/denglt拷贝程序包cp /u/wrf_xp/src/wrfpost20210829.tar释放程序包tar xvf wrfpost20210829.tar2）配置编译环境进入到后处理程序目录cd wrfpost编辑编译规则文件makefile，进行编译环境配置vi makefile3）编译程序在src目录下，编译MICAPS处理程序make micaps在src目录下，编译GrADS处理程序make grads2处理MICAPS数据1)编辑参数配置文件namelist.input进入MICA