Changeset 1011 for trunk/NEMO/TOP_SRC/trcini.F90

Timestamp:

2008-05-30T11:57:12+02:00 (16 years ago)

Author:

cetlod

Message:

re-organization of TOP initialization and outputs phases see ticket 171

File:

• trunk/NEMO/TOP_SRC/trcini.F90 (modified) (1 diff)

trunk/NEMO/TOP_SRC/trcini.F90

@@ -4 +4 @@
    !! TOP :   Manage the passive tracer initialization
    !!======================================================================
-   !! History :    -   !  2000-04  (O. Aumont, M.A. Foujols)  original code
-   !!             1.0  !  2005-03  (O. Aumont, A. El Moussaoui) F90
-   !!              -   !  2005-10  (C. Ethe) Module
-   !!             2.0  !  2007-12  (C. Ethe, G. Madec) re-writing
+   !! History :    -   !  1991-03  ()  original code
+   !!             1.0  !  2005-03 (O. Aumont, A. El Moussaoui) F90
+   !!              -   !  2005-10 (C. Ethe) print control
+   !!             2.0  !  2005-10 (C. Ethe, G. Madec) revised architecture
+   !!----------------------------------------------------------------------
+#if defined key_top
+   !!----------------------------------------------------------------------
+   !!   'key_top'                                                TOP models
+   !!----------------------------------------------------------------------
    !!----------------------------------------------------------------------
    !!   trc_ini :   Initialization for passive tracer
    !!----------------------------------------------------------------------
-   USE par_trc         ! TOP parameters
+   USE oce_trc
+   USE trc
+   USE trcrst
+   USE trcctl
+   USE trclec
    USE trcini_cfc      ! CFC     initialisation
    USE trcini_lobster  ! LOBSTER initialisation
    USE trcini_pisces   ! PISCES  initialisation
    USE trcini_my_trc   ! MY_TRC  initialisation
+   USE trcdta          ! 
    USE in_out_manager  ! I/O manager
-
+   USE prtctl_trc      ! Print control passive tracers (prt_ctl_trc_init routine)
+   USE lib_mpp         ! distributed memory computing library
+   
    IMPLICIT NONE
    PRIVATE
+   
+   PUBLIC   trc_ini   ! called by opa
 
-   PUBLIC   trc_ini   ! called by ???
-
+    !! * Substitutions
+#  include "domzgr_substitute.h90"
    !!----------------------------------------------------------------------
-   !! NEMO/TOP 2.0 , LOCEAN-IPSL (2007) 
-   !! $Header:$ 
+   !! NEMO/TOP 1.0 , LOCEAN-IPSL (2005) 
+   !! $Id$ 
    !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt)
    !!----------------------------------------------------------------------
+  
+CONTAINS
+   
+   SUBROUTINE trc_ini
+      !!---------------------------------------------------------------------
+      !!                     ***  ROUTINE trc_ini  ***
+      !!
+      !! ** Purpose :   Initialization of the passive tracer fields 
+      !!
+      !! ** Method  : - read namelist
+      !!              - control the consistancy 
+      !!              - compute specific initialisations
+      !!              - set initial tracer fields (either read restart 
+      !!                or read data or analytical formulation
+      !!---------------------------------------------------------------------
+      INTEGER ::   jk, jn    ! dummy loop indices
+      REAL(wp), DIMENSION(jpi,jpj,jpk) ::   zbt   ! workspace: masked grid volume
+      CHARACTER (len=25) :: charout
 
-CONTAINS
+      !!---------------------------------------------------------------------
 
-   SUBROUTINE trc_ini
-      !!-------------------------------------------------------------------
-      !!                    ***  ROUTINE trc_ini  ***
-      !!              
-      !! ** Purpose :   Initialization of passive tracer to zero
-      !!
-      !! ** Method  :   call the initialisation of each defined tracer
-      !!             model (LOBSTER, PISCES, CFC, MY_TRC)
-      !!-------------------------------------------------------------------
-      !
       IF(lwp) WRITE(numout,*)
-      IF(lwp) WRITE(numout,*) 'trc_ini : initial set up of the tracer models'
+      IF(lwp) WRITE(numout,*) 'trc_ini : initial set up of the passive tracers'
       IF(lwp) WRITE(numout,*) '~~~~~~~'
-      !
+
+      !                 ! masked grid volume
+      DO jk = 1, jpk
+         zbt(:,:,jk) = e1t(:,:) * e2t(:,:) * fse3t(:,:,jk) * tmask(:,:,jk) * tmask_i(:,:)
+      END DO
+#if defined key_off_degrad
+      DO jk = 1, jpk
+         zbt(:,:,:) = zbt(:,:,:) * facvol(:,:,jk)      ! degrad option: reduction by facvol
+      ENDDO
+#endif
+
+      !                 ! total volume of the ocean
+      areatot = SUM( zbt(:,:,:) )
+      IF( lk_mpp )   CALL mpp_sum( areatot )     ! sum over the global domain  
+
+      CALL trc_lec      ! READ passive tracers namelists
+
+      CALL trc_ctl      ! control consistency between parameters, cpp key and namelists
+
       IF( lk_lobster ) THEN   ;   CALL trc_ini_lobster      ! LOBSTER bio-model
       ELSE                    ;   IF(lwp) WRITE(numout,*) '          LOBSTER not used'
       ENDIF
-
+      
       IF( lk_pisces  ) THEN   ;   CALL trc_ini_pisces       ! PISCES  bio-model
       ELSE                    ;   IF(lwp) WRITE(numout,*) '          PISCES not used'
       ENDIF
-
+      
       IF( lk_cfc     ) THEN   ;   CALL trc_ini_cfc          ! CFC     tracers
       ELSE                    ;   IF(lwp) WRITE(numout,*) '          CFC not used'
       ENDIF
-
+      
       IF( lk_my_trc  ) THEN   ;   CALL trc_ini_my_trc       ! MY_TRC  tracers
       ELSE                    ;   IF(lwp) WRITE(numout,*) '          MY_TRC not used'
       ENDIF
+
+      IF( .NOT. lrsttr ) THEN 
+           
+# if defined key_dtatrc
+         ! Initialization of tracer from a file that may also be used for damping
+         CALL trc_dta( nittrc000 )
+         DO jn = 1, jptra
+            IF( lutini(jn) )   trn(:,:,:,jn) = trdta(:,:,:,jn) * tmask(:,:,:)   ! initialisation from file if required
+         END DO
+# endif
+         trb(:,:,:,:) = trn(:,:,:,:)
+
+      ELSE
+
+         CALL trc_rst_read      ! restart from a file
+
+      ENDIF
+
+      !                 ! Computation content of all tracers
+      trai = 0.e0
+      DO jn = 1, jptra
+         trai = trai + SUM( trn(:,:,:,jn) * zbt(:,:,:) )
+      END DO      
+      IF( lk_mpp )   CALL mpp_sum( trai )     ! sum over the global domain  
+
+
+      !                 ! control print
+      IF(lwp) WRITE(numout,*) '          *** Total number of passive tracer jptra = ', jptra
+      IF(lwp) WRITE(numout,*) '          *** Total volume of ocean                = ', areatot
+      IF(lwp) WRITE(numout,*) '          *** Total inital content of all tracers  = ', trai
+      IF(lwp) WRITE(numout,*)
+
+      IF( ln_ctl )   CALL prt_ctl_trc_init      ! control print
       !
+
+      IF(ln_ctl) THEN      ! print mean trends (used for debugging)
+         WRITE(charout, FMT="('ini ')")
+         CALL prt_ctl_trc_info( charout )
+         CALL prt_ctl_trc( tab4d=trn, mask=tmask, clinfo=ctrcnm )
+      ENDIF
+
    END SUBROUTINE trc_ini
+
+#else
+   !!----------------------------------------------------------------------
+   !!  Empty module :                                     No passive tracer
+   !!----------------------------------------------------------------------
+CONTAINS
+   SUBROUTINE trc_ini                      ! Dummy routine   
+   END SUBROUTINE trc_ini
+#endif
 
    !!======================================================================