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 858 for branches/dev_001_GM/NEMO/TOP_SRC/PISCES/p4zsed.F90 – NEMO

Ignore:
Timestamp:
2008-03-13T15:17:04+01:00 (16 years ago)
Author:
cetlod
Message:

include the new version of PISCES model , see ticket:91

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/dev_001_GM/NEMO/TOP_SRC/PISCES/p4zsed.F90

    r775 r858  
    1717   USE sms 
    1818   USE lib_mpp 
     19   USE prtctl_trc 
     20 
    1921 
    2022   IMPLICIT NONE 
     
    5052      REAL(wp), DIMENSION(jpi,jpj)     ::   zsidep 
    5153      REAL(wp), DIMENSION(jpi,jpj,jpk) ::   znitrpot, zirondep 
     54      CHARACTER (len=25) :: charout 
    5255      !!--------------------------------------------------------------------- 
    5356 
     
    221224            DO ji = 2, jpim1 
    222225               zdenitot = zdenitot + denitr(ji,jj,jk) * rdenit * e1t(ji,jj) * e2t(ji,jj)   & 
    223                   &    *fse3t(ji,jj,jk) * tmask(ji,jj,jk) * tmask_i(ji,jj) * znegtr(ji,jj,jk) 
     226                  &    *fse3t(ji,jj,jk) * tmask(ji,jj,jk) * tmask_i(ji,jj) * xnegtr(ji,jj,jk) 
    224227            END DO 
    225228         END DO 
     
    231234      ! ------------------------------------------------------------- 
    232235 
     236!CDIR NOVERRCHK 
    233237      DO jk = 1, jpk 
    234          DO jj = 1, jpj 
     238!CDIR NOVERRCHK 
     239         DO jj = 1, jpj 
     240!CDIR NOVERRCHK 
    235241            DO ji = 1, jpi 
    236242               zlim = ( 1.- xnanono3(ji,jj,jk) - xnanonh4(ji,jj,jk) ) 
     
    265271            DO ji = 1, jpi 
    266272# if ! defined key_cfg_1d && ( defined key_orca_r4 || defined key_orca_r2 || defined key_orca_r05 || defined key_orca_r025 ) 
    267                zfact = znitrpot(ji,jj,jk) * zdenitot / znitrpottot 
     273!!               zfact = znitrpot(ji,jj,jk) * zdenitot / znitrpottot 
     274               zfact = znitrpot(ji,jj,jk) * 1.e-7 
    268275# else 
    269276               zfact = znitrpot(ji,jj,jk) * 1.e-7 
     
    285292# endif 
    286293      ! 
     294       IF(ln_ctl)   THEN  ! print mean trends (used for debugging) 
     295         WRITE(charout, FMT="('sed ')") 
     296         CALL prt_ctl_trc_info(charout) 
     297         CALL prt_ctl_trc(tab4d=trn, mask=tmask, clinfo=ctrcnm) 
     298       ENDIF 
     299 
    287300   END SUBROUTINE p4z_sed 
    288301 
Note: See TracChangeset for help on using the changeset viewer.