ROMS+SWAN+SEDIMENTS (One-way nesting)

Discussion of how to use ROMS on different regional and basin scale applications.

Moderators: arango, robertson

Post Reply
Message
Author
rgon
Posts: 12
Joined: Thu Jan 29, 2009 4:00 pm
Location: University of Plymouth

ROMS+SWAN+SEDIMENTS (One-way nesting)

#1 Unread post by rgon »

Hello,
I've been studying waves and tidal currents with the coupled system ROMS+SWAN. Three nested levels (from L1 to L3) were run in SWAN, and in the finest level (L3) the ROMS+SWAN system was coupled. Now I am in a new test case on a similar way, the difference is I am coupling both models at L2, then, a one-way nested (interpolation) to L3 with both models, and again coupling them at L3; finally the sediment transport model has been incorporated. BTW: I am waiting desperately for the new ROMS-NESTED version.

As you can see in my output.log , the Input open boundary data file (BRYNAME) is read one record before ROMS starts the coupled running, I mean, It starts with a one record-lag.
In order to correct the one record-lag;
Q1: About the boundary conditions, do they look reasonable? because the resulting sea surface is lower than the input (nested) one.
Q2: With regards to my cpp.h file/options, do you see anything unusual for my test case? In terms of Boundary conditions, Vertical mixing, Sediment transport.
I appreciate your comments.
rgs.

CPP OPTIONS:

Code: Select all

** Options for L3_WH Case, waves-ocean (SWAN/ROMS) two-way coupling.
**
** Application flag:   WH_COAST
** Input script:       ocean_whcoast.in
**                     coupling_whcoast.in
**                     sediment_whcoast.in
*/

#define MASKING
#define UV_ADV
#define UV_COR
#define TS_MPDATA
#define WET_DRY
#define SOLVE3D
#define SPLINES

/* ANALYTICAL */
#define ANA_INITIAL
#define ANA_SMFLUX
#define ANA_STFLUX
#define ANA_SSFLUX
#define ANA_BPFLUX
#define ANA_BTFLUX
#define ANA_BSFLUX
#define ANA_SPFLUX
#define ANA_SRFLUX

/* BOUNDARIES */
#define NORTH_FSCHAPMAN
#define NORTH_M2FLATHER
#define NORTH_TRADIATION
#define NORTH_M3RADIATION

#define WEST_FSCHAPMAN
#define WEST_M2FLATHER
#define WEST_TRADIATION
#define WEST_M3RADIATION

#define SOUTH_FSCHAPMAN
#define SOUTH_M2FLATHER
#define SOUTH_TRADIATION
#define SOUTH_M3RADIATION

#define EAST_FSCHAPMAN
#define EAST_M2FLATHER
#define EAST_TRADIATION
#define EAST_M3RADIATION

/* TIDES FROM ONE-WAY OFF NESTING */
#define ADD_FSOBC
#define ADD_M2OBC

/* COUPLING  */ 
#define SWAN_COUPLING
#define NEARSHORE_MELLOR05
#ifdef SWAN_COUPLING
# define MCT_LIB
#endif

/* BOTTOM FRICTION */
#define SSW_BBL
#if defined SSW_BBL
#define SSW_CALC_ZNOT
#endif

#define SEDIMENT
#ifdef SEDIMENT
# define SUSPLOAD
# define BEDLOAD_SOULSBY
# undef  SED_MORPH
#endif

#if defined SEDIMENT || defined SSW_BBL
#define ANA_SEDIMENT
#endif

/* VERTICAL MIXING */
#define MY25_MIXING
#undef GLS_MIXING
#if defined GLS_MIXING || defined MY25_MIXING
# define KANTHA_CLAYSON
# define N2S2_HORAVG
# undef CRAIG_BANNER
# undef CHARNOK
#endif
OUTPUT.LOG FILE

Code: Select all

 Resolution, Grid 01: 0298x0144x010,  Parallel Nodes:   4,  Tiling: 002x002


 Physical Parameters, Grid: 01
 =============================

     178560  ntimes          Number of timesteps for 3-D equations.
     30.000  dt              Timestep size (s) for 3-D equations.
         30  ndtfast         Number of timesteps for 2-D equations between
                               each 3D timestep.
          1  ERstr           Starting ensemble/perturbation run number.
          1  ERend           Ending ensemble/perturbation run number.
          0  nrrec           Number of restart records to read from disk.
          T  LcycleRST       Switch to recycle time-records in restart file.
        120  nRST            Number of timesteps between the writing of data
                               into restart fields.
          1  ninfo           Number of timesteps between print of information
                               to standard output.
          T  ldefout         Switch to create a new output NetCDF file(s).
        120  nHIS            Number of timesteps between the writing fields
                               into history file.
       2880  ndefHIS         Number of timesteps between creation of new
                               history files.
 5.0000E-06  Akt_bak(01)     Background vertical mixing coefficient (m2/s)
                               for tracer 01: temp
 5.0000E-06  Akt_bak(02)     Background vertical mixing coefficient (m2/s)
                               for tracer 02: salt
 5.0000E-05  Akv_bak         Background vertical mixing coefficient (m2/s)
                               for momentum.
 5.0000E-06  Akk_bak         Background vertical mixing coefficient (m2/s)
                               for turbulent energy.
 5.0000E-06  Akp_bak         Background vertical mixing coefficient (m2/s)
                               for turbulent generic statistical field.
 3.0000E-04  rdrg            Linear bottom drag coefficient (m/s).
 2.5000E-02  rdrg2           Quadratic bottom drag coefficient.
 1.5000E-02  Zob             Bottom roughness (m).
 1.0000E-01  Dcrit           Minimum depth for wetting and drying (m).
          2  Vtransform      S-coordinate transformation equation.
          2  Vstretching     S-coordinate stretching function.
 1.0000E+00  theta_s         S-coordinate surface control parameter.
 1.0000E+00  theta_b         S-coordinate bottom  control parameter.
     20.000  Tcline          S-coordinate surface/bottom layer width (m) used
                               in vertical coordinate stretching.
   1025.000  rho0            Mean density (kg/m3) for Boussinesq approximation.
      0.000  dstart          Time-stamp assigned to model initialization (days).
20101020.00  time_ref        Reference time for units attribute (yyyymmdd.dd)
 0.0000E+00  Tnudg(01)       Nudging/relaxation time scale (days)
                               for tracer 01: temp
 0.0000E+00  Tnudg(02)       Nudging/relaxation time scale (days)
                               for tracer 02: salt
 0.0000E+00  Znudg           Nudging/relaxation time scale (days)
                               for free-surface.
 0.0000E+00  M2nudg          Nudging/relaxation time scale (days)
                               for 2D momentum.
 0.0000E+00  M3nudg          Nudging/relaxation time scale (days)
                               for 3D momentum.
 0.0000E+00  obcfac          Factor between passive and active
                               open boundary conditions.
     10.000  T0              Background potential temperature (C) constant.
     30.000  S0              Background salinity (PSU) constant.
   1025.000  R0              Background density (kg/m3) used in linear Equation
                               of State.
 1.7000E-04  Tcoef           Thermal expansion coefficient (1/Celsius).
 7.6000E-04  Scoef           Saline contraction coefficient (1/PSU).
      1.000  gamma2          Slipperiness variable: free-slip (1.0) or 
                                                    no-slip (-1.0).
          T  Hout(idFsur)    Write out free-surface.
          T  Hout(idUbar)    Write out 2D U-momentum component.
          T  Hout(idVbar)    Write out 2D V-momentum component.
          T  Hout(idUvel)    Write out 3D U-momentum component.
          T  Hout(idVvel)    Write out 3D V-momentum component.
          T  Hout(idWvel)    Write out W-momentum component.
          T  Hout(idOvel)    Write out omega vertical velocity.
          T  Hout(idTvar)    Write out tracer 01: temp
          T  Hout(idTvar)    Write out tracer 02: salt
          T  Hout(idUsms)    Write out surface U-momentum stress.
          T  Hout(idVsms)    Write out surface V-momentum stress.
          T  Hout(idUbms)    Write out bottom U-momentum stress.
          T  Hout(idVbms)    Write out bottom V-momentum stress.
          T  Hout(idUbrs)    Write out bottom U-current stress.
          T  Hout(idVbrs)    Write out bottom V-current stress.
          T  Hout(idUbws)    Write out wind-induced, bottom U-wave stress.
          T  Hout(idVbws)    Write out wind-induced, bottom V-wave stress.
          T  Hout(idUbcs)    Write out max wind + current, bottom U-wave stress.
          T  Hout(idVbcs)    Write out max wind + current, bottom V-wave stress.
          T  Hout(idUbot)    Write out bed wave orbital U-velocity.
          T  Hout(idVbot)    Write out bed wave orbital V-velocity.
          T  Hout(idUbur)    Write out bottom U-momentum above bed.
          T  Hout(idVbvr)    Write out bottom V-momentum above bed.
          T  Hout(idW3xx)    Write out 3D radiation stress, Sxx.
          T  Hout(idW3xy)    Write out 3D radiation stress, Sxy.
          T  Hout(idW3yy)    Write out 3D radiation stress, Syy.
          T  Hout(idW3zx)    Write out 3D radiation stress, Szx.
          T  Hout(idW3zy)    Write out 3D radiation stress, Szy.
          T  Hout(idU3rs)    Write out total 3D u-radiation stress.
          T  Hout(idV3rs)    Write out total 3D v-radiation stress.
          T  Hout(idU3Sd)    Write out 3D u-momentum stokes velocity.
          T  Hout(idV3Sd)    Write out 3D v-momentum stokes velocity.
          T  Hout(idWamp)    Write out wave height.
          T  Hout(idWlen)    Write out wavelength.
          T  Hout(idWdir)    Write out wave direction.
          T  Hout(idBott)    Write out bottom property 01: grain_diameter
          T  Hout(idBott)    Write out bottom property 02: grain_density
          T  Hout(idBott)    Write out bottom property 03: settling_vel
          T  Hout(idBott)    Write out bottom property 04: erosion_stress
          T  Hout(idBott)    Write out bottom property 05: ripple_length
          T  Hout(idBott)    Write out bottom property 06: ripple_height
          T  Hout(idBott)    Write out bottom property 07: bed_wave_amp
          T  Hout(idBott)    Write out bottom property 08: Zo_def
          T  Hout(idBott)    Write out bottom property 09: Zo_app
          T  Hout(idMtke)    Write out turbulent kinetic energy.

 Output/Input Files:

             Output Restart File:  wh_rst.nc
        Prefix for History Files:  wh_his
        Physical parameters File:  ocean_wh.in
                 Input Grid File:  seagridL3_146x300.nc
             Input Boundary File:  WHsma_his_frc.nc

 Tile partition information for Grid 01:  0298x0144x0010  tiling: 002x002

     tile     Istr     Iend     Jstr     Jend     Npts

        0        1      149        1       72   107280
        1      150      298        1       72   107280
        2        1      149       73      144   107280
        3      150      298       73      144   107280

 Tile minimum and maximum fractional grid coordinates:
   (interior points only)

     tile     Xmin     Xmax     Ymin     Ymax     grid

        0     0.50   149.50     0.50    72.50  RHO-points
        1   149.50   298.50     0.50    72.50  RHO-points
        2     0.50   149.50    72.50   144.50  RHO-points
        3   149.50   298.50    72.50   144.50  RHO-points

        0     1.00   149.50     0.50    72.50    U-points
        1   149.50   298.00     0.50    72.50    U-points
        2     1.00   149.50    72.50   144.50    U-points
        3   149.50   298.00    72.50   144.50    U-points

        0     0.50   149.50     1.00    72.50    V-points
        1   149.50   298.50     1.00    72.50    V-points
        2     0.50   149.50    72.50   144.00    V-points
        3   149.50   298.50    72.50   144.00    V-points

 Maximum halo size in XI and ETA directions:

               HaloSizeI(1) =     486
               HaloSizeJ(1) =     255
                TileSide(1) =     156
                TileSize(1) =   12324



 Sediment Parameters, Grid: 01
 =============================


 Size     Sd50        Csed        Srho        Wsed        Erate       poros
 Class    (mm)       (kg/m3)     (kg/m3)     (mm/s)     (kg/m2/s)    (nondim)

   1    1.0000E+00  0.0000E+00  2.6500E+03  1.0000E+00  1.0000E-05  5.0000E-01

         tau_ce      tau_cd      nl_tnu2     nl_tnu4     Akt_bak      Tnudg
         (N/m2)      (N/m2)      (m2/s)      (m4/s)       (m2/s)      (day)

   1    7.0000E-02  2.0000E-02  0.0000E+00  0.0000E+00  5.0000E-06  0.0000E+00

         morph_fac
         (nondim)

   1    1.0000E+00

 New bed layer formed when deposition exceeds 0.10000E-01 (m).
 Two first layers are combined when 2nd layer smaller than 0.00000E+00 (m).
 Rate coefficient for bed load transport = 0.50000E-01

          T  Hout(idTvar)   Write out sediment01: sand_01
          T  Hout(idfrac)   Write out bed fraction, sediment 01: sandfrac_01
          T  Hout(idfrac)   Write out mass, sediment 01: sandmass_01
          T  Hout(idUbld)   Write out U-bedload, sediment 01: bedload_Usand_01
          T  Hout(idVbld)   Write out V-bedload, sediment 01: bedload_Vsand_01
          T  Hout(idSbed)   Write out BED property 01: bed_thickness
          T  Hout(idSbed)   Write out BED property 02: bed_age
          T  Hout(idSbed)   Write out BED property 03: bed_porosity

 Activated C-preprocessing Options:

 WH_COAST            Nested Level3 Test Case: Waves+Currents+Sediments
 ANA_BPFLUX          Analytical bottom passive tracers fluxes.
 ANA_BSFLUX          Analytical kinematic bottom salinity flux.
 ANA_BTFLUX          Analytical kinematic bottom temperature flux.
 ANA_INITIAL         Analytical initial conditions.
 ANA_SEDIMENT        Analytical sediment initial conditions.
 ANA_SMFLUX          Analytical kinematic surface momentum flux.
 ANA_SPFLUX          Analytical surface passive tracer fluxes.
 ANA_SSFLUX          Analytical kinematic surface salinity flux.
 ANA_STFLUX          Analytical kinematic surface temperature flux.
 ASSUMED_SHAPE       Using assumed-shape arrays.
 BEDLOAD_SOULSBY     Activate bed load sediment transport Soulsby formula.
 DOUBLE_PRECISION    Double precision arithmetic.
 EAST_FSCHAPMAN      Eastern edge, free-surface, Chapman condition.
 EAST_M2FLATHER      Eastern edge, 2D momentum, Flather condition.
 EAST_M3RADIATION    Eastern edge, 3D momentum, radiation condition.
 EAST_TRADIATION     Eastern edge, tracers, radiation condition.
 KANTHA_CLAYSON      Kantha and Clayson stability function formulation.
 MASKING             Land/Sea masking.
 MCT_LIB             Using Model Coupling Toolkit library.
 MPI                 MPI distributed-memory configuration.
 MY25_MIXING         Mellor/Yamada Level-2.5 mixing closure.
 NEARSHORE_MELLOR05  Nearshore Radiation Stress Terms (Mellor 2005).
 NONLINEAR           Nonlinear Model.
 !NONLIN_EOS         Linear Equation of State for seawater.
 NORTH_FSCHAPMAN     Northern edge, free-surface, Chapman condition.
 NORTH_M2FLATHER     Northern edge, 2D momentum, Flather condition.
 NORTH_M3RADIATION   Northern edge, 3D momentum, radiation condition.
 NORTH_TRADIATION    Northern edge, tracers, radiation condition.
 N2S2_HORAVG         Horizontal smoothing of buoyancy and shear.
 POWER_LAW           Power-law shape time-averaging barotropic filter.
 PRSGRD31            Standard density Jacobian formulation (Song, 1998).
 PROFILE             Time profiling activated .
 K_GSCHEME           Third-order upstream advection of TKE fields.
 RHO_SURF            Include difference between rho0 and surface density.
 !RST_SINGLE         Double precision fields in restart NetCDF file.
 SEDIMENT            Cohesive and noncohesive sediments.
 SUSPLOAD            Activate suspended sediment transport.
 SOLVE3D             Solving 3D Primitive Equations.
 SOUTH_FSCHAPMAN     Southern edge, free-surface, Chapman condition.
 SOUTH_M2FLATHER     Southern edge, 2D momentum, Flather condition.
 SOUTH_M3RADIATION   Southern edge, 3D momentum, radiation condition.
 SOUTH_TRADIATION    Southern edge, tracers, radiation condition.
 SPLINES             Conservative parabolic spline reconstruction.
 SSW_BBL             Styles and Glenn Bottom Boundary Layer - modified.
 SSW_CALC_ZNOT       Internal computation of bottom roughness.
 SWAN_COUPLING       Two-way SWAN/ROMS coupling.
 THREE_GHOST         Using three Ghost Points in halo regions.
 TS_MPDATA           Recursive flux corrected MPDATA 3D advection of tracers.
 UV_ADV              Advection of momentum.
 UV_COR              Coriolis term.
 UV_U3HADVECTION     Third-order upstream horizontal advection of 3D momentum.
 UV_C4VADVECTION     Fourth-order centered vertical advection of momentum.
 VAR_RHO_2D          Variable density barotropic mode.
 WAVES_OCEAN         Two-way wave-ocean models coupling.
 WEST_FSCHAPMAN      Western edge, free-surface, Chapman condition.
 WEST_M2FLATHER      Western edge, 2D momentum, Flather condition.
 WEST_M3RADIATION    Western edge, 3D momentum, radiation condition.
 WEST_TRADIATION     Western edge, tracers, radiation condition.
 WET_DRY             Wetting and drying activated.

 INITIAL: Configuring and initializing forward nonlinear model ...


 Vertical S-coordinate System: 

 level   S-coord     Cs-curve          at_hmin  over_slope     at_hmax

    10   0.0000000   0.0000000           0.000       0.000       0.000
     9  -0.1000000  -0.0103874          -0.410      -1.793      -2.458
     8  -0.2000000  -0.0452437          -0.845      -4.452      -6.818
     7  -0.3000000  -0.1078761          -1.308      -8.093     -13.338
     6  -0.4000000  -0.1987188          -1.799     -12.732     -22.053
     5  -0.5000000  -0.3154004          -2.315     -18.285     -32.776
     4  -0.6000000  -0.4527378          -2.853     -24.569     -45.106
     3  -0.7000000  -0.6026562          -3.403     -31.298     -58.414
     2  -0.8000000  -0.7540355          -3.954     -38.079     -71.835
     1  -0.9000000  -0.8924786          -4.492     -44.402     -84.251
     0  -1.0000000  -1.0000000          -5.000     -49.631     -94.262

 Time Splitting Weights: ndtfast =  30    nfast =  48

    Primary            Secondary            Accumulated to Current Step

  1-0.0004593313153502 0.0333333333333333-0.0004593313153502 0.0333333333333333
  2-0.0006450812691645 0.0333486443771783-0.0011044125845148 0.0666819777105117
  3-0.0005737703233078 0.0333701470861505-0.0016781829078225 0.1000521247966621
  4-0.0002619189396448 0.0333892727635941-0.0019401018474674 0.1334413975602562
  5 0.0002739524199594 0.0333980033949156-0.0016661494275080 0.1668394009551718
  6 0.0010173232936399 0.0333888716475836-0.0006488261338681 0.2002282726027554
  7 0.0019516732195320 0.0333549608711289 0.0013028470856639 0.2335832334738843
  8 0.0030604817357706 0.0332899050971445 0.0043633288214344 0.2668731385710288
  9 0.0043272283804909 0.0331878890392855 0.0086905572019253 0.3000610276103143
 10 0.0057353926918280 0.0330436480932691 0.0144259498937533 0.3331046757035835
 11 0.0072684542079170 0.0328524683368749 0.0216944041016703 0.3659571440404584
 12 0.0089098924668930 0.0326101865299443 0.0306042965685633 0.3985673305704027
 13 0.0106431870068911 0.0323131901143812 0.0412474835754543 0.4308805206847839
 14 0.0124518173660464 0.0319584172141515 0.0536993009415007 0.4628389378989354
 15 0.0143192630824940 0.0315433566352833 0.0680185640239948 0.4943822945342187
 16 0.0162290036943691 0.0310660478658668 0.0842475677183638 0.5254483424000855
 17 0.0181645187398067 0.0305250810760545 0.1024120864581705 0.5559734234761401
 18 0.0201092877569419 0.0299195971180610 0.1225213742151124 0.5858930205942011
 19 0.0220467902839098 0.0292492875261629 0.1445681644990222 0.6151423081203641
 20 0.0239605058588456 0.0285143945166993 0.1685286703578678 0.6436567026370633
 21 0.0258339140198844 0.0277157109880711 0.1943625843777522 0.6713724136251343
 22 0.0276504943051612 0.0268545805207416 0.2220130786829133 0.6982269941458760
 23 0.0293937262528111 0.0259328973772362 0.2514068049357244 0.7241598915231122
 24 0.0310470894009693 0.0249531065021425 0.2824538943366938 0.7491129980252548
 25 0.0325940632877709 0.0239182035221102 0.3150479576244647 0.7730312015473649
 26 0.0340181274513510 0.0228317347458512 0.3490660850758157 0.7958629362932161
 27 0.0353027614298446 0.0216977971641395 0.3843688465056603 0.8175607334573556
 28 0.0364314447613869 0.0205210384498113 0.4208002912670472 0.8380817719071669
 29 0.0373876569841130 0.0193066569577651 0.4581879482511602 0.8573884288649321
 30 0.0381548776361580 0.0180604017249613 0.4963428258873183 0.8754488305898933
 31 0.0387165862556570 0.0167885724704227 0.5350594121429753 0.8922374030603161
 32 0.0390562623807451 0.0154980195952342 0.5741156745237204 0.9077354226555502
 33 0.0391573855495574 0.0141961441825426 0.6132730600732778 0.9219315668380929
 34 0.0390034353002291 0.0128908979975574 0.6522764953735068 0.9348224648356502
 35 0.0385778911708951 0.0115907834875498 0.6908543865444019 0.9464132483232000
 36 0.0378642326996907 0.0103048537818533 0.7287186192440925 0.9567181021050533
 37 0.0368459394247509 0.0090427126918636 0.7655645586688434 0.9657608147969168
 38 0.0355064908842108 0.0078145147110385 0.8010710495530542 0.9735753295079553
 39 0.0338293666162056 0.0066309650148982 0.8349004161692597 0.9802062945228536
 40 0.0317980461588703 0.0055033194610247 0.8666984623281300 0.9857096139838782
 41 0.0293960090503400 0.0044433845890623 0.8960944713784700 0.9901529985729406
 42 0.0266067348287499 0.0034635176207177 0.9227012062072198 0.9936165161936582
 43 0.0234137030322352 0.0025766264597593 0.9461149092394550 0.9961931426534175
 44 0.0198003931989306 0.0017961696920182 0.9659153024383856 0.9979893123454356
 45 0.0157502848669716 0.0011361565853871 0.9816655873053572 0.9991254689308228
 46 0.0112468575744933 0.0006111470898214 0.9929124448798505 0.9997366160206442
 47 0.0062735908596305 0.0002362518373383 0.9991860357394811 0.9999728678579826
 48 0.0008139642605187 0.0000271321420173 0.9999999999999998 0.9999999999999999

 ndtfast, nfast =   30  48   nfast/ndtfast = 1.60000

 Centers of gravity and integrals (values must be 1, 1, approx 1/2, 1, 1):

    1.000000000000 1.089661200636 0.544830600318 1.000000000000 1.000000000000

 Power filter parameters, Fgamma, gamma =  0.11810   0.07873

 Minimum X-grid spacing, DXmin =  3.50876449E-01 km
 Maximum X-grid spacing, DXmax =  3.56312991E-01 km
 Minimum Y-grid spacing, DYmin =  5.51709086E-01 km
 Maximum Y-grid spacing, DYmax =  5.60247247E-01 km
 Minimum Z-grid spacing, DZmin =  4.10387450E-01 m
 Maximum Z-grid spacing, DZmax =  1.34215486E+01 m

 Minimum barotropic Courant Number =  2.36019280E-02
 Maximum barotropic Courant Number =  1.02707004E-01
 Maximum Coriolis   Courant Number =  3.39286263E-03

                                                           [color=#FF0000]Here the nested input file is read:[/color]

GET_NGFLD   - free-surface western boundary condition,   t =     0 00:00:00
                   (Rec=0001, Index=2, File: WHsma_his_frc.nc)
                   (Tmin=          0.0000 Tmax=         62.0000)
                   (Min =  0.00000000E+00 Max =  0.00000000E+00)
    GET_NGFLD   - free-surface eastern boundary condition,   t =     0 00:00:00
                   (Rec=0001, Index=2, File: WHsma_his_frc.nc)
                   (Tmin=          0.0000 Tmax=         62.0000)
                   (Min =  0.00000000E+00 Max =  0.00000000E+00)
    GET_NGFLD   - free-surface southern boundary condition,  t =     0 00:00:00
                   (Rec=0001, Index=2, File: WHsma_his_frc.nc)
                   (Tmin=          0.0000 Tmax=         62.0000)
                   (Min =  0.00000000E+00 Max =  0.00000000E+00)
    GET_NGFLD   - free-surface northern boundary condition,  t =     0 00:00:00
                   (Rec=0001, Index=2, File: WHsma_his_frc.nc)
                   (Tmin=          0.0000 Tmax=         62.0000)
                   (Min =  0.00000000E+00 Max =  0.00000000E+00)
    GET_NGFLD   - 2D u-momentum western boundary condition,  t =     0 00:00:00
                   (Rec=0001, Index=2, File: WHsma_his_frc.nc)
                   (Tmin=          0.0000 Tmax=         62.0000)
                   (Min =  9.99997199E-01 Max =  1.00000250E+00)
    GET_NGFLD   - 2D v-momentum western boundary condition,  t =     0 00:00:00
                   (Rec=0001, Index=2, File: WHsma_his_frc.nc)
                   (Tmin=          0.0000 Tmax=         62.0000)
                   (Min =  9.99997556E-01 Max =  1.00000286E+00)
    GET_NGFLD   - 2D u-momentum eastern boundary condition,  t =     0 00:00:00
                   (Rec=0001, Index=2, File: WHsma_his_frc.nc)
                   (Tmin=          0.0000 Tmax=         62.0000)
                   (Min =  0.00000000E+00 Max =  9.99997914E-01)
    GET_NGFLD   - 2D v-momentum eastern boundary condition,  t =     0 00:00:00
                   (Rec=0001, Index=2, File: WHsma_his_frc.nc)
                   (Tmin=          0.0000 Tmax=         62.0000)
                   (Min =  0.00000000E+00 Max =  1.00000238E+00)
    GET_NGFLD   - 2D u-momentum southern boundary condition, t =     0 00:00:00
                   (Rec=0001, Index=2, File: WHsma_his_frc.nc)
                   (Tmin=          0.0000 Tmax=         62.0000)
                   (Min =  0.00000000E+00 Max =  1.00000226E+00)
    GET_NGFLD   - 2D v-momentum southern boundary condition, t =     0 00:00:00
                   (Rec=0001, Index=2, File: WHsma_his_frc.nc)
                   (Tmin=          0.0000 Tmax=         62.0000)
                   (Min =  0.00000000E+00 Max =  9.99997735E-01)
    GET_NGFLD   - 2D u-momentum northern boundary condition, t =     0 00:00:00
                   (Rec=0001, Index=2, File: WHsma_his_frc.nc)
                   (Tmin=          0.0000 Tmax=         62.0000)
                   (Min =  9.99997377E-01 Max =  9.99997497E-01)
    GET_NGFLD   - 2D v-momentum northern boundary condition, t =     0 00:00:00
                   (Rec=0001, Index=2, File: WHsma_his_frc.nc)
                   (Tmin=          0.0000 Tmax=         62.0000)
                   (Min =  1.00000238E+00 Max =  1.00000262E+00)

 Maximum grid stiffness ratios:  rx0 =   3.781332E-01 (Beckmann and Haidvogel)
                                 rx1 =   7.008454E+00 (Haney)


 Initial basin volumes: TotVolume =  3.5413929774E+11 m3
                        MinVolume =  7.9800308972E+04 m3
                        MaxVolume =  2.5979619182E+06 m3
                          Max/Min =  3.2555787712E+01

      OCN2WAV   - (06) imported and (05) exported fields,    t =     0 00:00:00
                - ROMS coupling exchanges wait clock (s):
                   (Recv= 1.08249543E+02 Send= 2.00000000E-03)
                - ROMS Import: wave direction
                   (Min=  0.00000000E+00 Max=  4.75629812E+00)
                - ROMS Import: significant wave height
                   (Min=  0.00000000E+00 Max=  1.55755445E-01)
                - ROMS Import: average wave length
                   (Min=  0.00000000E+00 Max=  1.10770714E+00)
                - ROMS Import: surface wave relative peak period
                   (Min=  0.00000000E+00 Max=  1.21969092E+00)
                - ROMS Import: bottom wave period
                   (Min=  0.00000000E+00 Max=  0.00000000E+00)
                - ROMS Import: wave bottom orbital velocity
                   (Min=  0.00000000E+00 Max=  0.00000000E+00)
                - ROMS Export: bathymetry
                   (Min=  5.00000000E+00 Max=  9.42617159E+01)
                - ROMS Export: free-surface
                   (Min=  0.00000000E+00 Max=  0.00000000E+00)
                - ROMS Export: vertically integrated u-momentum component
                   (Min=  0.00000000E+00 Max=  0.00000000E+00)
                - ROMS Export: vertically integrated v-momentum component
                   (Min=  0.00000000E+00 Max=  0.00000000E+00)
                - ROMS Export: bottom roughness
                   (Min=  1.00000000E-04 Max=  1.00000000E-04)


NL ROMS/TOMS: started time-stepping: (Grid: 01 TimeSteps: 00000001 - 00178560)
      WAV2OCN   - (05) imported and (06) exported fields,    t = 20101020.000000
                - SWAN coupling exchanges wait clock (s):
                   (Recv= 2.18967000E-01 Send= 5.99900000E-03)
    GET_NGFLD   - free-surface western boundary condition,   t =     0 01:00:00
                   (Rec=0002, Index=1, File: WHsma_his_frc.nc)
                   (Tmin=          0.0000 Tmax=         62.0000)
                   (Min =  2.36198846E-02 Max =  4.32650894E-02)
    GET_NGFLD   - free-surface eastern boundary condition,   t =     0 01:00:00
                   (Rec=0002, Index=1, File: WHsma_his_frc.nc)
                   (Tmin=          0.0000 Tmax=         62.0000)
                   (Min = -1.47665814E-02 Max =  0.00000000E+00)
    GET_NGFLD   - free-surface southern boundary condition,  t =     0 01:00:00
                   (Rec=0002, Index=1, File: WHsma_his_frc.nc)
                   (Tmin=          0.0000 Tmax=         62.0000)
                   (Min =  0.00000000E+00 Max =  4.32650894E-02)
    GET_NGFLD   - free-surface northern boundary condition,  t =     0 01:00:00
                   (Rec=0002, Index=1, File: WHsma_his_frc.nc)
                   (Tmin=          0.0000 Tmax=         62.0000)
                   (Min = -1.47665814E-02 Max =  2.36198846E-02)
    GET_NGFLD   - 2D u-momentum western boundary condition,  t =     0 01:00:00
                   (Rec=0002, Index=1, File: WHsma_his_frc.nc)
                   (Tmin=          0.0000 Tmax=         62.0000)
                   (Min = -1.02874478E-02 Max =  2.26140977E-03)
    GET_NGFLD   - 2D v-momentum western boundary condition,  t =     0 01:00:00
                   (Rec=0002, Index=1, File: WHsma_his_frc.nc)
                   (Tmin=          0.0000 Tmax=         62.0000)
                   (Min = -1.00057228E-02 Max =  6.36333227E-03)
    GET_NGFLD   - 2D u-momentum eastern boundary condition,  t =     0 01:00:00
                   (Rec=0002, Index=1, File: WHsma_his_frc.nc)
                   (Tmin=          0.0000 Tmax=         62.0000)
                   (Min =  0.00000000E+00 Max =  6.83444180E-03)
    GET_NGFLD   - 2D v-momentum eastern boundary condition,  t =     0 01:00:00
                   (Rec=0002, Index=1, File: WHsma_his_frc.nc)
                   (Tmin=          0.0000 Tmax=         62.0000)
                   (Min =  0.00000000E+00 Max =  4.28247824E-03)
    GET_NGFLD   - 2D u-momentum southern boundary condition, t =     0 01:00:00
                   (Rec=0002, Index=1, File: WHsma_his_frc.nc)
                   (Tmin=          0.0000 Tmax=         62.0000)
                   (Min = -1.17228618E-02 Max =  0.00000000E+00)
    GET_NGFLD   - 2D v-momentum southern boundary condition, t =     0 01:00:00
                   (Rec=0002, Index=1, File: WHsma_his_frc.nc)
                   (Tmin=          0.0000 Tmax=         62.0000)
                   (Min =  0.00000000E+00 Max =  7.22288387E-03)
    GET_NGFLD   - 2D u-momentum northern boundary condition, t =     0 01:00:00
                   (Rec=0002, Index=1, File: WHsma_his_frc.nc)
                   (Tmin=          0.0000 Tmax=         62.0000)
                   (Min = -2.00623134E-03 Max =  8.42359103E-03)
    GET_NGFLD   - 2D v-momentum northern boundary condition, t =     0 01:00:00
                   (Rec=0002, Index=1, File: WHsma_his_frc.nc)
                   (Tmin=          0.0000 Tmax=         62.0000)
                   (Min = -6.89323479E-03 Max =  5.61088603E-03)

   STEP   Day HH:MM:SS  KINETIC_ENRG   POTEN_ENRG    TOTAL_ENRG    NET_VOLUME

      0     0 00:00:00  0.000000E+00  3.036525E+02  3.036525E+02  3.657815E+11
      DEF_HIS   - creating history file: wh_his_0001.nc
      WRT_HIS   - wrote history  fields (Index=1,1) into time record = 0000001
      1     0 00:00:30  4.906146E-03  3.134780E+02  3.134829E+02  3.543763E+11
      2     0 00:01:00  1.042724E-02  3.135361E+02  3.135466E+02  3.543998E+11
      3     0 00:01:30  1.306311E-02  3.135891E+02  3.136021E+02  3.544185E+11
      4     0 00:02:00  1.714217E-02  3.136400E+02  3.136571E+02  3.544371E+11
      5     0 00:02:30  2.098172E-02  3.136881E+02  3.137091E+02  3.544549E+11
      6     0 00:03:00  2.477972E-02  3.137332E+02  3.137580E+02  3.544719E+11
      7     0 00:03:30  2.851435E-02  3.137755E+02  3.138040E+02  3.544882E+11
      8     0 00:04:00  3.218923E-02  3.138149E+02  3.138470E+02  3.545038E+11
      9     0 00:04:30  3.579667E-02  3.138516E+02  3.138874E+02  3.545187E+11
     10     0 00:05:00  3.934106E-02  3.138856E+02  3.139249E+02  3.545328E+11
     11     0 00:05:30  4.282086E-02  3.139169E+02  3.139597E+02  3.545463E+11
     12     0 00:06:00  4.624059E-02  3.139456E+02  3.139919E+02  3.545591E+11
     13     0 00:06:30  4.960309E-02  3.139722E+02  3.140218E+02  3.545713E+11
     14     0 00:07:00  5.291256E-02  3.139967E+02  3.140497E+02  3.545827E+11
     15     0 00:07:30  5.617092E-02  3.140193E+02  3.140755E+02  3.545935E+11
     16     0 00:08:00  5.937834E-02  3.140398E+02  3.140991E+02  3.546036E+11
     17     0 00:08:30  6.253123E-02  3.140581E+02  3.141206E+02  3.546132E+11
     18     0 00:09:00  6.562480E-02  3.140743E+02  3.141400E+02  3.546221E+11
     19     0 00:09:30  6.865456E-02  3.140884E+02  3.141571E+02  3.546306E+11
     20     0 00:10:00  7.161919E-02  3.141005E+02  3.141721E+02  3.546385E+11
     21     0 00:10:30  7.452147E-02  3.141107E+02  3.141852E+02  3.546459E+11
     22     0 00:11:00  7.736695E-02  3.141194E+02  3.141967E+02  3.546527E+11
     23     0 00:11:30  8.016096E-02  3.141267E+02  3.142068E+02  3.546591E+11
     24     0 00:12:00  8.290590E-02  3.141326E+02  3.142155E+02  3.546650E+11
     25     0 00:12:30  8.560064E-02  3.141372E+02  3.142228E+02  3.546704E+11
     26     0 00:13:00  8.824311E-02  3.141406E+02  3.142288E+02  3.546754E+11
     27     0 00:13:30  9.083395E-02  3.141430E+02  3.142338E+02  3.546799E+11
     28     0 00:14:00  9.337683E-02  3.141449E+02  3.142382E+02  3.546839E+11
     29     0 00:14:30  9.587534E-02  3.141468E+02  3.142427E+02  3.546874E+11
     30     0 00:15:00  9.833086E-02  3.141492E+02  3.142475E+02  3.546905E+11
     31     0 00:15:30  1.007422E-01  3.141519E+02  3.142527E+02  3.546932E+11
     32     0 00:16:00  1.031061E-01  3.141549E+02  3.142580E+02  3.546955E+11
     33     0 00:16:30  1.054188E-01  3.141580E+02  3.142635E+02  3.546973E+11
     34     0 00:17:00  1.076762E-01  3.141612E+02  3.142689E+02  3.546987E+11
     35     0 00:17:30  1.098716E-01  3.141642E+02  3.142741E+02  3.546996E+11
     36     0 00:18:00  1.119955E-01  3.141669E+02  3.142789E+02  3.547001E+11
     37     0 00:18:30  1.140382E-01  3.141691E+02  3.142831E+02  3.547001E+11
     38     0 00:19:00  1.159954E-01  3.141704E+02  3.142864E+02  3.546997E+11
     39     0 00:19:30  1.178718E-01  3.141706E+02  3.142885E+02  3.546988E+11
     40     0 00:20:00  1.196777E-01  3.141698E+02  3.142895E+02  3.546975E+11
      OCN2WAV   - (06) imported and (05) exported fields,    t =     0 00:20:00
                - ROMS coupling exchanges wait clock (s):
                   (Recv= 5.92559940E+01 Send= 3.99900000E-03)
                - ROMS Import: wave direction
                   (Min=  4.47097449E+00 Max=  5.39240799E+00)
                - ROMS Import: significant wave height
                   (Min=  4.95239049E-02 Max=  2.21756533E-01)
                - ROMS Import: average wave length
                   (Min=  1.15064120E+00 Max=  2.99546432E+00)
                - ROMS Import: surface wave relative peak period
                   (Min=  1.21969092E+00 Max=  2.13405609E+00)
                - ROMS Import: bottom wave period
                   (Min=  1.26561022E+00 Max=  8.73337078E+00)
                - ROMS Import: wave bottom orbital velocity
                   (Min=  7.08658257E-16 Max=  3.94234030E-06)
                - ROMS Export: bathymetry
                   (Min=  5.00000000E+00 Max=  9.42617159E+01)
                - ROMS Export: free-surface
                   (Min= -1.86578775E+01 Max=  2.67778400E+00)
                - ROMS Export: vertically integrated u-momentum component
                   (Min= -6.65802467E-01 Max=  8.85190554E-01)
                - ROMS Export: vertically integrated v-momentum component
                   (Min= -1.31424837E-01 Max=  8.62409079E-01)
                - ROMS Export: bottom roughness
                   (Min=  2.50000000E-03 Max=  2.50000000E-03)
      WAV2OCN   - (05) imported and (06) exported fields,    t = 20101020.002000
                - SWAN coupling exchanges wait clock (s):
                   (Recv= 1.84871900E+00 Send= 5.00000000E-03)
     41     0 00:20:30  1.211807E-01  3.141683E+02  3.142895E+02  3.546962E+11
     42     0 00:21:00  1.227867E-01  3.141662E+02  3.142890E+02  3.546944E+11
     43     0 00:21:30  1.244604E-01  3.141639E+02  3.142883E+02  3.546923E+11
     44     0 00:22:00  1.260628E-01  3.141616E+02  3.142877E+02  3.546900E+11
     45     0 00:22:30  1.276106E-01  3.141597E+02  3.142873E+02  3.546875E+11
     46     0 00:23:00  1.290974E-01  3.141583E+02  3.142874E+02  3.546847E+11
     47     0 00:23:30  1.305306E-01  3.141574E+02  3.142879E+02  3.546817E+11
     48     0 00:24:00  1.319321E-01  3.141569E+02  3.142889E+02  3.546785E+11
     49     0 00:24:30  1.333212E-01  3.141567E+02  3.142900E+02  3.546751E+11
     50     0 00:25:00  1.347026E-01  3.141562E+02  3.142909E+02  3.546715E+11
     51     0 00:25:30  1.360710E-01  3.141553E+02  3.142913E+02  3.546679E+11
     52     0 00:26:00  1.374231E-01  3.141534E+02  3.142908E+02  3.546642E+11
     53     0 00:26:30  1.387579E-01  3.141505E+02  3.142892E+02  3.546605E+11
     54     0 00:27:00  1.400645E-01  3.141464E+02  3.142865E+02  3.546569E+11
     55     0 00:27:30  1.413155E-01  3.141413E+02  3.142827E+02  3.546534E+11
     56     0 00:28:00  1.424802E-01  3.141355E+02  3.142780E+02  3.546502E+11
     57     0 00:28:30  1.435472E-01  3.141291E+02  3.142726E+02  3.546473E+11
     58     0 00:29:00  1.445305E-01  3.141222E+02  3.142667E+02  3.546446E+11
     59     0 00:29:30  1.454461E-01  3.141149E+02  3.142604E+02  3.546423E+11
     60     0 00:30:00  1.462924E-01  3.141073E+02  3.142536E+02  3.546402E+11
     61     0 00:30:30  1.470650E-01  3.140994E+02  3.142465E+02  3.546384E+11
     62     0 00:31:00  1.477793E-01  3.140914E+02  3.142391E+02  3.546368E+11
     63     0 00:31:30  1.484604E-01  3.140832E+02  3.142317E+02  3.546354E+11
     64     0 00:32:00  1.491213E-01  3.140751E+02  3.142242E+02  3.546341E+11
     65     0 00:32:30  1.497617E-01  3.140670E+02  3.142168E+02  3.546330E+11
     66     0 00:33:00  1.503763E-01  3.140590E+02  3.142094E+02  3.546321E+11
     67     0 00:33:30  1.509594E-01  3.140512E+02  3.142021E+02  3.546315E+11
     68     0 00:34:00  1.515083E-01  3.140434E+02  3.141949E+02  3.546312E+11
     69     0 00:34:30  1.520320E-01  3.140356E+02  3.141877E+02  3.546312E+11
     70     0 00:35:00  1.525473E-01  3.140279E+02  3.141805E+02  3.546315E+11
     71     0 00:35:30  1.530568E-01  3.140202E+02  3.141733E+02  3.546321E+11
     72     0 00:36:00  1.535527E-01  3.140126E+02  3.141661E+02  3.546331E+11
     73     0 00:36:30  1.540267E-01  3.140050E+02  3.141590E+02  3.546343E+11
     74     0 00:37:00  1.544694E-01  3.139976E+02  3.141521E+02  3.546359E+11
     75     0 00:37:30  1.548757E-01  3.139906E+02  3.141455E+02  3.546380E+11
     76     0 00:38:00  1.552489E-01  3.139841E+02  3.141393E+02  3.546404E+11
     77     0 00:38:30  1.555944E-01  3.139779E+02  3.141335E+02  3.546433E+11
     78     0 00:39:00  1.559167E-01  3.139721E+02  3.141280E+02  3.546467E+11
     79     0 00:39:30  1.562221E-01  3.139664E+02  3.141227E+02  3.546504E+11
     80     0 00:40:00  1.565208E-01  3.139607E+02  3.141172E+02  3.546544E+11

jcwarner
Posts: 1200
Joined: Wed Dec 31, 2003 6:16 pm
Location: USGS, USA

Re: ROMS+SWAN+SEDIMENTS (One-way nesting)

#2 Unread post by jcwarner »

I have continued to advance this modeling system to include one-way grid refinement in roms (ie on-line nesting), one-way grid refinement for swan, coupling on the refined levels, ability to run swan and roms on different grids, coupling to WRF, etc. So this would take care of a lot of your problems. If you want to work with our modeling system, send me an email:
jcwarner@usgs.gov

Post Reply