$ -------------------------------------------------------------------- $ $ WAVEWATCH III Grid preprocessor input file $ $ -------------------------------------------------------------------- $ $ Grid name (C*30, in quotes) $ 'TEST GRID (GULF OF NOWHERE) ' $ $ Frequency increment factor and first frequency (Hz) ---------------- $ $ number of frequencies (wavenumbers) and directions, relative offset $ of first direction in terms of the directional increment [-0.5,0.5]. $ In versions 1.18 and 2.22 of the model this value was by definiton 0, $ it is added to mitigate the GSE for a first order scheme. Note that $ this factor is IGNORED in the print plots in ww3_outp. $ 1.1 0.04118 25 24 0. $ $ Set model flags ---------------------------------------------------- $ $ - FLDRY Dry run (input/output only, no calculation). $ - FLCX, FLCY Activate X and Y component of propagation. $ - FLCTH, FLCK Activate direction and wavenumber shifts. $ - FLSOU Activate source terms. $ F T T T F T $ $ Set time steps ----------------------------------------------------- $ $ - Time step information (this information is always read) $ maximum global time step, maximum CFL time step for x-y and $ k-theta, minimum source term time step (all in seconds). $ 900. 950. 900. 300. $ $ Start of namelist input section ------------------------------------ $ $ Starting with WAVEWATCH III version 2.00, the tunable parameters $ for source terms, propagation schemes, and numerics are read using $ namelists. Any namelist found in the folowing sections up to the $ end-of-section identifier string (see below) is temporarily written $ to ww3_grid.scratch, and read from there if necessary. Namelists $ not needed for the given switch settings will be skipped $ automatically, and the order of the namelists is immaterial. $ As an example, namelist input to change SWELLF and ZWND in the $ Tolman and Chalikov input would be $ $ &SIN2 SWELLF = 0.1, ZWND = 15. / $ $ Define constants in source terms ----------------------------------- $ $ $ Stresses - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - $ TC 1996 with cap : Namelist FLX3 $ CDMAX : Maximum allowed CD (cap) $ CTYPE : Cap type : $ 0: Discontinuous (default). $ 1: Hyperbolic tangent. $ Hwang 2011 : Namelist FLX4 $ CDFAC : re-scaling of drag $ $ Linear input - - - - - - - - - - - - - - - - - - - - - - - - - - - - $ Cavaleri and M-R : Namelist SLN1 $ CLIN : Proportionality constant. $ RFPM : Factor for fPM in filter. $ RFHF : Factor for fh in filter. $ $ Exponential input - - - - - - - - - - - - - - - - - - - - - - - - - $ WAM-3 : Namelist SIN1 $ CINP : Proportionality constant. $ $ Tolman and Chalikov : Namelist SIN2 $ ZWND : Height of wind (m). $ SWELLF : swell factor in (n.nn). $ STABSH, STABOF, CNEG, CPOS, FNEG : $ c0, ST0, c1, c2 and f1 in . (n.nn) $ through (2.65) for definition of $ effective wind speed (!/STAB2). $ WAM4 and variants : Namelist SIN3 $ ZWND : Height of wind (m). $ ALPHA0 : minimum value of Charnock coefficient $ Z0MAX : maximum value of air-side roughness z0 $ BETAMAX : maximum value of wind-wave coupling $ SINTHP : power of cosine in wind input $ ZALP : wave age shift to account for gustiness $ TAUWSHELTER : sheltering of short waves to reduce u_star $ SWELLFPAR : choice of swell attenuation formulation $ (1: TC 1996, 3: ACC 2008) $ SWELLF : swell attenuation factor $ Extra parameters for SWELLFPAR=3 only $ SWELLF2, SWELLF3 : swell attenuation factors $ SWELLF4 : Threshold Reynolds number for ACC2008 $ SWELLF5 : Relative viscous decay below threshold $ Z0RAT : roughness for oscil. flow / mean flow $ BYDRZ input : Namelist SIN6 $ SINA0 : factor for negative input $ SINU10 : wind speed scaling option $ $ Nonlinear interactions - - - - - - - - - - - - - - - - - - - - - - - $ Discrete I.A. : Namelist SNL1 $ LAMBDA : Lambda in source term. $ NLPROP : C in sourc term. NOTE : default $ value depends on other source $ terms selected. $ KDCONV : Factor before kd in Eq. (n.nn). $ KDMIN, SNLCS1, SNLCS2, SNLCS3 : $ Minimum kd, and constants c1-3 $ in depth scaling function. $ Exact interactions : Namelist SNL2 $ IQTYPE : Type of depth treatment $ 1 : Deep water $ 2 : Deep water / WAM scaling $ 3 : Shallow water $ TAILNL : Parametric tail power. $ NDEPTH : Number of depths in for which $ integration space is established. $ Used for IQTYPE = 3 only $ Namelist ANL2 $ DEPTHS : Array with depths for NDEPTH = 3 $ Gen. Multiple DIA : Namelist SNL3 $ NQDEF : Number of quadruplets. $ MSC : Scaling constant 'm'. $ NSC : Scaling constant 'N'. $ KDFD : Deep water relative filter depth, $ KDFS : Shallow water relative filter depth, $ Namelist ANL3 $ QPARMS : 5 x NQDEF paramaters describing the $ quadruplets, repeating LAMBDA, MU, DT12. $ Cdeep and Cshal. See examples below. $ Two Scale Approx. : Namelist SNL4 $ INDTSA : Index for TSA/FBI computations $ (0 = FBI ; 1 = TSA) $ ALTLP : Index for alternate looping $ (1 = no ; 2 = yes) $ $ Traditional DIA setup (default): $ $ &SNL3 NQDEF = 1, MSC = 0.00, NSC = -3.50 / $ &ANL3 QPARMS = 0.250, 0.000, -1.0, 0.1000E+08, 0.0000E+00 / $ $ GMD3 from 2010 report (G13d in later paper) : $ $ &SNL3 NQDEF = 3, MSC = 0.00, NSC = -3.50 / $ &ANL3 QPARMS = 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 / $ $ G35d from 2010 report: $ $ &SNL3 NQDEF = 5, MSC = 0.00, NSC = -3.50 / $ &ANL3 QPARMS = 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 / $ $ Nonlinear filter based on DIA - - - - - - - - - - - - - - - - - - - $ Namelist SNLS $ A34 : Relative offset in quadruplet $ FHFC : Proportionality constants. $ DMN : Maximum relative change. $ FC1-3 : Constants in frequency filter. $ $ Whitecapping dissipation - - - - - - - - - - - - - - - - - - - - - $ WAM-3 : Namelist SDS1 $ CDIS, APM : As in source term. $ $ Tolman and Chalikov : Namelist SDS2 $ SDSA0, SDSA1, SDSA2, SDSB0, SDSB1, PHIMIN : $ Constants a0, a1, a2, b0, b1 and $ PHImin. $ $ WAM4 and variants : Namelist SDS3 $ SDSC1 : WAM4 Cds coeffient $ MNMEANP, WNMEANPTAIL : power of wavenumber $ for mean definitions in Sds and tail $ SDSDELTA1, SDSDELTA2 : relative weights $ of k and k^2 parts of WAM4 dissipation $ SDSLF, SDSHF : coefficient for activation of $ WAM4 dissipation for unsaturated (SDSLF) and $ saturated (SDSHF) parts of the spectrum $ SDSC2 : Saturation dissipation coefficient $ SDSC4 : Value of B0=B/Br for wich Sds is zero $ SDSBR : Threshold Br for saturation $ SDSP : power of (B/Br-B0) in Sds $ SDSBR2 : Threshold Br2 for the separation of $ WAM4 dissipation in saturated and non-saturated $ SDSC5 : coefficient for turbulence dissipation $ SDSC6 : Weight for the istropic part of Sds_SAT $ SDSDTH: Angular half-width for integration of B $ $ BYDRZ : Namelist SDS6 $ SDSET : Select threshold normalization spectra $ SDSA1, SDSA2, SDSP1, SDSP2 : $ Coefficients for dissipation terms T1 and T2 $ : Namelist SWL6 $ SWLB1 : Coefficient for swell dissipation $ $ Bottom friction - - - - - - - - - - - - - - - - - - - - - - - - - - $ JONSWAP : Namelist SBT1 $ GAMMA : Bottom friction emprical constant $ $ $ Surf breaking - - - - - - - - - - - - - - - - - - - - - - - - - - - $ Battjes and Janssen : Namelist SDB1 $ BJALFA : Dissipation constant (default = 1) $ BJGAM : Breaking threshold (default = 0.73) $ BJFLAG : TRUE - Use Hmax/d ratio only (default) $ FALSE - Use Hmax/d in Miche formulation $ $ Dissipation in the ice - - - - - - - - - - - - - - - - - - - - - - $ Generalization of Liu et al. : Namelist SIC2 $ IC2DISPER : If true uses Liu formulation with eddy viscosity $ If false, uses the generalization with turbulent $ to laminar transition $ IC2TURB : empirical factor for the turbulent part $ IC2ROUGH : under-ice roughness length $ IC2REYNOLDS: Re number for laminar to turbulent transition $ IC2SMOOTH : smoothing of transition reprensenting random waves $ IC2VISC : empirical factor for viscous part $ $ $ Scattering in the ice & creep dissipations- - - - - - - - - - - - - $ Generalization of Wiliams et al. : Namelist SIS2 $ ISC1 : scattering coefficient (default = 1) $ IS2BACKSCAT : fraction of energy back-scattered (default = 1 ) $ IS2BREAK : TRUE - changes floe max diameter $ : FALSE - does not change floe max diameter $ IS2C1 : scattering in pack ice $ IS2C2 : frequency dependance of scattering in pack ice $ IS2C3 : frequency dependance of scattering in pack ice $ ISBACKSCAT : fraction of scattered energy actualy redistributed $ IS2DISP : use of ice-specific dispersion relation (T/F) $ FRAGILITY : parameter between 0 and 1 that gives the shape of FSD $ IS2DMIN : minimum floe diameter in meters $ IS2DAMP : multiplicative coefficient for dissipation term from RP $ IS2UPDATE : TRUE - updates the max floe diameter with forcing only $ : FALSE - updates the max floe diameter at each time step $ $ Triad nonlinear interactions - - - - - - - - - - - - - - - - - - - - $ Lumped Triad Interaction (LTA) : Namelist STR1 (To be implemented) $ PTRIAD1 : Proportionality coefficient (default 0.05) $ PTRIAD2 : Multiple of Tm01 up to which interaction $ is computed (2.5) $ PTRIAD3 : Ursell upper limit for computing $ interactions (not used, default 10.) $ PTRIAD4 : Shape parameter for biphase $ computation (0.2) $ PTRIAD5 : Ursell number treshold for computing $ interactions (0.01) $ $ Shoreline reflections - - - - - - - - - - - - - - - - - - - - - - - - $ ref. parameters : Namelist REF1 $ REFCOAST : Reflection coefficient at shoreline $ REFFREQ : Activation of freq-dependent ref. $ REFMAP : Scale factor for bottom slope map $ REFRMAX : maximum ref. coeffient (default 0.8) $ REFFREQPOW: power of frequency $ REFICEBERG: Reflection coefficient for icebergs $ REFSUBGRID: Reflection coefficient for islands $ REFCOSP_STRAIGHT: power of cosine used for $ straight shoreline $ $ Bound 2nd order spectrum and free IG - - - - - - - - - - - - - - - - - $ IG1 parameters : Namelist SIG1 $ IGMETHOD : 1: Hasselmann, 2: Krasitskii-Janssen $ IGADDOUTP : activation of bound wave correction $ in ww3_outp / ww3_ounp $ IGSOURCE : 1: uses bound waves, 2: empirical $ IGSTERMS : > 0 : no source term in IG band $ IGMAXFREQ : maximum frequency of IG band $ IGEMPIRICAL: constant in empirical free IG source $ IGBCOVERWRITE: T: Replaces IG spectrum, does not add $ IGSWELLMAX: T: activates free IG sources for all freq. $ $ $ Propagation schemes ------------------------------------------------ $ $ First order : Namelist PRO1 $ CFLTM : Maximum CFL number for refraction. $ $ UQ/UNO with diffusion : Namelist PRO2 $ CFLTM : Maximum CFL number for refraction. $ DTIME : Swell age (s) in garden sprinkler $ correction. If 0., all diffusion $ switched off. If small non-zero $ (DEFAULT !!!) only wave growth $ diffusion. $ LATMIN : Maximum latitude used in calc. of $ strength of diffusion for prop. $ $ UQ/UNO with averaging : Namelist PRO3 $ CFLTM : Maximum CFL number for refraction. $ WDTHCG : Tuning factor propag. direction. $ WDTHTH : Tuning factor normal direction. $ $ Note that UQ and UNO schemes have no tunable parameters. $ All tuneable parameters are associated with the refraction $ limitation and the GSE alleviation. $ $ Unstructured grids ------------------------------------------------ $ $ UNST parameters : Namelist UNST $ UGOBCAUTO : TRUE: OBC points are taken from type 15 elements $ FALSE: OBC points must be listed in ww3_grid.inp $ UGOBCDEPTH: Threshold ( < 0) depth for OBC points if UGOBCAUTO is TRUE $ EXPFSN : Activation of N scheme $ EXPFSPSI : Activation of PSI scheme $ EXPFSFCT : Activation of FCT scheme $ IMPFSN : Activation of N implicit scheme $ $ SMC grid propagation : Namelist PSMC and default values $ CFLTM : Maximum CFL no. for propagation, 0.7 $ DTIME : Swell age for diffusion term (s), 0.0 $ LATMIN : Maximum latitude (deg) for GCT, 86.0 $ RFMAXD : Maximum refraction turning (deg), 80.0 $ LvSMC : No. of refinement level, default 1 $ ISHFT : Shift number of i-index, default 0 $ JEQT : Shift number of j-index, default 0 $ NBISMC : No. of input boundary points, 0 $ UNO3 : Use 3rd order advection scheme, .FALSE. $ UNO3 : Add extra spatial averaging, .FALSE. $ &PSMC DTIME = 39600.0, LATMIN=85.0, RFMAXD = 36.0, LvSMC=3, JEQT=1344 / $ $ Output of 3D arrays------------------------------------------------- $ $ In order to limit the use of memory, arrays for 3D output fiels (i.e. $ variables that are a function of both space and frequency, are not $ declared, and thus cannot be used, unless specified by namelists. $ NB: Output of 'first 5' moments E, th1m, sth1m, th2, sth2m allows to estimate the full $ directional spectrum using, e.g. MEM (Lygre&Krogstad 1986). $ $ Parameters (integers) : Namelist OUTS $ For the frequency spectrum E(f) $ E3D : <=0: not declared, > 0: declared $ I1E3D : First frequency index of output (default is 1) $ I2E3D : Last frequency index of output (default is NK) $ For the mean direction th1m(f), and spread sth1m(f) $ TH1MF, STH1MF : <=0: not declared, > 0: declared $ I1TH1MF, I1STH1MF: First frequency index of output (default is 1) $ I2TH1MF, I2STH1MF: First frequency index of output (default is 1) $ For the mean direction th2m(f), and spread sth2m(f) $ TH2MF, STH2MF : <=0: not declared, > 0: declared $ I1TH2MF, I1STH2MF: First frequency index of output (default is 1) $ I2TH2MF, I2STH2MF: First frequency index of output (default is 1) $ For 2nd order pressure at K=0 (source of microseisms & microbaroms) $ P2SF : <=0: not declared, > 0: declared $ I1P2SF : First frequency index of output (default is 1) $ I2P2SF : Last frequency index of output (default is NK) $ $ Miscellaneous ------------------------------------------------------ $ $ Misc. parameters : Namelist MISC $ CICE0 : Ice concentration cut-off. $ CICEN : Ice concentration cut-off. $ PMOVE : Power p in GSE aleviation for $ moving grids in Eq. (D.4). $ XSEED : Xseed in seeding alg. (!/SEED). $ FLAGTR : Indicating presence and type of $ subgrid information : $ 0 : No subgrid information. $ 1 : Transparancies at cell boun- $ daries between grid points. $ 2 : Transp. at cell centers. $ 3 : Like 1 with cont. ice. $ 4 : Like 2 with cont. ice. $ XP, XR, XFILT $ Xp, Xr and Xf for the dynamic $ integration scheme. $ IHMAX : Number of discrete levels in part. $ HSPMIN : Minimum Hs in partitioning. $ WSM : Wind speed multiplier in part. $ WSC : Cut of wind sea fraction for $ identifying wind sea in part. $ FLC : Flag for combining wind seas in $ partitioning. $ NOSW : Number of partitioned swell fields $ in field output. $ FMICHE : Constant in Miche limiter. $ STDX : Space-Time Extremes X-Length $ STDY : Space-Time Extremes Y-Length $ STDT : Space-Time Extremes Duration $ P2SF : ...... $ $ Diagnostic Sea-state Dependent Stress- - - - - - - - - - - - - - - - - $ Reichl et al. 2014 : Namelist FLD1 $ TAILTYPE : High Frequency Tail Method $ 0: Constant value (prescribed) $ 1: Wind speed dependent $ (Based on GFDL Hurricane $ Model Z0 relationship) $ TAILLEV : Level of high frequency tail $ (if TAILTYPE==0) $ Valid choices: $ Capped min: 0.001, max: 0.02 $ TAILT1 : Tail transition ratio 1 $ TAILT1*peak input frequency $ is the first transition point of $ the saturation specturm $ Default is 1.25 $ TAILT1 : Tail transition ratio 2 $ TAILT2*peak input frequency $ is the second transition point of $ the saturation specturm $ Default is 3.00 $ Donelan et al. 2012 : Namelist FLD2 $ TAILTYPE : See above (FLD1) $ TAILLEV : See above (FLD1) $ TAILT1 : See above (FLD1) $ TAILT2 : See above (FLD1) $ $ In the 'Out of the box' test setup we run with sub-grid obstacles $ and with continuous ice treatment. $ &MISC CICE0 = 0.25, CICEN = 0.75, FLAGTR = 4 / &FLX3 CDMAX = 3.5E-3 , CTYPE = 0 / $ &SDB1 BJGAM = 1.26, BJFLAG = .FALSE. / $ $ Mandatory string to identify end of namelist input section. $ END OF NAMELISTS $ $ Define grid -------------------------------------------------------- $ $ $ Five records containing : $ $ 1 Type of grid, coordinate system and type of closure: GSTRG, FLAGLL, $ CSTRG. Grid closure can only be applied in spherical coordinates. $ GSTRG : String indicating type of grid : $ 'RECT' : rectilinear $ 'CURV' : curvilinear $ 'UNST' : unstructured (triangle-based) $ FLAGLL : Flag to indicate coordinate system : $ T : Spherical (lon/lat in degrees) $ F : Cartesian (meters) $ CSTRG : String indicating the type of grid index space closure : $ 'NONE' : No closure is applied $ 'SMPL' : Simple grid closure : Grid is periodic in the $ : i-index and wraps at i=NX+1. In other words, $ : (NX+1,J) => (1,J). A grid with simple closure $ : may be rectilinear or curvilinear. $ 'TRPL' : Tripole grid closure : Grid is periodic in the $ : i-index and wraps at i=NX+1 and has closure at $ : j=NY+1. In other words, (NX+1,J<=NY) => (1,J) $ : and (I,NY+1) => (NX-I+1,NY). Tripole $ : grid closure requires that NX be even. A grid $ : with tripole closure must be curvilinear. $ 2 NX, NY. As the outer grid lines are always defined as land $ points, the minimum size is 3x3. $ $ Branch here based on grid type $ $ IF ( RECTILINEAR GRID ) THEN $ $ 3 Grid increments SX, SY (degr.or m) and scaling (division) factor. $ If CSTRG='SMPL', then SX is set to 360/NX. $ 4 Coordinates of (1,1) (degr.) and scaling (division) factor. $ $ ELSE IF ( CURVILINEAR GRID ) THEN $ $ 3 Unit number of file with x-coordinate. $ Scale factor and add offset: x <= scale_fac * x_read + add_offset. $ IDLA, IDFM, format for formatted read, FROM and filename. $ IDLA : Layout indicator : $ 1 : Read line-by-line bottom to top. $ 2 : Like 1, single read statement. $ 3 : Read line-by-line top to bottom. $ 4 : Like 3, single read statement. $ IDFM : format indicator : $ 1 : Free format. $ 2 : Fixed format with above format descriptor. $ 3 : Unformatted. $ FROM : file type parameter $ 'UNIT' : open file by unit number only. $ 'NAME' : open file by name and assign to unit. $ $ If the above unit number equals 10, then the x-coord is read from this $ file. The x-coord must follow the above record. No comment lines are $ allowed within the x-coord input. $ $ 4 Unit number of file with y-coordinate. $ Scale factor and add offset: y <= scale_fac * y_read + add_offset. $ IDLA, IDFM, format for formatted read, FROM and filename. $ IDLA : Layout indicator : $ 1 : Read line-by-line bottom to top. $ 2 : Like 1, single read statement. $ 3 : Read line-by-line top to bottom. $ 4 : Like 3, single read statement. $ IDFM : format indicator : $ 1 : Free format. $ 2 : Fixed format with above format descriptor. $ 3 : Unformatted. $ FROM : file type parameter $ 'UNIT' : open file by unit number only. $ 'NAME' : open file by name and assign to unit. $ $ If the above unit number equals 10, then the y-coord is read from this $ file. The y-coord must follow the above record. No comment lines are $ allowed within the y-coord input. $ $ ELSE IF ( UNSTRUCTURED GRID ) THEN $ Nothing to declare: all the data will be read from the GMESH file $ END IF ( CURVILINEAR GRID ) $ $ 5 Limiting bottom depth (m) to discriminate between land and sea $ points, minimum water depth (m) as allowed in model, unit number $ of file with bottom depths, scale factor for bottom depths (mult.), $ IDLA, IDFM, format for formatted read, FROM and filename. $ IDLA : Layout indicator : $ 1 : Read line-by-line bottom to top. $ 2 : Like 1, single read statement. $ 3 : Read line-by-line top to bottom. $ 4 : Like 3, single read statement. $ IDFM : format indicator : $ 1 : Free format. $ 2 : Fixed format with above format descriptor. $ 3 : Unformatted. $ FROM : file type parameter $ 'UNIT' : open file by unit number only. $ 'NAME' : open file by name and assign to unit. $ $ If the above unit number equals 10, then the bottom depths are read from $ this file. The depths must follow the above record. No comment lines are $ allowed within the depth input. In the case of unstructured grids, the file $ is expected to be a GMESH grid file containing node and element lists. $ $ ------------------------------------------------------------------------ $ Example for rectilinear grid with spherical (lon/lat) coordinate system. $ Note that for Cartesian coordinates the unit is meters (NOT km). $ 'RECT' T 'NONE' 12 12 1. 1. 4. -1. -1. 4. -0.1 2.50 10 -10. 3 1 '(....)' 'NAME' 'bottom.inp' $ 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 5 4 2 0 2 4 5 6 6 6 6 6 5 4 2 0 2 4 5 6 6 6 6 6 5 4 2 0 2 4 5 6 6 6 6 6 5 4 2 0 0 4 5 6 6 6 6 6 5 4 4 2 2 4 5 6 6 6 6 6 6 5 5 4 4 5 6 6 6 6 6 6 6 6 6 5 5 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 $ $ ------------------------------------------------------------------------ $ Example for curvilinear grid with spherical (lon/lat) coordinate system. $ Same spatial grid as preceding rectilinear example. $ Note that for Cartesian coordinates the unit is meters (NOT km). $ $ 'CURV' T 'NONE' $ 12 12 $ $ 10 0.25 -0.5 3 1 '(....)' 'NAME' 'x.inp' $ $ 1 2 3 4 5 6 7 8 9 10 11 12 $ 1 2 3 4 5 6 7 8 9 10 11 12 $ 1 2 3 4 5 6 7 8 9 10 11 12 $ 1 2 3 4 5 6 7 8 9 10 11 12 $ 1 2 3 4 5 6 7 8 9 10 11 12 $ 1 2 3 4 5 6 7 8 9 10 11 12 $ 1 2 3 4 5 6 7 8 9 10 11 12 $ 1 2 3 4 5 6 7 8 9 10 11 12 $ 1 2 3 4 5 6 7 8 9 10 11 12 $ 1 2 3 4 5 6 7 8 9 10 11 12 $ 1 2 3 4 5 6 7 8 9 10 11 12 $ 1 2 3 4 5 6 7 8 9 10 11 12 $ $ 10 0.25 0.5 3 1 '(....)' 'NAME' 'y.inp' $ $ 1 1 1 1 1 1 1 1 1 1 1 1 $ 2 2 2 2 2 2 2 2 2 2 2 2 $ 3 3 3 3 3 3 3 3 3 3 3 3 $ 4 4 4 4 4 4 4 4 4 4 4 4 $ 5 5 5 5 5 5 5 5 5 5 5 5 $ 6 6 6 6 6 6 6 6 6 6 6 6 $ 7 7 7 7 7 7 7 7 7 7 7 7 $ 8 8 8 8 8 8 8 8 8 8 8 8 $ 9 9 9 9 9 9 9 9 9 9 9 9 $ 10 10 10 10 10 10 10 10 10 10 10 10 $ 11 11 11 11 11 11 11 11 11 11 11 11 $ 12 12 12 12 12 12 12 12 12 12 12 12 $ $ -0.1 2.50 10 -10. 3 1 '(....)' 'NAME' 'bottom.inp' $ $ 6 6 6 6 6 6 6 6 6 6 6 6 $ 6 6 6 5 4 2 0 2 4 5 6 6 $ 6 6 6 5 4 2 0 2 4 5 6 6 $ 6 6 6 5 4 2 0 2 4 5 6 6 $ 6 6 6 5 4 2 0 0 4 5 6 6 $ 6 6 6 5 4 4 2 2 4 5 6 6 $ 6 6 6 6 5 5 4 4 5 6 6 6 $ 6 6 6 6 6 6 5 5 6 6 6 6 $ 6 6 6 6 6 6 6 6 6 6 6 6 $ 6 6 6 6 6 6 6 6 6 6 6 6 $ 6 6 6 6 6 6 6 6 6 6 6 6 $ 6 6 6 6 6 6 6 6 6 6 6 6 $ $ ------------------------------------------------------------- $ SMC grid use the same spherical lat-lon grid parameters $ 'RECT' T 'SMPL' $ 1024 704 $ SMC grid base level resolution dlon dlat and start lon lat $ 0.35156250 0.23437500 1. $ 0.17578125 -78.6328125 1. $ $ Normal depth input line is used to passing the minimum depth $ though the depth file is not read for SMC grid. $ -0.1 10.0 30 -1. 1 1 '(....)' 'NAME' 'SMC25Depth.dat' $ SMC cell and face arrays and obstruction ratio: $ 32 1 1 '(....)' 'S6125MCels.dat' $ 33 1 1 '(....)' 'S6125ISide.dat' $ 34 1 1 '(....)' 'S6125JSide.dat' $ 31 1.0 1 1 '(...)' 'NAME' 'SMC25Subtr.dat' $ The input boundary cell file is only needed when NBISMC > 0. $ 35 1 1 '(....)' 'S6125Bundy.dat' $ Extra cell and face arrays for Arctic part if ARC is selected. $ 36 1 1 '(....)' 'S6125MBArc.dat' $ 37 1 1 '(....)' 'S6125AISid.dat' $ 38 1 1 '(....)' 'S6125AJSid.dat' $ Normal land-sea mask file input line is kept but file is not used. $ 39 1 1 '(....)' 'NAME' 'S6125Masks.dat' $ Boundary cell id list file (unit 35) is only required if boundary $ cell number entered above is non-zero. The cell id number should be $ the sequential number in the cell array (unit 32) S625MCels.dat. $ $ If sub-grid information is available as indicated by FLAGTR above, $ additional input to define this is needed below. In such cases a $ field of fractional obstructions at or between grid points needs to $ be supplied. First the location and format of the data is defined $ by (as above) : $ - Unit number of file (can be 10, and/or identical to bottom depth $ unit), scale factor for fractional obstruction, IDLA, IDFM, $ format for formatted read, FROM and filename $ 10 0.2 3 1 '(....)' 'NAME' 'obstr.inp' $ $ *** NOTE if this unit number is the same as the previous bottom $ depth unit number, it is assumed that this is the same file $ without further checks. *** $ $ If the above unit number equals 10, the bottom data is read from $ this file and follows below (no intermediate comment lines allowed, $ except between the two fields). $ 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 5 0 0 0 0 0 0 0 0 0 0 0 5 0 0 0 0 0 0 0 0 0 0 0 4 0 0 0 0 0 0 0 0 0 0 0 4 0 0 0 0 0 0 0 0 0 0 0 5 0 0 0 0 0 0 0 0 0 0 0 5 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 $ 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 5 5 5 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 $ $ *** NOTE size of fields is always NX * NY *** $ $ Input boundary points and excluded points -------------------------- $ $ The first line identifies where to get the map data, by unit number $ IDLA and IDFM, format for formatted read, FROM and filename $ if FROM = 'PART', then segmented data is read from below, else $ the data is read from file as with the other inputs (as INTEGER) $ 10 3 1 '(....)' 'PART' 'mapsta.inp' $ $ Read the status map from file ( FROM != PART ) --------------------- $ $ $ 3 3 3 3 3 3 3 3 3 3 3 3 $ 3 2 1 1 1 1 0 1 1 1 1 3 $ 3 2 1 1 1 1 0 1 1 1 1 3 $ 3 2 1 1 1 1 0 1 1 1 1 3 $ 3 2 1 1 1 1 0 0 1 1 1 3 $ 3 2 1 1 1 1 1 1 1 1 1 3 $ 3 2 1 1 1 1 1 1 1 1 1 3 $ 3 2 1 1 1 1 1 1 1 1 1 3 $ 3 2 1 1 1 1 1 1 1 1 1 3 $ 3 2 1 1 1 1 1 1 1 1 1 3 $ 3 2 1 1 1 1 1 1 1 1 1 3 $ 3 3 3 3 3 3 3 3 3 3 3 3 $ $ The legend for the input map is : $ $ 0 : Land point. $ 1 : Regular sea point. $ 2 : Active boundary point. $ 3 : Point excluded from grid. $ $ Input boundary points from segment data ( FROM = PART ) ------------ $ $ An unlimited number of lines identifying points at which input $ boundary conditions are to be defined. If the actual input data is $ not defined in the actual wave model run, the initial conditions $ will be applied as constant boundary conditions. Each line contains: $ Discrete grid counters (IX,IY) of the active point and a $ connect flag. If this flag is true, and the present and previous $ point are on a grid line or diagonal, all intermediate points $ are also defined as boundary points. $ 2 2 F 2 11 T $ $ Close list by defining point (0,0) (mandatory) $ 0 0 F $ $ Excluded grid points from segment data ( FROM != PART ) $ First defined as lines, identical to the definition of the input $ boundary points, and closed the same way. $ 0 0 F $ $ Second, define a point in a closed body of sea points to remove $ the entire body of sea points. Also close by point (0,0) $ 0 0 $ $ Sedimentary bottom map if namelist &SBT4 SEDMAPD50 = T $ $ 22 1. 1 1 '(f10.6)' 'NAME' 'SED.txt' $ $ Output boundary points --------------------------------------------- $ $ Output boundary points are defined as a number of straight lines, $ defined by its starting point (X0,Y0), increments (DX,DY) and number $ of points. A negative number of points starts a new output file. $ Note that this data is only generated if requested by the actual $ program. Example again for spherical grid in degrees. Note, these do $ not need to be defined for data transfer between grids in the multi $ grid driver. $ 1.75 1.50 0.25 -0.10 3 2.25 1.50 -0.10 0.00 -6 0.10 0.10 0.10 0.00 -10 $ $ Close list by defining line with 0 points (mandatory) $ 0. 0. 0. 0. 0 $ $ -------------------------------------------------------------------- $ $ End of input file $ $ -------------------------------------------------------------------- $