yasmi_waq.f

Go to the documentation of this file.

!                       ***********************
                        SUBROUTINE yasmi_waq
!                       ***********************
     &  (yasmi)
!
!
!***********************************************************************
! TELEMAC2D   V7P0
!***********************************************************************
!
!brieF tells which tracers will have implicit source terms
!
!
!history  R.ATA
!+        12/02/2016
!+        V7P2
!+        CREATION
!
!         R. ATA
!+        07/07/2016
!+        V7P3
!+        ADAPTATION FOR THE NEW MANAGEMENT OF TRACERS
!
!history  S.E. BOURBAN (HRW)
!+        21/09/2017
!+        V7P3
!+        WAQPROCESS is now a prime number, so that multiple processes
!+        can be called by multiplication of the prime numbers.
!
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!| YASMI          |<--| LOGICS FOR IMPLICIT SOURCE TERMS
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!
      USE bief
      USE declarations_special
      USE declarations_waqtel
      USE interface_waqtel, ex_yasmi_waq => yasmi_waq
!
      IMPLICIT NONE
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
!
      LOGICAL          , INTENT(INOUT)::  YASMI(*)
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      INTEGER :: I
      LOGICAL :: FOUND
      found = .false.
!
!-----------------------------------------------------------------------
!
!     O2 MODULE
!
      IF( 2*int(waqprocess/2).EQ.waqprocess ) THEN
!
!       DISSOLVED O2
        yasmi(ind_o2) = yasmi(ind_o2) .OR. .false.
!       ORGANIC LOAD
        yasmi(ind_ol) = .true.
!       NH4 LOAD
        yasmi(ind_nh4) = .true.
!
        found = .true.
      ENDIF
!
!-----------------------------------------------------------------------
!
!     BIOMASS MODULE
!
      IF( 3*int(waqprocess/3).EQ.waqprocess ) THEN
!
!       PHYTO BIOMASS
        yasmi(ind_phy) = yasmi(ind_phy) .OR. .false.
!       DISSOLVED PO4
        yasmi(ind_po4) = yasmi(ind_po4) .OR. .false.
!       POR NON ASSIM
        yasmi(ind_por) = yasmi(ind_por) .OR. .false.
!       DISSOLVED NO3
        yasmi(ind_no3) = yasmi(ind_no3) .OR. .false.
!       NO3 NON ASSIM
        yasmi(ind_nor) = yasmi(ind_nor) .OR. .false.
!
        found = .true.
      ENDIF
!
!-----------------------------------------------------------------------
!
!     EUTRO MODULE
!
      IF( 5*int(waqprocess/5).EQ.waqprocess ) THEN
!
!       PHYTO BIOMASS
        yasmi(ind_phy) = yasmi(ind_phy) .OR. .false.
!       DISSOLVED PO4
        yasmi(ind_po4) = yasmi(ind_po4) .OR. .false.
!       POR NON ASSIM
        yasmi(ind_por) = yasmi(ind_por) .OR. .false.
!       DISSOLVED NO3
        yasmi(ind_no3) = yasmi(ind_no3) .OR. .false.
!       NOR NON ASSIM
        yasmi(ind_nor) = yasmi(ind_nor) .OR. .false.
!       CHARGE NH4
        yasmi(ind_nh4) = .true.
!       ORGANIC LOAD
        yasmi(ind_ol) = .true.
!       DISSOLVED O2
        yasmi(ind_o2) = yasmi(ind_o2) .OR. .false.
!
        found = .true.
      ENDIF
!
!-----------------------------------------------------------------------
!
!     MICROPOL MODULE
!
      IF( 7*int(waqprocess/7).EQ.waqprocess ) THEN
!
!       SUSPENDED LOAD
        yasmi(ind_ss) = yasmi(ind_ss) .OR. .false.
!       BED SEDIMENTS
        yasmi(ind_sf) = yasmi(ind_sf) .OR. .false.
!       MICRO POLLUTANT
        yasmi(ind_c) = .true.
!       ADSORBED SUSPENSION LOAD
        yasmi(ind_css) = .true.
!       ADSORBED BED SED
        yasmi(ind_csf) = .true.
!
        found = .true.
      ENDIF
!
!-----------------------------------------------------------------------
!
!     THERMIC MODULE
!
      IF( 11*int(waqprocess/11).EQ.waqprocess ) THEN
!
!       CONSIDERED IMPLICIT BECAUSE ATABOS AND BTABOS ARE IMPLICIT
        yasmi(ind_t) = .true.
!
        found = .true.
      ENDIF
!
!-----------------------------------------------------------------------
!
!     AED2 COUPLING
!
      IF( 13*int(waqprocess/13).EQ.waqprocess ) THEN
        found = .true.
      ENDIF
!
!-----------------------------------------------------------------------
!
!     DEGRADATION LAW
!
      IF( 17*int(waqprocess/17).EQ.waqprocess ) THEN
        found = .true.
!
        DO i = 1,nwaq_degra
          yasmi(rank_degra(i)) = .true.
        ENDDO
!
      ENDIF
!
!-----------------------------------------------------------------------
!
!     GHOST PROCESS IN WAITING FOR THE MERGE WITH ICE PROCESS
!
      IF( 19*int(waqprocess/19).EQ.waqprocess ) THEN
        found = .true.
      ENDIF
!
!-----------------------------------------------------------------------
!
!      UNKNOWN PROCESS
!
      IF( .NOT.found ) THEN
        WRITE(lu,20) waqprocess
20    FORMAT(1x,'YASMI_WAQ: UNKNOWN WAQ MODULE: ',i4)
        CALL plante(1)
        stop
      ENDIF
!
!-----------------------------------------------------------------------
!
      RETURN
      END