algae_growth.f

Go to the documentation of this file.

!                   ************************
                    SUBROUTINE algae_growth
!                   ************************
!
     &(alg,cmax,ray,gt,nutr,tox,npoin )
!
!***********************************************************************
! TELEMAC2D   V7P1
!***********************************************************************
!
!brief    COMPUTES THE GROWTH RATE OF ALGAE
!
!history  R. ATA (LNHE)
!+        02/09/2015
!+        V7P1
!+
!history  R. ATA (LNHE)
!+        02/04/2016
!+        V7P2
!+        VARIABLE TEMPERATURE
!
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!| ALG            |<--| ALGAE GROWTH
!| CMAX           |-->| MAXIMAL GROWTH RATE AT 20C
!| RAY            |-->| EFFECT OF SUNSHINE in [0,1]
!| GT             |-->| EFFECT OF OF TEMPERATURE ON ALGAE GROWTH
!|                |   | GT=T/20,  T: WATER TEMPERATURE
!| NUTR           |-->| EFFECT OF PHOSPHORED AND NITROGENIOUS NUTRIMENTS
!| TOX            |-->| COEFFICIENT OF WATER TOXICITY
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!
      USE bief
      USE interface_waqtel, ex_algae_growth => algae_growth
      USE declarations_special
!
      IMPLICIT NONE
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      INTEGER         , INTENT(IN   ) :: npoin
      DOUBLE PRECISION, INTENT(IN   ) :: cmax,ray(npoin),nutr(npoin),tox
      DOUBLE PRECISION, INTENT(INOUT) :: alg(npoin)
      TYPE(bief_obj)   , INTENT(IN  ) :: gt
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      INTEGER i
!
      DO i=1,npoin
        alg(i)=cmax*ray(i)*gt%R(i)*nutr(i)*tox
      ENDDO
!
!-----------------------------------------------------------------------
!
      RETURN
      END