Changeset 858 for branches/dev_001_GM/NEMO/TOP_SRC/PISCES/p4zsink_kriest.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_kriest.F90 (modified) (6 diffs)

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

@@ -20 +20 @@
    USE sms
    USE p4zsink2
+   USE prtctl_trc
 
    IMPLICIT NONE
@@ -45 +46 @@
       !! ** Method  : - ???
       !!---------------------------------------------------------------------
-      INTEGER  ::   ji, jj, jk
-      INTEGER  ::   iksed
-      REAL(wp) ::   zagg1, zagg2, zagg3, zagg4, zagg5, zaggsi, zaggsh
-      REAL(wp) ::   znum , zeps, zfm, zgm, zsm
-      REAL(wp) ::   zdiv , zdiv1, zdiv2, zdiv3, zdiv4, zdiv5
-      REAL(wp) ::   zval1, zval2, zval3, zval4
-      REAL(wp) ::   zstep
+      INTEGER  :: ji, jj, jk
+      INTEGER  :: iksed
+      REAL(wp) :: zagg1, zagg2, zagg3, zagg4, zagg5, zaggsi, zaggsh
+      REAL(wp) :: zagg , zaggdoc, znumdoc
+      REAL(wp) :: znum , zeps, zfm, zgm, zsm
+      REAL(wp) :: zdiv , zdiv1, zdiv2, zdiv3, zdiv4, zdiv5
+      REAL(wp) :: zval1, zval2, zval3, zval4
+      REAL(wp) :: zstep
 #if defined key_trc_dia3d
       REAL(wp) ::   zrfact2
@@ -59 +61 @@
       REAL(wp), DIMENSION(jpi,jpj,jpk) ::   sinkfer
       REAL(wp), DIMENSION(jpi,jpj,jpk) ::   sinkcal, sinksil
+      CHARACTER (len=25) :: charout
+
       !!---------------------------------------------------------------------
 
@@ -184 +188 @@
                      &    )
 
-                  zaggsh = ( zagg1 + zagg2 + zagg3 ) * rfact2 * zdiss(ji,jj,jk) / 1000.
+                  zaggsh = ( zagg1 + zagg2 + zagg3 ) * rfact2 * xdiss(ji,jj,jk) / 1000.
 
 !    Aggregation of small into large particles
@@ -213 +217 @@
                   zaggsi = ( zagg4 + zagg5 ) * zstep / 10.
 
-                  xagg(ji,jj,jk) = 0.5 * xkr_stick * ( zaggsh + zaggsi )
+                  zagg = 0.5 * xkr_stick * ( zaggsh + zaggsi )
 
 !     Aggregation of DOC to small particles
 !     --------------------------------------
 
-                  xaggdoc(ji,jj,jk) = (     0.4 * trn(ji,jj,jk,jpdoc)               &
-                     &                 + 1018.  * trn(ji,jj,jk,jppoc)  ) * zstep    &
-# if defined key_off_degrad
-                     &              * facvol(ji,jj,jk)                              &
-# endif
-                     &              * zdiss(ji,jj,jk) * trn(ji,jj,jk,jpdoc)
+                  zaggdoc = ( 0.4 * trn(ji,jj,jk,jpdoc)               &
+                     &        + 1018.  * trn(ji,jj,jk,jppoc)  ) * zstep    &
+# if defined key_off_degrad
+                     &        * facvol(ji,jj,jk)                              &
+# endif
+                     &        * xdiss(ji,jj,jk) * trn(ji,jj,jk,jpdoc)
+                      
+                  znumdoc = trn(ji,jj,jk,jpnum) / ( trn(ji,jj,jk,jppoc) + rtrn )
+                  tra(ji,jj,jk,jppoc) = tra(ji,jj,jk,jppoc) + zaggdoc
+                  tra(ji,jj,jk,jpnum) = tra(ji,jj,jk,jpnum) + zaggdoc * znumdoc - zagg
+                  tra(ji,jj,jk,jpdoc) = tra(ji,jj,jk,jpdoc) - zaggdoc
 
                ENDIF
@@ -246 +255 @@
 #    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_kriest