Changeset 2819 for branches/2011/dev_r2787_LOCEAN3_TRA_TRP/NEMOGCM/NEMO/TOP_SRC/trcini.F90

Timestamp:

2011-08-09T10:29:53+02:00 (13 years ago)

Author:

cetlod

Message:

Improvment of branch dev_r2787_LOCEAN3_TRA_TRP

File:

• branches/2011/dev_r2787_LOCEAN3_TRA_TRP/NEMOGCM/NEMO/TOP_SRC/trcini.F90 (modified) (11 diffs)

branches/2011/dev_r2787_LOCEAN3_TRA_TRP/NEMOGCM/NEMO/TOP_SRC/trcini.F90

@@ -16 +16 @@
    !!   top_alloc :   allocate the TOP arrays
    !!----------------------------------------------------------------------
-   USE oce_trc
-   USE trc
-   USE trcrst
+   USE oce_trc         ! shared variables between ocean and passive tracers
+   USE trc             ! passive tracers common variables
+   USE trcrst          ! passive tracers restart
    USE trcnam          ! Namelist read
    USE trcini_cfc      ! CFC      initialisation
@@ -25 +25 @@
    USE trcini_c14b     ! C14 bomb initialisation
    USE trcini_my_trc   ! MY_TRC   initialisation
-   USE trcdta   
-   USE daymod
+   USE trcdta          ! initialisation form files
+   USE daymod          ! calendar manager
    USE zpshde          ! partial step: hor. derivative   (zps_hde routine)
    USE prtctl_trc      ! Print control passive tracers (prt_ctl_trc_init routine)
@@ -56 +56 @@
       !!                or read data or analytical formulation
       !!---------------------------------------------------------------------
-      INTEGER ::   jk, jn    ! dummy loop indices
+      INTEGER ::   jk, jn, jl    ! dummy loop indices
+      INTEGER ::   ierr          ! local integer
       CHARACTER (len=25) :: charout
+      REAL(wp), ALLOCATABLE, DIMENSION(:,:,:,:) ::  ztrcdta   ! 4D  workspace
       !!---------------------------------------------------------------------
 
@@ -65 +67 @@
 
       CALL top_alloc()              ! allocate TOP arrays
+
+      IF( ln_dm2dc .AND. ( lk_pisces .OR. lk_lobster ) )    &
+         &  CALL ctl_stop( ' The diurnal cycle is not compatible with PISCES or LOBSTER  ' )
+
+      IF( nn_cla == 1 )   &
+         &  CALL ctl_stop( ' Cross Land Advection not yet implemented with passive tracer ; nn_cla must be 0' )
 
       !                             ! masked grid volume
@@ -71 +79 @@
       END DO
 
-      !                             ! total volume of the ocean
-#if ! defined key_degrad
-      areatot = glob_sum( cvol(:,:,:) )
-#else
-      areatot = glob_sum( cvol(:,:,:) * facvol(:,:,:) )  ! degrad option: reduction by facvol
-#endif
+      !                           ! total volume of the ocean ( degrad option: reduction by facvol )
+      IF( .NOT.lk_degrad )  THEN  ;  areatot = glob_sum( cvol(:,:,:) )
+      ELSE                        ;  areatot = glob_sum( cvol(:,:,:) * facvol(:,:,:) )
+      ENDIF
 
       CALL trc_nam                  ! read passive tracers namelists
-
-      !                             ! restart for passive tracer (input)
-      IF( ln_rsttr ) THEN
-         IF(lwp) WRITE(numout,*) '       read a restart file for passive tracer : ', cn_trcrst_in
-         IF(lwp) WRITE(numout,*) ' '
-      ELSE
-         IF( lwp .AND. lk_dtatrc ) THEN
-            DO jn = 1, jptra
-               IF( lutini(jn) )  &                  ! open input FILE only IF lutini(jn) is true
-                  &  WRITE(numout,*) ' read an initial file for passive tracer number :', jn, ' traceur : ', ctrcnm(jn) 
-             END DO
-          ENDIF
-          IF( lwp ) WRITE(numout,*)
-      ENDIF
-
-      IF( ln_dm2dc .AND. ( lk_pisces .OR. lk_lobster ) )    &
-         &       CALL ctl_stop( ' The diurnal cycle is not compatible with PISCES or LOBSTER  ' )
-
-      IF( nn_cla == 1 )   &
-         &       CALL ctl_stop( ' Cross Land Advection not yet implemented with passive tracer ; nn_cla must be 0' )
 
       IF( lk_lobster ) THEN   ;   CALL trc_ini_lobster      ! LOBSTER bio-model
@@ -119 +105 @@
       ELSE                    ;   IF(lwp) WRITE(numout,*) '          MY_TRC not used'
       ENDIF
+
+      IF( ln_trcdta )             CALL trc_dta_init
 
       IF( ln_rsttr ) THEN
@@ -130 +118 @@
            CALL day_init               ! set calendar
         ENDIF
-#if defined key_dtatrc
-        CALL trc_dta( nit000 )      ! Initialization of tracer from a file that may also be used for damping
-        DO jn = 1, jptra
-           IF( lutini(jn) )   trn(:,:,:,jn) = trdta(:,:,:,jn) * tmask(:,:,:)   ! initialisation from file if required
-        END DO
-#endif
+        IF( ln_trcdta .AND. nb_trcdta > 0 ) THEN  ! Initialisation of tracer from a file that may also be used for damping
+            ALLOCATE( ztrcdta(jpi,jpj,jpk,nb_trcdta), STAT=ierr )
+            IF( ierr > 0 ) THEN
+               CALL ctl_stop( 'trc_ini: unable to allocate ztrcdta array' )   ;   RETURN
+            ENDIF
+            !
+            CALL trc_dta( nit000, ztrcdta )   ! read tracer data at nit000
+            !
+            DO jn = 1, jptra
+               IF( ln_trc_ini(jn) ) THEN      ! update passive tracers arrays with input data read from file
+                  jl = n_trc_index(jn) 
+                  trn(:,:,:,jn) = ztrcdta(:,:,:,jl) * tmask(:,:,:)  
+               ENDIF
+            ENDDO
+            DEALLOCATE( ztrcdta ) 
+        ENDIF
+        !
         trb(:,:,:,:) = trn(:,:,:,:)
         ! 
@@ -147 +146 @@
 
       !           
-      trai = 0._wp         ! Computation content of all tracers
+      trai(:) = 0._wp         ! Computation content of all tracers ( degrad option: reduction by facvol )
       DO jn = 1, jptra
-#if ! defined key_degrad
-         trai = trai + glob_sum( trn(:,:,:,jn) * cvol(:,:,:) )
-#else
-         trai = trai + glob_sum( trn(:,:,:,jn) * cvol(:,:,:) * facvol(:,:,:) ) ! degrad option: reduction by facvol
-#endif
+         IF( .NOT.lk_degrad )  THEN  ;  trai(jn) =  glob_sum( trn(:,:,:,jn) * cvol(:,:,:) ) 
+         ELSE                        ;  trai(jn) =  glob_sum( trn(:,:,:,jn) * cvol(:,:,:) * facvol(:,:,:) )
+         ENDIF
       END DO      
 
@@ -161 +158 @@
          WRITE(numout,*) '          *** Total number of passive tracer jptra = ', jptra
          WRITE(numout,*) '          *** Total volume of ocean                = ', areatot
-         WRITE(numout,*) '          *** Total inital content of all tracers  = ', trai
+         WRITE(numout,*) '          *** Total inital content of all tracers '
+         DO jn = 1, jptra
+            WRITE(numout,*) ' tracer nb : ', jn, '  name : ', ctrcnm(jn), ' initial content :', trai(jn)
+         ENDDO
          WRITE(numout,*)
       ENDIF
@@ -186 +186 @@
       USE trczdf        , ONLY:   trc_zdf_alloc
       USE trdmod_trc_oce, ONLY:   trd_mod_trc_oce_alloc
-#if ! defined key_iomput
-      USE trcdia        , ONLY:   trc_dia_alloc
-#endif
-#if defined key_trcdmp 
-      USE trcdmp        , ONLY:   trc_dmp_alloc
-#endif
-#if defined key_dtatrc
-      USE trcdta        , ONLY:   trc_dta_alloc
-#endif
-#if defined key_trdmld_trc   ||   defined key_esopa
+#if defined key_trdmld_trc 
       USE trdmld_trc    , ONLY:   trd_mld_trc_alloc
 #endif
@@ -207 +198 @@
       ierr = ierr + trc_zdf_alloc()
       ierr = ierr + trd_mod_trc_oce_alloc()
-#if ! defined key_iomput
-      ierr = ierr + trc_dia_alloc()
-#endif
-#if defined key_trcdmp 
-      ierr = ierr + trc_dmp_alloc()
-#endif
-#if defined key_dtatrc
-      ierr = ierr + trc_dta_alloc()
-#endif
-#if defined key_trdmld_trc   ||   defined key_esopa
+#if defined key_trdmld_trc 
       ierr = ierr + trd_mld_trc_alloc()
 #endif