wavewatch3.v4.18安装使用教学教程(浙江海洋大学修改于2016年度0322)

1、Wavewatch3.v4.18在服务器上安装使用教程浙江海洋学院盛叶新浙江海洋学院盛叶新11k11鉴于wavewatch3.v4.18的安装仅限于一本全英文的mannul手册，网络上没有任何有帮助的安装教程以及使用教程。因此，我于2015年10月在浙江海洋大学学习wavewatch3.v4.18时，在中国海大的两位师兄的帮助下成功安装了模式，仅针对我们学校的某个服务器安装的一些问题和步骤进行下列总结一、环境变量设置安装ww3之前需要加载编译器并设置环境变量，不同的服务器有不同的加载方式，一般情况下使用modulelist查看当前已加载编译器，moduleavail查看可加载编译器，modul

2、eloadname加载name编译器。安装ww3需要加载ifort、icc和openmpi-Intel（注意：使用openmpi-Intel是因为我们的集群上安装的是这个，其他mpi有关的编译器没有使用过，所以不太清楚）三个编译器。但我使用的服务器上不可使用module命令（具体原因和服务器配置有关，我没有权限深究，此处不做详解）,因此我们设置环境变量文件.bashrc文件。在修改之前，使用whichifort、whichicc和whichmpirun查看三个编译的安装目录。根目录下输入vi.bashrc添加如下代码：exportPATH=$PATH:/使用查找出来的ifort和icc安装目录

3、）/bin/（可能发生ifort和icc是同一个目录，此处填写该目录）exportPATH=$PATH:/使用查找出来的mpi安装目录）/binexportLD_LIBRARY_PATH=/（吏用查找出来的mpi安装目录）/lib/:$LD_LIBRARY_PATH此时环境变量则修改完毕，记得使用source.bashrc加载设置好的环境变量，然后继续完成安装。二ww3安装和编译安装个人建议直接上传ww3的压缩文件，这样可以节约时间，也可先在windows下解压后在上传，最好是建立一个新的文件夹，这样不会发生文件找不到的情况。我提供的是4.18版本，文件夹示例：/public/home/sha

4、o002113/wwatch3.v4.18上传完毕后，按如下步骤进行：（注意空格）tarxvfwwatch3.v4.18-C/public/home/shao002113/wwatch3.v4.18cdwwatch3.v4.18chmod+xinstall_ww3_tar./install_ww3_tar接下来进入安装了，他会问你一些问题，让你选择（y/n）的时候基本都是y有一个是选择安装模式，一个是L一个是G,安装选择G,然后有些设置需要修改，如下：Printerforlistingsprinter:printerCompilerforaux.f77:ifortCompilerforaux.

5、cc:iccScratchspace/tmp:这是一个关于tmp文件夹的设置，第一次安装不用修改Savesourcecodefiles:yesSavelistingfiles:yes之后出现的一律选择yes，会显示安装结束end。添加环境变量在根目录下输入命令：vi.bashrc在文本中添加:exportPATH=$PATH:/（ww3完整安装目录）/binexportPATH=$PATH:/（ww3完整安装目录）/exe保存并退出。根目录下输入命令：source.bashrc（这是修改.bashrc后的必要步骤）。验证路径是否添加成功，根目录下输入whichw3_make.若成功修改，则会出

6、现ww3完整目录。修改comp、link、switch进入bin文件夹，将comp.Intel和link.Intel复制成comp和link即：cpcomp.Intelcompcplink.Intellink使用vi编辑器修改comp和link。修改如下：comp中（大约96行处）修改成如下if$mpi_mod=yesthencomp=mpiifortwhichmpiifort1/dev/null2/dev/nullOK=$?if$OK!=0thencomp=mpif90fielsecomp=ifortfilink中（大约104行）修改成如下if$mpi_mod=yesthencomp=mpi

7、ifortwhichmpiifort1/dev/null2/dev/nullOK=$?if$OK!=0thencomp=mpif90fiopt=-O3-o$progelsecomp=ifortfiswitch中把SHRD改为MPIDIST（注意此处空格，否则会出错）。三个文件修改完毕后给予执行权限。chmod+ucomplinkswitch运行make_MPI直接输入命令make_MPI即可。如果没有错误显示则编译成功，一般会出现两个错误，但不影响模式的运行计算，到此ww3安装编译成功。据提取WW3对于输入数据文件的类型没有固定要求，本人常用的是风场为dat文件，也形为ASCII文件。那么输入

8、数据的格式有着一定的要求。对于风场来说，数据排列方式如下：Q1111111m丨11111111ViQ1I111AM丨111115口111111.&U1111111/12-1.96-2.-05-2.L4-2.21-2.27-2.18-2.39-2.34-2303-2.09-2.19-2,29-2.3&-2.40-2,45-2.47-2.4&-2454-2,10-2.22-2,34-2.42-2.46-2,se-2.54-2.55-2575-1.99-2.14-2.占-2.4-2.55-i.&l-267G-1.S8-2M-2.24-2.,42-2.51-2”晞-2.-2.73-27772.702.

9、&22.542.432.452-&52.422.552.6452.532.472.412.362.332.302.322.412.5092.352.30.381E722.052-052.es2.18-2.28111.971.9-31.,901.S61.891.911.932.18122-01506010600013-1.34-1.3-1,药-1.25-1,25-144-1.32-1,44-152-1.2B-1.23-1.L9-1.17-1.19-1.2-1.2.4-1.31-13B15-1.1-1.051&0-0.99-1.00

10、-1.-01-L.&3-1.09-1151S77-0.7370-0.63-.67-0.74B317-8.54-0.4541-0.37-0.33-Q.4&-e521S1.52L.B-31.741.671.591.881.451.591.&9191.721.6-31.541.461.391.321.31.371.43201.611.511.,421.341.2&1.20211.561471,371,291.191.080.9SB-97221.521斗孑1,331,241.W核吕20.79-74第一行是时间，如2015年6月1日0点0时0分，则写为：2015060100000

11、0（注意年月日和时分秒之间存在空格）接下来是两个矩阵，第一个是风场u，第二个是风场V。接下来是也形，说是也形，其实就是水深，大于等于0的设为无效值就好。当地行存为ASCII文件时，务必知道其绘出的地形图是什么形式的，在后面的修改脚本出需要用到。四、脚本修改1.ww3_grid.inp$WAVEWATCHIIIGridpreprocessorinputfile$Gridname(C*30,inquotes)$TESTGRID(GULFOFNOWHERE)(此处是使用的网格名称)$Frequencyincrementfactorandfirstfrequency(Hz)$numberoffrequ

12、encies(wavenumbers)anddirections,relativeoffset$offirstdirectionintermsofthedirectionalincrement-0.5,0.5.$Inversions1.18and2.22ofthemodelthisvaluewasbydefiniton0,$itisaddedtomitigatetheGSEforafirstorderscheme.Notethat$thisfactorisIGNOREDintheprintplotsinww3_outp.$1.10.0411825240.(频率间隔初始频率波速波向)$namel

13、ists.Anynamelistfoundinthefolowingsectionsuptothe$namelists.Anynamelistfoundinthefolowingsectionsuptothe$Setmodelflags$-FLDRYDryrun(input/outputonly,nocalculation).$-FLCX,FLCY$-$-FLCX,FLCY$-FLCTH,FLCK$-FLSOUActivatesourceterms.$-FLSOUActivatesourceterms.FTTTFT（上方六项的开关）$Settimesteps$-Timestepinformat

14、ion(thisinformationisalwaysread)$maximumglobaltimestep,maximumCFLtimestepforx-yand$k-theta,minimumsourcetermtimestep(allinseconds).$900.950.900.300.(第一个数和第三个数一般是第二个数的两倍，第二个数使用程序计算，最后一个一般为30)$Startofnamelistinputsection$StartingwithWAVEWATCHIIIversion2.00,thetunableparameters$forsourceterms,propagati

15、onschemes,andnumericsarereadusing$end-of-sectionidentifierstring(seebelow)istemporarilywritten$toww3_grid.scratch,andreadfromthereifnecessary.Namelists$notneededforthegivenswitchsettingswillbeskipped$automatically,andtheorderofthenamelistsisimmaterial.$Asanexample,namelistinputtochangeSWELLFandZWNDi

16、nthe$TolmanandChalikovinputwouldbe$&SIN2SWELLF=0.1,ZWND=15./(此处均属于对源函数的调整，手册有相关介绍)$Defineconstantsinsourceterms$Stresses$TC1996withcap:NamelistFLX3$CDMAX:MaximumallowedCD(cap)$CTYPE:Captype:$0:Discontinuous(default).$1:Hyperbolictangent.$Hwang2011:NamelistFLX4$(n.nn)$of$(!/STAB2).$CDFAC:re-scalingof

17、drag$Linearinput$CavaleriandM-R:NamelistSLN1$CLIN:Proportionalityconstant.$RFPM:FactorforfPMinfilter.$RFHF:Factorforfhinfilter.$Exponentialinput$WAM-3:NamelistSIN1$CINP:Proportionalityconstant.$TolmanandChalikov:NamelistSIN2ZWND:Heightofwind(m).SWELLF:swellfactorin(n.nn).STABSH,STABOF,CNEG,CPOS,FNEG

18、:c0,ST0,c1,c2andf1in.through(2.65)fordefinitioneffectivewindspeed$WAM4andvariants:NamelistSIN3$Charnockcoefficient$roughnessz0$Charnockcoefficient$roughnessz0$wind-wavecoupling$input$forgustiness$toreduceu_star$formulationALPHA0:minimumvalueofZ0MAX:maximumvalueofair-sideBETAMAX:maximumvalueofSINTHP:

19、powerofcosineinwindZALP:waveageshifttoaccountTAUWSHELTER:shelteringofshortwavesSWELLFPAR:choiceofswellattenuation(1:TC1996,3:ACC2008)SWELLF:swellattenuationfactor$ExtraparametersforSWELLFPAR=3only$SWELLF2,SWELLF3:swellattenuationfactors$SWELLF4:ThresholdReynoldsnumberforACC2008$SWELLF5:Relativevisco

20、usdecay$SWELLF5:Relativeviscousdecaybelowthreshold$meanflow$belowthreshold$meanflow$BYDRZinput$DiscreteI.A.:NamelistSNL1Z0RAT:roughnessforoscil.flow/:NamelistSIN6SINA0:factorfornegativeinputSINU10:windspeedscalingoption$Nonlinearinteractions$LAMBDA:Lambdainsourceterm.$NLPROP:Cinsourcterm.NOTE:defaul

21、t$valuedependsonothersource$termsselected.$KDCONV:FactorbeforekdinEq.(n.nn).$KDMIN,SNLCS1,SNLCS2,SNLCS3:$Minimumkd,andconstantsc1-3$indepthscalingfunction.$KDFS:Shallowwaterrelativefilter$KDFS:Shallowwaterrelativefilter$scaling$which$established.$Exactinteractions:NamelistSNL2IQTYPE:Typeofdepthtreat

22、ment:Deepwater:Deepwater/WAM:ShallowwaterTAILNL:Parametrictailpower.NDEPTH:NumberofdepthsinforintegrationspaceisUsedforIQTYPE=3only$DEPTHS:ArraywithdepthsforNDEPTH=3$Gen.MultipleDIA:NamelistSNL3$NQDEF:Numberofquadruplets.$MSC:Scalingconstantm.$NSC:ScalingconstantN.$KDFD:Deepwaterrelativefilter$Namel

23、istANL2depth,depth,NamelistANL3QPARMS:5xNQDEFparamatersdescribingthequadruplets,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=

24、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

25、,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

26、$CDIS,APM:Asinsourceterm.$TolmanandChalikov:NamelistSDS2$SDSA0,SDSA1,SDSA2,SDSB0,SDSB1,PHIMIN:$Constantsa0,a1,a2,b0,b1$Constantsa0,a1,a2,b0,b1$and$PHImin.$WAM4andvariants:NamelistSDS3$wavenumber$andtail$weights$dissipation$activationof$(SDSLF)and$spectrumSDSC1:WAM4CdscoeffientMNMEANP,WNMEANPTAIL:pow

27、erofformeandefinitionsinSdsSDSDELTA1,SDSDELTA2:relativeofkandkA2partsofWAM4SDSLF,SDSHF:coefficientforWAM4dissipationforunsaturatedsaturated(SDSHF)partsoftheSDSC2:Saturationdissipationcoefficient$SDSC4:ValueofB0=B/BrforwichSdsiszero$Sdsiszero$separationof$non-saturated$dissipation$Sds_SAT$integration

28、ofB$BYDRZ$normalizationspectra$T1andT2SDSBR:ThresholdBrforsaturationSDSP:powerof(B/Br-B0)inSdsSDSBR2:ThresholdBr2fortheWAM4dissipationinsaturatedandSDSC5:coefficientforturbulenceSDSC6:WeightfortheistropicpartofSDSDTH:Angularhalf-widthfor:NamelistSDS6SDSET:SelectthresholdSDSA1,SDSA2,SDSP1,SDSP2:Coeff

29、icientsfordissipationterms:NamelistSWL6SWLB1:Coefficientforswelldissipation$Bottomfriction$JONSWAP:NamelistSBT1$GAMMA:Asitsays.$Surfbreaking$BattjesandJanssen:NamelistSDB1$=1)$=0.73)$=1)$=0.73)$only(default)$MicheformulationBJALFA:Dissipationconstant(defaultBJGAM:Breakingthreshold(defaultBJFLAG:TRUE

30、-UseHmax/dratioFALSE-UseHmax/din$Triadnonlinearinteractions$LumpedTriadInteraction(LTA):NamelistSTR1(Tobeimplemented)$(default0.05)$whichinteractionimplemented)$(default0.05)$whichinteraction$computing$default10.)$biphase$computing$PTRIAD1:ProportionalitycoefficientPTRIAD2:MultipleofTm01uptoiscomput

31、ed(2.5)PTRIAD3:Ursellupperlimitforinteractions(notused,PTRIAD4:Shapeparameterforcomputation(0.2)PTRIAD5:Ursellnumbertresholdforinteractions(0.01)$Shorelinereflections$ref.parameters:NamelistREF1REFCOAST:ReflectioncoefficientatshorelineREFFREQ:Activationoffreq-dependentref.$REFMAP:Scalefactorforbotto

32、mslopemap$REFRMAX:maximumref.coeffient(default0.8)$REFFREQPOW:poweroffrequency$REFICEBERG:Reflectioncoefficientforicebergs$REFSUBGRID:Reflectioncoefficientforislands$REFCOSP_STRAIGHT:powerofcosineusedfor$straightshoreline$Bound2ndorderspectrumandfreeIG$IG1parameters$IG1parameters$Krasitskii-Janssen$

33、wavecorrection$IGMETHOD:1:Hasselmann,2:IGADDOUTP:activationofboundinww3_outp/ww3_ounpIGSOURCE:1:usesboundwaves,2:empirical$IGSTERMS:0:nosourceterminIGband$IGMAXFREQ:maximumfrequencyofIGband$IGEMPIRICAL:constantinempiricalfreeIGsource$IGSWELLMAX:activatesfreeIGsourcesforallfreq.$Propagationschemes$SM

34、Cgridpropagation: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.$

35、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

36、:MaximumCFLnumberforrefraction.$WDTHCG:Tuningfactorpropag.direction.$WDTHTH:Tuningfactornormaldirection.$NotethatUQandUNOschemeshavenotunableparameters.$Alltuneableparametersareassociatedwiththerefraction$limitationandtheGSEalleviation.$Unstructuredgrids$UNSTparameters:NamelistUNST$UGOBCAUTO:Automat

37、icdetectionofOBCpoints$UGOBCDEPTH:Threshold(1,J).Agridwithsimpleclosure:mayberectilinearorcurvilinear.$:i-indexandwrapsati=NX+1andhasclosureat$:j=NY+1.Inotherwords,(NX+1,J(1,J)$:and(I,NY+1)=(MOD(NX-I+1,NX)+1,NY).Tripole$:gridclosurerequiresthatNXbeeven.Agrid$:withtripoleclosuremustbecurvilinear.$2NX

38、,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$3U

39、nitnumberoffilewithx-coordinate.$Scalefactorandaddoffset:x=scale_fac*x_read+add_offset.IDLA,IDFM,formatforformattedread,FROMandfilename.IDLA:Layoutindicator:Readline-by-linebottomtotop.:Like1,singlereadstatement.:Readline-by-linetoptobottom.:Like3,singlereadstatement.IDFM:formatindicator:Freeformat.

40、:Fixedformatwithaboveformatdescriptor.:Unformatted.FROM:filetypeparameterUNIT:openfilebyunitnumberonly.NAME:openfilebynameandassigntounit.Iftheaboveunitnumberequals10,thenthex-coordisreadfromthis$file.Thex-coordmustfollowtheaboverecord.Nocommentlinesare$allowedwithinthex-coordinput.$4Unitnumberoffil

41、ewithy-coordinate.$Scalefactorandaddoffset:y2maybemoreefficient$formassivelyparallelcomputations.OnlyIOSTYP=0requiresatrue$parallelfilesystemlikeGPFS.$IOSTYP=0:Nodataserverprocesses,directaccessoutputfromeachprocess(requiresetrueparallelfilesystem).$1:Nodataserverprocess.Alloutputforeachtype$computa

42、tionstoo.$2:Lastprocessisreservedforalloutput,anddoesno$computing.$3:Multiplededicatedoutputprocesses.$2般改为1$Fiveoutputtypesareavailable(seebelow).Alloutputtypesshare$asimilarformatforthefirstinputline:$-firsttimeinyyyymmddhhmmssformat,outputinterval(s),and$lasttimeinyyyymmddhhmmssformat(allintegers

43、).$Outputisdisabledbysettingtheoutputintervalto0.$Type1:Fieldsofmeanwaveparameters$Standardlineandlinewithlogicalflagstoactivateoutput$fieldsasdefinedinsection2.4ofthemanual.Thelogical$flagsarenotsuppliedifnooutputisrequested.Thelogical$flagscanbeplacedonmultipleconsecutivelines.However,thetotalnumb

44、erandorderofthelogicalflagsisfixed.Therawdatafileisout_grd.ww3,seew3iogo.ftnforadditionaldoc.19680606000000360019680608000000修改为你的起始和终止时间，3600是时间步长，即相邻的两个数据时间间隔转换成秒$Outputrequestflagsidentifyingfields.$Thetablebelowprovidesafulldefinitionoffieldoutputparameters$aswellasflagsindicatingiftheyareavaila

45、bleindifferentfield$outputoutputfiletypes(ASCII,grib,NetCDF).$Furtherdefinitionsarefoundinsection2.4ofthemanual.$Selectionoffieldoutputsmaybemadeintwoways:$F/Tflags:firstflagissettoF,requestsmadepergroup(1stline)$followedbyparameterflags(totalof10groups).$Namelists:firstlineissettoN,nextlinecontains

46、parameter$symbolaspertablebelow.$ExampleofF/Tflaguseisgiveninthissampleww3_shel.inp,below.$Fornamelistusage,seethesampleww3_ounf.inpforanexample.$Outputfieldparameterdefinitionstable$Allparameterslistedbelowareavailableinoutputfileofthetypes$ASCIIandNetCDF.Ifselectedoutputfiletypesaregradsorgrib,$so

47、meparametersmaynotbeavailable.Thefirsttwocolumnsinthe$tablebelowidentifysuchcasesbyflags,cols1(GRB)and2(GXO)$refertogrib(ww3_grib)andgrads(gx_outf),respectively.$Columns3and4providegroupandparameternumberspergroup.$Columns5,6and7provide:$5-codename(internal)$6-outputtags(namesusedisASCIIfileextensio

48、ns,NetCDFvariablenamesandnamelist-basedselection(seeww3_ounf.inp)$7-Longparametername/definition$GG$RXGrpParamCodeOutputParameter/Group$BONumbNumbrNameTagDefinition$1ForcingFields$TT11DWDPTWaterdepth.$TT12CX,YCURCurrentvelocity.$TT13UAWNDWindspeed.$TT14ASASTAir-seatemperaturedifference.$TT15WLVWLVWa

49、terlevels.$TT16ICEICEIceconcentration.$TT17IBGIBGIceberg-induceddamping$TT18D50D50Mediansedimentgrainsize$2Standardmeanwave$2Parameters$-3F31Ef$F$FF32th1ma1,b2$FF33sth1ma1,b2$FF34th2ma2,b2$FF35sth2m$TT21HSHSWaveheight.$TT22WLMLMMeanwavelength.$TT23T02T02Meanwaveperiod(Tm02).$TT24T0M1T0M1Meanwaveperi

50、od(Tm0,-1).$TT25T01T01Meanwaveperiod(Tm01).$TT26FP0FPPeakfrequency.$TT27THMDIRMeanwavedirection.$TT28THSSPRMeandirectionalspread.$TT29THP0DPPeakdirection.$TT210HIGHIGInfragravityheight$SpectralParameters(first5)EFWavefrequencyspectrumTH1MMeanwavedirectionfromSTH1MDirectionalspreadingfromMeanwavedire

51、ctionfromDirectionalspreadingfroma2,b2$FF36WNWNWavenumberarray$4SpectralPartitionParameters$TT41PHSPHSPartitionedwaveheights.$TT42PTPPTPPartitionedpeakperiod.$TT43PLPPLPPartitionedpeakwavelength.$TT44PDIRPDIRPartitionedmeandirection.$TT45PSIPSPRPartitionedmeandirectionalspread.$TT46PWSPWSPartitionne

52、dwindseafraction.$TT47PWSTTWSTotalwindseafraction.$TT48PNRPNRNumberofpartitions.$5Atmosphere-waveslayer$TT51USTUSTFrictionvelocity.$FT52CHARNCHACharnockparameter$FT53CGECGEEnergyflux$FT54PHIAWFAWAir-seaenergyflux$FT55TAUWIX,YTAWNetwave-supportedstress$FT56TAUWNX,YTWANegativepartofthewave-supportedst

53、ress$FF57WHITECAPWCCWhitecapcoverage$FF58WHITECAPWCFWhitecapthickness$FF59WHITECAPWCHMeanbreakingheight$FF510WHITECAPWCMWhitecapmoment$6Wave-oceanlayer$FF61SXX,YY,XYSXYRadiationstresses.$FF62TAUOX,YTWOWavetooceanmomentumflux$FF63BHDBHDBernoullihead(Jterm)$FF64PHIOCFOCWavetooceanenergyflux$FF65TUSX,Y

54、TUSStokestransport$FF66USSX,YUSSSurfaceStokesdrift$FF67PR,TPMSP2SSecond-ordersumpressure$FF68US3DUSFSpectrumofsurfaceStokesdrift$FF69P2SMSP2LMicroseismsourceterm$7Wave-bottomlayer$FF71ABAABRNearbottomrmsamplitides.$FF72UBAUBRNearbottomrmsvelocities.$FF73BEDFORMSBEDBedforms$FF74PHIBBLFBBEnergyfluxdue

55、tobottomfriction$FF75TAUBBLTBBMomentumfluxduetobottomfriction$8Spectrumparameters$FF81MSSX,YMSSMeansquareslopes$FF82MSCX,YMSCSpectrallevelathighfrequencytail$9Numericaldiagnostics$TT91DTDYNDTDAveragetimestepinintegration.$TT92FCUTFCCut-offfrequency.$TT93CFLXYMAXCFXMax.CFLnumberforspatialadvection.$T

56、T94CFLTHMAXCFDMax.CFLnumberfortheta-advection.$FF95CFLKMAXCFKMax.CFLnumberfork-advection.$10Userdefined$FF101U1Userdefined#1.(requirescoding.)$FF102U2Userdefined#1.(requirescoding.)$Section4consistofasetoffields,index0=windsea,index$1:NOSWLLarefirstNOSWLLswellfields.$Actualactiveparameterselectionse

57、ction$以下均为上述介绍的输出参数，自选开关$(1)ForcingFields$DPTCURWNDASTWLVICEIBGD50TTTTTFFF$(2)StandardmeanwaveParametersT$HSLMT02T0M1T01FPDIRSPRDPTTTTTTTTT$(3)Frequency-dependentparametersT$EFTH1MSTH1MTH2MSTH2MWNTTTFFF$(4)SpectralPartitionParametersT$PHSPTPPLPPDIRPSPRPWSTWSPNRTTTTTTTT$(5)Atmosphere-waveslayerT$USTC

58、HACGEFAWTAWTWAWCCWCFWCHWCMTTTTTTTTTT$(6)Wave-OceanlayerT$SXYTWOBHDFOCTUSUSSP2SUSFP2LTTTTTTTFFF$(7)Wave-bottomlayerT$ABRUBRBEDFBBTBBTTTTT$(8)SpectrumparametersT$MSSMSCTT$(9)NumericaldiagnosticsT$DTDFCCFXCFDCFKTTTTT$(10)Userdefined(NOEXTRflagsneeded)F$U1U2$TT$Type2:PointoutputStandardlineandanumberofl

59、inesidentifyingthe$longitude,latitudeandname(C*10)ofoutputpoints.$longitude,latitudeandname(C*10)ofoutputpoints.$Thelistisclosedbydefiningapointwiththename$STOPSTRING.Nopointinforeadifnopointoutputis$requested(i.e.,noSTOPSTRINGneeded).$Exampleforsphericalgrid.$Therawdatafileisout_pnt.ww3,$seew3iogo.

60、ftnforadditionaldoc.$NOTE:Spacesmaybeincludedinthename,butthisisnot$advised,becauseitwillbreaktheGrADSutilityto$plotsspectraandsourceterms,andwillmakeitmore$diffuculttousepointnamesindatafiles.$1968060600000090019680608000000此处与上一个时间相同$(此处针对点输出，如若不需要，则加上$即可，两个数据各位xy坐标)-0.25-0.25Land0.00.0Point_12.01