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/p4zlys.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/p4zlys.F90

    r1152 r1329  
    6666      IF( kt == nittrc000  )   CALL p4z_lys_init      ! Initialization (first time-step only) 
    6767 
    68  
     68      zco3(:,:,:) = 0. 
    6969 
    7070      !     ------------------------------------------- 
     
    150150      END DO 
    151151 
    152 # if defined key_trc_dia3d 
    153          trc3d(:,:,:,jp_pcs0_3d    ) = hi(:,:,:) 
    154          trc3d(:,:,:,jp_pcs0_3d + 1) = zco3(:,:,:) 
    155          trc3d(:,:,:,jp_pcs0_3d + 2) = aksp(:,:,:) / calcon 
     152# if defined key_trc_diaadd &&  defined key_trc_dia3d 
     153      DO jk = 1, jpk 
     154         DO jj = 1, jpj 
     155            DO ji = 1, jpi 
     156               trc3d(ji,jj,jk,jp_pcs0_3d    ) = hi  (ji,jj,jk)          * tmask(ji,jj,jk) 
     157               trc3d(ji,jj,jk,jp_pcs0_3d + 1) = zco3(ji,jj,jk)          * tmask(ji,jj,jk) 
     158               trc3d(ji,jj,jk,jp_pcs0_3d + 2) = aksp(ji,jj,jk) / calcon * tmask(ji,jj,jk) 
     159            ENDDO 
     160         ENDDO 
     161      ENDDO 
    156162# endif 
    157163      ! 
Note: See TracChangeset for help on using the changeset viewer.