Ignore:
Timestamp:
2018-12-01T12:33:31+01:00 (22 months ago)
Author:
smasson
Message:

dev_r10164_HPC09_ESIWACE_PREP_MERGE: merge with dev_r9866_HPC_03_globcom, see #2133

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2018/dev_r10164_HPC09_ESIWACE_PREP_MERGE/src/ICE/icedyn_rhg_evp.F90

    r10345 r10365  
    143143      REAL(wp), DIMENSION(jpi,jpj) ::   zds                             ! shear 
    144144      REAL(wp), DIMENSION(jpi,jpj) ::   zs1, zs2, zs12                  ! stress tensor components 
    145       REAL(wp), DIMENSION(jpi,jpj) ::   zu_ice, zv_ice, zresr           ! check convergence 
     145!!$      REAL(wp), DIMENSION(jpi,jpj) ::   zu_ice, zv_ice, zresr           ! check convergence 
    146146      REAL(wp), DIMENSION(jpi,jpj) ::   zssh_lead_m                     ! array used for the calculation of ice surface slope: 
    147147      !                                                                 !    ocean surface (ssh_m) if ice is not embedded 
     
    324324      DO jter = 1 , nn_nevp                           !    loop over jter    ! 
    325325         !                                            !----------------------!         
    326          IF(ln_ctl) THEN   ! Convergence test 
    327             DO jj = 1, jpjm1 
    328                zu_ice(:,jj) = u_ice(:,jj) ! velocity at previous time step 
    329                zv_ice(:,jj) = v_ice(:,jj) 
    330             END DO 
    331          ENDIF 
     326!!$         IF(ln_ctl) THEN   ! Convergence test 
     327!!$            DO jj = 1, jpjm1 
     328!!$               zu_ice(:,jj) = u_ice(:,jj) ! velocity at previous time step 
     329!!$               zv_ice(:,jj) = v_ice(:,jj) 
     330!!$            END DO 
     331!!$         ENDIF 
    332332 
    333333         ! --- divergence, tension & shear (Appendix B of Hunke & Dukowicz, 2002) --- ! 
     
    636636            ! 
    637637         ENDIF 
    638           
    639          IF(ln_ctl) THEN   ! Convergence test 
    640             DO jj = 2 , jpjm1 
    641                zresr(:,jj) = MAX( ABS( u_ice(:,jj) - zu_ice(:,jj) ), ABS( v_ice(:,jj) - zv_ice(:,jj) ) ) 
    642             END DO 
    643             zresm = MAXVAL( zresr( 1:jpi, 2:jpjm1 ) ) 
    644             IF( lk_mpp )   CALL mpp_max( 'icedyn_rhg_evp', zresm )   ! max over the global domain 
    645          ENDIF 
     638 
     639!!$         IF(ln_ctl) THEN   ! Convergence test 
     640!!$            DO jj = 2 , jpjm1 
     641!!$               zresr(:,jj) = MAX( ABS( u_ice(:,jj) - zu_ice(:,jj) ), ABS( v_ice(:,jj) - zv_ice(:,jj) ) ) 
     642!!$            END DO 
     643!!$            zresm = MAXVAL( zresr( 1:jpi, 2:jpjm1 ) ) 
     644!!$            IF( lk_mpp )   CALL mpp_max( 'icedyn_rhg_evp', zresm )   ! max over the global domain 
     645!!$         ENDIF 
    646646         ! 
    647647         !                                                ! ==================== ! 
Note: See TracChangeset for help on using the changeset viewer.