nemuro.in

From WikiROMS
Jump to navigationJump to search
NEMURO Model Input Script - nemuro.in

The nemuro.in file sets the parameters for the NEMURO model. The name of this file is set by the BPARNAM keyword in the ocean.in file. A default nemuro.in standard input ASCII file can be found in the User/External subdirectory of the ROMS source code. In order to include the NEMURO model in ROMS you must set BPARNAM correctly and activate the NEMURO CPP option.


NEMURO equations and representative parameters may be found in:


Note Notice: Detailed information about ROMS input script file syntax can be found here.



Biological Model Parameters

  • This switch to control the computation of NEMURO within nested and/or multiple connected grids. By default this switch is set to TRUE in mod_scalars.F for all grids. Ngrids values are expected. The user has the option, for example, to compute the biology in just one of the nested grids. If so, this switch needs to be consistent with the dimension parameter NBT in mod_param.F. In order to make the model more efficient in memory usage, NBT(:) should be zero in such grids.
Lbiology == T
  • Maximum number of iterations to achieve convergence of the nonlinear solution.
BioIter == 1


Light Parameters

Input parameter units are specified within brackets ([ ]). Unless specified, [1:Ngrids] values expected for each parameter.

  • Light attenuation due to seawater [1/m].
AttSW == 0.04d0
  • Fraction of shortwave radiation that is photosynthetically active, [nondimensional].
PARfrac == 0.43d0
  • Phytoplankton photochemical reaction coefficient, initial slope of the P-I curve [1/(W/m2) 1/day].
AlphaPS == 0.01d0  ! small biomass
AlphaPL == 0.01d0  ! large biomass
  • Phytoplankton photoinhibition coefficient, [1/(W/m2) 1/day].
BetaPS == 4.5d-4  ! small biomass
BetaPL == 4.5d-4  ! large biomass


Phytoplankton Parameters

Input parameter units are specified within brackets ([ ]). Unless specified, [1:Ngrids] values expected for each parameter.

  • Phytoplankton maximum photosynthetic rate at 0 Celsius [1/day].
VmaxS == 0.4d0  ! small biomass
VmaxL == 0.8d0  ! large biomass
  • Phytoplankton half saturation constant for Nitrate [millimole_N/m3].
KNO3S == 1.0d0  ! small biomass
KNO3L == 3.0d0  ! large biomass
  • Phytoplankton half saturation constant for Ammonium [millimole_N/m3].
KNH4S == 0.1d0  ! small biomass
KNH4L == 0.3d0  ! large biomass
  • Phytoplankton half saturation constant for Silicate [millimole_Si/m3].
KSiL == 6.0d0  ! large biomass
  • Phytoplankton Ammonium inhibition coefficient [m3/millimole_N].
PusaiS == 1.5d0  ! small biomass
PusaiL == 1.5d0  ! large biomass
  • Phytoplankton temperature coefficient for photosynthetic rate [1/Celsius].
KGppS == 6.93d-2  ! small biomass
KGppL == 6.93d-2  ! large biomass
  • Phytoplankton respiration rate at 0 Celsius [1/day].
ResPS0 == 0.03d0  ! small biomass
ResPL0 == 0.03d0  ! large biomass
  • Phytoplankton temperature coefficient for respiration [1/Celsius].
KResPS == 0.0519d0  ! small biomass
KResPL == 0.0519d0  ! large biomass
  • Phytoplankton ratio of extracellular excretion to photosynthesis [nondimensional].
GammaS == 0.135d0  ! small biomass
GammaL == 0.135d0  ! large biomass
  • Phytoplankton mortality rate at 0 Celsius [m3/millimole_N 1/day].
MorPS0 == 58.5d-3  ! small biomass
MorPL0 == 29.0d-3  ! large biomass
  • Phytoplankton temperature coefficient for mortality [1/Celsius].
KMorPS == 6.93d-2  ! small biomass
KMorPL == 6.93d-2  ! large biomass


Zooplankton Parameters

Input parameter units are specified within brackets ([ ]). Unless specified, [1:Ngrids] values expected for each parameter.

  • Zooplankton maximum grazing rate at 0 Celsius [1/day].
GRmaxSps == 0.40d0  ! small Zoo on small phy
GRmaxLps == 0.10d0  ! large Zoo on small Phy
GRmaxLpl == 0.40d0  ! large Zoo on large Phy
GRmaxLzs == 0.40d0  ! large Zoo on small Zoo
GRmaxPpl == 0.20d0  ! predator Zoo on large Phy
GRmaxPzs == 0.20d0  ! predator Zoo on small Zoo
GRmaxPzl == 0.20d0  ! predator Zoo on large Zoo
  • Zooplankton temperature coefficient for grazing [1/Celsius].
KGraS == 6.93d-2  ! small biomass
KGraL == 6.93d-2  ! large biomass
KGraP == 6.93d-2  ! predator biomass

Zooplankton Ivlev constant [m3/millimole_N].

LamS == 1.4d0  ! small biomass
LamL == 1.4d0  ! large biomass
LamP == 1.4d0  ! predator biomass
  • Zooplankton half-saturation coefficient (squared) for ingestion used only when the Holling-type grazing formulation is activated [millimole_N/m3]2.
KPS2ZS == 0.16d0  ! small Zoo on small Phy
KPS2ZL == 0.16d0  ! large Zoo on small Phy
KPL2ZL == 0.16d0  ! large Zoo on large Phy
KZS2ZL == 0.16d0  ! large Zoo on small Zoo
KPL2ZP == 0.16d0  ! predator Zoo on large Phy
KZS2ZP == 0.16d0  ! predator Zoo on small Zoo
KZL2ZP == 0.16d0  ! predator Zoo on large Zoo
  • Zooplankton threshold value for grazing [millimole_N/m3].
PS2ZSstar == 4.3d-2  ! small Zoo on small Phy
PS2ZLstar == 4.0d-2  ! large Zoo on small Phy
PL2ZLstar == 4.0d-2  ! large Zoo on large Phy
ZS2ZLstar == 4.0d-2  ! large Zoo on small Zoo
PL2ZPstar == 4.0d-2  ! predator Zoo on large Phy
ZS2ZPstar == 4.0d-2  ! predator Zoo on small Zoo
ZL2ZPstar == 4.0d-2  ! predator Zoo on large Zoo
  • Zooplankton grazing inhibition coefficient [m3/millimole_N].
PusaiPL == 4.605d0  ! predator Zoo on large Phy
PusaiZS == 3.010d0  ! predator Zoo on small Zoo
  • Zooplankton mortality rate at 0 Celsius [m3/millimole_N 1/day].
MorZS0 == 58.5d-3  ! small biomass
MorZL0 == 58.5d-3  ! large biomass
MorZP0 == 58.5d-3  ! predator biomass
  • Zooplankton temperature coefficient for mortality [1/Celsius].
KMorZS == 0.0693d0  ! small biomass
KMorZL == 0.0693d0  ! large biomass
KMorZP == 0.0693d0  ! predator biomass
  • Zooplankton assimilation efficiency [nondimemsional].
AlphaZS == 0.70d0  ! small biomass
AlphaZL == 0.70d0  ! large biomass
AlphaZP == 0.70d0  ! predator biomass
  • Zooplankton growth efficiency [nondimensional].
BetaZS == 0.30d0  ! small biomass
BetaZL == 0.30d0  ! large biomass
BetaZP == 0.30d0  ! predator biomass

Nutrient Parameters

Input parameter units are specified within brackets ([ ]). Unless specified, [1:Ngrids] values expected for each parameter.

  • Decomposition rates at 0 Celsius [1/day].
Nit0 == 0.03d0  ! NH4 nitrification
VP2N0 == 0.10d0  ! PON to NH4
VP2D0 == 0.10d0  ! PON to DON
VD2N0 == 0.20d0  ! DON to NH4
VO2S0 == 0.10d0  ! Opal to Silicate
  • Temperature coefficients for decomposition [1/Celsius]
KNit == 6.93d-2  ! NH4 nitrification
KP2D == 6.93d-2  ! PON to DON
KP2N == 6.93d-2  ! PON to NH4
KD2N == 6.93d-2  ! DON to NH4
KO2S == 6.93d-2  ! Opal to Silicate
  • Si:N ratio [millimole_Si/millimole_N].
RSiN == 2.0d0
  • Settling (sinking) velocities [m/day].
setVPON == 40.0d0  ! PON
setVOpal == 40.0d0  ! Opal

Other Parameters

Unless specified, [1:NBT,Ngrids] values expected for each parameter.

  • Lateral, constant, harmonic/biharmonic horizontal diffusion of biological tracer for nonlinear model and adjoint-based algorithms: [1:NBT,Ngrids] values expected.
TNU2 == 11*0.0d0  ! m2/s
TNU4 == 11*0.0d0  ! m4/s

ad_TNU2 == 11*0.0d0  ! m2/s
ad_TNU4 == 11*0.0d0  ! m4/s
  • Vertical mixing coefficients for biological tracers: [1:NBT,Ngrids] values expected.
AKT_BAK == 11*1.0d-6  ! m2/s

ad_AKT_fac == 11*1.0d0  ! nondimensional
  • Nudging/relaxation time scale. Inverse scale will be computed internally: [1:NBT,Ngrids] values expected.
TNUDG == 11*0.0d0  ! days
  • Logical switches to specify which variables to consider on tracers point Sources/Sinks (like river runoff): [1:NBT,Ngrids] values expected.
LtracerSrc == 11*F
  • Logical switches to write out biological data into output History/Average NetCDF files. Currently, NEMURO has NBT=11 ecosystem compartments.
Hout(idTvar) == 11*T  ! biological tracers
Hout(idTsur) == 11*F  ! biological tracers surface flux

! idTvar(idbio( 1))=iPhyL Large Phytoplankton biomass
! idTvar(idbio( 2))=iPhyS Small Phytoplankton biomass
! idTvar(idbio( 3))=iZooL Large Zooplankton biomass
! idTvar(idbio( 4))=iZooS Small Zooplankton biomass
! idTvar(idbio( 5))=iZooP Predator Zooplankton biomass
! idTvar(idbio( 6))=iNO3_ Nitrate concentration
! idTvar(idbio( 7))=iNH4_ Ammonium concentration
! idTvar(idbio( 8))=iPON_ Particulate Organic Nitrogen
! idTvar(idbio( 9))=iDON_ Dissolved Organic Nitrogen
! idTvar(idbio(10))=iSiOH Silicate concentration
! idTvar(idbio(11))=iopal Particulate organic silica
As an example, if you wanted to exclude large zooplankton biomass and silicate concentration from writing to the history and averages files you would enter the following:
Hout(idTvar) == 2*T F 6*T F T
Retrieved from "https://www.myroms.org/wiki/index.php?title=nemuro.in&oldid=3774"