Changeset 1180 for trunk/NEMO/TOP_SRC/PISCES/p4zsed.F90

Timestamp:

2008-09-11T18:41:47+02:00 (16 years ago)

Author:

cetlod

Message:

update PISCES modules to couple with the sediment model, see ticket:249

File:

• trunk/NEMO/TOP_SRC/PISCES/p4zsed.F90 (modified) (8 diffs)

trunk/NEMO/TOP_SRC/PISCES/p4zsed.F90

@@ -63 +63 @@
    !!----------------------------------------------------------------------
    !! NEMO/TOP 2.0 , LOCEAN-IPSL (2007) 
-   !! $Id$ 
+   !! $Header:$ 
    !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt)
    !!----------------------------------------------------------------------
@@ -82 +82 @@
       INTEGER  ::   ji, jj, jk
       INTEGER  ::   ikt
+#if ! defined key_sed
       REAL(wp) ::   zsumsedsi, zsumsedpo4, zsumsedcal
+#endif
       REAL(wp) ::   zconctmp , zdenitot  , znitrpottot
       REAL(wp) ::   zlim, zconctmp2, zstep, zfact
@@ -137 +139 @@
          &       + zirondep(:,:,jk) + ironsed(:,:,jk) * rfact2
       END DO
+
+
+#if ! defined key_sed
 
       ! Initialisation of variables used to compute Sinking Speed
@@ -174 +179 @@
          CALL mpp_sum( zsumsedpo4 )   ! sums over the global domain
       ENDIF
+
+#endif
 
       ! Then this loss is scaled at each bottom grid cell for
@@ -192 +199 @@
 
             trn(ji,jj,ikt,jpdsi) = trn(ji,jj,ikt,jpdsi) - zconctmp
+
+#if ! defined key_sed
             trn(ji,jj,ikt,jpsil) = trn(ji,jj,ikt,jpsil) + zconctmp   &
             &      * ( 1.- ( sumdepsi + rivalkinput / raass / 6. ) / zsumsedsi )
+#endif
          END DO
       END DO
@@ -202 +212 @@
             zconctmp = trn(ji,jj,ikt,jpcal) * wscal(ji,jj,ikt) * zstep / fse3t(ji,jj,ikt)
             trn(ji,jj,ikt,jpcal) = trn(ji,jj,ikt,jpcal) - zconctmp
+
+#if ! defined key_sed
             trn(ji,jj,ikt,jptal) = trn(ji,jj,ikt,jptal) + zconctmp   &
                &   * ( 1.- ( rivalkinput / raass ) / zsumsedcal ) * 2.e0
             trn(ji,jj,ikt,jpdic) = trn(ji,jj,ikt,jpdic) + zconctmp   &
                &   * ( 1.- ( rivalkinput / raass ) / zsumsedcal )
+#endif
          END DO
       END DO
@@ -219 +232 @@
             trn(ji,jj,ikt,jpgoc) = trn(ji,jj,ikt,jpgoc) - zconctmp  * wsbio4(ji,jj,ikt) * zfact
             trn(ji,jj,ikt,jppoc) = trn(ji,jj,ikt,jppoc) - zconctmp2 * wsbio3(ji,jj,ikt) * zfact
+#if ! defined key_sed
             trn(ji,jj,ikt,jpdoc) = trn(ji,jj,ikt,jpdoc)    &
             &      + ( zconctmp  * wsbio4(ji,jj,ikt) + zconctmp2 * wsbio3(ji,jj,ikt) ) * zfact   &
             &      * ( 1.- rivpo4input / (raass * zsumsedpo4 ) )
+#endif
             trn(ji,jj,ikt,jpbfe) = trn(ji,jj,ikt,jpbfe) - trn(ji,jj,ikt,jpbfe) * wsbio4(ji,jj,ikt) * zfact
             trn(ji,jj,ikt,jpsfe) = trn(ji,jj,ikt,jpsfe) - trn(ji,jj,ikt,jpsfe) * wsbio3(ji,jj,ikt) * zfact
@@ -233 +248 @@
             trn(ji,jj,ikt,jppoc) = trn(ji,jj,ikt,jppoc)   &
             &      - zconctmp2 * wsbio3(ji,jj,ikt) * zfact
+#if ! defined key_sed
             trn(ji,jj,ikt,jpdoc) = trn(ji,jj,ikt,jpdoc)    &
             &      + ( zconctmp2 * wsbio3(ji,jj,ikt) )   &
             &      * zfact * ( 1.- rivpo4input / ( raass * zsumsedpo4 ) )
+#endif
             trn(ji,jj,ikt,jpsfe) = trn(ji,jj,ikt,jpsfe)   &
             &      - trn(ji,jj,ikt,jpsfe) * wsbio3(ji,jj,ikt) * zfact