Changeset 1735 for trunk/NEMO/TOP_SRC/PISCES/p4zsed.F90

Timestamp:

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

Author:

cetlod

Message:

suppress useless variables in TOP modules, see ticket:602

File:

• trunk/NEMO/TOP_SRC/PISCES/p4zsed.F90 (modified) (15 diffs)

trunk/NEMO/TOP_SRC/PISCES/p4zsed.F90

@@ -47 +47 @@
 
    !! * Module variables
-   INTEGER ::          &
-      numdust,          &  ! logical unit for surface fluxes data
-      nflx1 , nflx2,   &  !  first and second record used
+   INTEGER ::                   &
+     ryyss,                     &  !: number of seconds per year
+     rmtss                         !: number of seconds per month
+
+   INTEGER ::                   &
+      numdust,                  &  !: logical unit for surface fluxes data
+      nflx1 , nflx2,            &  !: first and second record used
       nflx11, nflx12      ! ???
    REAL(wp), DIMENSION(jpi,jpj,2) ::    &  !:
@@ -103 +107 @@
       IF( (jnt == 1) .and. ( ln_dustfer ) )  CALL p4z_sbc( kt )
 
-      zstep = rfact2 / rjjss      ! Time step duration for the biology
+      zstep = rfact2 / rday      ! Time step duration for the biology
 
       zirondep(:,:,:) = 0.e0          ! Initialisation of variables used to compute deposition
@@ -113 +117 @@
       DO jj = 1, jpj
          DO ji = 1, jpi
-            zirondep(ji,jj,1) = ( dustsolub * dust(ji,jj) / ( 55.85 * rmoss ) + 3.e-10 / raass )   &
+            zirondep(ji,jj,1) = ( dustsolub * dust(ji,jj) / ( 55.85 * rmtss ) + 3.e-10 / ryyss )   &
                &             * rfact2 / fse3t(ji,jj,1)
-            zsidep  (ji,jj)   = 8.8 * 0.075 * dust(ji,jj) * rfact2 / ( fse3t(ji,jj,1) * 28.1 * rmoss )
+            zsidep  (ji,jj)   = 8.8 * 0.075 * dust(ji,jj) * rfact2 / ( fse3t(ji,jj,1) * 28.1 * rmtss )
          END DO
       END DO
@@ -123 +127 @@
 
       DO jk = 2, jpkm1
-         zirondep(:,:,jk) = dust(:,:) / ( 10. * 55.85 * rmoss ) * rfact2 * 1.e-4
+         zirondep(:,:,jk) = dust(:,:) / ( 10. * 55.85 * rmtss ) * rfact2 * 1.e-4
       END DO
 
@@ -159 +163 @@
          DO ji = 1, jpi
             ikt = MAX( mbathy(ji,jj)-1, 1 )
-            zfact = e1t(ji,jj) * e2t(ji,jj) / rjjss * tmask_i(ji,jj)
+            zfact = e1t(ji,jj) * e2t(ji,jj) / rday * tmask_i(ji,jj)
 # if defined key_kriest
             zsumsedsi  = zsumsedsi  + zfact *  trn(ji,jj,ikt,jpdsi) * wscal (ji,jj,ikt)
@@ -199 +203 @@
 #if ! defined key_sed
             trn(ji,jj,ikt,jpsil) = trn(ji,jj,ikt,jpsil) + zconctmp   &
-            &      * ( 1.- ( sumdepsi + rivalkinput / raass / 6. ) / zsumsedsi )
+            &      * ( 1.- ( sumdepsi + rivalkinput / ryyss / 6. ) / zsumsedsi )
 #endif
          END DO
@@ -212 +216 @@
 #if ! defined key_sed
             trn(ji,jj,ikt,jptal) = trn(ji,jj,ikt,jptal) + zconctmp   &
-               &   * ( 1.- ( rivalkinput / raass ) / zsumsedcal ) * 2.e0
+               &   * ( 1.- ( rivalkinput / ryyss ) / zsumsedcal ) * 2.e0
             trn(ji,jj,ikt,jpdic) = trn(ji,jj,ikt,jpdic) + zconctmp   &
-               &   * ( 1.- ( rivalkinput / raass ) / zsumsedcal )
+               &   * ( 1.- ( rivalkinput / ryyss ) / zsumsedcal )
 #endif
          END DO
@@ -231 +235 @@
             trn(ji,jj,ikt,jpdoc) = trn(ji,jj,ikt,jpdoc)    &
             &      + ( zconctmp  * wsbio4(ji,jj,ikt) + zconctmp2 * wsbio3(ji,jj,ikt) ) * zfact   &
-            &      * ( 1.- rivpo4input / (raass * zsumsedpo4 ) )
+            &      * ( 1.- rivpo4input / (ryyss * zsumsedpo4 ) )
 #endif
             trn(ji,jj,ikt,jpbfe) = trn(ji,jj,ikt,jpbfe) - trn(ji,jj,ikt,jpbfe) * wsbio4(ji,jj,ikt) * zfact
@@ -246 +250 @@
             trn(ji,jj,ikt,jpdoc) = trn(ji,jj,ikt,jpdoc)    &
             &      + ( zconctmp2 * wsbio3(ji,jj,ikt) )   &
-            &      * zfact * ( 1.- rivpo4input / ( raass * zsumsedpo4 ) )
+            &      * zfact * ( 1.- rivpo4input / ( ryyss * zsumsedpo4 ) )
 #endif
             trn(ji,jj,ikt,jpsfe) = trn(ji,jj,ikt,jpsfe)   &
@@ -282 +286 @@
                zlim = ( 1.- xnanono3(ji,jj,jk) - xnanonh4(ji,jj,jk) )
                IF( zlim <= 0.2 )   zlim = 0.01
-               znitrpot(ji,jj,jk) = MAX( 0.e0, ( 0.6 * tgfunc(ji,jj,jk) - 2.15 ) / rjjss )   &
+               znitrpot(ji,jj,jk) = MAX( 0.e0, ( 0.6 * tgfunc(ji,jj,jk) - 2.15 ) / rday )   &
 # if defined key_off_degrad
                &                  * facvol(ji,jj,jk)   &
@@ -565 +569 @@
 
 
+      ! Number of seconds per year and per month
+      ryyss = nyear_len(1) * rday
+      rmtss = ryyss / raamo
+
       ! total atmospheric supply of Si
       ! ------------------------------
@@ -571 +579 @@
          DO jj = 2, jpjm1
             DO ji = fs_2, fs_jpim1
-               sumdepsi = sumdepsi + zdustmo(ji,jj,jm) / (12.*rmoss) * 8.8        &
+               sumdepsi = sumdepsi + zdustmo(ji,jj,jm) / (12.*rmtss) * 8.8        &
                   &     * 0.075/28.1 * e1t(ji,jj) * e2t(ji,jj) * tmask(ji,jj,1) * tmask_i(ji,jj)
             END DO
@@ -582 +590 @@
       DO jj = 1, jpj
          DO ji = 1, jpi
-            zcoef = raass * e1t(ji,jj) * e2t(ji,jj) * fse3t(ji,jj,1) * tmask(ji,jj,1) * tmask_i(ji,jj)
+            zcoef = ryyss * e1t(ji,jj) * e2t(ji,jj) * fse3t(ji,jj,1) * tmask(ji,jj,1) * tmask_i(ji,jj)
             cotdep(ji,jj) =  river(ji,jj)                  *1E9 / ( 12. * zcoef + rtrn )
             rivinp(ji,jj) = (river(ji,jj)+riverdoc(ji,jj)) *1E9 / ( 31.6* zcoef + rtrn )
-            nitdep(ji,jj) = 7.6 * ndepo(ji,jj)                  / ( 14E6*raass*fse3t(ji,jj,1) + rtrn )
+            nitdep(ji,jj) = 7.6 * ndepo(ji,jj)                  / ( 14E6*ryyss*fse3t(ji,jj,1) + rtrn )
          END DO
       END DO
@@ -596 +604 @@
       DO jj = 2 , jpjm1
          DO ji = fs_2, fs_jpim1
-            zcoef = cvol(ji,jj,1) * raass
+            zcoef = cvol(ji,jj,1) * ryyss
             rivpo4input = rivpo4input + rivinp(ji,jj) * zcoef
             rivalkinput = rivalkinput + cotdep(ji,jj) * zcoef
@@ -612 +620 @@
       ! -------------------------
       DO jk = 1, jpkm1
-         ironsed(:,:,jk) = sedfeinput * cmask(:,:,jk) / ( fse3t(:,:,jk) * rjjss )
+         ironsed(:,:,jk) = sedfeinput * cmask(:,:,jk) / ( fse3t(:,:,jk) * rday )
       END DO
       CALL lbc_lnk( ironsed , 'T', 1. )      ! Lateral boundary conditions on ( ironsed )   (sign unchanged)