教程plots這個示例展示如何GRIB文件的信息并從中讀取數(shù)據(jù)使用

Exmaple4— line這個示例展示如何獲得GRIB文件的信息并從中數(shù)據(jù)，如何使用“if-endifresources值，如何創(chuàng)建并繪制流線圖，如何把數(shù)據(jù)寫入netCDF文件，如何用“do”來進行循這個示例從GRIB文件數(shù)據(jù)并且創(chuàng)建了三個流線圖，每個圖都使用不同的數(shù)據(jù)。設(shè)置resources來改變每個流線圖的部分。為了得到打開的GRIB文件信息，這個示例使用了NCL函數(shù)isatt,getfilevarnames,getfilevaratts,getfilevardims，來確定文件中有哪些變量和屬NCLReferenceManual中的“Built-inNCLfunctionsandprocedures”GRIB數(shù)據(jù)寫入到netCDF運行這個示例，必須以下文件：gsun04n.ncl，然后鍵入：ncl示例4代碼及解load載入本示例使用的包含函數(shù)和程序（以gsn開頭的）的NCL。NCL中的load語句的作用和C和Fortran90程序中include作用一樣。開始NCLdata_dir=grb_file=addfile(data_dir+用addfile打開GRIBGRIBnames= ;Getthevariablenamesin ;GRIBfileandprintthemgetfilevarnamesGRIB文件中變量的名字，并把它們打印出來。Names是一個字符串數(shù)組，每個字符串都是變量名字。這個功能在你不知道用addfile打開的文件中有什么變量的時候特別有用。atts=getfilevaratts(grb_file,names(0));Getthevariableattributesdims=getfilevardims(grb_file,names(0));dimensionnamesfrom

;fileandprintthemgetfilevarattsgetfilevardimsGRIB文件中第一個變量的wks=gsn_open_wks("x11","gsun04n");OpenanX11X11 Beginfirst =準備設(shè)置一系列的resources1的關(guān)于如何設(shè)置resources的解釋。流線圖resources是“StreamlinePlot”組的一部分，以“ectorFielduvar=grb_file->U_GRD_6_ISBL(0,:,:);StoretwoGRIBvariablesvvar=grb_file->V_GRD_6_ISBL(0,:,:);localNCLGRIBU_GRD_6_ISBLV_GRD_6_ISBL。這些變（0if(isatt(uvar,"units").eq.True);Checkif"uvar"hasanattributeresources@tiMainString="GRD_6_ISBL(u,v"+uvar@units+resources@tiMainString=endif Thespaceisrequiredbetweenend"andif".NCLC和Fortran包含“if-else-endif”語句，但有以下例外：NCL中的if語句，不需要thenelseif在end和if之間必須有空格（NCL不能識別“endif.eq..lt.Fotran語法一樣。對于if語句的信息，請參閱NCLReferenceManual的“NCLStatements”一節(jié)。性，否則是False。在上述if語句中，如果isatt是True的話，那么“units”屬性返回的是一個邏輯值，如果沒有.eq.True的話，if語句的部分就應(yīng)該寫成： ="Times-resources@tiXAxisString=設(shè)置一些resources來改變標題的默認字體和X軸的。請記住：事先定義的字符串，就像字體表fonttable中列出的那些，是不區(qū)分大小寫的，你可以plot=gsn_streamline(wks,uvar(::2,::2),vvar(::2,::2),resources);

;drawastreamline創(chuàng)建并繪制二維數(shù)組uvar和vvar的流線圖。gsn_streamline參數(shù)是前面調(diào)用gsn_open_wks的返回工作站變量。接下來的兩個參數(shù)代表了要被畫的兩個二維矢量場（float、double或者integer，順序是Y×X。最后一個參數(shù)是一個邏輯值表示你是否設(shè)置了任何resources。這里設(shè)置的resources僅僅是為了流線圖的標題和X軸的。（2X/Y0到n-1，n;-----------Beginsecondplot----------------------------------------用不同的數(shù)據(jù)創(chuàng)建一個流線圖，給流線，改變標題和的字體delete(uvar);Deleteuvarandvvarsinceyouaregoingtodelete(vvar);themagain,butwithdifferentsized因為又要使用這些變量，所以刪除變量uvar和vvar。這只有在當你打算給uvar和vvar分配和初始化這兩個變量的時候的不同類型、不同維度的數(shù)據(jù)時才uvar=grb_file->U_GRD_6_TRO;StoretwodifferentGRIBfilevvar=grb_file->V_GRD_6_TRO;touvarand從GRIB文件中獲得變量U_GRD_6_TRO和V_GRD_6_TRO。這個量變量resources@tiMainString="GRD_6_TRO(u,v"+uvar@units+resources@tiMainString=end ="Times-Roman";Changethedefault

resources@tmXBLabelFont="Times-resources@tmYLLabelFont="Times-設(shè)置一些改變X軸和刻度字體的resources?？潭萺esources是“TickMark”組的一部分，以“tm”開頭。在tmXBLabelFonts中的“XB”表示“XBottomtmYLLabelFont中的“YL”表示“Yleft8resources@stLineColor="green";Changestreamlinesto使用NCL名顏色的功能，將流線的顏色改成了綠色。有650中有效顏色名稱可以使用，但是你的顏中必須包含你想要用名的顏色，否NCL將會使用最近的顏色代替。例如，如果你的顏中包含下列的“ConrflowerBlue，.,.,./（/0.,.,./，名的顏色，然后為你的顏指定他們。GSUN示例7將會展示如何用命名的來設(shè)置顏有關(guān)命名的顏色，請參見“Basics”一章中的“CreatingyourowncolormapusingnamedcolorsandRGBvalues”部分。plot=gsn_streamline(wks,uvar(::2,::2),vvar(::2,::2),resources);Drawstreamlineplot.用新設(shè)置的resources Beginthirdgetvalues ;Retrievesomeresource :"stMinLineSpacingF":end正如第一個示例中所示，resources都有默認值。一些resources是根據(jù)一些resources的值，你可以改變它們或者使用它們的值去改變其他的resources。StArrowLengthF和stMinLineSpacingFresources是兩個動態(tài)resources的例子，你可以獲得它們的值NCL中resources的值通過使用getvalues來獲得。Getvalues語句就像一個塊語句，因為它有一個開始語句（getvaluesplotid）和一個結(jié)束語句（endgetvaluesresources的列表（每個resource應(yīng)當在雙引號中，緊接著一個冒號，然后是一個要寫入數(shù)據(jù)的變量。變量plotid是gsn_*繪圖進程的一個返回值，并且resources要用到的（換句話說，如果已經(jīng)調(diào)用了gsn_vector來創(chuàng)建矢量圖的話，你就不應(yīng)該從這個圖中獲得流線圖resourcesgsn_streamline，變量名稱是plot。的getvalues的信息可以在NCLReferenceManual的“NCLresources@stMinLineSpacingF=spacing*2.0 ;Setsomeresourcesbased =arrowlength*2.0;onresourcesyou = ;Changelinecolorsback = ;既然已經(jīng)得到了stMinLineSpaingF和stArrowLengthF的值并且分別在變量spacing和arrowlength中，你可以使用這些值來改變這些resources。流線之間的距離以4倍增長，流線之間的箭頭長度翻倍。最后兩個resources將曲線50%的線條的粗細（1uvar=grb_file->U_GRD_6_GPML(0,:,:);GetnewGRIBvvar=grb_file-GRIBU_GRD_6_GPMLV_GRD_6_GPML的第一層resources@tiMainString="GRD_6_GPML(u,v"+uvar@units+resources@tiMainString=endplot=gsn_streamline(wks,uvar(::2,::2),vvar(::2,::2),resources);Drawstreamlineplot.用新建的resourcesdelete(uvar);Removesomevariableswedon'tneedany;-----------WriteGRIBdatatonetCDFfile-------------------------GRIB變量到netCDF文件中以“PRES”開始。這里展示的寫文件的方GRIB文件到netCDF中，請參閱“Basics”章中的“Efficiencyissuescdf_filename=system("/bin/rm-f"+cdf_file=addfile(cdf_filename,"c");CreateanewnetCDFGRIB數(shù)據(jù)的netCDF文件。首先刪除這個文件，如果他已經(jīng)存在的話，然后使用addfile來創(chuàng)建一個文件。通過“.ncNCL知道要創(chuàng)建的是netCDF文件。cdf_file@title="datafromaGRIBfile";Addsomeglobalattributescdf_file@date= ;thenetCDF添加叫做“title”和“date”的全局變量到netCDF文件。systemfunc函數(shù)doi=0,dimsizes(names)-NCL中的Do循環(huán)和Fortran中的比較類似。它們都有一個整型的起始值、enddo（enddo沒有用）結(jié)束。有關(guān)do循環(huán)的信息，請參見NCLReferenceManual的“NCLStatements”部分。在上面的doGRIB文件中獲取的變量的名字進行循環(huán)。請記住，dimsizes返回的一個數(shù)組的維度大小，因為名字是一個一維數(shù)組，所以dimsizes返回的是它的長度，是一個標量整數(shù)。names_char=組（NCL中的字符串被認為是一個單一的實體stringtocharif(names_char(0:3).eq."PRES") ;Onlywritevariablesthat ;withcdf_file->$names(i)$=grb_file-end上述4行中，names_char