Changeset 617 for trunk/NEMO/TOP_SRC/SMS/p4zmeso.F

Timestamp:

2007-02-21T14:03:10+01:00 (17 years ago)

Author:

opalod

Message:

* empty log message *

File:

• trunk/NEMO/TOP_SRC/SMS/p4zmeso.F (modified) (4 diffs)

trunk/NEMO/TOP_SRC/SMS/p4zmeso.F

@@ -146 +146 @@
 C    ----------------------------------
 C
+#if ! defined key_trc_kriest
         grazffe(ji,jj,jk) = 5E3*zstep*wsbio4(ji,jj,jk)
      &    *tgfunc2(ji,jj,jk)*trn(ji,jj,jk,jpgoc)*trn(ji,jj,jk,jpmes)
@@ -154 +155 @@
         grazfff(ji,jj,jk) = grazffe(ji,jj,jk)
      &    *trn(ji,jj,jk,jpbfe)/(trn(ji,jj,jk,jpgoc)+rtrn)
+#else
+C KRIEST3
+        grazffe(ji,jj,jk) = 0.5*1.3E-2/5.5E-7*0.3*zstep*wsbio3(ji,jj,jk)
+     &    *tgfunc(ji,jj,jk)*trn(ji,jj,jk,jppoc)*trn(ji,jj,jk,jpmes)
+     &    /(trn(ji,jj,jk,jppoc)*1E7+0.1)
+#    if defined key_off_degrad
+     &    *facvol(ji,jj,jk)
+#    endif
+
+
+C        grazffe(ji,jj,jk) = 5E3*zstep*wsbio3(ji,jj,jk)
+C     &    *tgfunc2(ji,jj,jk)*trn(ji,jj,jk,jppoc)*trn(ji,jj,jk,jpmes)
+C#    if defined key_off_degrad
+C     &    *facvol(ji,jj,jk)
+C#    endif
+        grazfff(ji,jj,jk) = grazffe(ji,jj,jk)
+     &    *trn(ji,jj,jk,jpsfe)/(trn(ji,jj,jk,jppoc)+rtrn)
+#endif
+
+
 C
             END DO
@@ -169 +190 @@
      &    +grazn(ji,jj,jk)+grazpoc(ji,jj,jk)+grazffe(ji,jj,jk))
      &    *(1.-epsher2-unass2)
-
+#if ! defined key_trc_kriest
         grafer2(ji,jj,jk)=(grazf(ji,jj,jk)+graznf(ji,jj,jk)
      &    +grazz(ji,jj,jk)*ferat3+grazpof(ji,jj,jk)
@@ -182 +203 @@
      &    (trn(ji,jj,jk,jpbfe)/(trn(ji,jj,jk,jpgoc)+rtrn)
      &    -ferat3),0.))
-
+#else
+        grafer2(ji,jj,jk)=(grazf(ji,jj,jk)+graznf(ji,jj,jk)
+     &    +grazz(ji,jj,jk)*ferat3+grazpof(ji,jj,jk)
+     &    +grazfff(ji,jj,jk))*(1.-epsher2-unass2)
+     &    +epsher2*(grazd(ji,jj,jk)*max(
+     &    (trn(ji,jj,jk,jpdfe)/(trn(ji,jj,jk,jpdia)+rtrn)
+     &    -ferat3),0.)+grazn(ji,jj,jk)*max(
+     &    (trn(ji,jj,jk,jpnfe)/(trn(ji,jj,jk,jpphy)+rtrn)
+     &    -ferat3),0.)+grazpoc(ji,jj,jk)*max(
+     &    (trn(ji,jj,jk,jpsfe)/(trn(ji,jj,jk,jppoc)+rtrn)
+     &    -ferat3),0.)+grazffe(ji,jj,jk)*max(
+     &    (trn(ji,jj,jk,jpsfe)/(trn(ji,jj,jk,jppoc)+rtrn)
+     &    -ferat3),0.))
+#endif
         grapoc2(ji,jj,jk)=(grazd(ji,jj,jk)+grazz(ji,jj,jk)
      &    +grazn(ji,jj,jk)+grazpoc(ji,jj,jk)+grazffe(ji,jj,jk))*unass2