Ignore:
Timestamp:
2017-09-12T20:46:13+02:00 (3 years ago)
Author:
clem
Message:

changes in style - part6 - one more round

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3/icerdgrft.F90

    r8515 r8517  
    127127      ! 
    128128      INTEGER, PARAMETER ::   nitermax = 20     
    129       ! 
    130       REAL(wp) :: zvi_b, zsmv_b, zei_b, zfs_b, zfw_b, zft_b  
    131129      !!----------------------------------------------------------------------------- 
    132       IF( nn_timing == 1 )  CALL timing_start('icerdgrft') 
     130      ! controls 
     131      IF( nn_timing == 1 )   CALL timing_start('icerdgrft')                                                             ! timing 
     132      IF( ln_icediachk   )   CALL ice_cons_hsm(0, 'icerdgrft', rdiag_v, rdiag_s, rdiag_t, rdiag_fv, rdiag_fs, rdiag_ft) ! conservation 
    133133 
    134134      IF( kt == nit000 ) THEN 
     
    140140         ! 
    141141      ENDIF       
    142       !                    ! conservation test 
    143       IF( ln_icediachk )   CALL ice_cons_hsm(0, 'icerdgrft', zvi_b, zsmv_b, zei_b, zfw_b, zfs_b, zft_b) 
    144142 
    145143      !-----------------------------------------------------------------------------! 
     
    284282      CALL ice_var_agg( 1 )  
    285283 
    286       !-----------------------------------------------------------------------------! 
    287       ! control prints 
    288       !-----------------------------------------------------------------------------! 
    289       !                    ! conservation test 
    290       IF( ln_icediachk )   CALL ice_cons_hsm(1, 'icerdgrft', zvi_b, zsmv_b, zei_b, zfw_b, zfs_b, zft_b) 
    291  
    292       !                    ! control prints 
    293       IF( ln_ctl       )   CALL ice_prt3D( 'icerdgrft' ) 
    294       ! 
    295       IF( nn_timing == 1 )  CALL timing_stop('icerdgrft') 
     284      ! controls 
     285      IF( ln_icediachk   )   CALL ice_cons_hsm(1, 'icerdgrft', rdiag_v, rdiag_s, rdiag_t, rdiag_fv, rdiag_fs, rdiag_ft) ! conservation 
     286      IF( ln_ctl         )   CALL ice_prt3D   ('icerdgrft')                                                             ! prints 
     287      IF( nn_timing == 1 )   CALL timing_stop ('icerdgrft')                                                             ! timing 
    296288      ! 
    297289   END SUBROUTINE ice_rdgrft 
     
    795787      ELSEIF( ln_str_H79 ) THEN      ! Ice strength => Hibler (1979) method             ! 
    796788      !                              !--------------------------------------------------! 
    797          strength(:,:) = rn_pstar * vt_i(:,:) * EXP( - rn_crhg * ( 1._wp - at_i(:,:) )  ) * tmask(:,:,1) 
     789         strength(:,:) = rn_pstar * SUM( v_i(:,:,:), dim=3 ) * EXP( -rn_crhg * ( 1._wp - SUM( a_i(:,:,:), dim=3 ) ) ) 
    798790         ! 
    799791         ismooth = 1 
Note: See TracChangeset for help on using the changeset viewer.