MODULE trcctl
   !!==========================================================================
   !!
   !!                       *** MODULE trcctl ***
   !!
   !! Only for passive tracer
   !! control the cpp options for the run and IF files are availables
   !! control also consistancy between options and namelist values
   !!  O.Aumont and A.El Moussaoui 03/05 F90 
   !!=========================================================================
   !!  TOP 1.0,  LOCEAN-IPSL (2005)
   !! $Header$
   !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt
   !!----------------------------------------------------------------------
      !!----------------------------------------------------------------------
      !! * Modules used
      !! ==============
      USE oce_trc
      USE trc
      USE sms
      USE trctrp_ctl

     IMPLICIT NONE
     PRIVATE

    !! * Accessibility
    PUBLIC trc_ctl

CONTAINS

#if defined key_passivetrc

SUBROUTINE trc_ctl
   !!===========================================================================================
   !!
   !!
   !!                       ROUTINE trcctl
   !!                     ******************
   !!
   !!      we use IF/ENDIF inside #IF defined option-cpp
   !!      FILE name must not exceed 21 characters
   !!      
   !!===========================================================================================
     
   !!----------------------------------------------------------------------
   !! local declarations
   !! ==================
   INTEGER  :: istop, jn

!!!---------------------------------------------------------------------
!!!  OPA.9    03/2005  
!!!---------------------------------------------------------------------

! 0. Parameter
! ------------
      istop = 0

! 1. LOGICAL UNIT initialization for specifi! files for passive tracer
! --------------------------------------------------------------------
!     nutwrs : OUTPUT for passive tracer restart UNIT (always used)
!     nutrst : restart FILE  INPUT  UNIT (lrsttr=.TRUE.)
!     nutini(jptra) : UNIT for initial FILE for tracer

      nutwrs = 72
      nutrst = 73

! 2. restart for passive tracer (input)
! -----------------------------

      IF(lwp) WRITE(numout,*) ' '
      IF(lwp) WRITE(numout,*) ' *** PASSIVE TRACER MODEL OPTIONS'
      IF(lwp) WRITE(numout,*) ' *** CONTROL'
      IF(lwp) WRITE(numout,*) ' '

      IF(lwp) WRITE(numout,*) ' '
      IF(lwp) WRITE(numout,*) ' *** restart option for passive tracer'
      IF(lwp) WRITE(numout,*) ' '

      IF(lrsttr) THEN
          IF(lwp) WRITE(numout,*) ' READ a restart FILE for passive tracer'
          IF(lwp) WRITE(numout,*) ' '
      ELSE
          IF(lwp) WRITE(numout,*) ' no restart FILE'
          IF(lwp) WRITE(numout,*) ' '

! 3. OPEN FILES for initial tracer value
! --------------------------------------
          DO jn=1,jptra

! OPEN input FILE only IF lutini(jn) is true
! ------------------------------------------
            IF (lutini(jn)) THEN  

! prepare input FILE name a
! -------------------------                        
                IF(lwp) WRITE(numout,*)                          &
                        ' READ an initial FILE :',               &
                        ' for passive tracer number :',jn        &
                        ,' traceur : ',ctrcnm(jn) 
                IF(lwp) WRITE(numout,*) ' '
            END IF 
          END DO   
      ENDIF

! 4. Don't USE non penetrative convective mixing option
!     it's not implemented for passive tracer
! -----------------------------------------------------

      IF( ln_zdfnpc) THEN
        IF(lwp) WRITE (numout,*) ' ===>>>> : w a r n i n g '
        IF(lwp) WRITE (numout,*) ' =======   ============= '
        IF(lwp) WRITE (numout,*) ' STOP, this sheme is not implemented'
        IF(lwp) WRITE (numout,*) ' in passive tracer model:'
        IF(lwp) WRITE (numout,*) ' non penetrative convect. mixing scheme'
        istop = istop + 1
      ENDIF

! 5. transport scheme option
! --------------------------

      IF(lwp) WRITE(numout,*) '  '
      CALL trc_trp_ctl
      IF(lwp) WRITE(numout,*) '  '

! 6. SMS model
! ---------------------------------------------

      IF(lwp) WRITE(numout,*) '  '
      IF(lwp) WRITE(numout,*) ' *** Source/Sink model option'
      IF(lwp) WRITE(numout,*) '  '

#if defined key_trc_npzd
      #include "trcctl.npzd.h90" 
#elif defined key_trc_lobster1
      #include "trcctl.lobster1.h90"
#elif defined key_trc_pisces
      #include "trcctl.pisces.h90"
#elif defined key_trc_hamocc3
#   if defined key_trc_p3zd
      #include "trcctl.p3zd.h90"
#   else
      #include "trcctl.hamocc3.h90"
#   endif
#else

      IF(lwp) WRITE (numout,*) ' No Source/Sink model '
      IF(lwp) WRITE (numout,*) ' '
#endif

! E r r o r  control
! ------------------

      IF ( istop > 0  ) THEN
          IF(lwp)WRITE(numout,*)
          IF(lwp)WRITE(numout,*) istop,' E R R O R found : we stop'
          IF(lwp)WRITE(numout,*) '**************************'
          IF(lwp)WRITE(numout,*)
          STOP 'trcctl'
      ENDIF

END SUBROUTINE trc_ctl 

#else

SUBROUTINE trc_ctl
  !!==============
         ! no passive tracers
  !!==============
END SUBROUTINE trc_ctl

#endif

END MODULE trcctl