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 12065 for NEMO/branches/2019/dev_r10742_ENHANCE-12_SimonM-Tides/src/TOP/TRP/trcsink.F90 – NEMO

Ignore:
Timestamp:
2019-12-05T12:06:36+01:00 (4 years ago)
Author:
smueller
Message:

Synchronizing with /NEMO/trunk@12055 (ticket #2194)

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2019/dev_r10742_ENHANCE-12_SimonM-Tides/src/TOP/TRP/trcsink.F90

    r10376 r12065  
    8181            END DO 
    8282         END DO 
    83          iiter(ji,jj) = MIN( iiter(ji,jj), nitermax ) 
     83         iiter(:,:) = MIN( iiter(:,:), nitermax ) 
    8484      ENDIF 
    8585 
     
    8787         DO jj = 1, jpj 
    8888            DO ji = 1, jpi 
    89                IF( tmask(ji,jj,jk) == 1 ) THEN 
     89               IF( tmask(ji,jj,jk) == 1.0 ) THEN 
    9090                 zwsmax = 0.5 * e3t_n(ji,jj,jk) * rday / rsfact 
    9191                 zwsink(ji,jj,jk) = MIN( pwsink(ji,jj,jk), zwsmax * REAL( iiter(ji,jj), wp ) ) 
     92               ELSE 
     93                 ! provide a default value so there is no use of undefinite value in trc_sink2 for zwsink2 initialization 
     94                 zwsink(ji,jj,jk) = 0. 
    9295               ENDIF 
    9396            END DO 
     
    215218      REWIND( numnat_ref )              ! namtrc_rad in reference namelist  
    216219      READ  ( numnat_ref, namtrc_snk, IOSTAT = ios, ERR = 907) 
    217 907   IF( ios /= 0 )   CALL ctl_nam ( ios , 'namtrc_snk in reference namelist', lwp ) 
     220907   IF( ios /= 0 )   CALL ctl_nam ( ios , 'namtrc_snk in reference namelist' ) 
    218221      REWIND( numnat_cfg )              ! namtrc_rad in configuration namelist  
    219222      READ  ( numnat_cfg, namtrc_snk, IOSTAT = ios, ERR = 908 ) 
    220 908   IF( ios > 0 )   CALL ctl_nam ( ios , 'namtrc_snk in configuration namelist', lwp ) 
     223908   IF( ios > 0 )   CALL ctl_nam ( ios , 'namtrc_snk in configuration namelist' ) 
    221224      IF(lwm) WRITE( numont, namtrc_snk ) 
    222225 
Note: See TracChangeset for help on using the changeset viewer.