Changeset 617 for trunk/NEMO/TOP_SRC/initrc.F90

Timestamp:

2007-02-21T14:03:10+01:00 (17 years ago)

Author:

opalod

Message:

* empty log message *

File:

• trunk/NEMO/TOP_SRC/initrc.F90 (modified) (5 diffs)

trunk/NEMO/TOP_SRC/initrc.F90

@@ -32 +32 @@
    !! * Accessibility
    PUBLIC ini_trc
-   
+
+    !! * Substitutions
+#  include "domzgr_substitute.h90"
+  
 CONTAINS
    
@@ -60 +63 @@
       !!  OPA.9, 03-2005
       !!---------------------------------------------------------------------
+      INTEGER :: ji, jj, jk, jn    !: dummy loop indices
 
       !! 0.b PRINT the number of tracer
@@ -83 +87 @@
       CALL trc_ini
 
-      ! 4. restart from a FILE (nutrst)
-      ! ----------------------
+ 
+      ! 4. total volume of the ocean
+      !-----------------------------
+
+      areatot = 0.
+      DO jk = 1, jpk
+         DO jj = 1, jpj
+            DO ji = 1, jpi
+               areatot = areatot + tmask(ji,jj,jk) * tmask_i(ji,jj)  &
+#if defined key_off_degrad
+                  &                * facvol(ji,jj,jk)    &
+#endif
+                  &                * e1t(ji,jj) * e2t(ji,jj) * fse3t(ji,jj,jk) 
+            END DO
+         END DO
+      END DO
+      IF( lk_mpp ) THEN 
+         CALL mpp_sum(areatot)     ! sum over the global domain  
+      END IF
+
+      IF(lwp) WRITE(numout,*) ' '
+      IF (lwp) WRITE(numout,*) 'Total volume of ocean =',areatot
+      IF(lwp) WRITE(numout,*) ' '
+
+      ! 5. Initialization of tracers
+      ! -----------------------------
 
       IF( lrsttr ) THEN
 
-
-         CALL trc_rst
+         ! 5.1 restart from a file
+         !------------------------
+         CALL trc_rst_read
 
       ELSE
 
-         ! start from anything ELSE
-
+         !  5.2 analytical formulation or global data
+         !-------------------------------------
          CALL trc_dtr
 
       ENDIF
 
-      ! 5. Print control
+
+      ! 6. Computation integral of all tracers
+      !------------------
+
+      trai = 0.
+      DO jn = 1, jptra
+         DO jk = 1, jpk
+            DO jj = 1, jpj
+               DO ji = 1, jpi
+                  trai = trai + trn(ji,jj,jk,jn) * tmask(ji,jj,jk) * tmask_i(ji,jj)    &
+#if defined key_off_degrad
+                     &              * facvol(ji,jj,jk)   &
+#endif
+
+                     &              * e1t(ji,jj) * e2t(ji,jj) * fse3t(ji,jj,jk)
+               END DO
+            END DO
+         END DO
+      ENDDO
+      
+      IF( lk_mpp ) THEN 
+         CALL mpp_sum(trai)         ! sum over the global domain  
+      END IF
+
+      IF(lwp) WRITE(numout,*) ' '     
+      IF(lwp) WRITE(numout,*) 'Integral of all tracers over the full domain at initial time =',trai
+      IF(lwp) WRITE(numout,*) ' '
+
+      ! 6. Print control
       !------------------
 
@@ -106 +163 @@
    END SUBROUTINE ini_trc
 
+
 #else
    !!======================================================================
@@ -111 +169 @@
    !!======================================================================
 CONTAINS
-   SUBROUTINE ini_trc
-      
+   SUBROUTINE ini_trc      
    END SUBROUTINE ini_trc
 #endif