New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 1255 for trunk/NEMO/TOP_SRC/PISCES/p4zprod.F90 – NEMO

Ignore:
Timestamp:
2009-01-13T11:20:17+01:00 (15 years ago)
Author:
cetlod
Message:

minor modifications in all top models, see ticket:299

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMO/TOP_SRC/PISCES/p4zprod.F90

    r1180 r1255  
    7373      REAL(wp) ::   zmxltst, zmxlday, zlim1 
    7474      REAL(wp) ::   zpislopen  , zpislope2n 
    75       REAL(wp) ::   zrum, zcodel, zargu 
     75      REAL(wp) ::   zrum, zcodel, zargu, zvol 
    7676      REAL(wp), DIMENSION(jpi,jpj)     ::   zmixnano   , zmixdiat, zstrn 
    7777      REAL(wp), DIMENSION(jpi,jpj,jpk) ::   zpislopead , zpislopead2 
     
    355355        DO jj = 1, jpj 
    356356          DO ji = 1, jpi 
    357              tpp  = tpp + ( zprorca(ji,jj,jk) + zprorcad(ji,jj,jk) )  & 
     357             zvol = cvol(ji,jj,jk) 
    358358#if defined key_off_degrad 
    359              &              * facvol(ji,jj,jk)   & 
     359             zvol = zvol * facvol(ji,jj,jk) 
    360360#endif 
    361              &              * e1t(ji,jj) * e2t(ji,jj) * fse3t(ji,jj,jk) * tmask_i(ji,jj) 
     361             tpp  = tpp + ( zprorca(ji,jj,jk) + zprorcad(ji,jj,jk) ) * zvol 
    362362          END DO 
    363363        END DO 
    364364      END DO 
     365 
     366      IF( lk_mpp ) CALL mpp_sum( tpp ) 
    365367 
    366368      IF( MOD( kt, nspyr ) == 0 ) THEN 
Note: See TracChangeset for help on using the changeset viewer.