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
- 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:
