algae_death.f

Go to the documentation of this file.

!                   ***********************
                    SUBROUTINE algae_death
!                   ***********************
!
     &(ald,mp,cmor,trr,tresp,gt,tox,npoin )
!
!***********************************************************************
! TELEMAC2D   V7P1
!***********************************************************************
!
!brief    COMPUTES THE DISAPPEARANCE RATE OF ALGAE
!
!history  R. ATA (LNHE)
!+        02/09/2015
!+        V7P1
!+
!
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!| ALD            |<--| ALGAE DEATH
!| CMOR           |-->| COEFFICIENTS OF ALGAE DEATH RATE AT 20C
!| RAY            |-->| EFFECT OF SUNSHINE in [0,1]
!| GT             |-->| EFFECT OF OF TEMPERATURE ON ALGAE GROWTH
!|                |   | GT=T/20,  T: WATER TEMPERATURE
!| MP             |<--| DISAPPEAANCE RATE OF ALGAL BIOMASS AT 20 DEG C
!| NPOIN          |-->| NUMBER OF NODES
!| TOX            |-->| COEFFICIENT OF WATER TOXICITY
!| TRR            |-->| TRACER (PHYTOPLANCTOPN BIOMASS)
!| TRESP          |-->| RESPIRATION RATE OF ALGAL BIOMASS
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!
      USE bief
      USE interface_waqtel, ex_algae_death=>algae_death
      USE declarations_special
!
      IMPLICIT NONE
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      INTEGER         , INTENT(IN   ) :: npoin
      DOUBLE PRECISION, INTENT(IN   ) :: cmor(2),trr(npoin),tox,tresp
      DOUBLE PRECISION, INTENT(INOUT) :: ald(npoin),mp(npoin)
      TYPE(bief_obj)  , INTENT(IN   ) :: gt
!     LOCAL VARIABLES
      INTEGER i
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      DO i=1,npoin
        mp(i) = cmor(1)+cmor(2)*trr(i)+tox
        ald(i)= gt%R(i)*(mp(i)+tresp)
      ENDDO
!
!-----------------------------------------------------------------------
!
      RETURN
      END