Changeset 1736 for trunk/NEMO/TOP_SRC/PISCES/p4zprod.F90

Timestamp:

2009-11-16T16:49:43+01:00 (14 years ago)

Author:

cetlod

Message:

suppress useless variables in TOP modules, see ticket:602

File:

• trunk/NEMO/TOP_SRC/PISCES/p4zprod.F90 (modified) (7 diffs)

trunk/NEMO/TOP_SRC/PISCES/p4zprod.F90

@@ -113 +113 @@
 
 # if defined key_off_degrad
-      prmax(:,:,:) = 0.6 / rjjss * tgfunc(:,:,:) * facvol(:,:,:)
+      prmax(:,:,:) = 0.6 / rday * tgfunc(:,:,:) * facvol(:,:,:)
 # else
-      prmax(:,:,:) = 0.6 / rjjss * tgfunc(:,:,:)
+      prmax(:,:,:) = 0.6 / rday * tgfunc(:,:,:)
 # endif
 
@@ -161 +161 @@
                    zpislopen = zpislopead(ji,jj,jk) * trn(ji,jj,jk,jpnch)                 &
                      &         / ( trn(ji,jj,jk,jpphy) * 12.                   + rtrn )   &
-                     &         / ( prmax(ji,jj,jk) * rjjss * xlimphy(ji,jj,jk) + rtrn )
+                     &         / ( prmax(ji,jj,jk) * rday * xlimphy(ji,jj,jk) + rtrn )
 
                    zpislope2n = zpislopead2(ji,jj,jk) * trn(ji,jj,jk,jpdch)                &
                      &          / ( trn(ji,jj,jk,jpdia) * 12.                   + rtrn )   &
-                     &          / ( prmax(ji,jj,jk) * rjjss * xlimdia(ji,jj,jk) + rtrn )
+                     &          / ( prmax(ji,jj,jk) * rday * xlimdia(ji,jj,jk) + rtrn )
 
                    ! Computation of production function
@@ -212 +212 @@
          DO ji = 1, jpi
             zmxltst = MAX( 0.e0, hmld(ji,jj) - heup(ji,jj) )
-            zmxlday = zmxltst**2 / rjjss
+            zmxlday = zmxltst**2 / rday
             zmixnano(ji,jj) = 1.- zmxlday / ( 1.+ zmxlday )
             zmixdiat(ji,jj) = 1.- zmxlday / ( 3.+ zmxlday )
@@ -247 +247 @@
                   zpislopen = zpislopead(ji,jj,jk)          &
                   &         * trn(ji,jj,jk,jpnch) / ( rtrn + trn(ji,jj,jk,jpphy) * 12.)         &
-                  &         / ( prmax(ji,jj,jk) * rjjss * zmax + rtrn )
+                  &         / ( prmax(ji,jj,jk) * rday * zmax + rtrn )
 
                   zprbiochl = prmax(ji,jj,jk) * (  1.- EXP( -zpislopen * zetot2 )  )
@@ -255 +255 @@
                   zpronew(ji,jj,jk) = zprorca(ji,jj,jk) * xnanono3(ji,jj,jk)    &
                   &             / ( xnanono3(ji,jj,jk) + xnanonh4(ji,jj,jk) + rtrn )
-                  zprod = rjjss * zprorca(ji,jj,jk) * zprbiochl * trn(ji,jj,jk,jpphy) *zmax
+                  zprod = rday * zprorca(ji,jj,jk) * zprbiochl * trn(ji,jj,jk,jpphy) *zmax
 
                   zprofen(ji,jj,jk) = (fecnm)**2 * zprod / chlcnm            &
@@ -279 +279 @@
                   zpislope2n = zpislopead2(ji,jj,jk) * trn(ji,jj,jk,jpdch)        &
                   &           / ( rtrn + trn(ji,jj,jk,jpdia) * 12.)        &
-                  &           / ( prmax(ji,jj,jk) * rjjss * zmax + rtrn )
+                  &           / ( prmax(ji,jj,jk) * rday * zmax + rtrn )
 
                   zprdiachl = prmax(ji,jj,jk) * (  1.- EXP( -zetot2 * zpislope2n )  )
@@ -288 +288 @@
                   &              / ( xdiatno3(ji,jj,jk) + xdiatnh4(ji,jj,jk) + rtrn )
 
-                  zprod = rjjss * zprorcad(ji,jj,jk) * zprdiachl * trn(ji,jj,jk,jpdia) * zmax
+                  zprod = rday * zprorcad(ji,jj,jk) * zprdiachl * trn(ji,jj,jk,jpdia) * zmax
 
                   zprofed(ji,jj,jk) = (fecdm)**2 * zprod / chlcdm                   &