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 1329 for trunk/NEMO/TOP_SRC/PISCES/p4zprod.F90 – NEMO

Ignore:
Timestamp:
2009-02-20T11:32:52+01:00 (15 years ago)
Author:
cetlod
Message:

update modules to take into account the mask land points in NetCDF outputs, see ticket:322

File:
1 edited

Legend:

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

    r1298 r1329  
    359359      ENDIF 
    360360 
    361 #if defined key_trc_dia3d 
     361#if defined key_trc_diaadd && defined key_trc_dia3d 
    362362      zrfact2 = 1.e3 * rfact2r 
    363363      !   Supplementary diagnostics 
    364       trc3d(:,:,:,jp_pcs0_3d + 3)  = etot(:,:,:) 
    365       trc3d(:,:,:,jp_pcs0_3d + 4)  = zprorca(:,:,:)  * zrfact2 
    366       trc3d(:,:,:,jp_pcs0_3d + 5)  = zprorcad(:,:,:) * zrfact2 
    367       trc3d(:,:,:,jp_pcs0_3d + 6)  = zpronew(:,:,:)  * zrfact2 
    368       trc3d(:,:,:,jp_pcs0_3d + 7)  = zpronewd(:,:,:) * zrfact2 
    369       trc3d(:,:,:,jp_pcs0_3d + 8)  = zprorcad(:,:,:) * zysopt(:,:,:) * zrfact2 
    370       trc3d(:,:,:,jp_pcs0_3d + 9) = zprofed(:,:,:) * zrfact2 
     364      DO jk = 1, jpk 
     365         DO jj = 1, jpj 
     366            DO ji = 1, jpi  
     367               trc3d(ji,jj,jk,jp_pcs0_3d + 3)  = etot(ji,jj,jk)               * tmask(ji,jj,jk) 
     368               trc3d(ji,jj,jk,jp_pcs0_3d + 4)  = zprorca (ji,jj,jk) * zrfact2 * tmask(ji,jj,jk) 
     369               trc3d(ji,jj,jk,jp_pcs0_3d + 5)  = zprorcad(ji,jj,jk) * zrfact2 * tmask(ji,jj,jk) 
     370               trc3d(ji,jj,jk,jp_pcs0_3d + 6)  = zpronew (ji,jj,jk) * zrfact2 * tmask(ji,jj,jk) 
     371               trc3d(ji,jj,jk,jp_pcs0_3d + 7)  = zpronewd(ji,jj,jk) * zrfact2 * tmask(ji,jj,jk) 
     372               trc3d(ji,jj,jk,jp_pcs0_3d + 8)  = zprorcad(ji,jj,jk) * zrfact2 * tmask(ji,jj,jk) & 
     373     &                                         * zysopt(ji,jj,jk) 
     374               trc3d(ji,jj,jk,jp_pcs0_3d + 9)  = zprofed (ji,jj,jk) * zrfact2 * tmask(ji,jj,jk) 
    371375#if ! defined key_kriest 
    372       trc3d(:,:,:,jp_pcs0_3d + 10) = zprofen(:,:,:) * zrfact2 
     376               trc3d(ji,jj,jk,jp_pcs0_3d + 10) = zprofen (ji,jj,jk) * zrfact2 * tmask(ji,jj,jk) 
    373377#endif 
     378            ENDDO 
     379         ENDDO 
     380      ENDDO 
    374381#endif 
    375382 
Note: See TracChangeset for help on using the changeset viewer.