Changeset 502

Timestamp:

2006-09-12T13:15:31+02:00 (18 years ago)

Author:

opalod

Message:

nemo_v1_update_069:CE:add a correction factor for degradation

Location:

trunk/NEMO/TOP_SRC/SMS

Files:

• p4zflx.F (modified) (3 diffs)
• p4zsed.F (modified) (5 diffs)

trunk/NEMO/TOP_SRC/SMS/p4zflx.F

@@ -120 +120 @@
           kgco2(ji,jj) = kgco2(ji,jj)/(100.*3600.)
      &      *(1-freeze(ji,jj))*tmask(ji,jj,1)
+#    if defined key_off_degrad
+     &        *facvol(ji,jj,1)
+#    endif
 C
          END DO
@@ -137 +140 @@
           kgo2(ji,jj) = (0.3*ws*ws + 2.5*(0.5246+ttc*(0.016256+
      &      ttc*0.00049946)))*sqrt(660./schmitto2)
+
 C
 C CONVERT TO m/s AND APPLY SEA ICE COVER
@@ -142 +146 @@
 C
           kgo2(ji,jj) = kgo2(ji,jj)/(100.*3600.)
-     $      *(1-freeze(ji,jj))*tmask(ji,jj,1)
+     &      *(1-freeze(ji,jj))*tmask(ji,jj,1)
+#    if defined key_off_degrad
+     &        *facvol(ji,jj,1)
+#    endif
 C
          ENDDO

trunk/NEMO/TOP_SRC/SMS/p4zsed.F

@@ -53 +53 @@
       REAL xlim,xconctmp2,zstep,zfact
       REAL irondep(jpi,jpj,jpk),sidep(jpi,jpj)
+      REAL zvol
 CC
 C
@@ -112 +113 @@
         END DO
 C
+
+C
 C     Add the external input of iron which is 3D distributed
 C     (dust, river and sediment mobilization)
@@ -234 +237 @@
         xlim=(1.-xnanono3(ji,jj,jk)-xnanonh4(ji,jj,jk))
         if (xlim.le.0.2) xlim=0.01
-        nitrpot(ji,jj,jk)=max(0.,(prmax(ji,jj,jk)-2.15/rjjss))
+        nitrpot(ji,jj,jk)=max(0.,(0.6*tgfunc(ji,jj,jk)-2.15)/rjjss)
+#if defined key_off_degrad
+     &    *facvol(ji,jj,jk)
+#endif
      &    *xlim*rfact2*trn(ji,jj,jk,jpfer)/(conc3
      &    +trn(ji,jj,jk,jpfer))*(1.-exp(-etot(ji,jj,jk)/50.))
@@ -241 +247 @@
        END DO
 C
-        nitrpottot=0.
+      nitrpottot=0.
       DO jk=1,jpkm1
         DO jj=2,jpj-1
@@ -259 +265 @@
         DO jj=1,jpj
           DO ji=1,jpi
+#if ! defined key_cfg_1d && ( defined key_orca_r4 || defined key_orca_r2 || defined key_orca_r05 || defined key_orca_r025 )
+        zfact=nitrpot(ji,jj,jk)*denitot/nitrpottot
+#else
         zfact=nitrpot(ji,jj,jk)*1.E-7
+#endif
         trn(ji,jj,jk,jpnh4)=trn(ji,jj,jk,jpnh4)+zfact
         trn(ji,jj,jk,jpoxy)=trn(ji,jj,jk,jpoxy)+zfact*o2nit