diff --git a/model/bin/switch_NCEP_rwps b/model/bin/switch_NCEP_rwps new file mode 100644 index 0000000000..2ce91c2aa7 --- /dev/null +++ b/model/bin/switch_NCEP_rwps @@ -0,0 +1 @@ +NCO PDLIB SCOTCH NOGRB BIN2NC DIST MPI PR3 UQ FLX0 SEED ST4 STAB0 NL1 BT1 DB1 MLIM TR1 BS0 RWND WNX1 WNT1 CRX1 CRT1 O0 O1 O2 O3 O4 O5 O6 O7 O14 O15 IC0 IS0 REF0 diff --git a/regtests/ww3_ufs1.1/info b/regtests/ww3_ufs1.1/info index 66250d12d9..a0550ce0c7 100644 --- a/regtests/ww3_ufs1.1/info +++ b/regtests/ww3_ufs1.1/info @@ -18,6 +18,8 @@ # # grid_a: Domain Decomposition (PDLIB) and Explicit solver # # # grid_b: Domain Decomposition (PDLIB) and Block Explicit solver # # # grid_c: Domain Decomposition (PDLIB) and Implicit solver # +# # grid_d: Domain Decomposition (PDLIB) and Block Explicit solver # +# # grid_e: Domain Decomposition (PDLIB) and Implicit solver # # --------------------------------------------------------------------------# # if "ufscoarse" option is selected in bin/matrix.base, a coarser grid with # # 5 deg. resolution is used for MODEL and INPUT grids. # @@ -34,11 +36,11 @@ # ww3_ounf and GRIB2 using ww3_grib programs. The grib2 files are with # # rectilinear projection(-o all option). # # # -# Model should be compiled with the switches : # +# Model is compiled with the switches : # # # -# NCO NOGRB DIST MPI SCRIP SCRIPNC WRST PR3 UQ FLX0 SEED ST4 STAB0 # -# NL1 BT1 DB1 MLIM FLD2 TR0 BS0 RWND WNX1 WNT1 CRX1 CRT1 O0 O1 O2 O3 O4 # -# O5 O6 O7 O14 O15 IC0 IS0 REF0. # +# NCO PDLIB SCOTCH NOGRB BIN2NC DIST MPI PR3 UQ FLX0 SEED ST4 STAB0 NL1 BT1 # +# DB1 MLIM FLD2 TR0 BS0 RWND WNX1 WNT1 CRX1 CRT1 O0 O1 O2 O3 O4 O5 O6 O7 # +# O14 O15 IC0 IS0 REF0 # # # # Remarks : # # # @@ -80,6 +82,6 @@ # --------------------------------------------------------------------------# # # # Ali Abdolali, April 2021 # -# Last Mod : Feb 2023 # +# Last Mod : Oct 2025 # # # ############################################################################# diff --git a/regtests/ww3_ufs1.1/input_unstr/namelists_b.nml b/regtests/ww3_ufs1.1/input_unstr/namelists_b.nml index 8d1c72eda1..09fa434eac 100644 --- a/regtests/ww3_ufs1.1/input_unstr/namelists_b.nml +++ b/regtests/ww3_ufs1.1/input_unstr/namelists_b.nml @@ -1,73 +1,57 @@ -&SLN1 CLIN = 80.0, RFPM = 1.00, RFHF = 0.50 / -&SIN4 ALPHA0=0.0095, -BETAMAX=1.33, -SINTHP=2.00, -Z0MAX=0.00, -ZALP=0.006, -ZWND=10.00, -TAUWSHELTER =1.00, -SWELLFPAR = 1, -SWELLF= 0.800, -SWELLF2=-0.018, -SWELLF3 =0.015, -SWELLF4 =100000.0, -SWELLF5 =1.200, -SWELLF6 =0.000, -SWELLF7 =230000.000, -Z0RAT =0.0400 / -&UNST UGOBCAUTO = F, -UGOBCDEPTH= -10., -UGBCCFL= F, -EXPFSN = F, -EXPFSPSI = F, -EXPFSFCT = F, -IMPFSN = F, -EXPTOTAL = T, -IMPTOTAL = F, -IMPREFRACTION = F, -IMPFREQSHIFT = F, -IMPSOURCE = F, -SETUP_APPLY_WLV = F, -SOLVERTHR_SETUP=1E-14, -CRIT_DEP_SETUP=0.1, -JGS_NLEVEL = 0, -JGS_USE_JACOBI = F, -JGS_BLOCK_GAUSS_SEIDEL = F, -JGS_TERMINATE_MAXITER = F, -JGS_MAXITER = 1000, -JGS_TERMINATE_NORM = F, -JGS_TERMINATE_DIFFERENCE = F, -JGS_DIFF_THR = 1.E-8, -JGS_PMIN = 3.0, -JGS_LIMITER = F, -JGS_NORM_THR = 1.E-6 / -&SNL1 LAMBDA = 0.250, NLPROP = 0.250E+08, KDCONV = 0.750, KDMIN = 0.500, -SNLCS1 = 5.500, SNLCS2 = 0.833, SNLCS3 = -1.250 / -&SDS4 SDSBCHOICE = 1, SDSC2 = -0.2200E-04, SDSCUM = -0.4034E+00, -SDSC4 = 0.1000E+01, SDSC5 = 0.0000E+00, SDSC6 = 0.3000E+00, -WNMEANP =0.50, FXPM3 =4.00,FXFM3 =9.90, -SDSBINT = 0.3000E+00, SDSBCK = 0.0000E+00, SDSABK = 1.500, SDSPBK = 4.000, -SDSHCK = 1.50, SDSBR = 0.9000E-03, SDSSTRAIN = 0.000, -SDSP = 2.00, SDSISO = 2, SDSCOS =2.0, SDSDTH = 80.0, -SDSBRF1 = 0.50, SDSBRFDF = 0, -SDSBM0 = 1.00, SDSBM1 = 0.00, SDSBM2 = 0.00, SDSBM3 = 0.00, SDSBM4 = 0.00, -WHITECAPWIDTH = 0.30/ -&SBT1 GAMMA = -0.6700E-01 / -&SDB1 BJALFA = 1.000, BJGAM = 0.730, BJFLAG = .TRUE. / -&PRO3 CFLTM = 0.70, WDTHCG = 1.50, WDTHTH = 1.50 / -&OUTS P2SF = 0, I1P2SF = 1, I2P2SF = 15, -US3D = 0, I1US3D = 1, I2US3D = 32, -E3D = 0, I1E3D = 1, I2E3D = 32, -TH1MF = 0, I1TH1M = 1, I2TH1M = 32, -STH1MF= 0, I1STH1M= 1, I2STH1M= 32, -TH2MF = 0, I1TH2M = 1, I2TH2M = 32, -STH2MF= 0, I1STH2M= 1, I2STH2M= 32, -E3D = 1, USSP = 1, IUSSP = 3, STK_WN = 0.04, 0.110, 0.3305 / -&REF1 REFCOAST=0.10, REFSLOPE=0.1, REFCOSP_STRAIGHT=4, REFFREQ=1., REFSUBGRID = 0.00 / -&SIN4 BETAMAX = 1.33 / -&MISC CICE0 = 0.250, CICEN = 0.750, LICE = 0.0, PMOVE = 0.500, -XSEED = 1.000, FLAGTR = 0, XP = 0.150, XR = 0.100, XFILT = 0.050 -IHM = 100, HSPM = 0.050, WSM = 1.700, WSC = 0.333, FLC = .TRUE. -NOSW = 3, FMICHE = 1.600, RWNDC = 1.000, -FACBERG = 1.0, GSHIFT = 0.000E+00 / +&OUTS + USSP = 1, + IUSSP = 3, + STK_WN = 0.04, 0.110, 0.3305 / +&SIN4 + BETAMAX = 1.315, + TAUWSHELTER = 1.0, + SWELLF = 0.798, + SWELLF2 = -0.0127, + SWELLF3 = 0.0151, + SWELLF4 = 100025.0, + SWELLF5 = 1.1999, + SWELLF7 = 235500.0 / +&SNL1 + NLPROP = 2.502E7 / +&SDS4 + FXFM3 = 2.501, + SDSC2 = -2.1975e-05, + SDSCUM = -0.4032, + SDSC6 = 0.2978, + SDSBR = 0.0009035 / +&MISC + CICE0 = 0.75, + CICEN = 0.75 / +&SBT1 + GAMMA = -0.038 / +&PRO3 + WDTHCG=1.5, + WDTHTH=1.5 / +&UNST + UGOBCAUTO = F, + UGOBCDEPTH = -10., + UGBCCFL = F, + EXPFSN = F, + EXPFSPSI = F, + EXPFSFCT = F, + IMPFSN = F, + EXPTOTAL = T, + IMPTOTAL = F, + IMPREFRACTION = F, + IMPFREQSHIFT = F, + IMPSOURCE = F, + SETUP_APPLY_WLV = F, + SOLVERTHR_SETUP = 1E-14, + CRIT_DEP_SETUP = 0.1, + JGS_NLEVEL = 0, + JGS_USE_JACOBI = F, + JGS_BLOCK_GAUSS_SEIDEL = F, + JGS_TERMINATE_MAXITER = F, + JGS_MAXITER = 1000, + JGS_TERMINATE_NORM = F, + JGS_TERMINATE_DIFFERENCE = F, + JGS_DIFF_THR = 1.E-8, + JGS_PMIN = 3.0, + JGS_LIMITER = F, + JGS_NORM_THR = 1.E-6 / END OF NAMELISTS diff --git a/regtests/ww3_ufs1.1/input_unstr/namelists_c.nml b/regtests/ww3_ufs1.1/input_unstr/namelists_c.nml index 013ea40af0..d4a540c6af 100644 --- a/regtests/ww3_ufs1.1/input_unstr/namelists_c.nml +++ b/regtests/ww3_ufs1.1/input_unstr/namelists_c.nml @@ -1,72 +1,56 @@ -&SLN1 CLIN = 80.0, RFPM = 1.00, RFHF = 0.50 / -&SIN4 ALPHA0=0.0095, -BETAMAX=1.33, -SINTHP=2.00, -Z0MAX=0.00, -ZALP=0.006, -ZWND=10.00, -TAUWSHELTER =1.00, -SWELLFPAR = 1, -SWELLF= 0.800, -SWELLF2=-0.018, -SWELLF3 =0.015, -SWELLF4 =100000.0, -SWELLF5 =1.200, -SWELLF6 =0.000, -SWELLF7 =230000.000, -Z0RAT =0.0400 / +&OUTS + USSP = 1, + IUSSP = 3, + STK_WN = 0.04, 0.110, 0.3305 / +&SIN4 + BETAMAX = 1.315, + TAUWSHELTER = 1.0, + SWELLF = 0.798, + SWELLF2 = -0.0127, + SWELLF3 = 0.0151, + SWELLF4 = 100025.0, + SWELLF5 = 1.1999, + SWELLF7 = 235500.0 / +&SNL1 + NLPROP = 2.502E7 / +&SDS4 + FXFM3 = 2.501, + SDSC2 = -2.1975e-05, + SDSCUM = -0.4032, + SDSC6 = 0.2978, + SDSBR = 0.0009035 / +&MISC + CICE0 = 0.75, + CICEN = 0.75 / +&SBT1 + GAMMA = -0.038 / +&PRO3 + WDTHCG=1.5, + WDTHTH=1.5 / &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-6, -JGS_PMIN = 3.0, -JGS_LIMITER = F, -JGS_NORM_THR = 1.E-6 / -&SNL1 LAMBDA = 0.250, NLPROP = 0.250E+08, KDCONV = 0.750, KDMIN = 0.500, -SNLCS1 = 5.500, SNLCS2 = 0.833, SNLCS3 = -1.250 / -&SDS4 SDSBCHOICE = 1, SDSC2 = -0.2200E-04, SDSCUM = -0.4034E+00, -SDSC4 = 0.1000E+01, SDSC5 = 0.0000E+00, SDSC6 = 0.3000E+00, -WNMEANP =0.50, FXPM3 =4.00,FXFM3 =9.90, -SDSBINT = 0.3000E+00, SDSBCK = 0.0000E+00, SDSABK = 1.500, SDSPBK = 4.000, -SDSHCK = 1.50, SDSBR = 0.9000E-03, SDSSTRAIN = 0.000, -SDSP = 2.00, SDSISO = 2, SDSCOS =2.0, SDSDTH = 80.0, -SDSBRF1 = 0.50, SDSBRFDF = 0, -SDSBM0 = 1.00, SDSBM1 = 0.00, SDSBM2 = 0.00, SDSBM3 = 0.00, SDSBM4 = 0.00, -, WHITECAPWIDTH = 0.30/ -&SBT1 GAMMA = -0.6700E-01 / -&SDB1 BJALFA = 1.000, BJGAM = 0.730, BJFLAG = .TRUE. / -&PRO3 CFLTM = 0.70, WDTHCG = 1.50, WDTHTH = 1.50 / -&OUTS P2SF = 0, I1P2SF = 1, I2P2SF = 15, -US3D = 0, I1US3D = 1, I2US3D = 32, -E3D = 0, I1E3D = 1, I2E3D = 32, -TH1MF = 0, I1TH1M = 1, I2TH1M = 32, -STH1MF= 0, I1STH1M= 1, I2STH1M= 32, -TH2MF = 0, I1TH2M = 1, I2TH2M = 32, -STH2MF= 0, I1STH2M= 1, I2STH2M= 32, -E3D = 1, USSP = 1, IUSSP = 3, STK_WN = 0.04, 0.110, 0.3305 / -&REF1 REFCOAST=0.10, REFSLOPE=0.1, REFCOSP_STRAIGHT=4, REFFREQ=1., REFSUBGRID = 0.00 / -&SIN4 BETAMAX = 1.33 / -&MISC CICE0 = 0.250, CICEN = 0.750, LICE = 0.0, PMOVE = 0.500, -XSEED = 1.000, FLAGTR = 0, XP = 0.150, XR = 0.100, XFILT = 0.050 -IHM = 100, HSPM = 0.050, WSM = 1.700, WSC = 0.333, FLC = .TRUE. -NOSW = 3, FMICHE = 1.600, RWNDC = 1.000, -FACBERG = 1.0, GSHIFT = 0.000E+00 / + 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_NLEVEL = 0, + JGS_USE_JACOBI = T, + JGS_BLOCK_GAUSS_SEIDEL = F, + JGS_TERMINATE_MAXITER = T, + JGS_MAXITER = 1000, + JGS_TERMINATE_NORM = F, + JGS_TERMINATE_DIFFERENCE = T, + JGS_DIFF_THR = 1.E-6, + JGS_PMIN = 3.0, + JGS_LIMITER = F, + JGS_NORM_THR = 1.E-6 / END OF NAMELISTS diff --git a/regtests/ww3_ufs1.1/input_unstr/switch_PDLIB b/regtests/ww3_ufs1.1/input_unstr/switch_PDLIB index 90d19f1c3a..1a4dd99a6b 100644 --- a/regtests/ww3_ufs1.1/input_unstr/switch_PDLIB +++ b/regtests/ww3_ufs1.1/input_unstr/switch_PDLIB @@ -1 +1 @@ -NCO PDLIB SCOTCH NOGRB DIST MPI PR3 UQ FLX0 SEED ST4 STAB0 NL1 BT1 DB1 MLIM FLD1 TR0 BS0 WNX1 WNT1 CRX1 CRT1 O0 O1 O2 O3 O4 O5 O6 O7 O14 O15 IC0 IS0 REF0 BIN2NC +NCO PDLIB SCOTCH NOGRB BIN2NC DIST MPI PR3 UQ FLX0 SEED ST4 STAB0 NL1 BT1 DB1 MLIM FLD2 TR0 BS0 RWND WNX1 WNT1 CRX1 CRT1 O0 O1 O2 O3 O4 O5 O6 O7 O14 O15 IC0 IS0 REF0 diff --git a/regtests/ww3_ufs1.1/input_unstr/ww3_grid_a.inp b/regtests/ww3_ufs1.1/input_unstr/ww3_grid_a.inp index d2f508ac2f..6ae42b9081 100644 --- a/regtests/ww3_ufs1.1/input_unstr/ww3_grid_a.inp +++ b/regtests/ww3_ufs1.1/input_unstr/ww3_grid_a.inp @@ -66,60 +66,6 @@ $ Implicit with ww3ifr code version EXPFSFCT = F, IMPFSN = F / $ -$ 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. -$ -$ &SNL1 LAMBDA = 0.250, NLPROP = 0.250E+08, KDCONV = 0.750, KDMIN = 0.500, SNLCS1 = 5.500, SNLCS2 = 0.833, SNLCS3 = -1.250 / &SDS4 SDSBCHOICE = 1, SDSC2 = -0.2200E-04, SDSCUM = -0.4034E+00, @@ -130,7 +76,7 @@ $ SDSP = 2.00, SDSISO = 2, SDSCOS =2.0, SDSDTH = 80.0, SDSBRF1 = 0.50, SDSBRFDF = 0, SDSBM0 = 1.00, SDSBM1 = 0.00, SDSBM2 = 0.00, SDSBM3 = 0.00, SDSBM4 = 0.00, -, WHITECAPWIDTH = 0.30/ + WHITECAPWIDTH = 0.30/ &SBT1 GAMMA = -0.6700E-01 / &SDB1 BJALFA = 1.000, BJGAM = 0.730, BJFLAG = .TRUE. / &PRO3 CFLTM = 0.70, WDTHCG = 1.50, WDTHTH = 1.50 / @@ -143,27 +89,12 @@ $ STH2MF= 0, I1STH2M= 1, I2STH2M= 32, E3D = 1, USSP = 1, IUSSP = 3, STK_WN = 0.04, 0.110, 0.3305 / $ -$ -$AW021317 &MISC P2SF = 1 ,I1P2SF = 2, I2P2SF = 16 / -$&REF1 REFCOAST=0.1 / -$&REF1 REFCOAST=0.10, REFSLOPE=0.1, REFCOSP_STRAIGHT=4, REFFREQ=1., REFSUBGRID = 0.00 / &REF1 REFCOAST=0.10, REFSLOPE=0.1, REFCOSP_STRAIGHT=4, REFFREQ=1., REFSUBGRID = 0.00 / -&SIN4 BETAMAX = 1.33 / &MISC CICE0 = 0.250, CICEN = 0.750, LICE = 0.0, PMOVE = 0.500, XSEED = 1.000, FLAGTR = 0, XP = 0.150, XR = 0.100, XFILT = 0.050 IHM = 100, HSPM = 0.050, WSM = 1.700, WSC = 0.333, FLC = .TRUE. NOSW = 3, FMICHE = 1.600, RWNDC = 1.000, FACBERG = 1.0, GSHIFT = 0.000E+00 / -$&MISC CICE0 = 0.250, CICEN = 0.750, LICE = 0.0, PMOVE = 0.500, -$ XSEED = 1.000, FLAGTR = 4, XP = 0.150, XR = 0.100, XFILT = 0.050 -$ IHM = 100, HSPM = 0.050, WSM = 1.700, WSC = 0.333, FLC = .TRUE. -$ NOSW = 3, FMICHE = 1.600, RWNDC = 1.000, WCOR1 = 99.00, WCOR2 = 0.00, -$ FACBERG = 1.0, GSHIFT = 0.000E+00, STDX = -1.00, STDY = -1.00, -$ STDT = -1.00, ICEHMIN = 0.20, ICEHFAC = 1.00, -$ ICEHINIT = 0.50, ICEDISP = F, ICEHDISP = 0.60, -$ ICESLN = 1.00, ICEWIND = 1.00, ICESNL = 1.00, ICESDS = 1.00, -$ ICEDDISP = 80.00, ICEFDISP = 2.00, CALTYPE = standard , TRCKCMPR = T, -$ BTBET = 1.20 / $ $ Mandatory string to identify end of namelist input section. $ diff --git a/regtests/ww3_ufs1.1/input_unstr/ww3_grid_b.inp b/regtests/ww3_ufs1.1/input_unstr/ww3_grid_b.inp index 3474c9a869..74b3cb0921 100644 --- a/regtests/ww3_ufs1.1/input_unstr/ww3_grid_b.inp +++ b/regtests/ww3_ufs1.1/input_unstr/ww3_grid_b.inp @@ -1,9 +1,8 @@ $ -------------------------------------------------------------------- $ $ WAVEWATCH III Grid preprocessor input file $ $ -------------------------------------------------------------------- $ -$ Grid name (C*30, in quotes) $ - 'GLOBAL' + 'Global Unstructured' $ $ Frequency increment factor and first frequency (Hz) ---------------- $ $ number of frequencies (wavenumbers) and directions, relative offset @@ -12,7 +11,9 @@ $ 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.07 0.035 33 36 0. +$gfs: 1.07 0.035 50 36 0.5 +$gefs: 1.1 0.035 33 36 0.5 + 1.1 0.035 33 36 0.5 $ $ Set model flags ---------------------------------------------------- $ $ - FLDRY Dry run (input/output only, no calculation). @@ -28,40 +29,43 @@ $ maximum global time step, maximum CFL time step for x-y and $ k-theta, minimum source term time step (all in seconds). $ $ - 720. 360. 360. 30. + 720. 180. 360. 30. $ $ 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. $ -&SLN1 CLIN = 80.0, RFPM = 1.00, RFHF = 0.50 / -$ -&SIN4 ALPHA0=0.0095, - BETAMAX=1.33, - SINTHP=2.00, - Z0MAX=0.00, - ZALP=0.006, - ZWND=10.00, - TAUWSHELTER =1.00, - SWELLFPAR = 1, - SWELLF= 0.800, - SWELLF2=-0.018, - SWELLF3 =0.015, - SWELLF4 =100000.0, - SWELLF5 =1.200, - SWELLF6 =0.000, - SWELLF7 =230000.000, - Z0RAT =0.0400 / -$ -$ Implicit with ww3ifr code version -&UNST UGOBCAUTO = F, - UGOBCDEPTH= -10., - UGBCCFL= F, +&OUTS + USSP = 1, + IUSSP = 3, + STK_WN = 0.04, 0.110, 0.3305 / +&SIN4 + BETAMAX = 1.315, + TAUWSHELTER = 1.0, + SWELLF = 0.798, + SWELLF2 = -0.0127, + SWELLF3 = 0.0151, + SWELLF4 = 100025.0, + SWELLF5 = 1.1999, + SWELLF7 = 235500.0 / +&SNL1 + NLPROP = 2.502E7 / +&SDS4 + FXFM3 = 2.501, + SDSC2 = -2.1975e-05, + SDSCUM = -0.4032, + SDSC6 = 0.2978, + SDSBR = 0.0009035 / +&MISC + CICE0 = 0.75, + CICEN = 0.75 / +&SBT1 + GAMMA = -0.038 / +&PRO3 + WDTHCG=1.5, + WDTHTH=1.5 / +&UNST + UGOBCAUTO = F, + UGOBCDEPTH = -10., + UGBCCFL = F, EXPFSN = F, EXPFSPSI = F, EXPFSFCT = F, @@ -72,8 +76,8 @@ $ Implicit with ww3ifr code version IMPFREQSHIFT = F, IMPSOURCE = F, SETUP_APPLY_WLV = F, - SOLVERTHR_SETUP=1E-14, - CRIT_DEP_SETUP=0.1, + SOLVERTHR_SETUP = 1E-14, + CRIT_DEP_SETUP = 0.1, JGS_NLEVEL = 0, JGS_USE_JACOBI = F, JGS_BLOCK_GAUSS_SEIDEL = F, @@ -86,105 +90,6 @@ $ Implicit with ww3ifr code version JGS_LIMITER = F, JGS_NORM_THR = 1.E-6 / $ -$ 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. -$ -$ -&SNL1 LAMBDA = 0.250, NLPROP = 0.250E+08, KDCONV = 0.750, KDMIN = 0.500, - SNLCS1 = 5.500, SNLCS2 = 0.833, SNLCS3 = -1.250 / -&SDS4 SDSBCHOICE = 1, SDSC2 = -0.2200E-04, SDSCUM = -0.4034E+00, - SDSC4 = 0.1000E+01, SDSC5 = 0.0000E+00, SDSC6 = 0.3000E+00, - WNMEANP =0.50, FXPM3 =4.00,FXFM3 =9.90, - SDSBINT = 0.3000E+00, SDSBCK = 0.0000E+00, SDSABK = 1.500, SDSPBK = 4.000, - SDSHCK = 1.50, SDSBR = 0.9000E-03, SDSSTRAIN = 0.000, - SDSP = 2.00, SDSISO = 2, SDSCOS =2.0, SDSDTH = 80.0, - SDSBRF1 = 0.50, SDSBRFDF = 0, - SDSBM0 = 1.00, SDSBM1 = 0.00, SDSBM2 = 0.00, SDSBM3 = 0.00, SDSBM4 = 0.00, -, WHITECAPWIDTH = 0.30/ -&SBT1 GAMMA = -0.6700E-01 / -&SDB1 BJALFA = 1.000, BJGAM = 0.730, BJFLAG = .TRUE. / -&PRO3 CFLTM = 0.70, WDTHCG = 1.50, WDTHTH = 1.50 / -&OUTS P2SF = 0, I1P2SF = 1, I2P2SF = 15, - US3D = 0, I1US3D = 1, I2US3D = 32, - E3D = 0, I1E3D = 1, I2E3D = 32, - TH1MF = 0, I1TH1M = 1, I2TH1M = 32, - STH1MF= 0, I1STH1M= 1, I2STH1M= 32, - TH2MF = 0, I1TH2M = 1, I2TH2M = 32, - STH2MF= 0, I1STH2M= 1, I2STH2M= 32, - E3D = 1, USSP = 1, IUSSP = 3, STK_WN = 0.04, 0.110, 0.3305 / -$ -$ -$AW021317 &MISC P2SF = 1 ,I1P2SF = 2, I2P2SF = 16 / -$&REF1 REFCOAST=0.1 / -$&REF1 REFCOAST=0.10, REFSLOPE=0.1, REFCOSP_STRAIGHT=4, REFFREQ=1., REFSUBGRID = 0.00 / -&REF1 REFCOAST=0.10, REFSLOPE=0.1, REFCOSP_STRAIGHT=4, REFFREQ=1., REFSUBGRID = 0.00 / -&SIN4 BETAMAX = 1.33 / -&MISC CICE0 = 0.250, CICEN = 0.750, LICE = 0.0, PMOVE = 0.500, - XSEED = 1.000, FLAGTR = 0, XP = 0.150, XR = 0.100, XFILT = 0.050 - IHM = 100, HSPM = 0.050, WSM = 1.700, WSC = 0.333, FLC = .TRUE. - NOSW = 3, FMICHE = 1.600, RWNDC = 1.000, - FACBERG = 1.0, GSHIFT = 0.000E+00 / -$&MISC CICE0 = 0.250, CICEN = 0.750, LICE = 0.0, PMOVE = 0.500, -$ XSEED = 1.000, FLAGTR = 4, XP = 0.150, XR = 0.100, XFILT = 0.050 -$ IHM = 100, HSPM = 0.050, WSM = 1.700, WSC = 0.333, FLC = .TRUE. -$ NOSW = 3, FMICHE = 1.600, RWNDC = 1.000, WCOR1 = 99.00, WCOR2 = 0.00, -$ FACBERG = 1.0, GSHIFT = 0.000E+00, STDX = -1.00, STDY = -1.00, -$ STDT = -1.00, ICEHMIN = 0.20, ICEHFAC = 1.00, -$ ICEHINIT = 0.50, ICEDISP = F, ICEHDISP = 0.60, -$ ICESLN = 1.00, ICEWIND = 1.00, ICESNL = 1.00, ICESDS = 1.00, -$ ICEDDISP = 80.00, ICEFDISP = 2.00, CALTYPE = standard , TRCKCMPR = T, -$ BTBET = 1.20 / -$ $ Mandatory string to identify end of namelist input section. $ END OF NAMELISTS @@ -226,77 +131,13 @@ $ $ 4.0 0.30 20 -1. 4 1 '(20f10.2)' 'NAME' '../input_unstr/global_1deg_unstr.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. -$ -$ are also defined as boundary points. -$ $ 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. @@ -309,14 +150,6 @@ $ 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) $ diff --git a/regtests/ww3_ufs1.1/input_unstr/ww3_grid_b.nml b/regtests/ww3_ufs1.1/input_unstr/ww3_grid_b.nml index e31a4eabbc..7f8ec4bc30 100644 --- a/regtests/ww3_ufs1.1/input_unstr/ww3_grid_b.nml +++ b/regtests/ww3_ufs1.1/input_unstr/ww3_grid_b.nml @@ -10,8 +10,10 @@ SPECTRUM%FREQ1 = 0.035 SPECTRUM%NK = 33 SPECTRUM%NTH = 36 + SPECTRUM%THOFF = 0.5 / + ! -------------------------------------------------------------------- ! ! Define the run parameterization via RUN_NML namelist ! -------------------------------------------------------------------- ! @@ -28,7 +30,7 @@ ! -------------------------------------------------------------------- ! &TIMESTEPS_NML TIMESTEPS%DTMAX = 720. - TIMESTEPS%DTXY = 360. + TIMESTEPS%DTXY = 180. TIMESTEPS%DTKTH = 360. TIMESTEPS%DTMIN = 30. / diff --git a/regtests/ww3_ufs1.1/input_unstr/ww3_grid_c.inp b/regtests/ww3_ufs1.1/input_unstr/ww3_grid_c.inp index 32db4238db..caaa4e78ac 100644 --- a/regtests/ww3_ufs1.1/input_unstr/ww3_grid_c.inp +++ b/regtests/ww3_ufs1.1/input_unstr/ww3_grid_c.inp @@ -1,9 +1,8 @@ $ -------------------------------------------------------------------- $ $ WAVEWATCH III Grid preprocessor input file $ $ -------------------------------------------------------------------- $ -$ Grid name (C*30, in quotes) $ - 'GLOBAL' + 'Global Unstructured' $ $ Frequency increment factor and first frequency (Hz) ---------------- $ $ number of frequencies (wavenumbers) and directions, relative offset @@ -12,7 +11,9 @@ $ 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.07 0.035 50 36 0. +$gfs: 1.07 0.035 50 36 0.5 +$gefs: 1.1 0.035 33 36 0.5 + 1.1 0.035 33 36 0.5 $ $ Set model flags ---------------------------------------------------- $ $ - FLDRY Dry run (input/output only, no calculation). @@ -27,162 +28,66 @@ $ - 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. 900. 900. 900. +$ + 900. 900. 900. 900. $ $ 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. -$ -&SLN1 CLIN = 80.0, RFPM = 1.00, RFHF = 0.50 / -$ -&SIN4 ALPHA0=0.0095, - BETAMAX=1.33, - SINTHP=2.00, - Z0MAX=0.00, - ZALP=0.006, - ZWND=10.00, - TAUWSHELTER =1.00, - SWELLFPAR = 1, - SWELLF= 0.800, - SWELLF2=-0.018, - SWELLF3 =0.015, - SWELLF4 =100000.0, - SWELLF5 =1.200, - SWELLF6 =0.000, - SWELLF7 =230000.000, - Z0RAT =0.0400 / -$ -$ -$ Implicit with ww3ifr code version -&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-6, -JGS_PMIN = 3.0, -JGS_LIMITER = F, -JGS_NORM_THR = 1.E-6 / -$ 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. -$ -$ -&SNL1 LAMBDA = 0.250, NLPROP = 0.250E+08, KDCONV = 0.750, KDMIN = 0.500, - SNLCS1 = 5.500, SNLCS2 = 0.833, SNLCS3 = -1.250 / -&SDS4 SDSBCHOICE = 1, SDSC2 = -0.2200E-04, SDSCUM = -0.4034E+00, - SDSC4 = 0.1000E+01, SDSC5 = 0.0000E+00, SDSC6 = 0.3000E+00, - WNMEANP =0.50, FXPM3 =4.00,FXFM3 =9.90, - SDSBINT = 0.3000E+00, SDSBCK = 0.0000E+00, SDSABK = 1.500, SDSPBK = 4.000, - SDSHCK = 1.50, SDSBR = 0.9000E-03, SDSSTRAIN = 0.000, - SDSP = 2.00, SDSISO = 2, SDSCOS =2.0, SDSDTH = 80.0, - SDSBRF1 = 0.50, SDSBRFDF = 0, - SDSBM0 = 1.00, SDSBM1 = 0.00, SDSBM2 = 0.00, SDSBM3 = 0.00, SDSBM4 = 0.00, -, WHITECAPWIDTH = 0.30/ -&SBT1 GAMMA = -0.6700E-01 / -&SDB1 BJALFA = 1.000, BJGAM = 0.730, BJFLAG = .TRUE. / -&PRO3 CFLTM = 0.70, WDTHCG = 1.50, WDTHTH = 1.50 / -&OUTS P2SF = 0, I1P2SF = 1, I2P2SF = 15, - US3D = 0, I1US3D = 1, I2US3D = 32, - E3D = 0, I1E3D = 1, I2E3D = 32, - TH1MF = 0, I1TH1M = 1, I2TH1M = 32, - STH1MF= 0, I1STH1M= 1, I2STH1M= 32, - TH2MF = 0, I1TH2M = 1, I2TH2M = 32, - STH2MF= 0, I1STH2M= 1, I2STH2M= 32, - E3D = 1, USSP = 1, IUSSP = 3, STK_WN = 0.04, 0.110, 0.3305 / -$ $ -$AW021317 &MISC P2SF = 1 ,I1P2SF = 2, I2P2SF = 16 / -$&REF1 REFCOAST=0.1 / -$&REF1 REFCOAST=0.10, REFSLOPE=0.1, REFCOSP_STRAIGHT=4, REFFREQ=1., REFSUBGRID = 0.00 / -&REF1 REFCOAST=0.10, REFSLOPE=0.1, REFCOSP_STRAIGHT=4, REFFREQ=1., REFSUBGRID = 0.00 / -&SIN4 BETAMAX = 1.33 / -&MISC CICE0 = 0.250, CICEN = 0.750, LICE = 0.0, PMOVE = 0.500, - XSEED = 1.000, FLAGTR = 0, XP = 0.150, XR = 0.100, XFILT = 0.050 - IHM = 100, HSPM = 0.050, WSM = 1.700, WSC = 0.333, FLC = .TRUE. - NOSW = 3, FMICHE = 1.600, RWNDC = 1.000, - FACBERG = 1.0, GSHIFT = 0.000E+00 / -$&MISC CICE0 = 0.250, CICEN = 0.750, LICE = 0.0, PMOVE = 0.500, -$ XSEED = 1.000, FLAGTR = 4, XP = 0.150, XR = 0.100, XFILT = 0.050 -$ IHM = 100, HSPM = 0.050, WSM = 1.700, WSC = 0.333, FLC = .TRUE. -$ NOSW = 3, FMICHE = 1.600, RWNDC = 1.000, WCOR1 = 99.00, WCOR2 = 0.00, -$ FACBERG = 1.0, GSHIFT = 0.000E+00, STDX = -1.00, STDY = -1.00, -$ STDT = -1.00, ICEHMIN = 0.20, ICEHFAC = 1.00, -$ ICEHINIT = 0.50, ICEDISP = F, ICEHDISP = 0.60, -$ ICESLN = 1.00, ICEWIND = 1.00, ICESNL = 1.00, ICESDS = 1.00, -$ ICEDDISP = 80.00, ICEFDISP = 2.00, CALTYPE = standard , TRCKCMPR = T, -$ BTBET = 1.20 / +&OUTS + USSP = 1, + IUSSP = 3, + STK_WN = 0.04, 0.110, 0.3305 / +&SIN4 + BETAMAX = 1.315, + TAUWSHELTER = 1.0, + SWELLF = 0.798, + SWELLF2 = -0.0127, + SWELLF3 = 0.0151, + SWELLF4 = 100025.0, + SWELLF5 = 1.1999, + SWELLF7 = 235500.0 / +&SNL1 + NLPROP = 2.502E7 / +&SDS4 + FXFM3 = 2.501, + SDSC2 = -2.1975e-05, + SDSCUM = -0.4032, + SDSC6 = 0.2978, + SDSBR = 0.0009035 / +&MISC + CICE0 = 0.75, + CICEN = 0.75 / +&SBT1 + GAMMA = -0.038 / +&PRO3 + WDTHCG=1.5, + WDTHTH=1.5 / +&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_NLEVEL = 0, + JGS_USE_JACOBI = T, + JGS_BLOCK_GAUSS_SEIDEL = F, + JGS_TERMINATE_MAXITER = T, + JGS_MAXITER = 1000, + JGS_TERMINATE_NORM = F, + JGS_TERMINATE_DIFFERENCE = T, + JGS_DIFF_THR = 1.E-6, + JGS_PMIN = 3.0, + JGS_LIMITER = F, + JGS_NORM_THR = 1.E-6 / $ $ Mandatory string to identify end of namelist input section. $ @@ -225,77 +130,13 @@ $ $ 4.0 0.30 20 -1. 4 1 '(20f10.2)' 'NAME' '../input_unstr/global_1deg_unstr.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. -$ -$ are also defined as boundary points. -$ $ 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. @@ -308,14 +149,6 @@ $ 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) $ diff --git a/regtests/ww3_ufs1.1/input_unstr/ww3_grid_c.nml b/regtests/ww3_ufs1.1/input_unstr/ww3_grid_c.nml index 9c68a4e317..00dcdf635b 100644 --- a/regtests/ww3_ufs1.1/input_unstr/ww3_grid_c.nml +++ b/regtests/ww3_ufs1.1/input_unstr/ww3_grid_c.nml @@ -6,10 +6,11 @@ ! Define the spectrum parameterization via SPECTRUM_NML namelist ! -------------------------------------------------------------------- ! &SPECTRUM_NML - SPECTRUM%XFR = 1.07 + SPECTRUM%XFR = 1.1 SPECTRUM%FREQ1 = 0.035 - SPECTRUM%NK = 50 + SPECTRUM%NK = 33 SPECTRUM%NTH = 36 + SPECTRUM%THOFF = 0.5 / ! -------------------------------------------------------------------- ! diff --git a/regtests/ww3_ufs1.1/input_unstr/ww3_ounf.inp b/regtests/ww3_ufs1.1/input_unstr/ww3_ounf.inp index c721aeffce..ccc3d4722b 100644 --- a/regtests/ww3_ufs1.1/input_unstr/ww3_ounf.inp +++ b/regtests/ww3_ufs1.1/input_unstr/ww3_ounf.inp @@ -19,7 +19,7 @@ $ WCF WCH WCM SXY TWO BHD FOC TUS USS P2S USF P2L TWI FIC ABR UBR BED $ FBB TBB MSS MSC DTD FC CFX CFD CFK U1 U2 $ N - WND CUR ICE HS T01 T02 DIR FP DP PHS PTP PDIR UST CHA + WND CUR ICE HS T01 T02 DIR FP DP PHS PTP PDIR $ $--------------------------------------------------------------------- $ $ netCDF version [3,4] diff --git a/regtests/ww3_ufs1.1/input_unstr/ww3_ounf.nml b/regtests/ww3_ufs1.1/input_unstr/ww3_ounf.nml index 2995d58833..ec1c8f0f86 100644 --- a/regtests/ww3_ufs1.1/input_unstr/ww3_ounf.nml +++ b/regtests/ww3_ufs1.1/input_unstr/ww3_ounf.nml @@ -28,8 +28,6 @@ FILE%VAR(10) = 'PHS' FILE%VAR(11) = 'PTP' FILE%VAR(12) = 'PDIR' - FILE%VAR(13) = 'UST' - FILE%VAR(14) = 'CHA' / ! -------------------------------------------------------------------- ! diff --git a/regtests/ww3_ufs1.1/input_unstr/ww3_shel.inp b/regtests/ww3_ufs1.1/input_unstr/ww3_shel.inp index c33d0a4d1c..4030e0593e 100644 --- a/regtests/ww3_ufs1.1/input_unstr/ww3_shel.inp +++ b/regtests/ww3_ufs1.1/input_unstr/ww3_shel.inp @@ -27,7 +27,7 @@ $ Fields of mean wave parameters $ 20210401 000000 3600 20210402 000000 N -CHA EF UST WND HS FP DP PHS PTP PDIR CUR ICE + WND CUR ICE HS T01 T02 DIR FP DP PHS PTP PDIR $ $ Point output $ diff --git a/regtests/ww3_ufs1.1/input_unstr/ww3_shel.nml b/regtests/ww3_ufs1.1/input_unstr/ww3_shel.nml index c3ceb1056e..0fa97ba304 100644 --- a/regtests/ww3_ufs1.1/input_unstr/ww3_shel.nml +++ b/regtests/ww3_ufs1.1/input_unstr/ww3_shel.nml @@ -25,7 +25,7 @@ ! Define the output types point parameters via OUTPUT_TYPE_NML namelist ! -------------------------------------------------------------------- ! &OUTPUT_TYPE_NML - TYPE%FIELD%LIST = 'CHA EF UST WND HS FP DP PHS PTP PDIR CUR ICE' + TYPE%FIELD%LIST = 'WND CUR ICE HS T01 T02 DIR FP DP PHS PTP PDIR' TYPE%POINT%FILE = '../input_unstr/ww3_points.list' /