Ignore:
Timestamp:
2015-05-29T18:13:37+02:00 (6 years ago)
Author:
mcastril
Message:

ticket #1523 Convergence Check

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2015/dev_r5302_CNRS18_HPC_scalability/NEMOGCM/NEMO/LIM_SRC_3/limhdf.F90

    r5123 r5320  
    121121         CALL lbc_lnk( zrlx, 'T', 1. )                   ! lateral boundary condition 
    122122         ! 
    123          zconv = 0._wp                                   ! convergence test 
    124          DO jj = 2, jpjm1 
    125             DO ji = fs_2, fs_jpim1 
    126                zconv = MAX( zconv, ABS( zrlx(ji,jj) - ptab(ji,jj) )  ) 
    127             END DO 
    128          END DO 
    129          IF( lk_mpp )   CALL mpp_max( zconv )            ! max over the global domain 
     123 
     124         !! The convergence check optimization can be activated in the namelist using  
     125         !! the parameters cc_opt to activate the optimization and cc_freq to set the frequency of  
     126         !! Convergence Checks. These parameters must go into namicedyn section. 
     127 
     128         IF ( cc_opt == .true. ) THEN 
     129            !Convergence test every cc_freq time-steps 
     130            IF ( MOD( iter-1 , cc_freq ) == 0 )  THEN 
     131               zconv = 0._wp                                   ! convergence test 
     132               DO jj = 2, jpjm1 
     133                  DO ji = fs_2, fs_jpim1 
     134                     zconv = MAX( zconv, ABS( zrlx(ji,jj) - ptab(ji,jj) )  ) 
     135                  END DO 
     136               END DO 
     137               IF( lk_mpp )   CALL mpp_max( zconv )            ! max over the global domain              
     138            ENDIF 
     139         ELSE 
     140             zconv = 0._wp                                   ! convergence test 
     141             DO jj = 2, jpjm1 
     142                DO ji = fs_2, fs_jpim1 
     143                   zconv = MAX( zconv, ABS( zrlx(ji,jj) - ptab(ji,jj) )  ) 
     144                END DO 
     145             END DO 
     146             IF( lk_mpp )   CALL mpp_max( zconv )            ! max over the global domain 
     147         END IF 
     148          
    130149         ! 
    131150         ptab(:,:) = zrlx(:,:) 
Note: See TracChangeset for help on using the changeset viewer.