Changeset 350 for trunk/NEMO/TOP_SRC/trcrst.F90

Timestamp:

2005-12-08T11:42:00+01:00 (18 years ago)

Author:

opalod

Message:

nemo_v1_update_032:RB: update restart routines ; update control for cfc

File:

• trunk/NEMO/TOP_SRC/trcrst.F90 (modified) (3 diffs)

trunk/NEMO/TOP_SRC/trcrst.F90

@@ -30 +30 @@
    CHARACTER (len=48) ::   &
       trestart = 'initial.trc.nc'   ! restart file name
+
+   !! * Substitutions
+#  include "passivetrc_substitute.h90"
    
 CONTAINS
@@ -85 +88 @@
 
       REAL(wp), DIMENSION(3) :: zinfo
+
+#if defined key_trc_pisces && ( defined key_orca_r4 || defined key_orca_r2 || defined key_orca_r05 || defined key_orca_r025 )
+      REAL(wp) , DIMENSION(jpi,jpj,jpk) ::   zvolk      
+      REAL(wp) ::   zareatot, zpo4tot
+#endif
 
       !!---------------------------------------------------------------------
@@ -238 +246 @@
 
 #if defined key_trc_pisces 
+
+#if defined key_orca_r4 || defined key_orca_r2 || defined key_orca_r05 || defined key_orca_r025
+
+      zvolk(:,:,:) = 0.
+      zareatot     = 0.
+      DO jk = 1, jpkm1
+         DO jj = 2, jpjm1
+            DO ji = 2, jpim1
+               zvolk(ji,jj,jk) = e1t(ji,jj) * e2t(ji,jj) * fse3t(ji,jj,jk) * tmask(ji,jj,jk)
+               zareatot        = zareatot + zvolk(ji,jj,jk)
+            ENDDO
+         ENDDO
+      ENDDO
+
+
+      zpo4tot = 0.
+      DO jk = 1, jpkm1
+         DO jj = 2, jpjm1
+            DO ji = 2, jpim1
+               zpo4tot = zpo4tot + trn(ji,jj,jk,jptal) * zvolk(ji,jj,jk)
+            END DO
+         END DO
+      END DO
+
+      WRITE(0,*) 'TALK moyen ', zpo4tot/zareatot*1E6
+      zpo4tot = zpo4tot/zareatot*1E6
+      trn(:,:,:,jptal) = trn(:,:,:,jptal)*2391./zpo4tot
+
+
+      zpo4tot = 0.
+      DO jk = 1, jpkm1
+         DO jj = 2, jpjm1
+            DO ji = 2, jpim1
+               zpo4tot = zpo4tot + trn(ji,jj,jk,jppo4) * zvolk(ji,jj,jk)
+            END DO
+         END DO
+      END DO
+
+      WRITE(0,*) 'PO4 moyen ', zpo4tot/zareatot*1E6/122.
+      zpo4tot = zpo4tot/zareatot*1E6/122.
+      trn(:,:,:,jppo4) = trn(:,:,:,jppo4)*2.165/zpo4tot
+
+
+      zpo4tot = 0.
+      DO jk = 1, jpkm1
+         DO jj = 2, jpjm1
+            DO ji = 2, jpim1
+               zpo4tot = zpo4tot + trn(ji,jj,jk,jpno3) * zvolk(ji,jj,jk)
+            END DO
+         END DO
+      END DO
+
+      WRITE(0,*) 'NO3 moyen ', zpo4tot/zareatot*1E6/7.6
+      zpo4tot = zpo4tot/zareatot*1E6/7.6
+      trn(:,:,:,jpno3) = trn(:,:,:,jpno3)*30.9/zpo4tot
+
+      zpo4tot = 0.
+      DO jk = 1, jpkm1
+         DO jj = 2, jpjm1
+            DO ji = 2, jpim1
+               zpo4tot = zpo4tot + trn(ji,jj,jk,jpsil) * zvolk(ji,jj,jk)
+            END DO
+         END DO
+      END DO
+
+      WRITE(0,*) 'SiO3 moyen ', zpo4tot/zareatot*1E6
+      zpo4tot = zpo4tot/zareatot*1E6
+      trn(:,:,:,jpsil) = MIN( 400E-6,trn(:,:,:,jpsil)*91.51/zpo4tot) 
+
+#endif
       !!  Initialization of chemical variables of the carbon cycle
       !!  --------------------------------------------------------
-
       DO jk = 1,jpk
          DO jj = 1,jpj