【正文】 ‘SI’, 或 ‘MM5’）寫成的，所以大部分用于看NetCDF文件的軟件都法使用了。盡管如此，一個以NCAR Graphics為基礎(chǔ)的應(yīng)用軟件—plotfmt（包含在WPS內(nèi)核中）可以解決這個問題。這個應(yīng)用程序制作出的要素場的等值線可以通過一個多媒體軟件 來查看。如果NCAR Graphics庫被正常安裝了，則當(dāng)編譯WPS時，plotfmt會和其它一些應(yīng)用程序一起被自動編譯。0 kamp。 Q) x ]amp。 ^, DS, oWPS應(yīng)用程序處理WPS的三個主要程序—geogrid，ungrib和metgrid—WPS還有很多其它的應(yīng)用程序會被編譯，并放在util目錄下。這些程序可以被用來檢查數(shù)據(jù)文件，看模擬區(qū)域的位置，計算氣壓場及計算平均的地表溫度場。 x% F, t4 Z39。 l1 J, Y: k [A． 程序是用來計算日平均地表溫度的，這些數(shù)據(jù)來自于過渡格式文件，“share”記錄中的相關(guān)參數(shù)定的， 并且同樣考慮到了文件之間的時間間隔，在計算平均值時，并用上盡可能多時次的 數(shù)據(jù)。如果一個中天沒有一個完整的時次，則不會有輸出文件生成，這個程序也會立即停止。類似的，如果一個文件的時次沒有包含在一個完整的24小時的周期 內(nèi)，則這個文件也將會被舍掉，例如，如果有5個過渡文件，之間的時間間隔為6個小時，則最后一個文件將被舍棄。計算出的平均場將以文件的形式輸出，名字為 TAVGSFC，格式同樣是過渡格式。日平均的地表溫度可以被metgrid所使用，方法是在“metgrid”記錄部分中設(shè)置 constants_name=“TAVGSFC”。B． ：press_pa = 201300 , 200100 , 100000 ,95000 , 90000 ,85000 , 80000 ,75000 , 70000 , c! x }$ T* e3 _65000 , 60000 ,45000 , 40000 ,35000 , 30000 ,25000 , 20000 ,! Z15000 , 10000 ,:5000 , 1000氣象,在amp。mod_levs記錄部分中的參數(shù)press_pa是用來保存層次列表的，這些確定的層次要與過渡格式文件中的xlvl（具體參考WPS intermediate format）的值相匹配。mod_levs程序要附帶兩個命令行參數(shù)作為它的輸入。第一個參數(shù)是要處理的過渡文件的名字，第二個是輸出文件的名字。從氣 象數(shù)據(jù)集中移除不需要的層次是很有用的。例如，當(dāng)一個數(shù)據(jù)集被用來作為模式的初始條件，另一個則作為邊界條件，可以通過在第一個時間周期中用 metgrid的插值來提供初始條件，而邊界條件的數(shù)據(jù)集則用于其它時次來實現(xiàn)。如果兩個數(shù)據(jù)集的垂直層次相同，則無需再做其它工作，但是當(dāng)它們的層次數(shù) 不同，則必須以最小的層次為標(biāo)準(zhǔn)，在有m層數(shù)據(jù)的文件中移除（mn）個層次，這里mn,并且m和n是兩個數(shù)據(jù)集中各自的層次數(shù)。必須有相同層次 ，它要求必須垂直層次必須是常值以便從中進(jìn)行插值。的應(yīng)用只 是處理兩個或更多具有不同垂直層次的數(shù)據(jù)集的一個臨時解決方法。如果用戶選擇使用mod_levs，則需要注意的是，盡管數(shù)據(jù)集之間垂直層次的設(shè)定可能不 同，但是所有數(shù)據(jù)都要有個地表的層次，參數(shù)p_top（）必須被設(shè)定成所有 數(shù)據(jù)集中最低層次以下的某一層。. o5 c4 |amp。 q) G B) s2 ^C． n M, c) v2 p/ Y, damp。 y39。 A, Y Q在垂直插值氣象要素場的過程中，real程序要求在每個相同的層次上和其它氣象要素一樣都有一個3D的氣壓場。 序可以被用于利用ECMWF的sigma坐標(biāo)數(shù)據(jù)集來生成這樣一個氣壓場?？紤]了地表氣壓場（地面氣壓場的記錄）及一些列系數(shù)A和 B，（i，j）點(diǎn)計算出氣壓，即Pijk = Ak Bk*Psfcij。這些系數(shù)可sigma高度系數(shù)的對應(yīng)表中找到（）。這個文件被寫在普通的txt文件，ecmwf_coeffs中。例如，16個sigma層的emcwf_coeffs文件會包含以下內(nèi)容：0 k* 1 2 3 4 5 ! w7 9 10 11 13 14 15 16 考 慮到ungrib生成的過渡格式文件和ecmwf_coeffs文件， 環(huán)，且每個時次都會生成一個額外的過渡格式文件PRES:YYYYMMDD_HH，文件里包含了完整的sigma層次和3D的相對濕度場。這個過渡 文件可以通過在參數(shù)fg_name中添加‘PRES’使之與ungrib生成的過渡文件一起被metgrid處理。D． 是一個以NCAR Graphics為基礎(chǔ)的應(yīng)用程序。，因此無需 運(yùn)行WPS的三個主程序。成功完成后，plotgrid生成了一個NCARG的元文件，gmeta，這個文件可以被idt所查看。粗糙區(qū)域（最外圍區(qū)） 占據(jù)整個畫面，一個地圖輪廓包圍了粗糙區(qū)域，并且如果有嵌套區(qū)域則會以矩形框的形式被標(biāo)注在最外圍區(qū)里（見下圖）。這個程序在最初設(shè)定模擬區(qū)域或每次調(diào)整 模擬區(qū)域位置時很有用。目前，這個程序還不能畫出ARW中用latlon投影方式（map_proj = 39。latlon39。）設(shè)置的區(qū)域。E． 程序會畫出文件中要素場，層次和數(shù)據(jù)的日期。下面是一個fnl文件的模板：Copen: File = fnl_080801_00_00 Fortran Unit = 0 I7 t3 @1 r2 Ob. u. V( X Jrec GRIB GRIBLvlLvlLvl Time Fcst氣Num Code nameCode onetwo hour 2 7 HGT 100975 020080801_00:00 00 3 7 HGT 100950 020080801_00:00 00 4 7 HGT 100925 020080801_00:00 00% \4 h, ^ 286 230 5WAVA100500 020080801_00:00 00Met***** EndOfFile on C unit 3. s( ]: A$ Y/ O7 c0 E. v，只是它針對的是GRIB2的數(shù)據(jù)文件。G． 是一個NCARG的程序，它畫的是過渡格式文件的內(nèi)容。程序唯一的命令行參數(shù)是要畫的文件名，并且生成NCARG的元文件，這個文件里包含了輸入文件 中各個要素場的等值線圖。graphics的源文件，gmeta，可以用idt來查看，并且用ctrans來轉(zhuǎn)換成其它格式。下面是一個模板：H．。================================================UNITS = Pa DESCRIPTION = Sealevel PressureDATE = 20080809_00:00:00 FCST = ]9SOURCE = NCEP GFS GDAS/FNL GRID 3LEVEL = ) a ^3 h39。 {3 h: J/ iI,J DIMS = 360, 181IPROJ = 0REF_X, REF_Y = , REF_LAT, REF_LON = , DLAT, DLON = , DATA(1,1)=! xamp。 p9 ^ ?) |$ b O! X================================================* T。 v0 xamp。 i5 t…………………………………[2 N7 I39。 m _: a/ \7 o39。 JFIELD = GHTDATE = 20080809_00:00:00 FCST = SOURCE = NCEP GFS GDAS/FNL GRID 38LEVEL = I,J DIMS = 360, 181IPROJ = 0REF_X, REF_Y = , DLAT, DLON = , DATA(1,1)=SUCCESSFUL COMPLETION OF PROGRAM RD_INTERMEDIATE把氣象要素場寫成過渡格式ungrib的作用是讀取GRIB數(shù)據(jù)集并寫成一個簡單的過渡格式，這個格式是可以被metgrid所識別的。如果某氣象要素并不存在于GRIB1/2格 式中，則用戶就有必要把這類數(shù)據(jù)寫入過渡格式文件中。幸運(yùn)的是，過渡格式是一種相對簡單的并且順序固定的，且是由Fortran寫成的無格式數(shù)據(jù)。需要注 意的是，這些無格式寫入是用得bigendian byte order（最高位字節(jié)在最前頭的排列順序），這種寫入方式可以由編譯所識別。下面來描述一下WPS的過渡格式；對SI或MM5的過渡格式熟悉的用戶會非 常容易熟悉WPS格式，因為它是非常類似的，而且WPS格式可以檢查Fortran的子程序，這些子程序（分別為metgrid/src /）可以讀寫所有的三種的過渡格式。當(dāng)在寫 成WPS過渡格式時，2D的要素場作為一個實型矩形數(shù)組來寫的。3D的數(shù)組則會以垂直維數(shù)為基礎(chǔ)分成多個2D數(shù)組，然后被獨(dú)立寫入。需要注意的是， 全球數(shù)據(jù)集—要么是高斯投影，要么是等距圓柱投影，而局地數(shù)據(jù)集—要么是麥卡托投影，要么是蘭伯托投影，再就是極射投影或等圓柱體投影。寫的順序用于在 WPS過渡格式中寫一個單獨(dú)的2D數(shù)組，正如下面所示（注意：并不是所有定義的變量都是用來做一個給定數(shù)據(jù)的投影）：) ? P. qinteger :: version ! Format version (must =5 for WPS format)Meinteger :: nx, ny ! x and ydimensions of 2d arrayinteger :: iproj ! Code for projection of data in array:c g9 n! {3 e) ] S/ K) }/ a! 0 = cylindrical equidistant! 1 = Mercator! 3 = Lambert conformal conic。! 4 = Gaussian (global only!)! 5 = Polar stereographic x(real :: nlats ! Number of latitudes north of equator! (for Gaussian grids)real :: xfcst ! Forecast hour of datareal :: xlvl ! Vertical level of data in 2d arrayreal :: startlat, startlon ! Lat/lon of point in array indicated by! startloc stringreal :: dx, dy ! Grid spacing, kmreal :: xlonc ! Standard longitude of projection6real :: truelat1, truelat2 ! True latitudes of projectionreal :: earth_radius ! Earth radius, kmreal, dimension(nx,ny) :: slab ! The 2d array holding the datalogical :: is_wind_grid_rel ! Flag indicating whether winds are! relative to source grid (TRUE) or, U8 V$ o! relative to earth (FALSE)character (len=8) :: startloc ! Which point in array is given by s: I, F/ ^/ p) n9 }! startlat/startlon。 set eitheramp。 f0 PF* J0 L! to 39。SWCORNER39。 or 39。CENTER 39。Mcharacter (len=9) :: field ! Name of the fieldcharacter (len=24) :: hdate ! Valid date for data YYYY:MM ![if !vml]![endif]D_HH:00:00character (len=25) :: units ! Units of datacharacter (len=32) :: map