【正文】 始時(shí)間）。這是因?yàn)?，嵌套是從它的父區(qū)域得到的邊界條件，因此，除了WRF中用到網(wǎng)格松弛，或者說松弛分析，嵌套只要在初始時(shí)間接受WPS的處理即可。但是很重要的一點(diǎn)是，當(dāng)運(yùn)行WRF時(shí)。剩下要做的就是對(duì)geogrid進(jìn)行設(shè)置了。在這塊記錄中，通過改變參數(shù)parent_id來設(shè)置每個(gè)區(qū)域的父區(qū)域，除了粗糙區(qū)域（最外圍區(qū)）的父區(qū)域是自己本身外，每個(gè)區(qū)域都必須是另一個(gè)區(qū)域的子區(qū)域。與每個(gè)區(qū)域的父區(qū)域的ID類似，每個(gè)父區(qū)域和其子區(qū)域的格點(diǎn)距離的比值是通過參數(shù)parent_grid_ratio來設(shè)定的。3 |$ q ?8 v d39。 h9 E1 e接下來，需要通過參數(shù)i_parent_start和j_parent_start來設(shè)定子區(qū)域左下角在父區(qū)域內(nèi)的位置（i，j），并且這個(gè)位置是在無參差格點(diǎn)的情況下確定的。每個(gè)區(qū)域的格點(diǎn)維數(shù)都要通過參數(shù)s_we, e_we,s_sn,和e_sn來設(shè)置。關(guān)于與格點(diǎn)有關(guān)的設(shè)置的模板在下圖已經(jīng)給出。其中，每個(gè)網(wǎng)格在南北（s_sn）向和東西（s_we）向的起始格點(diǎn)都要設(shè)置成1，而相應(yīng)的終止（e_sn和e_we）格點(diǎn)則是區(qū)域的總維數(shù)，為了保證每個(gè)嵌套區(qū)域的右上角與父區(qū)域的無參差格點(diǎn)想重合， e_sn和e_we的大小必須要比parent_grid_ratio所設(shè)數(shù)值的N倍再大一個(gè)。同樣的，每個(gè)區(qū)域所要插值的源數(shù)據(jù)的分辨率（可參考description of namelist variables）需要通過參數(shù)geog_data_res來設(shè)定。，獲得一個(gè)對(duì)namelist各個(gè)參數(shù)的完整描述。 USGS和基于MODIS的地形數(shù)據(jù)之間的選擇默認(rèn)情況下，geogrid方案將插值從美國(guó)地質(zhì)勘探局（USGS）的土地用途類別的24類數(shù)據(jù)。但是，用戶可以選擇基于國(guó)際地圈生物圈計(jì)劃MODIS遙感土地覆蓋分級(jí)的土地用途類別和修改Noah土壤表面模型來替代。雖然基于MODIS衛(wèi)星的數(shù)據(jù)包含20個(gè)土地使用類別，這些類別不是24個(gè)美國(guó)地質(zhì)調(diào)查局類別的一個(gè)子集。必須強(qiáng)調(diào)的是，在WRF模式中，基于MODIS衛(wèi)星為基礎(chǔ)的分類應(yīng)該只用于與Noah土壤表面模型。
在運(yùn)行時(shí)，通過“geogrid”名單記錄中的geog_data_res變量，20類基于MODIS衛(wèi)星的土地利用數(shù)據(jù)可被選來代替美國(guó)地質(zhì)調(diào)查局的數(shù)據(jù)。這是由添加每個(gè)靜態(tài)資料的分辨率前綴為字符串“modis_30s +”來完成的。例如，在一個(gè)三個(gè)區(qū)域的配置中， geog_data_res變量將通常被指定為Geog_data_res=’10m’,’20m’,’30’用戶應(yīng)改為指定geog_data_res =39。modis_30s十10m39。，39。modis_30s+2m39。，39。modis_30s +30 s39。這一變化的作用，39。modis_30s39。的靜態(tài)數(shù)據(jù)，若找不到分辨率為39。modis_30s39。的，則找分辨率字符串后帶39。+39。的靜態(tài)數(shù)據(jù)。因此，在此例中，由字符串39。modis_30s確定的基于MODIS衛(wèi)星的土地利用數(shù)據(jù)，將用來代替39。10m39。，39。2m39。，39。30s 39。分辨率的美國(guó)地質(zhì)調(diào)查局的數(shù)據(jù)。對(duì)所有其它要素場(chǎng)，39。10m39。，39。2m39。，39。30s 39。分辨率將分別用于第一，第二和第三區(qū)域。順便說一句，那么分辨率將會(huì)設(shè)為默認(rèn)值。 用于重力波拖曳方案的靜態(tài)數(shù)據(jù)在NMM（）中的重力波拖曳（GWD）方案要求14個(gè)來自WPS的靜態(tài)文件。事實(shí)上，不管GWD方案是否用于模式，這些文件都將被geogrid程序插值。當(dāng)GWD方案將不被使用時(shí)，該要素場(chǎng)將完全被WRF模式忽略。使用多種源氣象數(shù)據(jù)資料! Z$ T5 M( b E* R9 i+ H) D) w6 yMeteorological Numerical Model Union of China (MNMUC) 中國(guó)氣象數(shù)值模式聯(lián)盟metgrid程序可以插值不隨時(shí)間變化的要素場(chǎng)，也能從多個(gè)源數(shù)據(jù)中獲得要素場(chǎng)并進(jìn)行插值。第一個(gè)功能可以通過設(shè)置amp。metgrid記錄中的參數(shù)constants_name來實(shí)現(xiàn)，這個(gè)參數(shù)可以用來設(shè)置一系列具有過渡格式的文件的文件名（如果需要還可包含路徑信息），這些文件里面包含了不隨時(shí)間改變的要素場(chǎng)，并且在每次metgrid進(jìn)行處理時(shí)都可以被用上，例如，短時(shí)模擬可能會(huì)用上常值SST（海表溫度）場(chǎng)—這個(gè)場(chǎng)只需有一個(gè)時(shí)次的值就可以了，用法如下： M5,Grapes,HYCOM,LICOM,同化,集合預(yù)報(bào),模式聯(lián)盟,模式論壇5 Qj: z4 _3 W39。 Hamp。metgrid v* p1 x8 o I2 R2 g8 vconstants_name = 39。/data/ungribbed/constants/SST_FILE:20060816_1239。9 u9 R, ]/ n: E7 h8 u4 q/ ?D. {A a4 M/ \( n1 b9 _或者/ H9 Y% b) p+ p) amp。metgrid1 R5 H$ m1 u: T2 h: W。 = 39。LANDSEA39。, 39。SOILHGT39。9 V /! c8 pamp。 h4 o6 w Qamp。 V3 Z第二個(gè)功能—從多個(gè)源數(shù)據(jù)里獲得資源進(jìn)行插值—在有兩個(gè)或兩個(gè)以上補(bǔ)充數(shù)據(jù)（隨時(shí)間變化）?？梢酝ㄟ^amp。metgrid記錄中的參數(shù)fg_name設(shè)置一系列具有過渡格式的數(shù)據(jù)的前綴（如果需要還可包含路徑信息）來實(shí)現(xiàn)此功能。當(dāng)給出多個(gè)路徑前綴的設(shè)置，如果相同的氣象要素場(chǎng)在多個(gè)數(shù)據(jù)文件中出現(xiàn)，則最后一次出現(xiàn)的數(shù)值將被使用，不論之前出現(xiàn)過多少次，采用靠后原則。因此，可以通過調(diào)整源數(shù)據(jù)的順序來安排它們的優(yōu)先級(jí)別。下面是一個(gè)關(guān)于使用插值地面要素場(chǎng)和高空風(fēng)場(chǎng)到地形網(wǎng)格上的模板：8 t. h5 `7 }3 i* p氣象,數(shù)值,模式,氣象數(shù)值模式,數(shù)值模式,氣象模式,海洋模式,氣候模式,WRF,MM5,Grapes,HYCOM,LICOM,同化,集合預(yù)報(bào),模式聯(lián)盟,模式論壇amp。metgrid! n ]7 Q9 j0 Q5 F* q2 y1 ]fg_name = 39。/data/ungribbed/SFC39。, 39。/data/ungribbed/UPPER_AIR39。 rY, `8 |3 q// u K: I0 l3 ~7 ] amp。ungrib記錄中的參數(shù)prefix來達(dá)到簡(jiǎn)化從GRIB文件中提取要素場(chǎng)的步驟。這個(gè)參數(shù)允許用戶控制ungrib程序所生成的過渡文件的路徑和文件名。下面將給出這種方法在應(yīng)用上的一個(gè)模板。North American Regional Reanalysis (NARR)—北美區(qū)域再分析資料集被分為了三部分，分別是三維的大氣數(shù)據(jù)，地面數(shù)據(jù)和常值場(chǎng)數(shù)據(jù)，如果用戶想使用這種數(shù)據(jù)， GRIB數(shù)據(jù)進(jìn)行鏈接，然后再鏈接NARR Vtable到文件Vtable（ 程序ungrib一節(jié)），然后在運(yùn)行ungrib之前通過設(shè)置amp。ungrib記錄里的prefix使得產(chǎn)生的過渡文件有一個(gè)合適的前綴，模板如下：: k( f\3 a2 Xc R。 ]amp。( h) S9 H9 S* E: qamp。 Q0 Kout_format = 39。WPS39。,* ]4 q) B* ^% g) u2 p/ X氣象,數(shù)值,模式,氣象數(shù)值模式,數(shù)值模式,氣象模式,海洋模式,氣候模式,WRF,MM5,Grapes,HYCOM,LICOM,同化,集合預(yù)報(bào),模式聯(lián)盟,模式論壇prefix = 39。NARR_3D39。,氣象,數(shù)值,模式,氣象數(shù)值模式,數(shù)值模式,氣象模式,海洋模式,氣候盟, /氣象,數(shù)值,模式,氣象數(shù)值模式,數(shù)值模式,氣象模式,海洋模式,氣候模，將會(huì)出現(xiàn)如下的文件（附帶正確的日期）: @。 m5 O NARR_3D:20080816_129 J! U: g2 O2 `7 f+ g/ ` w氣象,數(shù)值,模式,氣象數(shù)值模式,數(shù)值模式,氣象模式,海洋模式,氣候模式,WRF,MM5,Grapes,HYCOM,LICOM,同化,集合預(yù)報(bào),模式聯(lián)盟,模式論壇NARR_3D:20080816_15氣象,數(shù)值,模式,氣象數(shù)值模式,數(shù)值模式,氣象模式,海洋模式,氣NARR_3D:20080816_18氣象,數(shù)值,模式,氣象數(shù)值模式,數(shù)值模式,氣象模式,海洋模式,氣.... J( P5 X$ c1 J+ Q4 _考慮到過渡文件是一個(gè)三維場(chǎng)的數(shù)據(jù)，因此當(dāng)處理地面數(shù)據(jù)時(shí)可以另外改一個(gè)prefix與之對(duì)應(yīng):amp。 S: W I8 U9 E氣象,數(shù)值,模式,氣象數(shù)值模式,數(shù)值模式,氣象模式,海洋模式,氣候模式,WRF,MM5,Grapes,HYCOM,LICOM,同化,集合預(yù)報(bào),模式聯(lián)盟,模式論壇amp。ungrib) n* Z/ F* ~8 l6 f6 w8 = 39。WPS39。,Meteorological Numerical Model Union of China (MNMUC) 中國(guó)氣prefix = 39。NARR_SFC39。,( f( Z! ?: u7 K6 _氣象,數(shù)值,模式,氣象數(shù)值模式,數(shù)值模式,氣象模式,海洋模式,氣候模式,WRF,MM5,Grapes,HYCOM,LICOM,同化,集合預(yù)報(bào),模式聯(lián)盟,模式論壇/* u m ?。 q5 {G1 q1 I9 C，會(huì)另外出現(xiàn)如下的數(shù)據(jù)文件: D2 |. e1 m/ t mNARR_SFC:20080816_129 m。 S J8 Q。 ] N2 fNARR_SFC:20080816_157 H H, K8 G o) M氣象,數(shù)值,模式,氣象數(shù)值模式,數(shù)值模式,氣象模式,海洋模式,氣候模式,WRF,MM5,Grapes,HYCOM,LICOM,同化,集合預(yù)報(bào),模式聯(lián)盟,模式論壇NARR_SFC:]) ~/ G/ C39。 r9 m0 w3 O Z...amp。 t0 `4 @* A7 W* S7 Z最后，當(dāng)處理常值數(shù)據(jù)時(shí)，同樣要給出合適的prefix:, v: Q8 w. o* w( B+ E4 G。 |$ \4 amp。ungrib, d。 BW6 M/ mMeteorological Numerical Model Union of China (MNMUC) 中國(guó)氣象數(shù)值模式聯(lián)盟out_format = 39。WPS39。,Meteorological Numerical Model Union of China (MNMUC) 中國(guó)氣prefix = 39。NARR_FIXED39。, N$ C+ F7 {1 `! y( V/7 t/ B5 R, D( N, V8 s，會(huì)得到如下數(shù)據(jù)文件:* B3 L8 P39。 k3 A% Mamp。 Bamp。 T6 C:19791108_00 8 p2 p mamp。 Y8 R) X! [6 F% ，常值數(shù)據(jù)最好省掉日期信息。如，因?yàn)閿?shù)據(jù)是靜態(tài)的，所以可以把上面提到的文件名改為NARR_FIXED。在這個(gè)例子中，NARR常值數(shù)據(jù)僅僅在一個(gè)特定的時(shí)次—1979年11月08日00時(shí)00分（UTC）可以使用，因此，當(dāng)運(yùn)行ungrib程序之前，用戶應(yīng)該設(shè)置好amp。share記錄中的起止日期；當(dāng)然，在運(yùn)行metgrid程序之前，這個(gè)日期還需再次被調(diào)整?？紤]到NARR數(shù)據(jù)集有三個(gè)過渡文件，amp。metgrid記錄里的相關(guān)參數(shù)—constants_name和fg_name，模板如下：amp。 X4 CC9 b, ^m氣象,數(shù)值,模式,氣象數(shù)值模式,數(shù)值模式,氣象模式,海洋模式,氣候模式,WRF,MM5,Grapes,HYCOM,LICOM,同化,集合預(yù)報(bào),模式聯(lián)盟,模式論壇amp。metgrid6 s8 D. S( [ V4 s6 ` Zconstants_name = 39。NARR_FIXED39。,氣象,數(shù)值,模式,氣象數(shù)值模式,數(shù)值模式,氣象模式,海洋fg_name = 39。NARR_3D39。, 39。NARR_SFC39。! V2 ~+ }39。 b0 F3 e* b6 o/氣象,數(shù)值,模式,氣象數(shù)值模式,數(shù)值模式,氣象模式,海洋模式,氣候模{( ]+ bamp。 O盡管可能性很小，但是還是會(huì)發(fā)生的一種情況是，當(dāng)一個(gè)來自區(qū)域模式的氣象源數(shù)據(jù)不足以覆蓋整個(gè)模擬區(qū)域時(shí)，則多個(gè)源數(shù)據(jù)將被使用，并且那些來著一個(gè)更大的區(qū)域模式或全球模式的數(shù)據(jù)將被插值到剩下的模擬網(wǎng)格上。例如，為了能在有可能的區(qū)域上應(yīng)用NAM的數(shù)據(jù)，剩下的區(qū)域應(yīng)用GFS的數(shù)據(jù)，需要設(shè)置namelist中的如下參數(shù)：2 p+ [!amp。metgridMeteorological Numerical Model Union of China (MNMUC) 中國(guó)氣象數(shù)值模式fg_name = 39。/data/ungribbed/GFS39。, 39。/data/ungribbed/NAM39。 Tlamp。 D9 b: Y/: z。 U C3 E* |$ Q5 p則插值后的結(jié)果如下圖：1 x。 k+ h。 F$ A4 q$ V, u: [如果要素場(chǎng)不會(huì)在其它源數(shù)據(jù)文件中出現(xiàn)，那就不會(huì)有優(yōu)先級(jí)一說，則metgrid進(jìn)行的插值就會(huì)像正常情況下一樣，每個(gè)源數(shù)據(jù)都要覆蓋整個(gè)模擬區(qū)域以免有的區(qū)域沒有數(shù)據(jù)。) a5 WPS中并行2 vF: c0 v7 i氣象,數(shù)值,模式,氣象數(shù)值模式,數(shù)值模式,氣象模式,海洋模式,氣候模式,WRF,MM5,Grapes,HYCOM,LICOM,同化,集合預(yù)報(bào),模式聯(lián)盟,模式論壇如果WPS中設(shè)置的模擬區(qū)域的維數(shù)過于龐大以至于單CPU無法勝任，則用分布式存儲(chǔ)器的方式來運(yùn)行g(shù)eogrid和metgrid程序是一個(gè)解決辦法。為了能把geogrid和metgrid程序以分布式存儲(chǔ)器的類型編譯，用戶必須在機(jī)器上安裝MPI庫(kù)，及在編譯時(shí)選中選項(xiàng)“DM parallel”。當(dāng)編譯成功后，geogrid和metgrid程序就可以在mpirun和mpiexec的命令，或通過批隊(duì)列系統(tǒng)下運(yùn)行了，具體的運(yùn)行方式則是根據(jù)機(jī)器的不同而定。正如之前提到的，ungrib程序是不能用并行運(yùn)行的，而且ungrib程序?qū)τ趦?nèi)存的需求是獨(dú)立于geogrid和metgrid程序之外的；因此，無論編譯時(shí)是否選擇“DM parallel”選項(xiàng)，u