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 15523 – NEMO

Changeset 15523


Ignore:
Timestamp:
2021-11-19T09:14:41+01:00 (6 months ago)
Author:
vancop
Message:

Some cleaning of lbclnks in VP routine

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2020/SI3-03_VP_rheology/src/ICE/icedyn_rhg_vp.F90

    r15516 r15523  
    88   !!            4.0  !  2020-09  (M. Vancoppenolle) Adaptation to SI3 
    99   !! 
     10   !! CLEAN-UP 2 — remove iomputs 
    1011   !!---------------------------------------------------------------------- 
    1112#if defined key_si3 
     
    174175      REAL(wp), DIMENSION(jpi,jpj) ::   zu_c, zv_c                      ! "current" ice velocity (m/s), average of previous two OL iterates 
    175176      REAL(wp), DIMENSION(jpi,jpj) ::   zu_b, zv_b                      ! velocity at previous sub-iterate 
    176       REAL(wp), DIMENSION(jpi,jpj) ::   zu_b_outer, zv_b_outer          ! velocity at previous outer-iterate 
    177177      REAL(wp), DIMENSION(jpi,jpj) ::   zuerr, zverr                    ! absolute U/Vvelocity difference between current and previous sub-iterates 
    178178      ! 
     
    217217                         
    218218      !!---------------------------------------------------------------------------------------------------------------------- 
    219       ! DEBUG put all forcing terms to zero 
    220       ! air-ice drag 
    221       !utau_ice(:,:) = 0._wp 
    222       !vtau_ice(:,:) = 0._wp 
    223  
    224       ! coriolis 
    225       !ff_t(:,:) = 0._wp 
    226  
    227       ! ice-ocean drag 
    228       !rn_cio = 0._wp 
    229  
    230219      ! sea surface height 
    231220      !    embedded sea ice: compute representative ice top surface 
    232221      !    non-embedded sea ice: use ocean surface for slope calculation 
    233       zsshdyn(:,:) = ice_var_sshdyn( ssh_m, snwice_mass, snwice_mass_b) ! DEBUG 
    234       !zsshdyn(:,:) = 0._wp  
    235  
    236       ! END DEBUG 
     222      zsshdyn(:,:) = ice_var_sshdyn( ssh_m, snwice_mass, snwice_mass_b) 
     223 
    237224      !!---------------------------------------------------------------------------------------------------------------------- 
    238225 
     
    315302         END DO 
    316303      END DO 
     304 
     305      !CALL lbc_lnk( 'icedyn_rhg_vp', zfmask, 'F', 1._wp ) ! OK to remove 
    317306 
    318307      ! Lateral boundary conditions on velocity (modify zfmask) 
     
    397386            ELSE                                                      ;   zmsk01y(ji,jj) = 1._wp   ;   ENDIF               
    398387 
    399 ! MV TEST DEBUG 
    400 !           IF ( ( zmt(ji,jj)   <= zmmin .OR. zmt(ji+1,jj)  <= zmmin )     .AND.  & 
    401 !              & ( at_i(ji,jj)  <= zamin .OR. at_i(ji+1,jj) <= zamin ) )    THEN   ;   zmsk01x(ji,jj) = 0._wp 
    402 !           ELSE                                                                   ;   zmsk01x(ji,jj) = 1._wp   ;   ENDIF 
    403  
    404 !           IF ( ( zmt(ji,jj)   <= zmmin .OR. zmt(ji,jj+1)  <= zmmin )     .AND.  & 
    405 !              & ( at_i(ji,jj)  <= zamin .OR. at_i(ji,jj+1) <= zamin ) )    THEN   ;   zmsk01y(ji,jj) = 0._wp 
    406 !           ELSE                                                                   ;   zmsk01y(ji,jj) = 1._wp   ;   ENDIF               
    407 ! END MV TEST DEBUG 
    408  
    409388         END DO 
    410389      END DO    
    411  
    412       CALL lbc_lnk_multi( 'icedyn_rhg_vp', zspgU   ,  'U', -1., zspgV   , 'V', -1. ) 
    413       CALL lbc_lnk_multi( 'icedyn_rhg_vp', ztaux_ai,  'U', -1., ztauy_ai, 'V', -1. ) 
    414  
    415       CALL iom_put( 'zmsk00x'    , zmsk00x  )   ! MV DEBUG 
    416       CALL iom_put( 'zmsk00y'    , zmsk00y  )   ! MV DEBUG 
    417       CALL iom_put( 'zmsk01x'    , zmsk01x  )   ! MV DEBUG 
    418       CALL iom_put( 'zmsk01y'    , zmsk01y  )   ! MV DEBUG 
    419       CALL iom_put( 'ztaux_ai'   , ztaux_ai )   ! MV DEBUG 
    420       CALL iom_put( 'ztauy_ai'   , ztauy_ai )   ! MV DEBUG 
    421       CALL iom_put( 'zspgU'      , zspgU    )   ! MV DEBUG 
    422       CALL iom_put( 'zspgV'      , zspgV    )   ! MV DEBUG 
    423390             
    424391      !------------------------------------------------------------------------------! 
     
    435402      DO i_out = 1, nn_nout_vp 
    436403 
    437          IF( lwp )   WRITE(numout,*) ' outer loop  i_out : ', i_out 
    438  
    439          zu_b_outer(:,:) = u_ice(:,:) ! velocity at previous outer-iterate 
    440          zv_b_outer(:,:) = v_ice(:,:) 
    441404 
    442405         ! Velocities used in the non linear terms are the average of the past two iterates 
    443          ! u_it = 0.5 * ( u_{it-1} + u_{it-2}) 
     406         ! u_it = 0.5 * ( u_{it-1} + u_{it-2} ) 
    444407         ! Also used in Hibler and Ackley (1983); Zhang and Hibler (1997); Lemieux and Tremblay (2009) 
    445  
    446408         zu_c(:,:)       = 0.5_wp * ( u_ice(:,:) + zu_c(:,:) ) 
    447409         zv_c(:,:)       = 0.5_wp * ( v_ice(:,:) + zv_c(:,:) ) 
     
    472434         END DO 
    473435 
    474          CALL lbc_lnk( 'icedyn_rhg_vp', zds, 'F', 1. ) ! MV TEST could be un-necessary according to Gurvan 
    475  
    476          IF ( i_out == nn_nout_vp ) CALL iom_put( 'zds'          , zds          )   ! MV DEBUG 
    477  
    478          IF( lwp )   WRITE(numout,*) ' outer loop  1a i_out : ', i_out 
    479  
    480          !DO jj = 2, jpj - 1    ! loop to jpi,jpj to avoid making a communication for zs1,zs2,zs12 
    481          !   DO ji = 2, jpi - 1 !  
    482  
    483 ! MV DEBUG 
     436         CALL lbc_lnk( 'icedyn_rhg_vp', zds, 'F', 1. ) ! necessary  
     437 
    484438         DO jj = 2, jpj        ! loop to jpi,jpj to avoid making a communication for zs1,zs2,zs12 
    485439            DO ji = 2, jpi     !  
    486 ! END MV DEBUG 
    487440 
    488441               ! shear**2 at T points (doc eq. A16) 
     
    518471            END DO 
    519472         END DO 
    520  
    521 !        zzt(:,:) = 0. ! test 
    522 !        zet(:,:) = 0. ! test 
    523           
    524          CALL lbc_lnk_multi( 'icedyn_rhg_vp', zp_deltastar_t , 'T', 1. , zzt , 'T', 1., zet, 'T', 1. , zdeltat, 'T', 1.) 
    525  
    526          IF ( i_out == nn_nout_vp ) CALL iom_put( 'zzt'        , zzt      )   ! MV DEBUG 
    527          IF ( i_out == nn_nout_vp ) CALL iom_put( 'zet'        , zet      )   ! MV DEBUG 
    528          IF ( i_out == nn_nout_vp ) CALL iom_put( 'zp_deltastar_t', zp_deltastar_t ) ! MV DEBUG 
    529          IF ( i_out == nn_nout_vp ) CALL iom_put( 'zdeltat'       , zdeltat        ) ! MV DEBUG 
    530  
    531          IF( lwp )   WRITE(numout,*) ' outer loop  1b i_out : ', i_out 
     473          
     474!         CALL lbc_lnk_multi( 'icedyn_rhg_vp', zp_deltastar_t , 'T', 1. , zzt , 'T', 1., zet, 'T', 1. , zdeltat, 'T', 1.) ! OK to 
     475!         remove zet and zzt 
     476         CALL lbc_lnk_multi( 'icedyn_rhg_vp', zp_deltastar_t , 'T', 1.,  zdeltat, 'T', 1.) 
    532477 
    533478         DO jj = 1, jpj - 1 
     
    543488         END DO 
    544489 
    545          CALL lbc_lnk( 'icedyn_rhg_vp', zef, 'F', 1. ) 
    546          IF ( i_out == nn_nout_vp ) CALL iom_put( 'zef'           , zef            ) ! MV DEBUG 
    547          IF( lwp )   WRITE(numout,*) ' outer loop  1c i_out : ', i_out 
     490!        CALL lbc_lnk( 'icedyn_rhg_vp', zef, 'F', 1. ) ! Ok to remove 
    548491 
    549492         !--------------------------------------------------- 
     
    581524         END DO 
    582525 
    583          IF( lwp )   WRITE(numout,*) ' outer loop  1d i_out : ', i_out 
    584           
    585          CALL lbc_lnk_multi( 'icedyn_rhg_vp', zCwU ,  'U', -1., zCwV, 'V', -1. ) 
    586          CALL lbc_lnk_multi( 'icedyn_rhg_vp', zCorU,  'U', -1., zCorV, 'V', -1. ) 
    587  
    588          IF ( i_out == nn_nout_vp ) CALL iom_put( 'zCwU'          , zCwU           ) ! MV DEBUG 
    589          IF ( i_out == nn_nout_vp ) CALL iom_put( 'zCwV'          , zCwV           ) ! MV DEBUG 
    590          IF ( i_out == nn_nout_vp ) CALL iom_put( 'zCorU'         , zCorU          ) ! MV DEBUG 
    591          IF ( i_out == nn_nout_vp ) CALL iom_put( 'zCorV'         , zCorV          ) ! MV DEBUG 
    592  
    593          IF( lwp )   WRITE(numout,*) ' outer loop  1f i_out : ', i_out 
    594           
     526!        CALL lbc_lnk_multi( 'icedyn_rhg_vp', zCwU ,  'U', -1., zCwV, 'V', -1. )! OK to remove 
     527!        CALL lbc_lnk_multi( 'icedyn_rhg_vp', zCorU,  'U', -1., zCorV, 'V', -1. )! OK to remove 
     528 
    595529         ! a priori, Coriolis and drag terms only affect diagonal or independent term of the linear system,  
    596530         ! so there is no need for lbclnk on drag and coriolis 
     
    621555         END DO 
    622556 
    623          IF ( i_out == nn_nout_vp ) CALL iom_put( 'zs1_rhsu'      , zs1_rhsu       ) ! MV DEBUG 
    624          IF ( i_out == nn_nout_vp ) CALL iom_put( 'zs2_rhsu'      , zs2_rhsu       ) ! MV DEBUG 
    625          IF ( i_out == nn_nout_vp ) CALL iom_put( 'zs1_rhsv'      , zs1_rhsv       ) ! MV DEBUG 
    626          IF ( i_out == nn_nout_vp ) CALL iom_put( 'zs2_rhsv'      , zs2_rhsv       ) ! MV DEBUG 
    627           
    628557         ! a priori, no lbclnk, because rhsu is only used in the inner domain 
    629558          
     
    631560         ! MV NOTE: I applied zfmask on zds, by mimetism on EVP, but without deep understanding of what I was doing 
    632561         ! My guess is that this is the way to enforce boundary conditions on strain rate tensor 
    633  
    634          IF( lwp )   WRITE(numout,*) ' outer loop 2 i_out : ', i_out 
    635562 
    636563         DO jj = 1, jpj - 1 
     
    646573         END DO 
    647574 
    648          CALL lbc_lnk( 'icedyn_rhg_vp', zs12_rhsu, 'F', 1. ) 
    649          CALL lbc_lnk( 'icedyn_rhg_vp', zs12_rhsv, 'F', 1. ) 
    650  
    651          IF ( i_out == nn_nout_vp ) CALL iom_put( 'zs12_rhsu'     , zs12_rhsu      ) ! MV DEBUG 
    652          IF ( i_out == nn_nout_vp ) CALL iom_put( 'zs12_rhsv'     , zs12_rhsv      ) ! MV DEBUG 
     575         ! CALL lbc_lnk( 'icedyn_rhg_vp', zs12_rhsu, 'F', 1. ) ! Ok to remove 
     576         ! CALL lbc_lnk( 'icedyn_rhg_vp', zs12_rhsv, 'F', 1. ) ! Ok to remove 
    653577 
    654578         ! a priori, no lbclnk, because rhsu are only used in the inner domain 
     
    674598            END DO 
    675599         END DO 
    676  
    677          IF ( i_out == nn_nout_vp ) CALL iom_put( 'zf_rhsu'       , zf_rhsu        ) ! MV DEBUG 
    678          IF ( i_out == nn_nout_vp ) CALL iom_put( 'zf_rhsv'       , zf_rhsv        ) ! MV DEBUG 
    679600          
    680601         !--------------------------- 
     
    692613         END DO 
    693614          
    694          CALL lbc_lnk_multi( 'icedyn_rhg_vp', zrhsu, 'U', -1., zrhsv, 'V',  -1.) 
    695          CALL lbc_lnk_multi( 'icedyn_rhg_vp', zmU_t, 'U', -1., zmV_t, 'V',  -1.) 
    696          CALL lbc_lnk_multi( 'icedyn_rhg_vp', ztaux_oi_rhsu, 'U', -1., ztauy_oi_rhsv, 'V',  -1.) 
    697  
    698          IF ( i_out == nn_nout_vp ) CALL iom_put( 'zmU_t'         , zmU_t          ) ! MV DEBUG 
    699          IF ( i_out == nn_nout_vp ) CALL iom_put( 'zmV_t'         , zmV_t          ) ! MV DEBUG 
    700          IF ( i_out == nn_nout_vp ) CALL iom_put( 'ztaux_oi_rhsu' , ztaux_oi_rhsu  ) ! MV DEBUG 
    701          IF ( i_out == nn_nout_vp ) CALL iom_put( 'ztauy_oi_rhsv' , ztauy_oi_rhsv  ) ! MV DEBUG 
    702          IF ( i_out == nn_nout_vp ) CALL iom_put( 'zrhsu'         , zrhsu          ) ! MV DEBUG 
    703          IF ( i_out == nn_nout_vp ) CALL iom_put( 'zrhsv'         , zrhsv          ) ! MV DEBUG 
    704           
    705          ! inner domain calculations -> no lbclnk 
    706  
    707          IF( lwp )   WRITE(numout,*) ' outer loop 4 i_out : ', i_out 
    708       
     615!        CALL lbc_lnk_multi( 'icedyn_rhg_vp', zrhsu, 'U', -1., zrhsv, 'V',  -1.)! OK to remove 
     616!        CALL lbc_lnk_multi( 'icedyn_rhg_vp', zmU_t, 'U', -1., zmV_t, 'V',  -1.)! OK to remove 
     617!        CALL lbc_lnk_multi( 'icedyn_rhg_vp', ztaux_oi_rhsu, 'U', -1., ztauy_oi_rhsv, 'V',  -1.)! OK to remove 
     618 
    709619         !---------------------------------------------------------------------------------------! 
    710620         ! 
     
    797707         END DO 
    798708 
    799          CALL lbc_lnk_multi( 'icedyn_rhg_vp', zAU  , 'U', 1., zAV  , 'V',  1. ) ! --> here normal that we use '1' and not '-1' ? 
    800          CALL lbc_lnk_multi( 'icedyn_rhg_vp', zBU  , 'U', 1., zBV  , 'V',  1. ) 
    801          CALL lbc_lnk_multi( 'icedyn_rhg_vp', zCU  , 'U', 1., zCV  , 'V',  1. ) 
    802          CALL lbc_lnk_multi( 'icedyn_rhg_vp', zDU  , 'U', 1., zDV  , 'V',  1. ) 
    803          CALL lbc_lnk_multi( 'icedyn_rhg_vp', zEU  , 'U', 1., zEV  , 'V',  1. ) 
    804                 
    805          IF ( i_out == nn_nout_vp ) CALL iom_put( 'zAU'           , zAU            ) ! MV DEBUG 
    806          IF ( i_out == nn_nout_vp ) CALL iom_put( 'zBU'           , zBU            ) ! MV DEBUG 
    807          IF ( i_out == nn_nout_vp ) CALL iom_put( 'zCU'           , zCU            ) ! MV DEBUG 
    808          IF ( i_out == nn_nout_vp ) CALL iom_put( 'zDU'           , zDU            ) ! MV DEBUG 
    809          IF ( i_out == nn_nout_vp ) CALL iom_put( 'zEU'           , zEU            ) ! MV DEBUG 
    810  
    811          IF ( i_out == nn_nout_vp ) CALL iom_put( 'zAV'           , zAV            ) ! MV DEBUG 
    812          IF ( i_out == nn_nout_vp ) CALL iom_put( 'zBV'           , zBV            ) ! MV DEBUG 
    813          IF ( i_out == nn_nout_vp ) CALL iom_put( 'zCV'           , zCV            ) ! MV DEBUG 
    814          IF ( i_out == nn_nout_vp ) CALL iom_put( 'zDV'           , zDV            ) ! MV DEBUG 
    815          IF ( i_out == nn_nout_vp ) CALL iom_put( 'zEV'           , zEV            ) ! MV DEBUG 
     709         ! OK to remove 
     710         !CALL lbc_lnk_multi( 'icedyn_rhg_vp', zAU  , 'U', 1., zAV  , 'V',  1. ) ! --> here normal that we use '1' and not '-1' ? 
     711         !CALL lbc_lnk_multi( 'icedyn_rhg_vp', zBU  , 'U', 1., zBV  , 'V',  1. ) 
     712         !CALL lbc_lnk_multi( 'icedyn_rhg_vp', zCU  , 'U', 1., zCV  , 'V',  1. ) 
     713         !CALL lbc_lnk_multi( 'icedyn_rhg_vp', zDU  , 'U', 1., zDV  , 'V',  1. ) 
     714         !CALL lbc_lnk_multi( 'icedyn_rhg_vp', zEU  , 'U', 1., zEV  , 'V',  1. ) 
    816715 
    817716      !------------------------------------------------------------------------------! 
     
    827726         DO i_inn = 1, nn_ninn_vp ! inner loop iterations 
    828727 
    829             IF( lwp )   WRITE(numout,*) ' inner loop 1 i_inn : ', i_inn 
    830           
    831728            !--- mitgcm computes initial value of residual here... 
    832729 
     
    857754                     ELSE                  ;   jj_min = 3 
    858755                     ENDIF 
    859  
    860                      IF ( lwp ) WRITE(numout,*) ' Into the U-zebra loop at step jn = ', jn, ', with jj_min = ', jj_min 
    861756 
    862757                     DO jj = jj_min, jpj - 1, nn_zebra_vp 
     
    883778                     END DO 
    884779                      
    885                      CALL lbc_lnk( 'icedyn_rhg_vp', zFU, 'U',  -1. ) ! Nov 5 bug test MV zFU is a vector 
     780                     CALL lbc_lnk( 'icedyn_rhg_vp', zFU, 'U',  -1. ) ! -1, zFU is a vector 
    886781                      
    887782                     !--------------- 
     
    904799 
    905800                     END DO 
    906  
    907                      CALL lbc_lnk_multi( 'icedyn_rhg_vp', zFU_prime, 'U',  -1., zBU_prime, 'U',  1. ) ! Nov 5 bug test zFU and zBU 
    908                                                                                                       ! are vectors 
    909   
     801                     ! CALL lbc_lnk_multi( 'icedyn_rhg_vp', zFU_prime, 'U',  -1., zBU_prime, 'U',  1. ) --> can be removed 
     802                                                                                                      
    910803                     !----------------------------- 
    911804                     ! Backward sweep & relaxation 
     
    963856                     ENDIF 
    964857 
    965                      IF ( lwp ) WRITE(numout,*) ' Into the V-zebra loop at step jn = ', jn, ', with ji_min = ', ji_min 
    966           
    967858                     DO ji = ji_min, jpi - 1, nn_zebra_vp  
    968859                      
     
    1011902                     END DO 
    1012903 
    1013                      CALL lbc_lnk_multi( 'icedyn_rhg_vp', zFV_prime, 'V',  -1., zBV_prime, 'V',  1. ) ! MV Nov 5 bug test 
     904                     ! new test 
     905                     ! CALL lbc_lnk_multi( 'icedyn_rhg_vp', zFV_prime, 'V',  -1., zBV_prime, 'V',  1. ) ! MV Nov 5 bug test 
    1014906                      
    1015907                     !----------------------------- 
     
    11281020 
    11291021               CALL rhg_cvg_vp( kt, i_out, i_inn, i_inn_tot, nn_nout_vp, nn_ninn_vp, nn_nvp,        & 
    1130                       &         u_ice, v_ice, zu_b, zv_b, zu_c, zv_c_outer                 ,        & 
     1022                      &         u_ice, v_ice, zu_b, zv_b, zu_c, zv_c,                               & 
    11311023                      &         zmt, za_iU, za_iV, zuerr_max, zverr_max, zglob_area,                & 
    11321024                      &         zrhsu, zAU, zBU, zCU, zDU, zEU, zFU,                                & 
     
    11361028            ENDIF 
    11371029 
    1138             IF ( lwp ) WRITE(numout,*) ' Done convergence tests ' 
    1139  
    11401030         END DO ! i_inn, end of inner loop 
    11411031 
    11421032      END DO ! End of outer loop (i_out) ============================================================================================= 
    11431033 
    1144       IF ( lwp ) WRITE(numout,*) ' We are out of outer loop ' 
    1145  
    1146       CALL iom_put( 'zFU'           , zFU            ) ! MV DEBUG 
    1147       CALL iom_put( 'zBU_prime'     , zBU_prime      ) ! MV DEBUG 
    1148       CALL iom_put( 'zCU_prime'     , zCU_prime      ) ! MV DEBUG 
    1149       CALL iom_put( 'zFU_prime'     , zFU_prime      ) ! MV DEBUG 
    1150  
    1151       CALL iom_put( 'zFV'           , zFV            ) ! MV DEBUG 
    1152       CALL iom_put( 'zBV_prime'     , zBV_prime      ) ! MV DEBUG 
    1153       CALL iom_put( 'zCV_prime'     , zCV_prime      ) ! MV DEBUG 
    1154       CALL iom_put( 'zFV_prime'     , zFV_prime      ) ! MV DEBUG 
    1155  
    11561034      CALL lbc_lnk_multi( 'icedyn_rhg_vp', u_ice, 'U', -1., v_ice, 'V', -1. ) 
    1157  
    11581035 
    11591036      IF( ln_rhg_chkcvg ) THEN 
     
    11671044         
    11681045      ENDIF ! ln_rhg_chkcvg 
    1169  
    1170       ! MV DEBUG test - replace ice velocity by ocean current to give the model the means to go ahead 
    1171       ! info: probably wrong because when doing du/dt = div ( sigma ) we have spurious waves instead of u = v = 0 
    1172  
    1173 !     DO jj = 2, jpj - 1 
    1174 !        DO ji = 2, jpi - 1    
    1175  
    1176 !            u_ice(ji,jj) =   zmsk00x(ji,jj)                               &  
    1177 !     &         * (           zmsk01x(ji,jj)   * u_oce(ji,jj) * 0.01_wp    & 
    1178 !                 + ( 1._wp - zmsk01x(ji,jj) ) * u_oce(ji,jj) * 0.01_wp    ) 
    1179  
    1180 !            v_ice(ji,jj) =   zmsk00y(ji,jj)                               &  
    1181 !     &         * (           zmsk01y(ji,jj)   * v_oce(ji,jj) * 0.01_wp    & 
    1182 !                 + ( 1._wp - zmsk01y(ji,jj) ) * v_oce(ji,jj) * 0.01_wp    ) 
    1183  
    1184 !        END DO 
    1185 !     END DO 
    1186  
    1187 !     CALL lbc_lnk_multi( 'icedyn_rhg_vp', u_ice, 'U', -1., v_ice, 'V', -1. ) 
    1188  
    1189 !     IF ( lwp ) WRITE(numout,*) ' Velocity replaced ' 
    1190  
    1191       ! END DEBUG 
    11921046 
    11931047      !------------------------------------------------------------------------------! 
     
    12421096            zdelta         = SQRT( zdiv2 + ( zdt2 + zds2 ) * z1_ecc2 ) 
    12431097            rswitch        = 1._wp - MAX( 0._wp, SIGN( 1._wp, -zdelta ) ) ! 0 if delta=0 
    1244              
    1245             !pdelta_i(ji,jj) = zdelta + rn_creepl * rswitch 
    1246             pdelta_i(ji,jj) = zdelta + rn_creepl 
     1098 
     1099            pdelta_i(ji,jj) = zdelta + rn_creepl ! * rswitch 
    12471100 
    12481101         END DO 
    12491102      END DO 
    1250  
    1251       IF ( lwp ) WRITE(numout,*) ' Deformation recalculated ' 
    12521103       
    12531104      CALL lbc_lnk_multi( 'icedyn_rhg_vp', pshear_i, 'T', 1., pdivu_i, 'T', 1., pdelta_i, 'T', 1. ) 
     
    13021153          
    13031154      ENDIF 
    1304  
    1305       IF ( lwp ) WRITE(numout,*) ' zs12f recalculated ' 
    13061155 
    13071156      ! 
     
    13581207      IF( iom_use('icestr') )   CALL iom_put( 'icestr' , strength * zmsk00 )   ! strength 
    13591208 
    1360       IF ( lwp ) WRITE(numout,*) 'Some terms recalculated ' 
    1361  
    13621209      ! --- Stress tensor invariants (SIMIP diags) --- ! 
    13631210      IF( iom_use('normstr') .OR. iom_use('sheastr') ) THEN 
     
    13871234          
    13881235      ENDIF 
    1389  
    1390       IF ( lwp ) WRITE(numout,*) 'SIMIP work done' 
    13911236 
    13921237      ! --- Normalized stress tensor principal components --- ! 
     
    14211266            END DO 
    14221267         END DO 
    1423          IF ( lwp ) WRITE(numout,*) 'Some shitty stress work done' 
    14241268         ! 
    14251269         CALL lbc_lnk_multi( 'icedyn_rhg_vp', zsig1_p, 'T', 1., zsig2_p, 'T', 1.) 
    1426          !       
    1427          IF ( lwp ) WRITE(numout,*) ' Beauaaaarflblbllll ' 
    14281270         ! 
    14291271         CALL iom_put( 'sig1_pnorm' , zsig1_p )  
     
    14311273 
    14321274         DEALLOCATE( zsig1_p , zsig2_p , zsig_I , zsig_II ) 
    1433  
    1434          IF ( lwp ) WRITE(numout,*) ' So what ??? ' 
    14351275          
    14361276      ENDIF 
     
    18831723         ENDIF 
    18841724 
    1885          ! spatially-dependent things (to be coded) 
    1886          ! velocity difference ? 
    1887          ! residual ? 
    1888 !        ! 
    1889 !        IF ( kiter == kitermax ) THEN 
    1890 !           WRITE(numout,*) ' Should plot the spatially dependent residual ' 
    1891 !           istatus = NF90_PUT_VAR( ncvgid, nvarid_ures_xy, (/zu_res/) )          ! U-residual, spatially dependent 
    1892 !           istatus = NF90_PUT_VAR( ncvgid, nvarid_vres_xy, (/zv_res/) )          ! V-residual, spatially dependent 
    1893 !        ENDIF 
    1894  
    1895          ! close file 
    1896          ! IF( kt == nitend )   istatus = NF90_CLOSE( ncvgid ) ! BUG Ed blockley 
    18971725         IF( kt == nitend - nn_fsbc + 1 .AND. kitinntot == kitinntotmax )    istatus = NF90_CLOSE( ncvgid ) 
    18981726      ENDIF 
Note: See TracChangeset for help on using the changeset viewer.