Changeset 339 for trunk/NEMO/TOP_SRC/SMS/p4zmicro.F

Timestamp:

2005-11-14T13:30:28+01:00 (18 years ago)

Author:

opalod

Message:

nemo_v1_update_027 : CE + RB + CT : update of SMS routines

File:

• trunk/NEMO/TOP_SRC/SMS/p4zmicro.F (modified) (6 diffs)

trunk/NEMO/TOP_SRC/SMS/p4zmicro.F

@@ -1 @@
-CCC$Header$
-CCC  TOP 1.0 , LOCEAN-IPSL (2005)
-C This software is governed by CeCILL licence see modipsl/doc/NEMO_CeCILL.txt
-C ---------------------------------------------------------------------------
 CDIR$ LIST
       SUBROUTINE p4zmicro
@@ -47 +43 @@
       INTEGER ji, jj, jk
       REAL compadi,compadi2,compaz,compaph,compapoc
-      REAL graze,zdenom
-      REAL zinano,zidiat,zipoc
-
+      REAL graze,zdenom,zdenom2
+      REAL zfact,zstep,zinano,zidiat,zipoc
+C
+C    Time step duration for biology
+C    ------------------------------
+C
+        zstep=rfact2/rjjss
 C
 
@@ -57 +57 @@
 C
         compaz = max((trn(ji,jj,jk,jpzoo)-1.E-9),0.)
+        zfact=zstep*tgfunc(ji,jj,jk)*compaz
+#    if defined key_off_degrad
+     &    *facvol(ji,jj,jk)
+#    endif
 C
 C     Respiration rates of both zooplankton
 C     -------------------------------------
 C
-        respz(ji,jj,jk) = resrat*rfact2/rjjss
+        respz(ji,jj,jk) = resrat*zfact
+     &    *(1.+3.*nitrfac(ji,jj,jk))
      &    *trn(ji,jj,jk,jpzoo)/(xkmort+trn(ji,jj,jk,jpzoo))
-     &    *compaz*tmask(ji,jj,jk)
-#    if defined key_off_degrad
-     &    *facvol(ji,jj,jk)
-#    endif
 C
 C     Zooplankton mortality. A square function has been selected with
@@ -73 +74 @@
 C     ---------------------------------------------------------------
 C
-          tortz(ji,jj,jk) = mzrat*rfact2*1E6/rjjss*tgfunc(ji,jj,jk)
-     &      *compaz*trn(ji,jj,jk,jpzoo)*tmask(ji,jj,jk)
-#    if defined key_off_degrad
-     &      *facvol(ji,jj,jk)
-#    endif
+          tortz(ji,jj,jk) = mzrat*1E6*zfact*trn(ji,jj,jk,jpzoo)
+C
             END DO
           END DO
@@ -87 +85 @@
 C
         compadi = max((trn(ji,jj,jk,jpdia)-1E-8),0.)
-        compadi2=min(compadi,2.E-6)
-        compaph = max((trn(ji,jj,jk,jpphy)-1E-8),0.)
-        compapoc=max((trn(ji,jj,jk,jppoc)-1E-9),0.)
+        compadi2=min(compadi,5.E-7)
+        compaph = max((trn(ji,jj,jk,jpphy)-2E-7),0.)
+        compapoc=max((trn(ji,jj,jk,jppoc)-1E-8),0.)
 C
 C     Microzooplankton grazing
 C     ------------------------
 C
-          graze = grazrat*rfact2/rjjss*tmask(ji,jj,jk)
-     &      *tgfunc(ji,jj,jk)
+          zdenom2 = 1./(zprefp*compaph
+     &      +zprefc*compapoc+zprefd*compadi2+rtrn)
+
+          graze = grazrat*zstep*tgfunc(ji,jj,jk)
+     &      *trn(ji,jj,jk,jpzoo)
 #    if defined key_off_degrad
      &      *facvol(ji,jj,jk)
 #    endif
-          zdenom = 1./(xkgraz+zprefp*compaph
-     &      +zprefc*trn(ji,jj,jk,jppoc)+zprefd*compadi2)
 
-          zinano=zprefp*trn(ji,jj,jk,jpphy)/
-     &      (zprefp*trn(ji,jj,jk,jpphy)
-     &      +zprefc*trn(ji,jj,jk,jppoc)
-     &      +zprefd*trn(ji,jj,jk,jpdia)+rtrn)
+          zinano=zprefp*compaph*zdenom2
+          zipoc=zprefc*compapoc*zdenom2
+          zidiat=zprefd*compadi2*zdenom2
 
-          zipoc=zprefc*trn(ji,jj,jk,jppoc)/
-     &      (zprefp*trn(ji,jj,jk,jpphy)
-     &      +zprefc*trn(ji,jj,jk,jppoc)
-     &      +zprefd*trn(ji,jj,jk,jpdia)+rtrn)
-
-          zidiat=zprefd*trn(ji,jj,jk,jpdia)/
-     &      (zprefp*trn(ji,jj,jk,jpphy)
-     &      +zprefc*trn(ji,jj,jk,jppoc)
-     &      +zprefd*trn(ji,jj,jk,jpdia)+rtrn)
+          zdenom = 1./(xkgraz+zinano*compaph
+     &      +zipoc*compapoc+zidiat*compadi2)
 
           grazp(ji,jj,jk) = graze*zinano*compaph*zdenom
-     &      *trn(ji,jj,jk,jpzoo)
+          grazm(ji,jj,jk) = graze*zipoc*compapoc*zdenom
+          grazsd(ji,jj,jk) = graze*zidiat*compadi2*zdenom
 
           grazpf(ji,jj,jk) = grazp(ji,jj,jk)*
@@ -126 +118 @@
      &      trn(ji,jj,jk,jpnch)/(trn(ji,jj,jk,jpphy)+rtrn)
 
-          grazm(ji,jj,jk) = graze*zipoc*compapoc
-     &      *zdenom*trn(ji,jj,jk,jpzoo)
-
           grazmf(ji,jj,jk) = grazm(ji,jj,jk)
      &      *trn(ji,jj,jk,jpsfe)/(trn(ji,jj,jk,jppoc)+rtrn)
-
-          grazsd(ji,jj,jk) = graze*zidiat*compadi2*zdenom
-     &      *trn(ji,jj,jk,jpzoo)
 
           grazsf(ji,jj,jk) = grazsd(ji,jj,jk)