source: branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC/IDTRA/trcini_idtra.F90 @ 6046

MODULE trcini_idtra
   !!======================================================================
   !!                         ***  MODULE trcini_idtra  ***
   !! TOP :   initialisation of the IDEAL-TRACER tracers
   !!======================================================================
   !! History :   2.0  !  2007-12  (C. Ethe, G. Madec) from trcini.idtra.h90
   !!----------------------------------------------------------------------
#if defined key_idtra
   !!----------------------------------------------------------------------
   !!   'key_idtra'                                               IDEAL-TRACER tracers
   !!----------------------------------------------------------------------
   !! trc_ini_idtra      : IDEAL-TRACER model initialisation
   !!----------------------------------------------------------------------
   USE oce_trc         ! Ocean variables
   USE par_trc         ! TOP parameters
   USE trc             ! TOP variables
   USE trcsms_idtra    ! IDEAL-TRACER sms trends
   ! USE par_idtra       ! IDEAL-TRACER parameters
   ! USE in_out_manager  ! I/O manager
   ! USE lib_mpp
   ! USE iom

   IMPLICIT NONE
   PRIVATE

   PUBLIC   trc_ini_idtra   ! called by trcini.F90 module

   INTEGER  ::   inum                   ! unit number

   !!----------------------------------------------------------------------
   !! NEMO/TOP 2.0 , LOCEAN-IPSL (2007)
   !! $Id$
   !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt)
   !!----------------------------------------------------------------------
CONTAINS

   SUBROUTINE trc_ini_idtra
      !!----------------------------------------------------------------------
      !!                     ***  trc_ini_idtra  ***
      !!
      !! ** Purpose :   initialization for idtra model
      !!
      !! ** Method  : - Read the namidtra namelist and check the parameter values
      !!----------------------------------------------------------------------
      INTEGER  ::    jn, jl 
      !!----------------------------------------------------------------------


      ! Initialization of trn in case of  no restart
      !----------------------------------------------
      IF( .NOT. ln_rsttr ) THEN
         IF(lwp) THEN
            WRITE(numout,*)
            WRITE(numout,*) 'Initialization de id-tracers ; No restart : '
            WRITE(numout,*) '                             ; Init field equal 1 at surface - zero elsewhere'
         ENDIF
         DO jn = jp_idtra0, jp_idtra1
             trn(:,:,:,jn) = 0.e0
             trn(:,:,1,jn) = 1.0
           IF(lwp) WRITE(numout,*) 'Idealise Tracer initialisation -- jn = ',jn
         END DO
      ENDIF


      !   Ideal traceur do not need any atmospheric concentration.
      ! We consider that sucface concentration is equal to 1,
      ! that it is advectied within the water circulation,
      ! and that it is regularly degraded as if it was a radiactive tracer (tricium for example)
      ! But we can play with tha caracteristic time of
      !--------------------------------------------------------------------



      IF(lwp) WRITE(numout,*) 'Initialization of IDEAL-TRACER tracers done'
      IF(lwp) WRITE(numout,*) ' '

   END SUBROUTINE trc_ini_idtra

#else
   !!----------------------------------------------------------------------
   !!   Dummy module                                         No IDEAL-TRACER tracers
   !!----------------------------------------------------------------------
CONTAINS
   SUBROUTINE trc_ini_idtra             ! Empty routine


   END SUBROUTINE trc_ini_idtra
#endif

   !!======================================================================
END MODULE trcini_idtra