$ -------------------------------------------------------------------- $
$ 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                                                    $
$ -------------------------------------------------------------------- $