wavewatch3.v4.18安裝使用教學(xué)教程(浙江海洋大學(xué)修改于2016年度0322)

.Wavewatch3.v4.18在服務(wù)器上安裝使用教程浙江海洋學(xué)院盛葉新鑒于wavewatch3.v4.18的安裝僅限于一本全英文的mannul手冊，網(wǎng)絡(luò)上沒有任何有幫助的安裝教程以及使用教程。因此，我于2015年10月在浙江海洋大學(xué)學(xué)習(xí)wavewatch3.v4.18時，在中國海大的兩位師兄的幫助下成功安裝了模式，僅針對我們學(xué)校的某個服務(wù)器安裝的一些問題和步驟進(jìn)行下列總結(jié)。一、環(huán)境變量設(shè)置安裝ww3之前需要加載編譯器并設(shè)置環(huán)境變量，不同的服務(wù)器有不同的加載方式，一般情況下使用modulelist查看當(dāng)前已加載編譯器，moduleavail查看可加載編譯器，moduleloadname加載name編譯器。安裝ww3需要加載ifort、icc和openmpi-Intel（注意：使用openmpi-Intel是因?yàn)槲覀兊募荷习惭b的是這個，其他mpi有關(guān)的編譯器沒有使用過，所以不太清楚）三個編譯器。但我使用的服務(wù)器上不可使用module命令（具體原因和服務(wù)器配置有關(guān)，我沒有權(quán)限深究，此處不做詳解），因此我們設(shè)置環(huán)境變量文件.bashrc文件。在修改之前，使用whichifort、whichicc和whichmpirun查看三個編譯的安裝目錄。根目錄下輸入vi.bashrc添加如下代碼：exportPATH=$PATH:/(使用查找出來的ifort和icc安裝目錄)/bin/(可能發(fā)生ifort和icc是同一個目錄，此處填寫該目錄)exportPATH=$PATH:/(使用查找出來的mpi安裝目錄)/binexportLD_LIBRARY_PATH=/(使用查找出來的mpi安裝目錄)/lib/:$LD_LIBRARY_PATH此時環(huán)境變量則修改完畢，記得使用source.bashrc加載設(shè)置好的環(huán)境變量，然后繼續(xù)完成安裝。二、ww3安裝和編譯1.安裝個人建議直接上傳ww3的壓縮文件，這樣可以節(jié)約時間，也可先在windows下解壓后在上傳，最好是建立一個新的文件夾，這樣不會發(fā)生文件找不到的情況。我提供的是4.18版本，文件夾示例：/public/home/shao002113/wwatch3.v4.18上傳完畢后，按如下步驟進(jìn)行：（注意空格）tarxvfwwatch3.v4.18-C/public/home/shao002113/wwatch3.v4.18cdwwatch3.v4.18chmod+xinstall_ww3_tar./install_ww3_tar接下來進(jìn)入安裝了，他會問你一些問題，讓你選擇（y/n）的時候基本都是y。有一個是選擇安裝模式，一個是L一個是G，安裝選擇G，然后有些設(shè)置需要修改，如下：Printerforlistings[printer]:printerCompilerforaux.[f77]:ifortCompilerforaux.[cc]:iccScratchspace[/tmp]：這是一個關(guān)于tmp文件夾的設(shè)置，第一次安裝不用修改Savesourcecodefiles:yesSavelistingfiles:yes之后出現(xiàn)的一律選擇yes，會顯示安裝結(jié)束end。2.添加環(huán)境變量在根目錄下輸入命令：vi.bashrc在文本中添加：exportPATH=$PATH:/(ww3完整安裝目錄)/binexportPATH=$PATH:/(ww3完整安裝目錄)/exe保存并退出。根目錄下輸入命令：source.bashrc（這是修改.bashrc后的必要步驟）。驗(yàn)證路徑是否添加成功，根目錄下輸入whichw3_make.若成功修改，則會出現(xiàn)ww3完整目錄。3.修改comp、link、switch進(jìn)入bin文件夾，將comp.Intel和link.Intel

復(fù)制成comp和link即：cpcomp.Intelcompcplink.Intellink使用vi編輯器修改comp和link。修改如下：comp中（大約96行處）修改成如下if["$mpi_mod"='yes']thencomp=mpiifortwhichmpiifort1>/dev/null2>/dev/nullOK=$?if[$OK!=0]thencomp=mpif90fielsecomp=ifortfilink中（大約104行）修改成如下if["$mpi_mod"='yes']thencomp=mpiifortwhichmpiifort1>/dev/null2>/dev/nullOK=$?if[$OK!=0]thencomp='mpif90'fiopt="-O3-o$prog"elsecomp=ifortfiswitch中把SHRD改為MPIDIST(注意此處空格，否則會出錯)。三個文件修改完畢后給予執(zhí)行權(quán)限。chmod+ucomplinkswitch4.運(yùn)行make_MPI直接輸入命令make_MPI即可。如果沒有錯誤顯示則編譯成功，一般會出現(xiàn)兩個錯誤，但不影響模式的運(yùn)行計算，到此ww3安裝編譯成功。三、數(shù)據(jù)提取WW3對于輸入數(shù)據(jù)文件的類型沒有固定要求，本人常用的是風(fēng)場為dat文件，地形為ASCII文件。那么輸入數(shù)據(jù)的格式有著一定的要求。對于風(fēng)場來說，數(shù)據(jù)排列方式如下：第一行是時間，如2015年6月1日0點(diǎn)0時0分，則寫為：20150601000000（注意年月日和時分秒之間存在空格）。接下來是兩個矩陣，第一個是風(fēng)場u，第二個是風(fēng)場v。接下來是地形，說是地形，其實(shí)就是水深，大于等于0的設(shè)為無效值就好。當(dāng)?shù)匦写鏋锳SCII文件時，務(wù)必知道其繪出的地形圖是什么形式的，在后面的修改腳本出需要用到。四、腳本修改1.ww3_grid.inp$--------------------------------------------------------------------$$WAVEWATCHIIIGridpreprocessorinputfile$$--------------------------------------------------------------------$$Gridname(C*30,inquotes)$'TESTGRID(GULFOFNOWHERE)'(此處是使用的網(wǎng)格名稱)$ $Frequencyincrementfactorandfirstfrequency(Hz)----------------$$numberoffrequencies(wavenumbers)anddirections,relativeoffset$offirstdirectionintermsofthedirectionalincrement[-0.5,0.5].$Inversions1.18and2.22ofthemodelthisvaluewasbydefiniton0,$itisaddedtomitigatetheGSEforafirstorderscheme.Notethat$thisfactorisIGNOREDintheprintplotsinww3_outp.$1.10.0411825240.（頻率間隔初始頻率波速波向）$$Setmodelflags----------------------------------------------------$$-FLDRYDryrun(input/outputonly,nocalculation).$-FLCX,FLCYActivateXandYcomponentofpropagation.$-FLCTH,FLCKActivatedirectionandwavenumbershifts.$-FLSOUActivatesourceterms.$FTTTFT（上方六項(xiàng)的開關(guān)）$$Settimesteps-----------------------------------------------------$$-Timestepinformation(thisinformationisalwaysread)$maximumglobaltimestep,maximumCFLtimestepforx-yand$k-theta,minimumsourcetermtimestep(allinseconds).$900.950.900.300.（第一個數(shù)和第三個數(shù)一般是第二個數(shù)的兩倍，第二個數(shù)使用程序計算，最后一個一般為30）$$Startofnamelistinputsection------------------------------------$$StartingwithWAVEWATCHIIIversion2.00,thetunableparameters$forsourceterms,propagationschemes,andnumericsarereadusing$namelists.Anynamelistfoundinthefolowingsectionsuptothe$end-of-sectionidentifierstring(seebelow)istemporarilywritten$toww3_grid.scratch,andreadfromthereifnecessary.Namelists$notneededforthegivenswitchsettingswillbeskipped$automatically,andtheorderofthenamelistsisimmaterial.$Asanexample,namelistinputtochangeSWELLFandZWNDinthe$TolmanandChalikovinputwouldbe$$&SIN2SWELLF=0.1,ZWND=15./（此處均屬于對源函數(shù)的調(diào)整，手冊有相關(guān)介紹）$$Defineconstantsinsourceterms-----------------------------------$$$Stresses------------------------------$TC1996withcap:NamelistFLX3$CDMAX:MaximumallowedCD(cap)$CTYPE:Captype:$0:Discontinuous(default).$1:Hyperbolictangent.$Hwang2011:NamelistFLX4$CDFAC:re-scalingofdrag$$Linearinput----------------------------$CavaleriandM-R:NamelistSLN1$CLIN:Proportionalityconstant.$RFPM:FactorforfPMinfilter.$RFHF:Factorforfhinfilter.$$Exponentialinput-------------------------$WAM-3:NamelistSIN1$CINP:Proportionalityconstant.$$TolmanandChalikov:NamelistSIN2$ZWND:Heightofwind(m).$SWELLF:swellfactorin(n.nn).$STABSH,STABOF,CNEG,CPOS,FNEG:$c0,ST0,c1,c2andf1in.(n.nn)$through(2.65)fordefinitionof$effectivewindspeed(!/STAB2).$WAM4andvariants:NamelistSIN3$ZWND:Heightofwind(m).$ALPHA0:minimumvalueofCharnockcoefficient$Z0MAX:maximumvalueofair-sideroughnessz0$BETAMAX:maximumvalueofwind-wavecoupling$SINTHP:powerofcosineinwindinput$ZALP:waveageshifttoaccountforgustiness$TAUWSHELTER:shelteringofshortwavestoreduceu_star$SWELLFPAR:choiceofswellattenuationformulation$(1:TC1996,3:ACC2008)$SWELLF:swellattenuationfactor$ExtraparametersforSWELLFPAR=3only$SWELLF2,SWELLF3:swellattenuationfactors$SWELLF4:ThresholdReynoldsnumberforACC2008$SWELLF5:Relativeviscousdecaybelowthreshold$Z0RAT:roughnessforoscil.flow/meanflow$BYDRZinput:NamelistSIN6$SINA0:factorfornegativeinput$SINU10:windspeedscalingoption$$Nonlinearinteractions-----------------------$DiscreteI.A.:NamelistSNL1$LAMBDA:Lambdainsourceterm.$NLPROP:Cinsourcterm.NOTE:default$valuedependsonothersource$termsselected.$KDCONV:FactorbeforekdinEq.(n.nn).$KDMIN,SNLCS1,SNLCS2,SNLCS3:$Minimumkd,andconstantsc1-3$indepthscalingfunction.$Exactinteractions:NamelistSNL2$IQTYPE:Typeofdepthtreatment$1:Deepwater$2:Deepwater/WAMscaling$3:Shallowwater$TAILNL:Parametrictailpower.$NDEPTH:Numberofdepthsinforwhich$integrationspaceisestablished.$UsedforIQTYPE=3only$NamelistANL2$DEPTHS:ArraywithdepthsforNDEPTH=3$Gen.MultipleDIA:NamelistSNL3$NQDEF:Numberofquadruplets.$MSC:Scalingconstant'm'.$NSC:Scalingconstant'N'.$KDFD:Deepwaterrelativefilterdepth,$KDFS:Shallowwaterrelativefilterdepth,$NamelistANL3$QPARMS:5xNQDEFparamatersdescribingthe$quadruplets,repeatingLAMBDA,MU,DT12.$CdeepandCshal.Seeexamplesbelow.$$TraditionalDIAsetup(default):$$&SNL3NQDEF=1,MSC=0.00,NSC=-3.50/$&ANL3QPARMS=0.250,0.000,-1.0,0.1000E+08,0.0000E+00/$$GMD3from2010report(G13dinlaterpaper):$$&SNL3NQDEF=3,MSC=0.00,NSC=-3.50/$&ANL3QPARMS=0.126,0.000,-1.0,0.4790E+08,0.0000E+00,$0.237,0.000,-1.0,0.2200E+08,0.0000E+00,$0.319,0.000,-1.0,0.1110E+08,0.0000E+00/$$G35dfrom2010report:$$&SNL3NQDEF=5,MSC=0.00,NSC=-3.50/$&ANL3QPARMS=0.066,0.018,21.4,0.170E+09,0.000E+00,$0.127,0.069,19.6,0.127E+09,0.000E+00,$0.228,0.065,2.0,0.443E+08,0.000E+00,$0.295,0.196,40.5,0.210E+08,0.000E+00,$0.369,0.226,11.5,0.118E+08,0.000E+00/$$NonlinearfilterbasedonDIA-------------------$NamelistSNLS$A34:Relativeoffsetinquadruplet$FHFC:Proportionalityconstants.$DMN:Maximumrelativechange.$FC1-3:Constantsinfrequencyfilter.$$Dissipation----------------------------$WAM-3:NamelistSDS1$CDIS,APM:Asinsourceterm.$$TolmanandChalikov:NamelistSDS2$SDSA0,SDSA1,SDSA2,SDSB0,SDSB1,PHIMIN:$Constantsa0,a1,a2,b0,b1and$PHImin.$$WAM4andvariants:NamelistSDS3$SDSC1:WAM4Cdscoeffient$MNMEANP,WNMEANPTAIL:powerofwavenumber$formeandefinitionsinSdsandtail$SDSDELTA1,SDSDELTA2:relativeweights$ofkandk^2partsofWAM4dissipation$SDSLF,SDSHF:coefficientforactivationof$WAM4dissipationforunsaturated(SDSLF)and$saturated(SDSHF)partsofthespectrum$SDSC2:Saturationdissipationcoefficient$SDSC4:ValueofB0=B/BrforwichSdsiszero$SDSBR:ThresholdBrforsaturation$SDSP:powerof(B/Br-B0)inSds$SDSBR2:ThresholdBr2fortheseparationof$WAM4dissipationinsaturatedandnon-saturated$SDSC5:coefficientforturbulencedissipation$SDSC6:WeightfortheistropicpartofSds_SAT$SDSDTH:Angularhalf-widthforintegrationofB$$BYDRZ:NamelistSDS6$SDSET:Selectthresholdnormalizationspectra$SDSA1,SDSA2,SDSP1,SDSP2:$CoefficientsfordissipationtermsT1andT2$:NamelistSWL6$SWLB1:Coefficientforswelldissipation$$Bottomfriction--------------------------$JONSWAP:NamelistSBT1$GAMMA:Asitsays.$$$Surfbreaking---------------------------$BattjesandJanssen:NamelistSDB1$BJALFA:Dissipationconstant(default=1)$BJGAM:Breakingthreshold(default=0.73)$BJFLAG:TRUE-UseHmax/dratioonly(default)$FALSE-UseHmax/dinMicheformulation$$Triadnonlinearinteractions--------------------$LumpedTriadInteraction(LTA):NamelistSTR1(Tobeimplemented)$PTRIAD1:Proportionalitycoefficient(default0.05)$PTRIAD2:MultipleofTm01uptowhichinteraction$iscomputed(2.5)$PTRIAD3:Ursellupperlimitforcomputing$interactions(notused,default10.)$PTRIAD4:Shapeparameterforbiphase$computation(0.2)$PTRIAD5:Ursellnumbertresholdforcomputing$interactions(0.01)$$Shorelinereflections------------------------$ref.parameters:NamelistREF1$REFCOAST:Reflectioncoefficientatshoreline$REFFREQ:Activationoffreq-dependentref.$REFMAP:Scalefactorforbottomslopemap$REFRMAX:maximumref.coeffient(default0.8)$REFFREQPOW:poweroffrequency$REFICEBERG:Reflectioncoefficientforicebergs$REFSUBGRID:Reflectioncoefficientforislands$REFCOSP_STRAIGHT:powerofcosineusedfor$straightshoreline$$Bound2ndorderspectrumandfreeIG-----------------$IG1parameters:NamelistIG1$IGMETHOD:1:Hasselmann,2:Krasitskii-Janssen$IGADDOUTP:activationofboundwavecorrection$inww3_outp/ww3_ounp$IGSOURCE:1:usesboundwaves,2:empirical$IGSTERMS:>0:nosourceterminIGband$IGMAXFREQ:maximumfrequencyofIGband$IGEMPIRICAL:constantinempiricalfreeIGsource$IGSWELLMAX:activatesfreeIGsourcesforallfreq.$$Propagationschemes------------------------------------------------$$SMCgridpropagation:NamelistPSMC$CFLTM:MaximumCFLnumberforpropagation.$DTIME:Swellage(s)fordiffusionterm.$LATMIN:Maximumlatitude(deg)forGCT.$RFMAXD:Maximumrefractionturning(deg).$&PSMCDTIME=39600.0,LATMIN=85.0,RFMAXD=36.0/$$Propagationschemes------------------------------------------------$$Firstorder:NamelistPRO1$CFLTM:MaximumCFLnumberforrefraction.$$UQ/UNOwithdiffusion:NamelistPRO2$CFLTM:MaximumCFLnumberforrefraction.$DTIME:Swellage(s)ingardensprinkler$correction.If0.,alldiffusion$switchedoff.Ifsmallnon-zero$(DEFAULT!!!)onlywavegrowth$diffusion.$LATMIN:Maximumlatitudeusedincalc.of$strengthofdiffusionforprop.$$UQ/UNOwithaveraging:NamelistPRO3$CFLTM:MaximumCFLnumberforrefraction.$WDTHCG:Tuningfactorpropag.direction.$WDTHTH:Tuningfactornormaldirection.$$NotethatUQandUNOschemeshavenotunableparameters.$Alltuneableparametersareassociatedwiththerefraction$limitationandtheGSEalleviation.$$Unstructuredgrids------------------------------------------------$$UNSTparameters:NamelistUNST$UGOBCAUTO:AutomaticdetectionofOBCpoints$UGOBCDEPTH:Threshold(<0)depthforOBCpoints$EXPFSN:ActivationofNscheme$EXPFSPSI:ActivationofPSIscheme$EXPFSFCT:ActivationofFCTscheme$IMPFSNIMP:ActivationofNimplicitscheme$$SMCgridpropagation:NamelistPSMC$CFLTM:MaximumCFLnumberforpropagation.$DTIME:Swellage(s)fordiffusionterm.$LATMIN:Maximumlatitude(deg)forGCT.$RFMAXD:Maximumrefractionturning(deg).$&PSMCDTIME=39600.0,LATMIN=85.0,RFMAXD=36.0/$$Miscellaneous------------------------------------------------------$$Misc.parameters:NamelistMISC$CICE0:Iceconcentrationcut-off.$CICEN:Iceconcentrationcut-off.$PMOVE:PowerpinGSEaleviationfor$movinggridsinEq.(D.4).$XSEED:Xseedinseedingalg.(!/SEED).$FLAGTR:Indicatingpresenceandtypeof$subgridinformation:$0:Nosubgridinformation.$1:Transparanciesatcellboun-$dariesbetweengridpoints.$2:Transp.atcellcenters.$3:Like1withcont.ice.$4:Like2withcont.ice.$XP,XR,XFILT$Xp,XrandXfforthedynamic$integrationscheme.$IHMAX:Numberofdiscretelevelsinpart.$HSPMIN:MinimumHsinpartitioning.$WSM:Windspeedmultiplierinpart.$WSC:Cutofwindseafractionfor$identifyingwindseainpart.$FLC:Flagforcombiningwindseasin$partitioning.$NOSW:Numberofpartitionedswellfields$infieldoutput.$FMICHE:ConstantinMichelimiter.$P2SF:......$$Inthe'Outofthebox'testsetupwerunwithsub-gridobstacles$andwithcontinuousicetreatment.$&MISCCICE0=0.25,CICEN=0.75,FLAGTR=4/（4改為0）&FLX3CDMAX=3.5E-3,CTYPE=0/$&SDB1BJGAM=1.26,BJFLAG=.FALSE./$$Mandatorystringtoidentifyendofnamelistinputsection.$ENDOFNAMELISTS$$Definegrid--------------------------------------------------------$$$Fiverecordscontaining:$$1Typeofgrid,coordinatesystemandtypeofclosure:GSTRG,FLAGLL,$CSTRG.Gridclosurecanonlybeappliedinsphericalcoordinates.$GSTRG:Stringindicatingtypeofgrid:$'RECT':rectilinear$'CURV':curvilinear$FLAGLL:Flagtoindicatecoordinatesystem:$T:Spherical(lon/latindegrees)$F:Cartesian(meters)$CSTRG:Stringindicatingthetypeofgridindexspaceclosure:$'NONE':Noclosureisapplied$'SMPL':Simplegridclosure:Gridisperiodicinthe$:i-indexandwrapsati=NX+1.Inotherwords,$:(NX+1,J)=>(1,J).Agridwithsimpleclosure$:mayberectilinearorcurvilinear.$'TRPL':Tripolegridclosure:Gridisperiodicinthe$:i-indexandwrapsati=NX+1andhasclosureat$:j=NY+1.Inotherwords,(NX+1,J<=NY)=>(1,J)$:and(I,NY+1)=>(MOD(NX-I+1,NX)+1,NY).Tripole$:gridclosurerequiresthatNXbeeven.Agrid$:withtripoleclosuremustbecurvilinear.$2NX,NY.Astheoutergridlinesarealwaysdefinedasland$points,theminimumsizeis3x3.$$Branchherebasedongridtype$$IF(RECTILINEARGRID)THEN$$3GridincrementsSX,SY(degr.orm)andscaling(division)factor.$IfCSTRG='SMPL',thenSXissetto360/NX.$4Coordinatesof(1,1)(degr.)andscaling(division)factor.$$ELSEIF(CURVILINEARGRID)THEN$$3Unitnumberoffilewithx-coordinate.$Scalefactorandaddoffset:x<=scale_fac*x_read+add_offset.$IDLA,IDFM,formatforformattedread,FROMandfilename.$IDLA:Layoutindicator:$1:Readline-by-linebottomtotop.$2:Like1,singlereadstatement.$3:Readline-by-linetoptobottom.$4:Like3,singlereadstatement.$IDFM:formatindicator:$1:Freeformat.$2:Fixedformatwithaboveformatdescriptor.$3:Unformatted.$FROM:filetypeparameter$'UNIT':openfilebyunitnumberonly.$'NAME':openfilebynameandassigntounit.$$Iftheaboveunitnumberequals10,thenthex-coordisreadfromthis$file.Thex-coordmustfollowtheaboverecord.Nocommentlinesare$allowedwithinthex-coordinput.$$4Unitnumberoffilewithy-coordinate.$Scalefactorandaddoffset:y<=scale_fac*y_read+add_offset.$IDLA,IDFM,formatforformattedread,FROMandfilename.$IDLA:Layoutindicator:$1:Readline-by-linebottomtotop.$2:Like1,singlereadstatement.$3:Readline-by-linetoptobottom.$4:Like3,singlereadstatement.$IDFM:formatindicator:$1:Freeformat.$2:Fixedformatwithaboveformatdescriptor.$3:Unformatted.$FROM:filetypeparameter$'UNIT':openfilebyunitnumberonly.$'NAME':openfilebynameandassigntounit.$$Iftheaboveunitnumberequals10,thenthey-coordisreadfromthis$file.They-coordmustfollowtheaboverecord.Nocommentlinesare$allowedwithinthey-coordinput.$$ELSEIF(UNSTRUCTUREDGRID)THEN$Nothingtodeclare:allthedatawillbereadfromtheGMESHfile$ENDIF(CURVILINEARGRID)$$5Limitingbottomdepth(m)todiscriminatebetweenlandandsea$points,minimumwaterdepth(m)asallowedinmodel,unitnumber$offilewithbottomdepths,scalefactorforbottomdepths(mult.),$IDLA,IDFM,formatforformattedread,FROMandfilename.$IDLA:Layoutindicator:$1:Readline-by-linebottomtotop.$2:Like1,singlereadstatement.$3:Readline-by-linetoptobottom.$4:Like3,singlereadstatement.$IDFM:formatindicator:$1:Freeformat.$2:Fixedformatwithaboveformatdescriptor.$3:Unformatted.$FROM:filetypeparameter$'UNIT':openfilebyunitnumberonly.$'NAME':openfilebynameandassigntounit.$$Iftheaboveunitnumberequals10,thenthebottomdepthsarereadfrom$thisfile.Thedepthsmustfollowtheaboverecord.Nocommentlinesare$allowedwithinthedepthinput.Inthecaseofunstructuredgrids,thefile$isexpectedtobeaGMESHgridfilecontainingnodeandelementlists.$$------------------------------------------------------------------------$Exampleforrectilineargridwithspherical(lon/lat)coordinatesystem.$NotethatforCartesiancoordinatestheunitismeters(NOTkm).$'RECT'T'NONE'（此處關(guān)于地形的設(shè)置）1212（網(wǎng)格大?。航?jīng)度格點(diǎn)數(shù)緯度格點(diǎn)數(shù)）1.1.4.（就是一度的方格里有多少格）-1.-1.4.（起始經(jīng)度起始緯度步長（一般為1））-0.12.5010（為通道號，最好使用其他數(shù)）-10.（地形數(shù)值乘以該數(shù)等于真實(shí)值）3（代表地形文件從上往下讀，1代表從下往上讀，根據(jù)地形文件務(wù)必保證數(shù)據(jù)讀取是從北緯到南緯）1'(....)''NAME''bottom.inp'（將‘bottom.inp’改為自己的地形文件名）$（下面加上$）666666666666666542024566666542024566666542024566666542004566666544224566666655445666666666556666666666666666666666666666666666666666666666666666$$------------------------------------------------------------------------$Exampleforcurvilineargridwithspherical(lon/lat)coordinatesystem.$Samespatialgridasprecedingrectilinearexample.$NotethatforCartesiancoordinatestheunitismeters(NOTkm).$$'CURV'T'NONE'$1212$$100.25-0.531'(....)''NAME''x.inp'$$123456789101112$123456789101112$123456789101112$123456789101112$123456789101112$123456789101112$123456789101112$123456789101112$123456789101112$123456789101112$123456789101112$123456789101112$$100.250.531'(....)''NAME''y.inp'$$111111111111$222222222222$333333333333$444444444444$555555555555$666666666666$777777777777$888888888888$999999999999$101010101010101010101010$111111111111111111111111$121212121212121212121212$$-0.12.5010-10.31'(....)''NAME''bottom.inp'$$666666666666$666542024566$666542024566$666542024566$666542004566$666544224566$6