Ignore:
Timestamp:
2019-03-20T18:53:44+01:00 (20 months ago)
Author:
aumont
Message:

bug fixes corresponding to tickets #2262

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/releases/release-4.0/src/TOP/PISCES/P4Z/p4zsed.F90

    r10425 r10780  
    5353      INTEGER  ::  ji, jj, jk, ikt 
    5454      REAL(wp) ::  zrivalk, zrivsil, zrivno3 
    55       REAL(wp) ::  zwflux, zfminus, zfplus 
    56       REAL(wp) ::  zlim, zfact, zfactcal 
     55      REAL(wp) ::  zwflux, zlim, zfact, zfactcal 
    5756      REAL(wp) ::  zo2, zno3, zflx, zpdenit, z1pdenit, zolimit 
    5857      REAL(wp) ::  zsiloss, zcaloss, zws3, zws4, zwsc, zdep 
     
    104103               zdep    = rfact2 / e3t_n(ji,jj,1) 
    105104               zwflux  = fmmflx(ji,jj) / 1000._wp 
    106                zfminus = MIN( 0._wp, -zwflux ) * trb(ji,jj,1,jpfer) * zdep 
    107                zfplus  = MAX( 0._wp, -zwflux ) * icefeinput * zdep 
    108                zironice(ji,jj) =  zfplus + zfminus 
     105               zironice(ji,jj) =  MAX( -0.99 * trb(ji,jj,1,jpfer), -zwflux * icefeinput * zdep ) 
    109106            END DO 
    110107         END DO 
Note: See TracChangeset for help on using the changeset viewer.