forked from NOAA-EMC/WW3
-
Notifications
You must be signed in to change notification settings - Fork 0
Commit
This commit does not belong to any branch on this repository, and may belong to a fork outside of the repository.
- Loading branch information
1 parent
ff0358a
commit f2081d2
Showing
14 changed files
with
6,342 additions
and
0 deletions.
There are no files selected for viewing
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| Original file line number | Diff line number | Diff line change |
|---|---|---|
| @@ -0,0 +1,29 @@ | ||
| ############################################################################# | ||
| # # | ||
| # ww3_tp2.6 Test script for WW-III, unstructured grid. # | ||
| # Example of Limon harbour, Costa Rica # | ||
| # # | ||
| # Model should be compiled with the switches : # | ||
| # # | ||
| # PR3 LN0 ST0 NL1 BT1 DB1 MLIM TR0 BS0 WNX1 WNT1 CRX1 CRT1 # | ||
| # O0 O1 O2 O3 O4 O5 O6 O7 NOGRB SHRD FLX0 # | ||
| # # | ||
| # The choice of propagation schemes for unstructured grids is done via # | ||
| # namelists. # | ||
| # # | ||
| # Remarks : # | ||
| # # | ||
| # - No other optional switches should be used. # | ||
| # # | ||
| # Sample run_test commands : # | ||
| # ./bin/run_test -s ST0 ../model ww3_tp2.6 # | ||
| # ./bin/run_test -s ST4 ../model ww3_tp2.6 # | ||
| # Domain Decomposition Implicit needs PDLIB and either of METIS or # | ||
| # SCOTCH switch # | ||
| # ./bin/run_test -s PDLIB ../model ww3_tp2.6 # | ||
| # # | ||
| # Aron Roland and Fabrice Ardhuin, Jan 2011 # | ||
| # Ali Abdolali Oct 2018 # | ||
| # Last Mod : Oct 2018 # | ||
| # # | ||
| ############################################################################# |
Large diffs are not rendered by default.
Oops, something went wrong.
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| Original file line number | Diff line number | Diff line change |
|---|---|---|
| @@ -0,0 +1,26 @@ | ||
| &SLN1 CLIN = 80.0, RFPM = 1.00, RFHF = 0.50 / | ||
| &UNST UGOBCAUTO = F, | ||
| UGOBCDEPTH= -10., | ||
| EXPFSN = F, | ||
| EXPFSPSI = F, | ||
| EXPFSFCT = F, | ||
| IMPFSN = F, | ||
| EXPTOTAL = F, | ||
| IMPTOTAL = T, | ||
| IMPREFRACTION = T, | ||
| IMPFREQSHIFT = T, | ||
| IMPSOURCE = T, | ||
| SETUP_APPLY_WLV = F, | ||
| SOLVERTHR_SETUP=1E-14, | ||
| CRIT_DEP_SETUP=0.1, | ||
| JGS_USE_JACOBI = T, | ||
| JGS_BLOCK_GAUSS_SEIDEL = T, | ||
| JGS_TERMINATE_MAXITER = T, | ||
| JGS_MAXITER = 1000, | ||
| JGS_TERMINATE_NORM = F, | ||
| JGS_TERMINATE_DIFFERENCE = T, | ||
| JGS_DIFF_THR = 1.E-8, | ||
| JGS_PMIN = 3.0, | ||
| JGS_LIMITER = F, | ||
| JGS_NORM_THR = 1.E-20 / | ||
| END OF NAMELISTS |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| Original file line number | Diff line number | Diff line change |
|---|---|---|
| @@ -0,0 +1,2 @@ | ||
| -83.02097840 9.98081977 'POINT1' | ||
| -83.01589326 9.98398516 'POINT2' |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| Original file line number | Diff line number | Diff line change |
|---|---|---|
| @@ -0,0 +1 @@ | ||
| ASCII NOGRB DIST MPI PDLIB SCOTCH PR3 UQ FLX2 LN0 ST0 NL1 BT1 DB1 MLIM TR0 BS0 IC0 IS0 REF0 WNT1 WNX1 CRT1 CRX1 O0 O1 O2 O3 O4 O5 O6 O7 O10 O11 |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| Original file line number | Diff line number | Diff line change |
|---|---|---|
| @@ -0,0 +1,235 @@ | ||
| $ -------------------------------------------------------------------- $ | ||
| $ WAVEWATCH III Grid preprocessor input file $ | ||
| $ -------------------------------------------------------------------- $ | ||
| $ Grid name (C*30, in quotes) | ||
| $ | ||
| 'LIMON' | ||
| $ | ||
| $ 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.05 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 T 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). | ||
| $ | ||
| $ | ||
| 10. 10. 10. 10. | ||
| $ | ||
| $ 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 | ||
| $ | ||
| &SLN1 CLIN = 80.0, RFPM = 1.00, RFHF = 0.50 / | ||
| $ | ||
| &UNST UGOBCAUTO = F, | ||
| UGOBCDEPTH= -10., | ||
| EXPFSN = F, | ||
| EXPFSPSI = F, | ||
| EXPFSFCT = F, | ||
| IMPFSN = F, | ||
| EXPTOTAL = F, | ||
| IMPTOTAL = T, | ||
| IMPREFRACTION = T, | ||
| IMPFREQSHIFT = T, | ||
| IMPSOURCE = T, | ||
| SETUP_APPLY_WLV = F, | ||
| SOLVERTHR_SETUP=1E-14, | ||
| CRIT_DEP_SETUP=0.1, | ||
| JGS_USE_JACOBI = T, | ||
| JGS_BLOCK_GAUSS_SEIDEL = T, | ||
| JGS_TERMINATE_MAXITER = T, | ||
| JGS_MAXITER = 1000, | ||
| JGS_TERMINATE_NORM = F, | ||
| JGS_TERMINATE_DIFFERENCE = T, | ||
| JGS_DIFF_THR = 1.E-8, | ||
| JGS_PMIN = 3.0, | ||
| JGS_LIMITER = F, | ||
| JGS_NORM_THR = 1.E-20 / | ||
| $ | ||
| $ Bottom friction - - - - - - - - - - - - - - - - - - - - - - - - - - | ||
| $ JONSWAP : Namelist SBT1 | ||
| $ GAMMA : As it says. | ||
| $ &SBT1 GAMMA = 0.15 / | ||
| $ | ||
| $ Propagation schemes ------------------------------------------------ $ | ||
| $ First order : Namelist PRO1 | ||
| $ CFLTM : Maximum CFL number for refraction. | ||
| $ | ||
| $ UQ with diffusion : Namelist PRO2 | ||
| $ CFLTM : Maximum CFL number for refraction. | ||
| $ FLSOFT : Flag for 'soft' land boundaries. | ||
| $ 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 with averaging : Namelist PRO3 | ||
| $ CFLTM : Maximum CFL number for refraction. | ||
| $ FLSOFT : Flag for 'soft' land boundaries. | ||
| $ WDTHCG : Tuning factor propag. direction. | ||
| $ WDTHTH : Tuning factor normal direction. | ||
| $ | ||
| $ UQ with divergence : Namelist PRO4 | ||
| $ CFLTM : Maximum CFL number for refraction. | ||
| $ FLSOFT : Flag for 'soft' land boundaries. | ||
| $ QTFAC : Tuning factor Eq. (3.41). | ||
| $ RSFAC : Tuning factor Eq. (3.42). | ||
| $ RNFAC : Tuning factor Eq. (3.43). | ||
| $ | ||
| $ Miscellaneous ------------------------------------------------------ $ | ||
| $ Misc. parameters : Namelist MISC | ||
| $ CICE0 : Ice concentration cut-off. | ||
| $ CICEN : Ice concentration cut-off. | ||
| $ 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. | ||
| $ | ||
| $ 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 / | ||
| $ | ||
| $ Mandatory string to identify end of namelist input section. | ||
| $ | ||
| END OF NAMELISTS | ||
| $ | ||
| $ Define grid -------------------------------------------------------- $ | ||
| $ | ||
| 'UNST' T 'NONE' | ||
| $ | ||
| 5.0 0.30 20 -1. 4 1 '(20f10.2)' 'NAME' './../input/limon_ll.msh' | ||
| $ | ||
| $ If the above unit number equals 10, the bottom data is read from | ||
| $ this file and follows below (no intermediate comment lines allowed). | ||
| $ | ||
| $ 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 | ||
| $ | ||
| $ If sub-grid information is avalaible 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 bottem 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 *** | ||
| $ | ||
| 10 3 1 '(....)' 'PART' 'mapsta.inp' | ||
| $ Input boundary points ---------------------------------------------- $ | ||
| $ 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. | ||
| $ | ||
| 66 1 F | ||
| 99 1 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 os sea points. Also close by point (0,0) | ||
| $ | ||
| 0 0 | ||
| $ | ||
| $ 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. | ||
| $ | ||
| $ -2.5312 48.5 0.00 0.008738 102 | ||
| $ -2.5312 49.3850 0.013554 0.00 51 | ||
| $ | ||
| $ Close list by defining line with 0 points (mandatory) | ||
| $ | ||
| 0. 0. 0. 0. 0 | ||
| $ | ||
| $ -------------------------------------------------------------------- $ | ||
| $ End of input file $ | ||
| $ -------------------------------------------------------------------- $ |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| Original file line number | Diff line number | Diff line change |
|---|---|---|
| @@ -0,0 +1,43 @@ | ||
| $ -------------------------------------------------------------------- $ | ||
| $ WAVEWATCH III Grid output post-processing $ | ||
| $--------------------------------------------------------------------- $ | ||
| $ First output time (yyyymmdd hhmmss), increment of output (s), | ||
| $ and number of output times. | ||
| $ | ||
| 20100801 000000 10 3600 | ||
| $ | ||
| $ Fields requested --------------------------------------------------- $ | ||
| $ | ||
| $ Output request flags identifying fields as in ww3_shel.inp. See that | ||
| $ file for a full documentation of field output options. Namelist type | ||
| $ selection is used here (for alternative F/T flags, see ww3_shel.inp). | ||
| $ | ||
| N | ||
| HS LM T02 T01 T0M1 UST CHA CGE DTD FC CFX CFD QP QKK | ||
| $ | ||
| $--------------------------------------------------------------------- $ | ||
| $ netCDF version [3,4] | ||
| $ and variable type 4 [2 = SHORT, 3 = it depends , 4 = REAL] | ||
| $ swell partitions [0 1 2 3 4 5] | ||
| $ variables in same file [T] or not [F] | ||
| $ | ||
| 3 4 | ||
| 0 1 2 | ||
| F | ||
| $ | ||
| $ -------------------------------------------------------------------- $ | ||
| $ File prefix | ||
| $ number of characters in date | ||
| $ IX, IY range | ||
| $ | ||
| ww3. | ||
| 6 | ||
| 1 1778 1 1 | ||
| $ | ||
| $ For each field and time a new file is generated with the file name | ||
| $ ww3.date_xxx.nc , where date is a conventional time indicator with S3 | ||
| $ characters, and xxx is a field identifier. | ||
| $ | ||
| $ -------------------------------------------------------------------- $ | ||
| $ End of input file $ | ||
| $ -------------------------------------------------------------------- $ |
Oops, something went wrong.