Incorrect water velocity with open boundaries

General scientific issues regarding ROMS

Moderators: arango, robertson

Post Reply
Message
Author
cloudy4568
Posts: 39
Joined: Thu Oct 17, 2019 8:56 pm
Location: Clarkson University

Incorrect water velocity with open boundaries

#1 Unread post by cloudy4568 »

Hello everyone:

I am simulating the Beaufort Sea with three open boundaries and ICE_MODEL activated, but I found the water velocities terrible.
There is an eddy in the middle and a stronger current along the east boundary, could it be the problem of boundary conditions? I have tried other boundary conditions for 3D velocities like 'RadNud' and 'Gra' but the calculation blows up with 'Diag speed problem'. Do you have any idea about this situation, any help will be appreciated! :D

Image

The ROMS version is Kate ROMS version.
The calculation grid is curvilinear grid with same spacing in kilometers.
The forcing data are from NCEP(CFsV2), the initial and boundary conditions are from HYcom.

Code: Select all

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

     434880  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.
       2880  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).
        720  nHIS              Number of timesteps between the writing fields
                                 into history file.
          0  nQCK              Number of timesteps between the writing fields
                                 into quicksave file.
          1  ntsAVG            Starting timestep for the accumulation of output
                                 time-averaged data.
       2880  nAVG              Number of timesteps between the writing of
                                 time-averaged data into averages file.
          1  ntsAVG2           Starting timestep for the accumulation of output
                                 secondary time-averaged data.
        120  nAVG2             Number of timesteps between the writing of
                                 secondary time-averaged data into averages file.
       2880  ndefAVG2          Number of timesteps between creation of new
                                 secondary time-averaged file.
 1.0000E+00  nl_tnu2(01)       NLM Horizontal, harmonic mixing coefficient
                                 (m2/s) for tracer 01: temp
 1.0000E+00  nl_tnu2(02)       NLM Horizontal, harmonic mixing coefficient
                                 (m2/s) for tracer 02: salt
 5.0000E+00  nl_visc2          NLM Horizontal, harmonic mixing coefficient
                                 (m2/s) for momentum.
          F  LuvSponge         Turning OFF sponge on horizontal momentum.
          F  LtracerSponge(01) Turning OFF sponge on tracer 01: temp
          F  LtracerSponge(02) Turning OFF sponge on tracer 02: salt
 1.0000E-06  Akt_bak(01)       Background vertical mixing coefficient (m2/s)
                                 for tracer 01: temp
 1.0000E-06  Akt_bak(02)       Background vertical mixing coefficient (m2/s)
                                 for tracer 02: salt
 1.0000E-05  Akv_bak           Background vertical mixing coefficient (m2/s)
                                 for momentum.
 3.0000E-04  rdrg              Linear bottom drag coefficient (m/s).
 3.0000E-03  rdrg2             Quadratic bottom drag coefficient.
 2.0000E-02  Zob               Bottom roughness (m).
 2.0000E+00  blk_ZQ            Height (m) of surface air humidity measurement.
 2.0000E+00  blk_ZT            Height (m) of surface air temperature measurement.
 2.0000E+00  blk_ZW            Height (m) of surface winds measurement.
          1  lmd_Jwt           Jerlov water type.
          2  Vtransform        S-coordinate transformation equation.
          4  Vstretching       S-coordinate stretching function.
 7.0000E+00  theta_s           S-coordinate surface control parameter.
 2.0000E+00  theta_b           S-coordinate bottom  control parameter.
    250.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).
     -1.000  tide_start        Reference time origin for tidal forcing (days).
20150802.00  time_ref          Reference time for units attribute (yyyymmdd.dd)
 3.6000E+02  Tnudg(01)         Nudging/relaxation time scale (days)
                                 for tracer 01: temp
 3.6000E+02  Tnudg(02)         Nudging/relaxation time scale (days)
                                 for tracer 02: salt
 9.0000E+01  Tnudg_SSS         Nudging/relaxation time scale (days)
                                 for sea surface salinity.
 3.6000E+02  Znudg             Nudging/relaxation time scale (days)
                                 for free-surface.
 3.6000E+02  M2nudg            Nudging/relaxation time scale (days)
                                 for 2D momentum.
 3.6000E+02  M3nudg            Nudging/relaxation time scale (days)
                                 for 3D momentum.
 1.2000E+02  obcfac            Factor between passive and active
                                 open boundary conditions.
          F  VolCons(1)        NLM western  edge boundary volume conservation.
          F  VolCons(2)        NLM southern edge boundary volume conservation.
          F  VolCons(3)        NLM eastern  edge boundary volume conservation.
          F  VolCons(4)        NLM northern edge boundary volume conservation.
     10.000  T0                Background potential temperature (C) constant.
     35.000  S0                Background salinity (PSU) constant.
      1.000  gamma2            Slipperiness variable: free-slip (1.0) or 
                                                      no-slip (-1.0).
          F  LuvSrc            Turning OFF momentum point Sources/Sinks.
          F  LwSrc             Turning OFF volume influx point Sources/Sinks.
          F  LtracerSrc(01)    Turning OFF point Sources/Sinks on tracer 01: temp
          F  LtracerSrc(02)    Turning OFF point Sources/Sinks on tracer 02: salt
          F  LsshCLM           Turning OFF processing of SSH climatology.
          F  Lm2CLM            Turning OFF processing of 2D momentum climatology.
          F  Lm3CLM            Turning OFF processing of 3D momentum climatology.
          F  LtracerCLM(01)    Turning OFF processing of climatology tracer 01: temp
          F  LtracerCLM(02)    Turning OFF processing of climatology tracer 02: salt
          F  LnudgeM2CLM       Turning OFF nudging of 2D momentum climatology.
          F  LnudgeM3CLM       Turning OFF nudging of 3D momentum climatology.
          F  LnudgeTCLM(01)    Turning OFF nudging of climatology tracer 01: temp
          F  LnudgeTCLM(02)    Turning OFF nudging of climatology tracer 02: salt
          F  LmiCLM            Turning OFF processing of ice momentum climatology.
          F  LaiCLM            Turning OFF processing of ice tracer climatology.
          F  LsiCLM            Turning OFF processing of ice stress climatology.
          F  LnudgeMICLM       Turning OFF nudging of ice momentum climatology.
          F  LnudgeAICLM       Turning OFF nudging of ice tracer climatology.
          F  LnudgeSICLM       Turning OFF nudging of ice stress climatology.

          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(idTsur)      Write out surface net heat flux.
          T  Hout(idTsur)      Write out surface net salt flux.
          T  Hout(idSrad)      Write out shortwave radiation flux.
          T  Hout(idLrad)      Write out longwave radiation flux.
          T  Hout(idLhea)      Write out latent heat flux.
          T  Hout(idShea)      Write out sensible heat flux.
          T  Hout(idEmPf)      Write out E-P flux.
          T  Hout(idevap)      Write out evaporation rate.
          T  Hout(idrain)      Write out rain rate.
          T  Hout(idDano)      Write out density anomaly.
          T  Hout(idVvis)      Write out vertical viscosity: AKv.
          T  Hout(idHsbl)      Write out depth of surface boundary layer.

          T  Aout(idFsur)      Write out averaged free-surface.
          T  Aout(idUbar)      Write out averaged 2D U-momentum component.
          T  Aout(idVbar)      Write out averaged 2D V-momentum component.
          T  Aout(idu2dE)      Write out averaged 2D U-eastward  at RHO-points.
          T  Aout(idv2dN)      Write out averaged 2D V-northward at RHO-points.
          T  Aout(idUvel)      Write out averaged 3D U-momentum component.
          T  Aout(idVvel)      Write out averaged 3D V-momentum component.
          T  Aout(idu3dE)      Write out averaged 3D U-eastward  at RHO-points.
          T  Aout(idv3dN)      Write out averaged 3D V-northward at RHO-points.
          T  Aout(idWvel)      Write out averaged W-momentum component.
          T  Aout(idOvel)      Write out averaged omega vertical velocity.
          T  Aout(idTvar)      Write out averaged tracer 01: temp
          T  Aout(idTvar)      Write out averaged tracer 02: salt
          T  Aout(idUsms)      Write out averaged surface U-momentum stress.
          T  Aout(idVsms)      Write out averaged surface V-momentum stress.
          T  Aout(idUbms)      Write out averaged bottom U-momentum stress.
          T  Aout(idVbms)      Write out averaged bottom V-momentum stress.
          T  Aout(idTsur)      Write out averaged surface net heat flux.
          T  Aout(idTsur)      Write out averaged surface net salt flux.
          T  Aout(idSrad)      Write out averaged shortwave radiation flux.
          T  Aout(idLrad)      Write out averaged longwave radiation flux.
          T  Aout(idLhea)      Write out averaged latent heat flux.
          T  Aout(idShea)      Write out averaged sensible heat flux.
          T  Aout(idevap)      Write out averaged evaporation rate.
          T  Aout(idrain)      Write out averaged rain rate.
          T  Aout(idDano)      Write out averaged density anomaly.
          T  Aout(idHsbl)      Write out averaged depth of surface boundary layer.

          T  Aout2(idFsur)     Write out averaged free-surface.
          T  Aout2(idUbar)     Write out averaged 2D U-momentum component.
          T  Aout2(idVbar)     Write out averaged 2D V-momentum component.
          T  Aout2(idu2dE)     Write out averaged 2D U-eastward  at RHO-points.
          T  Aout2(idv2dN)     Write out averaged 2D V-northward at RHO-points.
          T  Aout2(idUvel)     Write out averaged 3D U-momentum component.
          T  Aout2(idVvel)     Write out averaged 3D V-momentum component.
          T  Aout2(idu3dE)     Write out averaged 3D U-eastward  at RHO-points.
          T  Aout2(idv3dN)     Write out averaged 3D V-northward at RHO-points.
          T  Aout2(idTvar)     Write out averaged tracer 01: temp
          T  Aout2(idTvar)     Write out averaged tracer 02: salt
          T  Aout2(idUsms)     Write out averaged surface U-momentum stress.
          T  Aout2(idVsms)     Write out averaged surface V-momentum stress.
          T  Aout2(idUair)     Write out averaged surface U-wind component.
          T  Aout2(idVair)     Write out averaged surface V-wind component.
          T  Aout2(idUairE)    Write out averaged Eastward surface U-wind component.
          T  Aout2(idVairN)    Write out averaged Northward surface V-wind component.
          T  Aout2(idTsur)     Write out averaged surface net heat flux.
          T  Aout2(idTsur)     Write out averaged surface net salt flux.
          T  Aout2(idSrad)     Write out averaged shortwave radiation flux.
          T  Aout2(idLrad)     Write out averaged longwave radiation flux.
          T  Aout2(idLhea)     Write out averaged latent heat flux.
          T  Aout2(idShea)     Write out averaged sensible heat flux.
          1  shuffle           NetCDF-4/HDF5 file format shuffle filer flag.
          1  deflate           NetCDF-4/HDF5 file format deflate filer flag.
          1  deflate_level     NetCDF-4/HDF5 file format deflate level parameter.

Code: Select all

Ice Parameters, Grid: 01
 ============================

 Ice time step = ocean time step
 3.0000E+01  DTICE         Ice model time step (s).
          1  nstrs         Number of iterations for nonlinear ice dynamics.
         60  nevp          Number of elastic steps per plastic step in EVP.
 9.1900E+02  rhoice        Density of sea ice (kg/m3).
 7.5000E-03  cdiw          Ice-water drag coefficient (nondimensional).
 3.0000E-03  cdai          Air-ice drag coefficient (nondimensional).
 1.2670E+00  rho_air       Air density (kg/m3).
 3.3000E+02  rhosnow_dry   Dry snow density (kg/m3).
 4.5000E+02  rhosnow_wet   Wet snow density (kg/m3).
 4.5000E+01  alphai        Mohr-Coulomb stress angle (degrees).
 0.0000E+00  min_h         Minimum average ice thickness (m).
 1.0000E-30  min_a         Minimum ice concentration (nondimensional).
 1.0000E+00  max_a         Maximum ice concentration (nondimensional).
 0.0000E+00  stressang     Turning angle for ice-water drag (degrees).
 1.0000E-16  tol           Numerical tolerance in rheology calculations .
 9.7000E-01  ice_emiss     Ice emissivity.
 1.0040E+03  spec_heat_air Specific heat of air.
 1.7500E-03  trans_coeff   Transfer coefficient.
 2.8340E+06  sublim_latent_Latent_heat of sublimation.
 2.7315E+02  t0deg         Zero degrees Celsius in degrees Kelvin.
 8.0000E+00  lf_k1         Landfast ice tunable k1.
 1.5000E+01  lf_k2         Landfast ice tunable k2.
 5.0000E-05  lf_u0         Landfast ice minimal velocity.
          T  Hout(idUice)    Write out U-component ice velocity.
          T  Hout(idVice)    Write out V-component ice velocity.
          T  Hout(idUiceE)   Write out East component ice velocity.
          T  Hout(idViceN)   Write out North component ice velocity.
          T  Hout(idAice)    Write out ice concentration.
          T  Hout(idHice)    Write out average ice thickness.
          T  Hout(idHsno)    Write out snow thickness.
          T  Hout(idTice)    Write out ice/snow surface temperature.
          T  Hout(idTimid)   Write out interior ice temperature.
          T  Hout(idT2ice)   Write out top of ice temperature.
          T  Hout(idAgeice)  Write out ice age.
          T  Hout(idIomflx)  Write out ice-ocean mass flux
          T  Hout(idSig11)   Write out internal ice stress component 11.
          T  Hout(idSig12)   Write out internal ice stress component 12.
          T  Hout(idSig22)   Write out internal ice stress component 22.
          T  Hout(idT0mk)    Write out temperature of molecular sublayer under ice.
          T  Hout(idS0mk)    Write out salinity of molecular sublayer under ice.
          T  Hout(idWfr)     Write out frazil ice growth rate.
          T  Hout(idWai)     Write out ice growth/melt rate.
          T  Hout(idWao)     Write out ice growth/melt rate.
          T  Hout(idWio)     Write out ice growth/melt rate.
          T  Hout(idWro)     Write out ice melt runoff rate.
          T  Hout(idWdiv)    Write out ice divergence rate.
          T  Hout(idAlbe)    Write out ocean surface albedo.
          T  Hout(idAlbI)    Write out ice albedo.
          T  Hout(idLWdn)    Write out downward longwave radiation.
          T  Hout(idSWdn)    Write out downward shortwave radiation.
          T  Hout(idLatI)    Write out latent heat flux over ice.
          T  Hout(idSenI)    Write out sensible heat flux over ice.
          T  Hout(idLWupI)   Write out upward longwave over ice.
          T  Hout(idSWupI)   Write out upward shortwave over ice.
          T  Hout(idSWinI)   Write out upward shortwave over ice.
          T  Hout(idSWthruI) Write out shortwave radiation through ice.
          T  Hout(idSaltI)   Write out salt flux under ice.
          T  Hout(idSaltL)   Write out salt flux in leads.
          T  Hout(idQao)     Write out heat flux in leads.
          T  Hout(idQai)     Write out heat flux on top of ice.
          T  Hout(idQio)     Write out heat flux into ice.
          T  Hout(idQi2)     Write out heat flux in ice.
          T  Hout(idSnoI)    Write out heat flux into ice.
          T  Aout(idUice)    Write out time-averaged U-component ice velocity.
          T  Aout(idVice)    Write out time-averaged V-component ice velocity.
          T  Aout(idUiceE)   Write out time-averaged East component ice velocity.
          T  Aout(idViceN)   Write out time-averaged North component ice velocity.
          T  Aout(idAice)    Write out time-averaged ice concentration.
          T  Aout(idHice)    Write out time-averaged average ice thickness.
          T  Aout(idHsno)    Write out time-averaged snow thickness.
          T  Aout(idTice)    Write out time-averaged ice/snow surface temperature.
          T  Aout(idTimid)   Write out time-averaged interior ice temperature.
          T  Aout(idT2ice)   Write out time-averaged top of ice temperature.
          T  Aout(idAgeice)  Write out time-averaged ice age.
          T  Aout(idIomflx)  Write out time-averaged ice-ocean mass flux
          T  Aout(idSig11)   Write out time-averaged internal ice stress component 11.
          T  Aout(idSig12)   Write out time-averaged internal ice stress component 12.
          T  Aout(idSig22)   Write out time-averaged internal ice stress component 22.
          T  Aout(idT0mk)    Write out time-averaged under ice temperature.
          T  Aout(idS0mk)    Write out time-averaged under ice salinity.
          T  Aout(idWfr)     Write out time-averaged frazil ice growth rate.
          T  Aout(idWai)     Write out time-averaged ice growth/melt rate.
          T  Aout(idWao)     Write out time-averaged ice growth/melt rate.
          T  Aout(idWio)     Write out time-averaged ice growth/melt rate.
          T  Aout(idWro)     Write out time-averaged ice melt runoff rate.
          T  Aout(idWdiv)    Write out time-averaged ice divergence rate.
          T  Aout(idAlbe)    Write out ocean surface albedo.
          T  Aout(idAlbI)    Write out ice albedo.
          T  Aout(idLWdn)    Write out downward longwave radiation.
          T  Aout(idSWdn)    Write out downward shortwave radiation.
          T  Aout(idLatI)    Write out latent heat flux over ice.
          T  Aout(idSenI)    Write out sensible heat flux over ice.
          T  Aout(idLWupI)   Write out upward longwave over ice.
          T  Aout(idSWupI)   Write out upward shortwave over ice.
          T  Aout(idSWinI)   Write out upward shortwave over ice.
          T  Aout(idSWthruI) Write out shortwave radiation through ice.
          T  Aout(idSaltI)   Write out salt flux under ice.
          T  Aout(idSaltL)   Write out salt flux in leads.
          T  Aout(idQao)     Write out heat flux in leads.
          T  Aout(idQai)     Write out heat flux on top of ice.
          T  Aout(idQio)     Write out heat flux into ice.
          T  Aout(idQi2)     Write out heat flux in ice.
          T  Aout(idSnoI)    Write out heat flux into ice.
          T  Aout2(idUice)   Write out U-component ice velocity.
          T  Aout2(idVice)   Write out V-component ice velocity.
          T  Aout2(idUiceE)  Write out East component ice velocity.
          T  Aout2(idViceN)  Write out North component ice velocity.
          T  Aout2(idAice)   Write out ice concentration.
          T  Aout2(idHice)   Write out average ice thickness.
          T  Aout2(idHsno)   Write out snow thickness.
          T  Aout2(idTice)   Write out ice/snow surface temperature.
          T  Aout2(idTimid)  Write out interior ice temperature.
          T  Aout2(idAgeice) Write out ice age.
          T  Aout2(idSig11)  Write out internal ice stress component 11.
          T  Aout2(idSig12)  Write out internal ice stress component 12.
          T  Aout2(idSig22)  Write out internal ice stress component 22.
          T  Aout2(idT0mk)   Write out temperature of molecular sublayer under ice.
          T  Aout2(idS0mk)   Write out salinity of molecular sublayer under ice.
          T  Qout(idUice)    Write out U-component ice velocity.
          T  Qout(idVice)    Write out V-component ice velocity.
          T  Qout(idUiceE)   Write out East component ice velocity.
          T  Qout(idViceN)   Write out North component ice velocity.
          T  Qout(idAice)    Write out ice concentration.
          T  Qout(idHice)    Write out average ice thickness.
          T  Qout(idHsno)    Write out snow thickness.
          T  Qout(idTice)    Write out ice/snow surface temperature.
          T  Qout(idTimid)   Write out interior ice temperature.
          T  Qout(idAgeice)  Write out ice age.

Output/Input Files:
               Output Restart File:  beaufort_rst.nc
               Output History File:  beaufort_his.nc
              Output Averages File:  beaufort_avg.nc
       Prefix for Averages Files:  beaufort_avg2
                   Input Grid File:  /mnt/home/cyue/roms/bohai/dataarc/beaufort2.nc
      Input Nonlinear Initial File:  /mnt/home/cyue/roms/bohai/dataarc/initial_beaufort3.nc
                Tidal Forcing File:  /mnt/home/cyue/roms/bohai/dataarc/beauforttide2.nc
             Input Forcing File 01:  /mnt/home/cyue/roms/bohai/dataarc/gom_Pair_era.nc
             Input Forcing File 02:  /mnt/home/cyue/roms/bohai/dataarc/gom_Qair_era.nc
             Input Forcing File 03:  /mnt/home/cyue/roms/bohai/dataarc/gom_Tair_era.nc
             Input Forcing File 04:  /mnt/home/cyue/roms/bohai/dataarc/gom_rain_era.nc
             Input Forcing File 05:  /mnt/home/cyue/roms/bohai/dataarc/gom_swrad_era.nc
             Input Forcing File 06:  /mnt/home/cyue/roms/bohai/dataarc/gom_windU_era01.nc
             Input Forcing File 07:  /mnt/home/cyue/roms/bohai/dataarc/gom_windV_era01.nc
             Input Forcing File 08:  /mnt/home/cyue/roms/bohai/dataarc/gom_lwrad_era.nc
             Input Forcing File 09:  /mnt/home/cyue/roms/bohai/dataarc/albedobea.nc
             Input Forcing File 10:  /mnt/home/cyue/roms/bohai/dataarc/sss.nc
         Input Boundary File 01:  /mnt/home/cyue/roms/bohai/dataarc/boundary-beaufort4.nc
         Input Boundary File 02:  /mnt/home/cyue/roms/bohai/dataarc/boundary-beaufort-ice.nc
       Ice Model Parameters File:  /mnt/home/cyue/roms/beau/ice.in

 Tile partition information for Grid 01:  199x99x30  tiling: 6x6

Code: Select all

 Lateral Boundary Conditions: NLM
 ============================

 Variable               Grid  West Edge    South Edge   East Edge    North Edge
 ---------              ----  ----------   ----------   ----------   ----------

 zeta                     1   Chapman Exp  Closed       Chapman Exp  Chapman Exp

 ubar                     1   Shchepetkin  Closed       Shchepetkin  Shchepetkin

 vbar                     1   Shchepetkin  Closed       Shchepetkin  Shchepetkin

 u                        1   Clamped      Closed       Clamped      Clamped

 v                        1   Clamped      Closed       Clamped      Clamped

 temp                     1   Rad + Nud    Closed       Rad + Nud    Rad + Nud

 salt                     1   Rad + Nud    Closed       Rad + Nud    Rad + Nud

 uice                     1   Gradient     Closed       Gradient     Gradient

 vice                     1   Gradient     Closed       Gradient     Gradient

 aice                     1   Chapman Exp  Closed       Chapman Exp  Chapman Exp

 hice                     1   Chapman Exp  Closed       Chapman Exp  Chapman Exp

 tisrf                    1   Gradient     Closed       Gradient     Gradient

 snow_thick               1   Gradient     Closed       Gradient     Gradient

 apond                    1   Closed       Closed       Closed       Closed

 hpond                    1   Closed       Closed       Closed       Closed

 sig11                    1   Gradient     Closed       Gradient     Gradient

 sig12                    1   Gradient     Closed       Gradient     Gradient

 sig22                    1   Gradient     Closed       Gradient     Gradient

 Activated C-preprocessing Options:

 BEAUFORT                BEAUFORT Seas
 ADD_FSOBC               Adding tidal elevation to processed OBC data
 ADD_M2OBC               Adding tidal currents to processed OBC data
 ALBEDO_CSIM             Ice albedo computed according to CSIM3
 ALBEDO_CURVE            Ocean shortwave radiation adjusted by albedo curve
 ANA_BSFLUX              Analytical kinematic bottom salinity flux
 ANA_BTFLUX              Analytical kinematic bottom temperature flux
 ANA_DRAG_GRID           Analytical spatially varying quadratic drag coefficient
 ASSUMED_SHAPE           Using assumed-shape arrays
 AVERAGES                Writing out time-averaged nonlinear model fields
 AVERAGES2               Writing out secondary time-averaged fields
 !BOUNDARY_ALLGATHER     Using mpi_allreduce in mp_boundary routine
 BULK_FLUXES             Surface bulk fluxes parameterization
 CCSM_FLUXES             Surface bulk fluxes parameterization
 !COLLECT_ALL...         Using mpi_isend/mpi_recv in mp_collect routine
 CORE_FORCING            Specific humidity from CORE forcing files
 CURVGRID                Orthogonal curvilinear grid
 DEFLATE                 Setting compression in output NetCDF-4/HDF5 files
 DIFF_GRID               Horizontal diffusion coefficient scaled by grid size
 DJ_GRADPS               Parabolic Splines density Jacobian (Shchepetkin, 2002)
 DOUBLE_PRECISION        Double precision arithmetic numerical kernel
 EMINUSP                 Compute Salt Flux using E-P
 FILTERED                Writing out time-filtered output
 FILTERED_RST            Reading and writing time-filtered restart files
 HDF5                    Creating NetCDF-4/HDF5 format files
 ICE_ADVECT              Advection of ice tracers
 ICE_BULK_FLUXES         Ice bulk fluxes from the atmosphere
 ICE_CONVSNOW            Conversion of flooded snow to ice
 ICE_DIAGS               Include extra diagnostic output from sea ice
 ICE_EVP                 Elastic-viscous-plastic ice rheology
 ICE_I_O                 Put shortwave heat into ice interior
 ICE_LANDFAST            Lemieux shorefast ice
 ICE_MK                  Mellor-Kantha ice thermodynamics
 ICE_MODEL               Include sea ice model
 ICE_MOMENTUM            Compute ice momentum equations
 ICE_SMOLAR              Advect ice tracers with MPDATA scheme
 ICE_STRENGTH_QUAD       Ice strength a quadratic function of thickness
 ICE_THERMO              Include ice thermodynamics
 LIMIT_BSTRESS           Limit bottom stress to maintain bottom velocity direction
 LMD_BKPP                KPP bottom boundary layer mixing
 LMD_CONVEC              LMD convective mixing due to shear instability
 LMD_DDMIX               LMD double-diffusive mixing
 LMD_MIXING              Large/McWilliams/Doney interior mixing
 LMD_NONLOCAL            LMD convective nonlocal transport
 LMD_RIMIX               LMD diffusivity due to shear instability
 LMD_SHAPIRO             Shapiro filtering boundary layer depth
 LMD_SKPP                KPP surface boundary layer mixing
 LONGWAVE_OUT            Compute outgoing longwave radiation internally
 MASKING                 Land/Sea masking
 MIX_GEO_TS              Mixing of tracers along geopotential surfaces
 MIX_S_UV                Mixing of momentum along constant S-surfaces
 MPI                     MPI distributed-memory configuration
 NONLINEAR               Nonlinear Model
 NONLIN_EOS              Nonlinear Equation of State for seawater
 NO_WRITE_GRID           Not Writing grid arrays into NetCDF ouput files
 POWER_LAW               Power-law shape time-averaging barotropic filter
 PROFILE                 Time profiling activated
 RADIATION_2D            Use tangential phase speed in radiation conditions
 REDUCE_ALLGATHER        Using mpi_allgather in mp_reduce routine
 RI_SPLINES              Parabolic Spline Reconstruction for Richardson Number
 RST_SINGLE              Single precision fields in restart NetCDF file
 SALINITY                Using salinity
 SCORRECTION             Surface salinity flux correction
 SOLAR_SOURCE            Solar Radiation Source Term
 SOLVE3D                 Solving 3D Primitive Equations
 SPLINES_VDIFF           Parabolic Spline Reconstruction for Vertical Diffusion
 SPLINES_VVISC           Parabolic Spline Reconstruction for Vertical Viscosity
 SSH_TIDES               Add tidal elevation to SSH climatology
 TIDES_ASTRO             Compute astronomical (long-period tide) corrections
 TS_U3HADVECTION         Third-order upstream horizontal advection of tracers
 TS_C4VADVECTION         Fourth-order centered vertical advection of tracers
 TS_DIF2                 Harmonic mixing of tracers
 TWO_D_TRACER_SOURCE     One value per tracer per source
 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
 UV_DRAG_GRID            Spatially varying quadratic drag coefficient
 UV_QDRAG                Quadratic bottom stress
 UV_TIDES                Add tidal currents to 2D momentum climatologies
 UV_VIS2                 Harmonic mixing of momentum
 VAR_RHO_2D              Variable density barotropic mode
 VISC_GRID               Horizontal viscosity coefficient scaled by grid size

Code: Select all

Vertical S-coordinate System, Grid 01:

 level   S-coord     Cs-curve   Z   at hmin       at hc    half way     at hmax

    30   0.0000000   0.0000000        0.000       0.000       0.000       0.000
    29  -0.0333333  -0.0001156       -0.066      -4.181      -7.569      -8.240
    28  -0.0666667  -0.0004685       -0.132      -8.392     -15.541     -17.334
    27  -0.1000000  -0.0010779       -0.198     -12.635     -23.948     -27.352
    26  -0.1333333  -0.0019766       -0.265     -16.914     -32.847     -38.411
    25  -0.1666667  -0.0032130       -0.331     -21.235     -42.318     -50.687
    24  -0.2000000  -0.0048536       -0.397     -25.607     -52.476     -64.418
    23  -0.2333333  -0.0069862       -0.463     -30.040     -63.469     -79.921
    22  -0.2666667  -0.0097241       -0.529     -34.549     -75.490     -97.603
    21  -0.3000000  -0.0132122       -0.595     -39.152     -88.784    -117.987
    20  -0.3333333  -0.0176335       -0.662     -43.871    -103.663    -141.731
    19  -0.3666667  -0.0232170       -0.728     -48.735    -120.514    -169.660
    18  -0.4000000  -0.0302484       -0.794     -53.781    -139.824    -202.803
    17  -0.4333333  -0.0390807       -0.860     -59.052    -162.190    -242.431
    16  -0.4666667  -0.0501474       -0.927     -64.602    -188.350    -290.105
    15  -0.5000000  -0.0639753       -0.993     -70.497    -219.198    -347.722
    14  -0.5333333  -0.0811979       -1.059     -76.816    -255.808    -417.563
    13  -0.5666667  -0.1025651       -1.126     -83.654    -299.455    -502.327
    12  -0.6000000  -0.1289464       -1.193     -91.118    -351.613    -605.148
    11  -0.6333333  -0.1613228       -1.259     -99.332    -413.949    -729.557
    10  -0.6666667  -0.2007569       -1.326    -108.428    -488.266    -879.379
     9  -0.7000000  -0.2483317       -1.393    -118.541    -576.404   -1058.516
     8  -0.7333333  -0.3050421       -1.460    -129.797    -680.050   -1270.550
     7  -0.7666667  -0.3716257       -1.527    -142.287    -800.457   -1518.136
     6  -0.8000000  -0.4483212       -1.594    -156.040    -938.031   -1802.134
     5  -0.8333333  -0.5345589       -1.662    -170.987   -1091.803   -2120.493
     4  -0.8666667  -0.6286208       -1.730    -186.911   -1258.859   -2467.027
     3  -0.9000000  -0.7273553       -1.797    -203.419   -1433.846   -2830.386
     2  -0.9333333  -0.8261019       -1.865    -219.929   -1608.854   -3193.789
     1  -0.9666667  -0.9190220       -1.933    -235.711   -1773.971   -3536.211
     0  -1.0000000  -1.0000000       -2.000    -250.000   -1918.815   -3835.630

 Time Splitting Weights for Grid 01:    ndtfast =  30    nfast =  42
 ==================================

    Primary            Secondary            Accumulated to Current Step

  1-0.0008094437383769 0.0333333333333333-0.0008094437383769 0.0333333333333333
  2-0.0014053566728197 0.0333603147912792-0.0022148004111966 0.0666936481246126
  3-0.0017877524645903 0.0334071600137066-0.0040025528757869 0.1001008081383191
  4-0.0019566842408176 0.0334667517625262-0.0059592371166046 0.1335675599008453
  5-0.0019122901320372 0.0335319745705535-0.0078715272486418 0.1670995344713988
  6-0.0016548570247459 0.0335957175749547-0.0095263842733877 0.2006952520463536
  7-0.0011849025289723 0.0336508794757796-0.0107112868023600 0.2343461315221331
  8-0.0005032751608631 0.0336903762267453-0.0112145619632232 0.2680365077488784
  9 0.0003887272597151 0.0337071520654408-0.0108258347035081 0.3017436598143192
 10 0.0014892209965583 0.0336941944901169-0.0093366137069498 0.3354378543044362
 11 0.0027955815694920 0.0336445537902317-0.0065410321374578 0.3690824080946679
 12 0.0043042707117221 0.0335513677379153-0.0022367614257357 0.4026337758325831
 13 0.0060106451121704 0.0334078920475245 0.0037738836864347 0.4360416678801076
 14 0.0079087469427945 0.0332075372104522 0.0116826306292293 0.4692492050905598
 15 0.0099910761708920 0.0329439123123590 0.0216737068001212 0.5021931174029188
 16 0.0122483446563884 0.0326108764399960 0.0339220514565096 0.5348039938429148
 17 0.0146692120341107 0.0322025982847830 0.0485912634906203 0.5670065921276978
 18 0.0172400033810439 0.0317136245503127 0.0658312668716642 0.5987202166780105
 19 0.0199444086685725 0.0311389577709445 0.0857756755402367 0.6298591744489550
 20 0.0227631639997064 0.0304741441486588 0.1085388395399431 0.6603333185976138
 21 0.0256737146312910 0.0297153720153352 0.1342125541712341 0.6900486906129490
 22 0.0286498597812016 0.0288595815276255 0.1628624139524357 0.7189082721405746
 23 0.0316613792205220 0.0279045862015855 0.1945237931729577 0.7468128583421600
 24 0.0346736416507075 0.0268492068942347 0.2291974348236651 0.7736620652363948
 25 0.0376471948657328 0.0256934188392112 0.2668446296893979 0.7993554840756060
 26 0.0405373376992232 0.0244385123436867 0.3073819673886211 0.8237939964192927
 27 0.0432936737565711 0.0230872677537126 0.3506756411451922 0.8468812641730054
 28 0.0458596469320356 0.0216441452951603 0.3965352880772278 0.8685254094681656
 29 0.0481720587108284 0.0201154903974257 0.4447073467880562 0.8886408998655914
 30 0.0501605672561820 0.0185097551070648 0.4948679140442381 0.9071506549726561
 31 0.0517471682814030 0.0168377361985254 0.5466150823256412 0.9239883911711815
 32 0.0528456577069106 0.0151128305891453 0.5994607400325517 0.9391012217603267
 33 0.0533610761022577 0.0133513086655816 0.6528218161348094 0.9524525304259084
 34 0.0531891349131379 0.0115726061288397 0.7060109510479473 0.9640251365547481
 35 0.0522156244733761 0.0097996349650684 0.7582265755213234 0.9738247715198165
 36 0.0503158038019031 0.0080591141492892 0.8085423793232265 0.9818838856691057
 37 0.0473537721847153 0.0063819206892258 0.8558961515079418 0.9882658063583315
 38 0.0431818225418188 0.0048034616164019 0.8990779740497606 0.9930692679747335
 39 0.0376397765791564 0.0033640675316746 0.9367177506289170 0.9964333355064081
 40 0.0305543017255206 0.0021094083123694 0.9672720523544377 0.9985427438187775
 41 0.0217382098544505 0.0010909315881854 0.9890102622088881 0.9996336754069629
 42 0.0109897377911119 0.0003663245930371 1.0000000000000000 1.0000000000000000

 ndtfast, nfast =   30  42   nfast/ndtfast =  1.40000

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

    1.000000000000 1.047601458608 0.523800729304 1.000000000000 1.000000000000

 Power filter parameters, Fgamma, gamma =  0.28400   0.18933

 Metrics information for Grid 01:
 ===============================

 Minimum X-grid spacing, DXmin =  3.65231759E+00 km
 Maximum X-grid spacing, DXmax =  4.16564098E+00 km
 Minimum Y-grid spacing, DYmin =  1.11176924E+01 km
 Maximum Y-grid spacing, DYmax =  1.11911327E+01 km
 Minimum Z-grid spacing, DZmin =  6.61394005E-02 m
 Maximum Z-grid spacing, DZmax =  3.63402678E+02 m

 Minimum barotropic Courant Number =  1.77333713E-03
 Maximum barotropic Courant Number =  5.51830481E-02
 Maximum Coriolis   Courant Number =  4.27841684E-03

 Horizontal mixing scaled by grid size, GRDMAX =  6.80603442E+00 km
  Maximum grid stiffness ratios:  rx0 =   7.132708E-01 (Beckmann and Haidvogel)
                                 rx1 =   2.629503E+01 (Haney)

 Initial domain volumes:  TotVolume =  1.2559256133E+15 m3
                         MinCellVol =  6.8905412947E+06 m3
                         MaxCellVol =  1.5705931091E+10 m3
                            Max/Min =  2.2793464865E+03

 NL ROMS/TOMS: started time-stepping: (Grid: 01 TimeSteps: 000000000001 - 000000434880)

    GET_2DFLD   - surface u-wind component,                        2015-08-01 02:00:00.00
                   (Grid=01, Rec=0000002, Index=2, File: gom_windU_era01.nc)
                   (Tmin=         -1.0000 Tmax=        152.0000)      t =         -0.9167
                   (Min = -4.76425726E-01 Max =  1.14701714E+00)      regrid = T
    GET_2DFLD   - surface v-wind component,                        2015-08-01 02:00:00.00
                   (Grid=01, Rec=0000002, Index=2, File: gom_windV_era01.nc)
                   (Tmin=         -1.0000 Tmax=        152.0000)      t =         -0.9167
                   (Min = -9.72492107E-01 Max =  9.75210293E-01)      regrid = T
    GET_2DFLD   - surface air pressure,                            2015-08-01 02:00:00.00
                   (Grid=01, Rec=0000002, Index=2, File: gom_Pair_era.nc)
                   (Tmin=         -1.0000 Tmax=        152.0000)      t =         -0.9167
                   (Min =  8.36178548E+02 Max =  1.01244477E+03)      regrid = T
    GET_2DFLD   - surface air temperature,                         2015-08-01 02:00:00.00
                   (Grid=01, Rec=0000002, Index=2, File: gom_Tair_era.nc)
                   (Tmin=         -1.0000 Tmax=        152.0000)      t =         -0.9167
                   (Min = -1.28274624E+00 Max =  1.99322102E+01)      regrid = T
    GET_2DFLD   - surface air relative humidity,                   2015-08-01 02:00:00.00
                   (Grid=01, Rec=0000002, Index=2, File: gom_Qair_era.nc)
                   (Tmin=         -1.0000 Tmax=        152.0000)      t =         -0.9167
                   (Min =  2.95494827E-03 Max =  7.15484283E-03)      regrid = T
    GET_2DFLD   - solar shortwave radiation flux,                  2015-08-01 02:00:00.00
                   (Grid=01, Rec=0000002, Index=2, File: gom_swrad_era.nc)
                   (Tmin=         -1.0000 Tmax=        152.0000)      t =         -0.9167
                   (Min =  5.21664204E-06 Max =  1.33325287E-04)      regrid = T
    GET_2DFLD   - downwelling longwave radiation flux,             2015-08-01 02:00:00.00
                   (Grid=01, Rec=0000002, Index=2, File: gom_lwrad_era.nc)
                   (Tmin=         -1.0000 Tmax=        152.0000)      t =         -0.9167
                   (Min =  5.65535392E-05 Max =  8.76457447E-05)      regrid = T
    GET_2DFLD   - rain fall rate,                                  2015-08-01 02:00:00.00
                   (Grid=01, Rec=0000002, Index=2, File: gom_rain_era.nc)
                   (Tmin=         -1.0000 Tmax=        152.0000)      t =         -0.9167
                   (Min =  0.00000000E+00 Max =  0.00000000E+00)      regrid = T
    GET_2DFLD   - sea surface salinity climatology,                2015-08-02 00:00:00.00
                   (Grid=01, Rec=0000002, Index=2, File: sss.nc)
                   (Tmin=         -1.0000 Tmax=        150.3750)      t =          0.0000
                   (Min =  1.79424982E+01 Max =  3.14597013E+01)      regrid = T
    GET_NGFLD   - free-surface western boundary condition,         2015-08-02 00:00:00.00
                   (Grid= 01, Rec=0000002, Index=1, File: boundary-beaufort4.nc)
                   (Tmin=         -1.0000 Tmax=        151.0000)      t =          0.0000
                   (Min = -4.68922317E-01 Max =  0.00000000E+00)
    GET_NGFLD   - free-surface eastern boundary condition,         2015-08-02 00:00:00.00
                   (Grid= 01, Rec=0000002, Index=1, File: boundary-beaufort4.nc)
                   (Tmin=         -1.0000 Tmax=        151.0000)      t =          0.0000
                   (Min = -3.05613518E-01 Max =  0.00000000E+00)
    GET_NGFLD   - free-surface northern boundary condition,        2015-08-02 00:00:00.00
                   (Grid= 01, Rec=0000002, Index=1, File: boundary-beaufort4.nc)
                   (Tmin=         -1.0000 Tmax=        151.0000)      t =          0.0000
                   (Min = -3.00482035E-01 Max = -2.36000016E-01)
    GET_NGFLD   - 2D u-momentum western boundary condition,        2015-08-02 00:00:00.00
                   (Grid= 01, Rec=0000002, Index=1, File: boundary-beaufort4.nc)
                   (Tmin=         -1.0000 Tmax=        151.0000)      t =          0.0000
                   (Min = -1.01783112E-01 Max =  7.59038180E-02)
    GET_NGFLD   - 2D v-momentum western boundary condition,        2015-08-02 00:00:00.00
                   (Grid= 01, Rec=0000002, Index=1, File: boundary-beaufort4.nc)
                   (Tmin=         -1.0000 Tmax=        151.0000)      t =          0.0000
                   (Min = -4.07351963E-02 Max =  9.38245356E-02)
    GET_NGFLD   - 2D u-momentum eastern boundary condition,        2015-08-02 00:00:00.00
                   (Grid= 01, Rec=0000002, Index=1, File: boundary-beaufort4.nc)
                   (Tmin=         -1.0000 Tmax=        151.0000)      t =          0.0000
                   (Min = -1.16973929E-01 Max =  1.15304545E-01)
    GET_NGFLD   - 2D v-momentum eastern boundary condition,        2015-08-02 00:00:00.00
                   (Grid= 01, Rec=0000002, Index=1, File: boundary-beaufort4.nc)
                   (Tmin=         -1.0000 Tmax=        151.0000)      t =          0.0000
                   (Min = -2.65568998E-02 Max =  2.80445777E-02)
    GET_NGFLD   - 2D u-momentum northern boundary condition,       2015-08-02 00:00:00.00
                   (Grid= 01, Rec=0000002, Index=1, File: boundary-beaufort4.nc)
                   (Tmin=         -1.0000 Tmax=        151.0000)      t =          0.0000
                   (Min = -4.17499617E-02 Max =  4.32500280E-02)
    GET_NGFLD   - 2D v-momentum northern boundary condition,       2015-08-02 00:00:00.00
                   (Grid= 01, Rec=0000002, Index=1, File: boundary-beaufort4.nc)
                   (Tmin=         -1.0000 Tmax=        151.0000)      t =          0.0000
                   (Min = -1.18500076E-01 Max = -8.34892225E-03)
    GET_NGFLD   - 3D u-momentum western boundary condition,        2015-08-02 00:00:00.00
                   (Grid= 01, Rec=0000002, Index=1, File: boundary-beaufort4.nc)
                   (Tmin=         -1.0000 Tmax=        151.0000)      t =          0.0000
                   (Min = -1.73000008E-01 Max =  9.70000029E-02)
    GET_NGFLD   - 3D v-momentum western boundary condition,        2015-08-02 00:00:00.00
                   (Grid= 01, Rec=0000002, Index=1, File: boundary-beaufort4.nc)
                   (Tmin=         -1.0000 Tmax=        151.0000)      t =          0.0000
                   (Min = -6.16486706E-02 Max =  1.14790864E-01)
    GET_NGFLD   - 3D u-momentum eastern boundary condition,        2015-08-02 00:00:00.00
                   (Grid= 01, Rec=0000002, Index=1, File: boundary-beaufort4.nc)
                   (Tmin=         -1.0000 Tmax=        151.0000)      t =          0.0000
                   (Min = -3.04206371E-01 Max =  1.92375675E-01)
    GET_NGFLD   - 3D v-momentum eastern boundary condition,        2015-08-02 00:00:00.00
                   (Grid= 01, Rec=0000002, Index=1, File: boundary-beaufort4.nc)
                   (Tmin=         -1.0000 Tmax=        151.0000)      t =          0.0000
                   (Min = -1.07987799E-01 Max =  6.88706189E-02)
    GET_NGFLD   - 3D u-momentum northern boundary condition,       2015-08-02 00:00:00.00
                   (Grid= 01, Rec=0000002, Index=1, File: boundary-beaufort4.nc)
                   (Tmin=         -1.0000 Tmax=        151.0000)      t =          0.0000
                   (Min = -8.26295018E-02 Max =  1.32801592E-01)
    GET_NGFLD   - 3D v-momentum northern boundary condition,       2015-08-02 00:00:00.00
                   (Grid= 01, Rec=0000002, Index=1, File: boundary-beaufort4.nc)
                   (Tmin=         -1.0000 Tmax=        151.0000)      t =          0.0000
                   (Min = -1.41603172E-01 Max =  1.18677234E-02)
    GET_NGFLD   - temperature western boundary condition,          2015-08-02 00:00:00.00
                   (Grid= 01, Rec=0000002, Index=1, File: boundary-beaufort4.nc)
                   (Tmin=         -1.0000 Tmax=        151.0000)      t =          0.0000
                   (Min = -1.54173160E+00 Max =  1.37320738E+01)
    GET_NGFLD   - salinity western boundary condition,             2015-08-02 00:00:00.00
                   (Grid= 01, Rec=0000002, Index=1, File: boundary-beaufort4.nc)
                   (Tmin=         -1.0000 Tmax=        151.0000)      t =          0.0000
                   (Min =  2.63740005E+01 Max =  3.49451599E+01)
    GET_NGFLD   - temperature eastern boundary condition,          2015-08-02 00:00:00.00
                   (Grid= 01, Rec=0000002, Index=1, File: boundary-beaufort4.nc)
                   (Tmin=         -1.0000 Tmax=        151.0000)      t =          0.0000
                   (Min = -1.58200562E+00 Max =  7.52923489E+00)
    GET_NGFLD   - salinity eastern boundary condition,             2015-08-02 00:00:00.00
                   (Grid= 01, Rec=0000002, Index=1, File: boundary-beaufort4.nc)
                   (Tmin=         -1.0000 Tmax=        151.0000)      t =          0.0000
                   (Min =  2.59456940E+01 Max =  3.49460564E+01)
    GET_NGFLD   - temperature northern boundary condition,         2015-08-02 00:00:00.00
                   (Grid= 01, Rec=0000002, Index=1, File: boundary-beaufort4.nc)
                   (Tmin=         -1.0000 Tmax=        151.0000)      t =          0.0000
                   (Min = -1.61400104E+00 Max = -4.66918051E-01)
    GET_NGFLD   - salinity northern boundary condition,            2015-08-02 00:00:00.00
                   (Grid= 01, Rec=0000002, Index=1, File: boundary-beaufort4.nc)
                   (Tmin=         -1.0000 Tmax=        151.0000)      t =          0.0000
                   (Min =  2.86680660E+01 Max =  2.99954510E+01)
 

User avatar
kate
Posts: 4088
Joined: Wed Jul 02, 2003 5:29 pm
Location: CFOS/UAF, USA

Re: Incorrect water velocity with open boundaries

#2 Unread post by kate »

I assume you know what I used for a regional Beaufort domain. I would go back to the RadNud option and also use nudging to climatology. If it's blowing up due to high speeds, where is that happening? Is that the usual vertical instability needing a short, short timestep or something else?

As for your "Beaufort gyre", is it wind-driven? Do you have something crazy in the ice velocities? What are you doing about the ice boundary conditions? 2015 is before you can download daily CICE fields from the HYCOM runs.

cloudy4568
Posts: 39
Joined: Thu Oct 17, 2019 8:56 pm
Location: Clarkson University

Re: Incorrect water velocity with open boundaries

#3 Unread post by cloudy4568 »

Thank you for replying, Kate!

Yes, I referred most settings to your regional Beaufort domain. I have tried the RadNud option, the blowing up of high speed happens at the north-east corner (north and east are both open boundaries) and middle in depth, the speed of surrounding points are also high. The timestep of 3D equations is 30 seconds, I think it is short enough.

I also have tried different experiments to solve the 'Beaufort gyre', like: turn off tides, reduce the wind velocities artificially, and try other 3D momentum boundary options, but the 'Beaufort gyre' always exists. The ice velocities are similar to the water velocities. The ice concentration and ice thickness boundary conditions are from NCEP-CFsV2 product, other ice boundary conditions are Gradient. The only thing I want to do now is to eliminate the 'Beaufort gyre'. :?

User avatar
kate
Posts: 4088
Joined: Wed Jul 02, 2003 5:29 pm
Location: CFOS/UAF, USA

Re: Incorrect water velocity with open boundaries

#4 Unread post by kate »

Is this a time when there should be full pack ice preventing the gyre? I had to have really strong nudging on the ice fields to make things sort of behave.

Post Reply