Changeset 858 for branches/dev_001_GM/NEMO/TOP_SRC/PISCES/p4zsink.F90

Timestamp:

2008-03-13T15:17:04+01:00 (16 years ago)

Author:

cetlod

Message:

include the new version of PISCES model , see ticket:91

File:

• branches/dev_001_GM/NEMO/TOP_SRC/PISCES/p4zsink.F90 (modified) (6 diffs)

branches/dev_001_GM/NEMO/TOP_SRC/PISCES/p4zsink.F90

@@ -18 +18 @@
    USE sms
    USE p4zsink2        !
+   USE prtctl_trc
+
 
    IMPLICIT NONE
@@ -46 +48 @@
       INTEGER  ::   iksed
       REAL(wp) ::   zagg1, zagg2, zagg3, zagg4
+      REAL(wp) ::   zagg , zaggfe, zaggdoc, zaggdoc2
       REAL(wp) ::   zfact, zstep, zwsmax
 #if defined key_trc_dia3d
@@ -53 +56 @@
       REAL(wp), DIMENSION(jpi,jpj,jpk) ::   zsinkfer, zsinkfer2
       REAL(wp), DIMENSION(jpi,jpj,jpk) ::   zsinkcal, zsinksil
+      CHARACTER (len=25) :: charout
       !!---------------------------------------------------------------------
 
@@ -120 +124 @@
             DO ji = 1, jpi
 
-               zfact = zstep * zdiss(ji,jj,jk)
+               zfact = zstep * xdiss(ji,jj,jk)
 
 !    Part I : Coagulation dependent on turbulence
@@ -155 +159 @@
 # endif
 
-               xagg  (ji,jj,jk) = zagg1 + zagg2 + zagg3 + zagg4
-               xaggfe(ji,jj,jk) = xagg(ji,jj,jk) * trn(ji,jj,jk,jpsfe) / ( trn(ji,jj,jk,jppoc) + rtrn )
+               zagg   = zagg1 + zagg2 + zagg3 + zagg4
+               zaggfe = zagg * trn(ji,jj,jk,jpsfe) / ( trn(ji,jj,jk,jppoc) + rtrn )
 
 !     Aggregation of DOC to small particles
 !     --------------------------------------
 
-               xaggdoc(ji,jj,jk) = ( 80.* trn(ji,jj,jk,jpdoc) + 698. * trn(ji,jj,jk,jppoc) )       &
-# if defined key_off_degrad
-                  &              * facvol(ji,jj,jk)                           &
-# endif
-                  &              * zfact * trn(ji,jj,jk,jpdoc)
-
-               xaggdoc2(ji,jj,jk) = 1.05e4 * zfact * trn(ji,jj,jk,jpgoc)   &
-# if defined key_off_degrad
-                  &               * facvol(ji,jj,jk)                            &
+               zaggdoc = ( 80.* trn(ji,jj,jk,jpdoc) + 698. * trn(ji,jj,jk,jppoc) )       &
+# if defined key_off_degrad
+                  &      * facvol(ji,jj,jk)                           &
+# endif
+                  &      * zfact * trn(ji,jj,jk,jpdoc)
+
+               zaggdoc2 = 1.05e4 * zfact * trn(ji,jj,jk,jpgoc)   &
+# if defined key_off_degrad
+                  &        * facvol(ji,jj,jk)                            &
 # endif      
-                  &               * trn(ji,jj,jk,jpdoc)
+                  &        * trn(ji,jj,jk,jpdoc)
+
+               tra(ji,jj,jk,jppoc) = tra(ji,jj,jk,jppoc) - zagg + zaggdoc
+               tra(ji,jj,jk,jpgoc) = tra(ji,jj,jk,jpgoc) + zagg + zaggdoc2
+               tra(ji,jj,jk,jpsfe) = tra(ji,jj,jk,jpsfe) - zaggfe
+               tra(ji,jj,jk,jpbfe) = tra(ji,jj,jk,jpbfe) + zaggfe
+               tra(ji,jj,jk,jpdoc) = tra(ji,jj,jk,jpdoc) - zaggdoc - zaggdoc2
 
             END DO
@@ -187 +197 @@
 # endif
       !
+       IF(ln_ctl)   THEN  ! print mean trends (used for debugging)
+         WRITE(charout, FMT="('sink')")
+         CALL prt_ctl_trc_info(charout)
+         CALL prt_ctl_trc(tab4d=tra, mask=tmask, clinfo=ctrcnm)
+       ENDIF
+
    END SUBROUTINE p4z_sink