MODULE sms_pisces !!---------------------------------------------------------------------- !! *** sms_pisces.F90 *** !! TOP : PISCES Source Minus Sink variables are declared and allocated !!---------------------------------------------------------------------- !! History : 1.0 ! 2000-02 (O. Aumont) original code !! 3.2 ! 2009-04 (C. Ethe & NEMO team) style !!---------------------------------------------------------------------- USE par_oce USE par_trc IMPLICIT NONE PUBLIC INTEGER :: numnatp_ref = -1 !: Logical units for namelist pisces INTEGER :: numnatp_cfg = -1 !: Logical units for namelist pisces INTEGER :: numonp = -1 !: Logical unit for namelist pisces output !!* Model used LOGICAL :: ln_p2z !: Flag to use LOBSTER model LOGICAL :: ln_p4z !: Flag to use PISCES model LOGICAL :: ln_p5z !: Flag to use PISCES quota model LOGICAL :: ln_ligand !: Flag to enable organic ligands LOGICAL :: ln_sediment !: Flag to enable sediment module !!* Time variables INTEGER :: nrdttrc !: ??? REAL(wp) :: rfact , rfactr !: time step duration (in seconds) REAL(wp) :: rfact2, rfact2r !: time step duration (in seconds) when timesplitting is activated for PISCES REAL(wp) :: xstep !: Time step duration for biology REAL(wp) :: ryyss !: number of seconds per year REAL(wp) :: r1_ryyss !: inverse number of seconds per year !!* Biological parameters REAL(wp) :: rno3 !: C/N stoichiometric ratio REAL(wp) :: o2ut !: O2/N stoichiometric ratio for ammonification REAL(wp) :: po4r !: C/P stoichiometric ratio REAL(wp) :: rdenit !: C/N ratio for denitrification REAL(wp) :: rdenita !: C/N ratio for denitrification REAL(wp) :: o2nit !: O2/N ratio for nitrification REAL(wp) :: wsbio, wsbio2 !: Sinking speeds of particles REAL(wp) :: wsbio2max !: Maximum sinking speed of the largest particles REAL(wp) :: wsbio2scale !: Length scale for the variations of wsbio2 REAL(wp) :: xkmort !: Mortality half-saturation constant REAL(wp) :: ferat3 !: Fe/C in heterotrophic organisms REAL(wp) :: ldocp !: Ligand production ratio during PP REAL(wp) :: ldocz !: Ligand production ratio by grazing REAL(wp) :: lthet !: Uptake of ligand by phytoplankton REAL(wp) :: no3rat3 !: C/N ratio of zooplankton REAL(wp) :: po4rat3 !: C/P ratio of zooplankton !!* diagnostic parameters REAL(wp) :: tpp !: total primary production REAL(wp) :: t_oce_co2_exp !: total carbon export REAL(wp) :: t_oce_co2_flx !: Total ocean carbon flux REAL(wp) :: t_oce_co2_flx_cum !: Cumulative Total ocean carbon flux REAL(wp) :: t_atm_co2_flx !: global mean of atmospheric pco2 !!* restoring LOGICAL :: ln_pisdmp !: restoring or not of nutrients to a mean value INTEGER :: nn_pisdmp !: frequency of relaxation or not of nutrients to a mean value !!* Mass conservation LOGICAL :: ln_check_mass !: Flag to check mass conservation LOGICAL, PUBLIC :: ln_ironice !: boolean for Fe input from sea ice !!* Biological fluxes for light : variables shared by pisces & lobster INTEGER , ALLOCATABLE, SAVE, DIMENSION(:,:) :: neln !: number of T-levels + 1 in the euphotic layer REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) :: heup !: euphotic layer depth REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: etot !: par (photosynthetic available radiation) REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: etot_ndcy !: PAR over 24h in case of diurnal cycle REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: enano, ediat !: PAR for phyto, nano and diat REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: enanom, ediatm !: mean PAR for phyto, nano and diat REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: epico !: PAR for pico REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: epicom !: mean PAR for pico REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: emoy !: averaged PAR in the mixed layer REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) :: heup_01 !: Absolute euphotic layer depth REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) :: xksi !: Half-saturation con,stant for diatoms !!* Biological fluxes for primary production REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) :: xksimax !: Maximum half-saturation constant over the year (Si) REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: biron !: bioavailable fraction of iron REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: plig !: proportion of iron organically complexed !!* Sinking speed REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: wsbio3 !: POC sinking speed REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: wsbio4 !: GOC sinking speed !!* SMS for the organic matter REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: xfracal !: Fraction of nanophytoplankton that are calcifying organisms REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: nitrfac !: OMZ REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: nitrfac2 !: N depleted indice REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: orem !: oxic remineralisation REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: xdiss !: Shear rate REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: prodcal !: Calcite production REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: prodpoc !: POC production REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: conspoc !: POC consumption REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: prodgoc !: GOC production REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: consgoc !: GOC consumption REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: blim !: bacterial production factor REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: sizen !: size of nanophyto REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: sizep !: size of picophyto REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: sized !: size of diatoms REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: sizena !: size of nanophytoplankton, after REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: sizepa !: size of picophyto, after REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: sizeda !: size of diatomss, after !!* Variable for chemistry of the CO2 cycle REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: ak13 !: Carbonate chemistry constant REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: ak23 !: Carbonate chemistry constant REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: aksp !: Solubility product of CaCO3 REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: hi !: Proton concentration REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: excess !: CO3 saturation REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: aphscale !: !!* Temperature dependancy of SMS terms REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: tgfunc !: Temp. dependancy of various biological rates REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: tgfunc2 !: Temp. dependancy of mesozooplankton rates REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: tgfunc3 !: Temp. dependancy of picophytoplankton rates !!---------------------------------------------------------------------- !! NEMO/TOP 4.0 , NEMO Consortium (2018) !! $Id$ !! Software governed by the CeCILL license (see ./LICENSE) !!---------------------------------------------------------------------- CONTAINS INTEGER FUNCTION sms_pisces_alloc() !!---------------------------------------------------------------------- !! *** ROUTINE sms_pisces_alloc *** !!---------------------------------------------------------------------- USE lib_mpp , ONLY: ctl_stop INTEGER :: ierr(11) ! Local variables !!---------------------------------------------------------------------- ierr(:) = 0 ALLOCATE( etot(jpi,jpj,jpk), neln(jpi,jpj), heup(jpi,jpj), & & heup_01(jpi,jpj) , xksi(jpi,jpj) , STAT=ierr(1) ) ! IF( ln_p4z .OR. ln_p5z ) THEN !* Optics ALLOCATE( enano(jpi,jpj,jpk) , ediat(jpi,jpj,jpk) , & & enanom(jpi,jpj,jpk) , ediatm(jpi,jpj,jpk), & & etot_ndcy(jpi,jpj,jpk), emoy(jpi,jpj,jpk) , STAT=ierr(2) ) !* Biological SMS ALLOCATE( xksimax(jpi,jpj) , biron(jpi,jpj,jpk) , STAT=ierr(3) ) ! Biological SMS ALLOCATE( xfracal (jpi,jpj,jpk), orem(jpi,jpj,jpk) , & & nitrfac(jpi,jpj,jpk) , nitrfac2(jpi,jpj,jpk), & & prodcal(jpi,jpj,jpk) , xdiss (jpi,jpj,jpk), & & prodpoc(jpi,jpj,jpk) , conspoc(jpi,jpj,jpk) , & & prodgoc(jpi,jpj,jpk) , consgoc(jpi,jpj,jpk) , & & blim (jpi,jpj,jpk) , STAT=ierr(4) ) !* Carbonate chemistry ALLOCATE( ak13 (jpi,jpj,jpk) , & & ak23(jpi,jpj,jpk) , aksp (jpi,jpj,jpk) , & & hi (jpi,jpj,jpk) , excess(jpi,jpj,jpk) , & & aphscale(jpi,jpj,jpk), STAT=ierr(5) ) ! !* Temperature dependency of SMS terms ALLOCATE( tgfunc (jpi,jpj,jpk) , tgfunc2(jpi,jpj,jpk), & & tgfunc3(jpi,jpj,jpk) , STAT=ierr(6) ) ! !* Sinking speed ALLOCATE( wsbio3 (jpi,jpj,jpk) , wsbio4 (jpi,jpj,jpk), STAT=ierr(7) ) !* Size of phytoplankton cells ALLOCATE( sizen (jpi,jpj,jpk), sized (jpi,jpj,jpk), & & sizena(jpi,jpj,jpk), sizeda(jpi,jpj,jpk), STAT=ierr(8) ) ! !* Prognostic ligand IF( ln_ligand ) THEN ALLOCATE( plig(jpi,jpj,jpk) , STAT=ierr(9) ) ENDIF ENDIF ! IF( ln_p5z ) THEN ! PISCES-QUOTA specific part ALLOCATE( epico(jpi,jpj,jpk) , epicom(jpi,jpj,jpk) , STAT=ierr(10) ) !* Size of phytoplankton cells ALLOCATE( sizep(jpi,jpj,jpk), sizepa(jpi,jpj,jpk), STAT=ierr(11) ) ENDIF ! sms_pisces_alloc = MAXVAL( ierr ) ! IF( sms_pisces_alloc /= 0 ) CALL ctl_stop( 'STOP', 'sms_pisces_alloc: failed to allocate arrays' ) ! END FUNCTION sms_pisces_alloc !!====================================================================== END MODULE sms_pisces