Sea-Ice Output

Discussion about modeling ice with ROMS

Moderators: arango, robertson

Post Reply
Message
Author
lushilima
Posts: 5
Joined: Tue Apr 17, 2018 4:59 pm
Location: INPE

Sea-Ice Output

#1 Unread post by lushilima »

Hi all,
I am using Kate's version of ROMS (https://github.com/kshedstrom/roms.git) - (github - Jun/2021), trying to run with ICE_MODEL and ICESHELF activated. And the sea-ice module that I am using is from Budgell.
And I have configured:

Code: Select all

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

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

 zeta                     1   Chapman Imp  Closed       Chapman Imp  Chapman Imp

 ubar                     1   Flather      Closed       Flather      Flather

 vbar                     1   Flather      Closed       Flather      Flather

 u                        1   Rad + Nud    Closed       Rad + Nud    Rad + Nud

 v                        1   Rad + Nud    Closed       Rad + Nud    Rad + Nud

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

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

 tke                      1   Radiation    Closed       Radiation    Radiation

 uice                     1   Rad + Nud    Closed       Rad + Nud    Rad + Nud

 vice                     1   Rad + Nud    Closed       Rad + Nud    Rad + Nud

 aice                     1   Mixed        Closed       Mixed        Mixed

 hice                     1   Mixed        Closed       Mixed        Mixed

 tisrf                    1   Gradient     Closed       Gradient     Gradient

 snow_thick               1   Mixed        Closed       Mixed        Mixed

 sig11                    1   Gradient     Closed       Gradient     Gradient

 sig12                    1   Gradient     Closed       Gradient     Gradient

 sig22                    1   Gradient     Closed       Gradient     Gradient

 Activated C-preprocessing Options:


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
 ASSUMED_SHAPE           Using assumed-shape arrays
 AVERAGES                Writing out time-averaged nonlinear model fields
 !BOUNDARY_ALLGATHER     Using mpi_allreduce in mp_boundary routine
 BULK_FLUXES             Surface bulk fluxes parameterization
 CCSM_FLUXES             Surface bulk fluxes parameterization
 CCSM_ICE_SHORTWAVE      Surface bulk fluxes parameterization
 !COLLECT_ALL...         Using mpi_isend/mpi_recv in mp_collect routine
 CHARNOK                 Charnok surface roughness from wind stress
 CRAIG_BANNER            Craig and Banner wave breaking surface flux
 CORE_FORCING            Specific humidity from CORE forcing files
 CURVGRID                Orthogonal curvilinear grid
 DEFLATE                 Setting compression in output NetCDF-4/HDF5 files
 DIAGNOSTICS_TS          Computing and writing tracer diagnostic terms
 DIAGNOSTICS_UV          Computing and writing momentum diagnostic terms
 DJ_GRADPS               Parabolic Splines density Jacobian (Shchepetkin, 2002)
 DOUBLE_PRECISION        Double precision arithmetic numerical kernel
 EMINUSP                 Compute Salt Flux using E-P
 GLS_MIXING              Generic Length-Scale turbulence closure
 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_EVP                 Elastic-viscous-plastic ice rheology
 ICE_I_O                 Put shortwave heat into ice interior
 ICE_MK                  Mellor-Kantha ice thermodynamics
 ICE_MODEL               Include sea ice model
 ICE_MOMENTUM            Compute ice momentum equations
 ICE_SHALLOW_LIMIT       Grounded ice doesn't grow thermodynamically.
 ICE_SMOLAR              Advect ice tracers with MPDATA scheme
 ICE_STRENGTH_QUAD       Ice strength a quadratic function of thickness
 ICE_THERMO              Include ice thermodynamics
 ICESHELF                Include Ice Shelf Cavities
 ICESHELF_3EQ            Include 3eq Ice Shelf Thermodynamics
 KANTHA_CLAYSON          Kantha and Clayson stability function formulation
LONGWAVE_OUT            Compute outgoing longwave radiation internally
 MASKING                 Land/Sea masking
 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_HIS                  Turn off writing of history files
 NO_LBC_ATT              Not checking NetCDF global attribute NLM_LBC during restart
 NO_SCORRECTION_ICE      Not including SCORRECTION under the ice
 NO_WRITE_GRID           Not Writing grid arrays into NetCDF ouput files
 N2S2_HORAVG             Horizontal smoothing of buoyancy and shear
 POWER_LAW               Power-law shape time-averaging barotropic filter
 PROFILE                 Time profiling activated
 K_GSCHEME               Third-order upstream advection of TKE fields
 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
 SNOWFALL                Surface snow accumulation on ice
 SOLAR_SOURCE            Solar Radiation Source Term
 SOLVE3D                 Solving 3D Primitive Equations
 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_QDRAG                Quadratic bottom stress
 UV_VIS2                 Harmonic mixing of momentum
 VAR_RHO_2D              Variable density barotropic mode
 VISC_GRID               Horizontal viscosity coefficient scaled by grid size
When we prepare the grid file:
- I noticed that the bathymetry values are positive, and the zice has negative values. Is this correct?
About the outputs:
- The data on ice-shelf regions I need to consider the s_rho starts under the ice (zice) or the surface (zeta).

Also, I am trying to understand how the melt rate from the ice shelf (?) and sea ice works. Also, the sea ice outputs units (time-averaged melt rates - in m3 s-1) seem a little confusing. The highest values in my monthly outputs appear in lower latitudes(wai,wio), where we, apparently, don't have sea ice:
Attachments
image_wro.png
image_aice.png

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

Re: Sea-Ice Output

#2 Unread post by kate »

You are using an unsupported code! I never ran with the ICESHELF, though there are those who did. Yes, the outputs will go from sea bottom to under the ice.

There are diagnostics which should be zero where there's no ice that aren't zero. I wouldn't swear by diagnostics ever being right. I abandoned this code
to move to something supported by a more careful scientist than I am (SIS2).

Post Reply