1、WRFV2.2.1,萬文龍（整理） 2008年11月10日 CUIT,目錄,1 預備知識 2 Linux操作系統(tǒng)介紹 3 模式簡介 4 編譯器安裝 5 WRF模塊編譯 6 WPS模塊編譯 7 模式運行 8 ARWPOST后處理編譯 9 WRFDomainWizard 10 相關源代碼及英文教程下載網(wǎng)業(yè),1 預備知識,常用linux、unix命令： 打開終端，可以輸入命令，如： cd directory 進入指定路徑的工作目錄 cd . 退到上一層目錄 ls 顯示目錄下所有子目錄與文件(包括隱藏文件) Vi file_name 瀏覽和編輯文件 gunzip file_name.tar.gz 解壓

2、文件（對于.tar.gz文件） tar xvf file_name.tar 解壓文件 （對于.tar文件） ./configure 連接 ./compile 編譯 ./install 安裝 Chmod +x file_name 文件權限設置,#echo $SHELL (查看所用計算機的shell) Bsh設置環(huán)境變量命令：例如： NETCDF=/usr/local;export NETCDF Csh設置環(huán)境變量命令：例如： setenv NETCDF /usr/local/netcdf,2 Linux操作系統(tǒng)介紹 （摘自于百度）,Linux操作系統(tǒng)核心最早是由芬蘭的Linus Torvalds

3、 1991年8月在芬蘭赫爾辛基大學上學時發(fā)布的那年Torvals 25歲，后來經(jīng)過眾多世界頂尖的軟件工程師的不斷修改和完善，Linux得以在全球普及開來，在服務器領域及個人桌面版得到越來越多的應用，在嵌入式開發(fā)方面更是具有其它操作系統(tǒng)無可比擬的優(yōu)勢，并以每年100%的用戶遞增數(shù)量顯示了Linux強大的力量。Linux的是一套免費的32位多人多工的操作系統(tǒng)，運行方式同UNIX系統(tǒng)很像，但Linux系統(tǒng)的穩(wěn)定性、多工能力與網(wǎng)絡功能已是許多商業(yè)操作系統(tǒng)無法比擬的，Linux還有一項最大的特色在于源代碼完全公開，在符合GNU GPL(General Public License)的原則下，任何人皆可自

4、由取得、散布、甚至修改源代碼。,與其它操作系統(tǒng)相比，Linux還具有以下特色：采用階層式目錄結(jié)構，文件歸類清楚、容易管理 支持多種文件系統(tǒng)，如Ext2FS,ISOFS以及Windows的文件系統(tǒng)FAT16,FAT32,NTFS等具有可移植性，系統(tǒng)核心只有小于10%的源代碼采用匯編語言編寫，其余均是采用C語言編寫，因此具備高度移植性可與其它的操作系統(tǒng)如Windows98/2000/xp等并存于同一臺計算機上,現(xiàn)在linux版本很多，應用比較廣泛的是redhat以及開源的fedora。本文模式的安裝是在redhat企業(yè)版5.0為例。 NOTES:經(jīng)陳功同學測試，fedora最新版本（指9版）在超線

5、程奔騰單機上要比redhat運算快1倍;PGI編譯使用f95比f90快25%。感謝陳功同學在模式編譯運算優(yōu)化上的試驗支持！,3 模式簡介,WRF（Weather Research and Forecasting Model）模式是由美國環(huán)境預測中心（NCEP）,美國國家大氣研究中心（NCAR）等美國科研機構中心著手開發(fā)的一種統(tǒng)一的中尺度天氣預報模式。 模式分為ARW(the Advanced Research WRF)和NMM(the Nonhydrostatic Mesoscale Model)兩種，即研究用和業(yè)務用兩種形式，本文介紹的是 ARW WRF。,3.1 模式特點,WRF模式系統(tǒng)具

6、有可移植，易維護，可擴充，高效率，方便等許多特點，將成為改進從云尺度到各種不同天氣尺度的重要天氣特征預報精度的工具。 WRF模式為完全可壓縮以及非靜力模式，采用F90語言編寫。水平方向采用Arakawa C(荒川C)網(wǎng)格點（重點考慮1-10km ），垂直方向則采用地形跟隨質(zhì)量坐標。WRF模式在時間積分方面采用三階或者四階的Runge-Kutta算法。WRF模式不僅可以用于真實天氣的個案模擬，也可以用其包含的模塊組作為基本物理過程探討的理論根據(jù)。此外，WRF模式還具有多重嵌套和方便的定位于不同地理位置的能力。,Arakawa-C 特點,運用Arakawa-C 交錯格點假設所有3維變量（U， V，

7、 和質(zhì)量）關于這些點是交錯格點。對于定義的非交錯格點，U格點向上交錯了0.5個格點，V格點向右交錯了0.5個格點，質(zhì)量網(wǎng)格分別向上向右交錯了0.5個格點。 為了便于說明，下面給出一個(XDIM，YDIM) = (4，4)的例子： (+) 為根據(jù)參數(shù)定義的點。(T)為由WRF預報模式提供和輸出的質(zhì)量變量的格點位置。(U)點為由WRF模式提供和輸出的U動量變量的格點位置。 (V)點為由WRF模式提供和輸出的V動量變量的格點位置。這樣，如果使用維數(shù)(XDIM， YDIM)，則模式輸出如下： (XDIM-1，YDIM-1)維的質(zhì)量變量 (XDIM，YDIM-1)維的U動量 (XDIM-1，YDIM)維

8、的V動量。,3.2 WRF模式流程,出處：Users Guide for Advanced Research WRF (ARW) Modeling System Version 2.2,4 編譯器安裝,Software requirements（軟件需求） - Fortran 90 or 95 and c compiler（模式編譯語言f90/f95）- perl 5.04 or better- If MPI and OpenMP compilation is desired, it requires MPI or OpenMP libraries- WRF I/O API supports

9、netCDF, PHD5, GriB 1 and GriB 2 formats, hence one of these libraries needs to be available on the computer where you compile and run WRF,安裝PGI 設置環(huán)境變量 安裝NETCDF,WRF模式系統(tǒng)是采用Fortran 90、Fortran 77及c+編譯語言進行編譯與測試的，所以，運行WRF必須先安裝編譯軟件， 目前運用的編譯軟件主要有PGI和INTEL。本文以PGI（6.0版本）編譯器、NETCDF3.6.1為例進行說明。,4.1 PGI6.0安裝,（1）

10、：把linux86-64文件復制到指定路徑下,如/usr/,用命令 # tar zxvf linux86-64.tar.gz 解壓文件 （2）：編譯并安裝 # cd /usr/linux86-64 進入安裝目錄 # ./install 執(zhí)行安裝Do you accept these terms? accept,declineacceptInstall the ACML? y/nnInstallation directory? /usr/pgi/usr/local/pgiCreate an evaluation license? y/nnDo you want the files in the

11、install directory to be read-only? y,nn （3）：用root用戶把license.dat（在使用期范圍的其他版本的license.dat也可以用）復制(cp)到/usr/local/pgi下$cp r /usr/linux86-64/license.dat /usr/local/pgi,4.2 環(huán)境變量設置,為了方便，我們可以將環(huán)境變量登錄到.bashrc（或者.bash_profile)里。下文將PGI以及NETCDF一起設置進去。在根目錄root下打開隱藏文件.bash_profile 或著在終端下輸入：#vi .bash_profile （一般我們修

12、改在/home/usr/cuit目錄里： #cd /home/usr/cuit） 出現(xiàn)如下內(nèi)容： # .bash_profile # Get the aliases and functions if -f /.bashrc ; then . /.bashrc,然后添加如下語句： # User specific environment and startup programs PATH=$PATH:$HOME/bin export PATH unset USERNAME,# set pgi # （設置PGI) export PGI=/usr/local/pgi export PATH=$PGI/

13、linux86/6.0/bin:$PGI/linux86/6.0/include:$PATH export MANPATH=$MANPATH:$PGI/linux86/6.0/man export LM_LICENSE_FILE=$PGI/license.dat # set netcdf # （設置NETCDF) export NETCDF=/usr/local/netcdf export PATH=$NETCDF/bin:$NETCDF/include:$NETCDF/lib:$NETCDF/man:$PATH export CC=/usr/bin/gcc export CXX=$PGI/l

14、inux86/6.0/bin/pgCC export FC=$PGI/linux86/6.0/bin/pgf90 export F90=$PGI/linux86/6.0/bin/pgf90,然后保存退出，用如下命令 #:wq! #source .bash_profile 保存完后，注銷一次。最后可以用命令whichpgf90來查看pgi是否安裝成功。 #which pgf90 /usr/local/pgi/linux86/6.0/bin/pgf90 #pgf90 pgf90-Warning-No files to process 顯示如上則pgi安裝成功,4.3 NETCDF3.6.1安裝,N

15、OTES:PGI安裝不成功，后面的NETCDF是編譯安裝不上的！ #cd /usr/local #mkdir netcdf #cd . 回到usr目錄下 #cd . 回到根目錄 #cd usr #tar zxvf netcdf.tar.gz #cd netcdf-3.6.1/src,#./configure prefix=/usr/local/netcdf 編譯并指示安裝路徑 #make check #make install 執(zhí)行安裝 此時NETCDF安裝完成，環(huán)境變量在前一個步驟設置好，注意前面的環(huán)境變量里的路徑必須與此安裝路徑相同。,5 WRF（主）模塊編譯安裝,首先在/home/use

16、r/cuit下解壓WRFV2.2.1.tar.gz文件： #tar zxvf WRFV2.2.1.tar.gz #ls #mv WRFV2.2.1 WRFV2 #cd WRFV2 下一步就要開始連接然后編譯WRF主模塊： #./configure(or：configure有的shell命令如此) 然后出現(xiàn)如下內(nèi)容：,checking for perl5. no checking for perl. found /usr/bin/perl (perl) Will use NETCDF in dir: /usr/local/netcdf-pgi PHDF5 not set in environme

17、nt. Will configure WRF for use without. - Please select from among the following supported platforms. 1. PC Linux i486 i586 i686,PGI compiler (Single-threaded, no nesting) 2. PC Linux i486 i586 i686, PGI compiler (single threaded, allows nesting using RSL without MPI) 3. PC Linux i486 i586 i686, PGI

18、 compiler SM-Parallel (OpenMP, no nesting) 4. PC Linux i486 i586 i686, PGI compiler SM-Parallel (OpenMP, allows nesting using RSL without MPI) 5. PC Linux i486 i586 i686, PGI compiler DM-Parallel (RSL, MPICH, Allows nesting) 6. PC Linux i486 i586 i686, PGI compiler DM-Parallel (RSL_LITE, MPICH, Allo

19、ws nesting) 7. Intel xeon i686 ia32 Xeon Linux, ifort compiler (single-threaded, no nesting) 8. Intel xeon i686 ia32 Xeon Linux, ifort compiler (single threaded, allows nesting using RSL without MPI) 9. Intel xeon i686 ia32 Xeon Linux, ifort compiler (OpenMP) 10. Intel xeon i686 ia32 Xeon Linux, ifo

20、rt compiler SM-Parallel (OpenMP, allows nesting using RSL without MPI) 11. Intel xeon i686 ia32 Xeon Linux, ifort+icc compiler DM-Parallel (RSL,MPICH, allows nesting) 12. Intel xeon i686 ia32 Xeon Linux, ifort+gcc compiler DM-Parallel (RSL,MPICH, allows nesting) 13. PC Linux i486 i586 i686, PGI comp

21、iler, ESMF (Single-threaded, ESMF coupling, no nesting) Enter selection1-13 :,我們選擇2，使用單機，允許嵌套，但不并行運算： Enter selection 1-13 : 2 - Configuration successful. To build the WRF, type: compile - #./compile 編譯安裝開始 Usage: compile wrf compile wrf in run dir (Note, no real.exe, ndown.exe or ideal.exe generate

22、d) test cases (see README_test_cases for details): compile em_b_wave compile em_grav2d_x compile em_hill2d_x compile em_quarter_ss compile em_real compile em_squall2d_x compile em_squall2d_y compile h help message,我們選擇real模擬，則輸入 #./compile em_real & compile.log 接下來就等待主模塊的編譯安裝 如果編譯安裝成功,在/home/user/WR

23、FV2/run路徑下生成real.exe和wrf.exe可執(zhí)行文件,在終端上顯示為藍色的即可(紅色表示沒有安裝成功，需重新編譯)。 若沒有編譯成功，則輸入clean命令再重新安裝（未成功，看compile.log里面的報錯） #clean a 再查找原因重新編譯安裝,6 WPS編譯安裝,首先在/home/user/cuit下解壓WPSV2.2.1.tar.gz文件： #tar zxvf WRFV2.2.1.tar.gz #ls #cd WPS 下一步就要開始連接然后編譯WPS模塊： #./configure(or：configure有的shell命令如此) 然后出現(xiàn)如下內(nèi)容：,Will use

24、 NETCDF in dir: /usr/local/netcdf - Please select from among the following supported platforms. 1. PC Linux i486 i586 i686, PGI compiler serial 2.PC Linux i486 i586 i686, PGI compiler serial, NO GRIB2 3. PC Linux i486 i586 i686, PGI compiler DM parallel 4. PC Linux i486 i586 i686, PGI compiler DM pa

25、rallel ,NO GRIB2 5. PC Linux i486 i586 i686, Intel compiler serial 6. PC Linux i486 i586 i686, Intel compiler serial, NO GRIB2 7. PC Linux i486 i586 i686, Intel compiler DM parallel 8. PC Linux i486 i586 i686, Intel compiler DM parallel , NO GRIB2 Enter selection 1-8 :,我們選擇2（不用grib2格式的資料，若做業(yè)務運行，即預報N

26、MM，而非ARW，必須要安裝可以使用grib2資料的，大部分業(yè)務運行要使用此類數(shù)據(jù)） Enter selection 1-8 : 2 - Configuration successful. To build the WPS, type: compile - 然后下一步就開始編譯安裝,#./compile & compile.log 安裝成功后，會在/wrf/WPS路徑下生成geogrid.exe，ungrib.exe，metgrid.exe 三個可執(zhí)行文件（藍色）。若沒有此三個文件，請查看ARWUsersGuideV2.pdf里面有詳細的trouble shooting。查看如下： #ls *

27、.exe geogrid.exe metgrid.exe ungrib.exe,GEOG安裝,到此WRF的前處理以及主模塊已經(jīng)安裝完成。下一步要運行還需要前處理模塊需要的地形資料文件geog，下一步解壓geog.tar.gz 首先在/home/user/cuit下解壓 #tar zxvf geog.tar.gz #ls #cd geog,#ls (出現(xiàn)一下內(nèi)容) albedo_ncep landuse_30s greenfrac landuse_5m Islope maxsnowalb landuse_10m soiltemp_1deg landuse_2m soiltype_bot_30s

28、如此地形文件解壓安裝完成，此步驟需要的時間比較長，因為此文件較大大約接近10G，請耐心的等待?。ㄓ嘘Pgeog的說明請參考英文教程，在此不解釋說明） 下一步我們可以真正的運行模式了,7 WRF模式運行,7.1 WPS運行,geogrid和ungrib 屬并列關系，運行不分先后。 geogrid 建立“靜態(tài)的”地面數(shù)據(jù)。 ungrib 解壓GRIB 氣象數(shù)據(jù)，并歸納成一個 intermediate 文件格式。 metgrid 把氣象數(shù)據(jù)水平插入模式領域內(nèi)。 metgrid的輸出文件將被用作WRFV2.2.1主模塊的輸入文件。,編輯namelist.wps 以2006年8月16日12時-8月16日1

29、8時NECP資料為背景場作為例子做一次模擬。 #cd /home/user/cuit/WPS #ls #vi namelist.wps (or gedit namelist.wps) 然后出現(xiàn)以下內(nèi)容,&share wrf_core = ARW, max_dom = 2, (最大嵌套數(shù)，2層) start_date = 2006-08-16_12:00:00,2006-08-16_12:00:00, end_date = 2006-08-16_18:00:00,2006-08-16_18:00:00, interval_seconds = 21600(前處理程序的兩次分析時間之間的時間間隔，以

30、秒為單位。也即模式的實時輸入數(shù)據(jù)的時間間隔，一般為輸入邊界條件的文件的時間間隔。) io_form_geogrid = 2, /,&geogrid parent_id = 1, 1, （嵌套區(qū)域的母區(qū)域的標號。注意MOAD 本身沒有母區(qū)域，因此PARENT_ID 的第一列總是設為1。第二列必須等于1?？偭袛?shù)必須等于NUM_DOMAINS） parent_grid_ratio = 1, 3, （嵌套時，母網(wǎng)格相對于嵌套網(wǎng)格的水平網(wǎng)格比例。在真實大氣方案中，此比例必須為奇數(shù)；在理想大氣方案中，如果將返饋選項feedback設置為0的話，則此比例也可以為偶數(shù)） i_parent_start = 1,

31、 31（嵌套網(wǎng)格的左下角（LLC）在上一級網(wǎng)格（母網(wǎng)格）中x方向的起始位置 ） j_parent_start = 1, 17（嵌套網(wǎng)格的左下角（LLC）在上一級網(wǎng)格（母網(wǎng)格）中y方向的起始位置） s_we = 1, 1,e_we = 74,112, s_sn = 1, 1, e_sn = 61,97, geog_data_res = 10m,2m, dx = 30000, dy = 30000, map_proj = lambert, ref_lat = 34.83 ref_lon = -81.03 truelat1 = 30.0, truelat2 = 60.0, stand_lon = -

32、98. geog_data_path = /home/user/cuit/geog opt_geogrid_tbl_path = geogrid/ /,e_we = 74, 112, （x方向(西-東方向)的終止格點值 (通常為x方向的格點數(shù))） e_sn = 61, 97, （y方向(南-北方向)的終止格點值 (通常為y方向的格點數(shù))） geog_data_res = 10m,2m, （區(qū)域?qū)x擇的地表面靜態(tài)數(shù)據(jù)） dx = 30000, （指定x方向的格距（單位為米）。在真實大氣方案中，此參數(shù)值必須與輸入數(shù)據(jù)中的x方向格距一致） dy = 30000, （指定y方向的格距（單位為米）。在

33、真實大氣方案中，此參數(shù)值必須與輸入數(shù)據(jù)中的y方向格距一致） map_proj = lambert, （地圖投影）polar - 極射投影 lambert - 蘭伯托等角投影（正割和正切） mercator - 麥卡托 ref_lat = 30.0, (中心緯度) ref_lon = 115.0, (中心經(jīng)度) truelat1 = 30.0, (真實緯度) truelat2 = 60.0, stand_lon = 115.0, (標準經(jīng)度) geog_data_path = /home/user/cuit/geot (放置地表面靜態(tài)數(shù)據(jù)路徑) /,&ungrib out_format = WP

34、S, prefix = FILE, / &metgrid fg_name = FILE io_form_metgrid = 2, / &mod_levs press_pa = 201300 , 200100 , 100000 , 95000 , 90000 , 85000 , 80000 , 75000 , 70000 , 65000 , 60000 , 55000 , 50000 , 45000 , 40000 , 35000 , 30000 , 25000 , 20000 , 15000 , 10000 , 5000 , 1000 /,計算公式,xdim = (uri-lli)*ratio

35、_to_parent+1 ydim = (urj-llj)*ratio_to_parent+1,修改完成后然后保存退出，以下為運行命令： #./geogrid.exe ! ! Successful completion of geogrid. ! ! #ls (如下圖) #ln sf ungrib/Variable_Tables/Vtable.GFS空格Vtable #ls #./link_grib.csh空格/home/user/cuit/data/*fnl #./ungrib.exe,! ! Successful completion of ungrib. ! ! #ls (如下圖) #.

36、/metgrid.exe ! ! Successful completion of metgrid.! ! #ls (如下圖),圖 .1,圖 .2,圖 .3,此過程前處理就處理完成了，下一步我們就可以初始化然后模式運行了 把wps里生成的met打頭的文件拷貝到wrfv2里面的run文件夾里： #cd /home/user/cuit/WPS #cp r met_em* /home/user/cuit/WRFV2/run 這樣就可以進行下一步了,7.2 WRF模塊運行,#cd /home/user/cuit/WRFV2/run #vi namelist.input (可以用gedit命令） 然后出

37、現(xiàn)以下內(nèi)容： &time_control run_days = 0, (運行時間（天) run_hours = 06, (時) run_minutes = 0, (分) run_seconds = 0, (秒) start_year = 2006, 2006, 2003, （起始年份） start_month = 08, 08, 07, （起始月份） start_day = 16, 16, 09, （起始日數(shù)） start_hour = 12, 12, 00, （起始小時） start_minute = 00, 00, 0 （起始分鐘） start_second = 00, 00, 00, （

38、起始秒數(shù)） end_year = 2006, 2006, 2003, end_month = 08, 08, 07, （對應的結(jié)束時間） end_day = 16, 16, 10, end_hour = 18, 18, 00, end_minute = 00, 00, 00, end_second = 00, 00, 00,interval_seconds = 21600 （前處理程序的兩次分析時間之間的時間間隔） input_from_file = .true.,.true.,.false., （嵌套初始場輸入選項。嵌套時，指定嵌套網(wǎng)格是否用不同的初始場文件） history_interva

39、l = 180, 60, 60, （指定模式結(jié)果輸出的時間間隔，以分鐘為單位） frames_per_outfile = 1000, 1000, 1000, restart = .false., （是否進行重行啟動 restart_interval = 1440, （重起時間間隔） io_form_history = 2 （ 2 = NetCDF ） io_form_restart = 2 （指定模式斷點重啟輸出的格式, 2為netCDF格式） io_form_input = 2 （2 = NetCDF） io_form_boundary = 2 （指定模式邊界條件數(shù)據(jù)的格式） 1 二進制格式

40、 2 NetCDF格式 4 PHD5格式 5 GRIB1格式,debug_level = 0 （此選項指定模式運行時的調(diào)試信息輸出等級。取值可為0,50,100,200,300，數(shù)值越大，調(diào)試信息輸出就越多，默認值為0） / &domains time_step = 180, （積分的時間步長，為整型數(shù)，單位為秒，在真實大氣中推薦值為dx公里數(shù)的6倍） time_step_fract_num = 0, （實數(shù)型時間步長的分子部分） time_step_fract_den = 1, （實數(shù)型時間步長的分母部分） max_dom = 2, （最大區(qū)域數(shù)） s_we = 1, 1, 1, （x方向(

41、西-東方向)的起始格點值 (通常為1) e_we = 74, 112, 94, （x方向(西-東方向)的終止格點值 (通常為x方向的格點數(shù)) s_sn = 1, 1, 1, (y方向(南-北方向)的起始格點值 (通常為1),e_sn = 61, 97, 91, (y方向(南-北方向)的終止格點值 (通常為y方向的格點數(shù)) s_vert = 1, 1, 1, (z方向(垂直方向)的起始格點值 ) e_vert = 28, 28, 28, (z方向(垂直方向)的終止格點值，即全垂直eta層的總層數(shù)。垂直層數(shù)在各嵌套網(wǎng)格中必須保持一致 ) dx = 30000, 10000, 3333, (指定x方

42、向的格距（單位為米）。在真實大氣方案中，此參數(shù)值必須與輸入數(shù)據(jù)中的x方向格距一致) dy = 30000, 10000, 3333, (指定y方向的格距（單位為米）。在真實大氣方案中，此參數(shù)值必須與輸入數(shù)據(jù)中的x方向格距一致) grid_id = 1, 2, 3, (計算區(qū)域的編號，一般是從1開始) parent_id = 0, 1, 2, (嵌套網(wǎng)格的上一級網(wǎng)格（母網(wǎng)格）的編號，一般是從0開始),i_parent_start = 0, 31, 30, (嵌套網(wǎng)格的左下角（LLC）在上一級網(wǎng)格（母網(wǎng)格）中x方向的起始位置) j_parent_start = 0, 17, 30, (嵌套網(wǎng)格的左

43、下角（LLC）在上一級網(wǎng)格（母網(wǎng)格）中y方向的起始位置) parent_grid_ratio = 1, 3, 3, (母網(wǎng)格相對于嵌套網(wǎng)格的水平網(wǎng)格比例) parent_time_step_ratio = 1, 3, 3, (嵌套時，母網(wǎng)格相對于嵌套網(wǎng)格的時間步長比例) feedback = 1, (嵌套時，嵌套網(wǎng)格向母網(wǎng)格得反饋作用。設置為0時，無反饋作用。而反饋作用也只有在母網(wǎng)格和子網(wǎng)格的網(wǎng)格比例(parent_grid_ratio)為奇數(shù)時才起作用) smooth_option = 0 (向上一級網(wǎng)格（母網(wǎng)格）反饋的平滑選項，只有設置了反饋選項為1時才起作用的) /,&physics m

44、p_physics = 3, 3, 3, 設置微物理過程方案，默認值為 0。 0 不采用微物理過程方案 1 Kessler 方案 (暖雨方案) 2 Lin 等的方案 (水汽、雨、雪、云水、冰、冰雹) 3 WSM 3類簡單冰方案 4 WSM 5類方案 5 Ferrier(new Eta)微物理方案(水汽、云水) 6 WSM 6類冰雹方案 8 新Thompson的冰雹方案 98 NCEP 3類簡冰方案 (水汽、云/冰和雨/雪) 99 NCEP 5類方案(水汽、雨、雪、云水和冰) ra_lw_physics = 1, 1, 1, 此選項指定長波輻射方案，默認值為 0。 0 不采用長波輻射方案 1 r

45、rtm 方案 99 GFDL (Eta) 長波方案 (semi-supported),ra_sw_physics = 1, 1, 1, 此選項指定短波輻射方案，默認值為 0。 0 不采用短波輻射方案 1 Dudhia 方案 2 Goddard 短波方案 99 GFDL (Eta) 短波方案 (semi-supported) radt = 10, 10, 10, （此參數(shù)指定調(diào)用輻散物理方案的時間間隔，默認值為 0, 單位為分鐘。建議與dx的公里數(shù)取同樣的值） sf_sfclay_physics = 1, 1, 1, 此選項指定近地面層(surface-layer)方案，默認值為 0。 0 不采

46、用近地面層方案 1 Monin-Obukhov 方案 2 MYJ Monin-Obukhov 方案 (僅用于MYJ 邊界層方案) sf_surface_physics = 1, 1, 1, 此選項指定陸面過程方案，默認值為 0。 0 不采用陸面過程方案 1 熱量擴散方案 2 Noah 陸面過程方案 3 RUC 陸面過程方案,bl_pbl_physics = 1, 1, 1, 此選項指定邊界層方案，默認值為 0 0 不采用邊界層方案 1 YSU 方案 2 Eta Mellor-Yamada-Janjic TKE(湍流動能) 方案 3 NCEP Global Forecast System方案 9

47、9 MRF 方案 bldt = 0, 0, 0, (此參數(shù)指定調(diào)用邊界層物理方案的時間間隔，默認值為 0，單位為分鐘。此參數(shù)指定調(diào)用邊界層物理方案的時間間隔，默認值為 0，單位為分鐘。0 (推薦值)表示每一個時間步長都調(diào)用邊界層物理方案) cu_physics = 1, 1, 0, 此選項指定積云參數(shù)化方案，默認值為 0。 0 不采用積云參數(shù)化方案 1 淺對流Kain-Fritsch (new Eta)方案 2 Betts-Miller-Janjic 方案 3 Grell-Devenyi 集合方案 4 Simplified Arakawa-Schubert方案 99 老Kain-Fritsch

48、 方案,cud = 5, 5, 5, (積云參數(shù)化方案的調(diào)用時間間隔，默認值為 0, 單位為分鐘。 一般的積云參數(shù)化方案是每一步都要調(diào)用，但如果是用Kain-Fritsch 方案(cu_physics=1)，則可以設cudt=5) isfflx = 1, 在選用擾動邊界層和陸面物理過程時(sf_sfclay_physics = 1)是否考慮地面熱量和水汽通量，默認值為1。 0 不考慮地面通量 1 考慮地面通量 ifsnow = 0, 是否考慮雪蓋效應?？紤]雪蓋效應時，必須要有雪蓋輸入場。默認值為0，只有在利用擾動邊界層PBL預報土壤溫度時才有效，即sf_surface_physics = 1。

49、 0 不考慮雪蓋效應 1 考慮雪蓋效應 icloud = 1, 輻射光學厚度中是否考慮云的影響，默認值為1。僅當ra_sw_physics = 1 和 ra_lw_physics = 1時有效。 0 不考慮云的影響 1 考慮云的影響,surface_input_source = 1, 1 SI/gridgen（由SI的gridgen_model.exe程序產(chǎn)生） 土地利用類型和土壤類型數(shù)據(jù)的來源格式，默認值為1。 2 其他模式產(chǎn)生的GRIB碼數(shù)據(jù)(VEGCAT/SOILCAT 數(shù)據(jù)) num_soil_layers = 5, 指定陸面模式中的土壤層數(shù)，默認值為5 5 熱量擴散方案 4 Noah

50、 陸面過程方案 6 RUC 陸面過程方案,修改完退出保存 然后運行WRF模塊，命令如下： #ulimit s unlimited (or:unlimited ,此命令防止溢出問題出現(xiàn)，當出現(xiàn)此問題為溢出： （5 points exceeded cfl=2 in domain 1 at time 4.200000 MAX AT i,j,k: 123 48 3 cfl,w,d(eta)= 4.165821 21 points exceeded cfl=2 in domain 1 at time 4.200000 MAX AT i,j,k: 123 49 4 cfl,w,d(eta)= 10.662

51、90），解決方案檢查時間積分步長，減少時間積分步長） #./real.exe 運行成功會提示：SUCCESS COMPLETE REAL_EM INIT,#ls 將會有這幾個文件生成 wrfinput_d01 （初始場文件） wrfinput_d02 （初始場文件） wrfbdy_d01 （邊界條件文件） 然后接下來運行主模塊 #./wrf.exe 運行成功后會提示： SUCCESS COMPLETE WRF #ls wrfout_d01_2006-08-16_12:00:00 wrfout_d02_2006-08-16_12:00:00,到此WRF模式運行就完畢了，接下來就是后處理，大家也可

52、以用NCVIEW軟件直接看wrfout文件里面的變量（圖形化界面） 不過大部分同學要用模式轉(zhuǎn)出的數(shù)據(jù)做分析，所以我們要轉(zhuǎn)換格式！ 接下來介紹大家比較常用的畫圖軟件grads，那么后處理就用ARWpost軟件,8 ARWPOST后處理編譯8.1 ARWpost安裝,首先在/home/user/cuit下解壓tar文件 #tar zxvf ARWpost.tar.gz #ls #mv ARWpost arwpost #cd arwpost #ls 接下來就連接編譯安裝,#./configure Will use NETCDF in dir: /usr/local/netcdf-pgi - Plea

53、se select from among the following supported platforms. 1. PC Linux i486 i586 i686, PGI compiler (no vis5d) 2. PC Linux i486 i586 i686, PGI compiler (vis5d) 3. PC Linux i486 i586 i686, Intel compiler (no vis5d) 4. PC Linux i486 i586 i686, Intel compiler (vis5d) Enter selection 1-41 (輸入1),然后核實NETCDF的

54、路徑是否正確 #./compile 編譯安裝成功即可產(chǎn)生ARWpost.exe #ls *.exe ARWpost.exe 到此，后處理軟件ARWpost完成安裝，接下來即可修改namelist.ARWpost，然后把wrfout文件轉(zhuǎn)換成.dat和.ctl文件（接下來的這兩個文件大家都應該會處理了）,8.2 ARWpost運行,#cd arwpost #vi namelist.ARWpost &datetime start_date = 2006-08-16_12:00:00, end_date = 2006-08-16_18:00:00, interval_seconds =3600, （資料的時長，以秒計算，此例為domain2的就是60m,3600s） tacc = 0, debug_level = 0, /,&io io_form_input = 2, (NETCDF=2,GRIB1=5)