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/icethd.F90

    r8515 r8517  
    8181      INTEGER  :: ji, jj, jk, jl   ! dummy loop indices 
    8282      REAL(wp) :: zfric_u, zqld, zqfr, zqfr_neg 
    83       REAL(wp) :: zvi_b, zsmv_b, zei_b, zfs_b, zfw_b, zft_b ! conservation check 
    8483      REAL(wp), PARAMETER :: zfric_umin = 0._wp           ! lower bound for the friction velocity (cice value=5.e-04) 
    8584      REAL(wp), PARAMETER :: zch        = 0.0057_wp       ! heat transfer coefficient 
     
    8786      ! 
    8887      !!------------------------------------------------------------------- 
    89  
    90       IF( nn_timing == 1 )   CALL timing_start('icethd') 
     88      ! controls 
     89      IF( nn_timing == 1 )   CALL timing_start('icethd')                                                             ! timing 
     90      IF( ln_icediachk   )   CALL ice_cons_hsm(0, 'icethd', rdiag_v, rdiag_s, rdiag_t, rdiag_fv, rdiag_fs, rdiag_ft) ! conservation 
    9191 
    9292      IF( kt == nit000 .AND. lwp ) THEN 
     
    9696      ENDIF 
    9797       
    98       ! conservation test 
    99       IF( ln_icediachk )   CALL ice_cons_hsm( 0, 'icethd', zvi_b, zsmv_b, zei_b, zfw_b, zfs_b, zft_b ) 
    100  
    10198      CALL ice_var_glo2eqv 
    10299 
     
    257254      oa_i(:,:,:) = o_i(:,:,:) * a_i(:,:,:) 
    258255 
    259       IF( ln_icediachk )   CALL ice_cons_hsm( 1, 'icethd', zvi_b, zsmv_b, zei_b, zfw_b, zfs_b, zft_b ) 
     256      IF( ln_icediachk )   CALL ice_cons_hsm(1, 'icethd', rdiag_v, rdiag_s, rdiag_t, rdiag_fv, rdiag_fs, rdiag_ft) 
    260257      ! 
    261258                           CALL ice_var_zapsmall           ! --- remove very small ice concentration (<1e-10) --- ! 
     
    266263      IF( ln_icedO )       CALL ice_thd_lac                ! --- frazil ice growing in leads --- ! 
    267264      ! 
    268       IF( ln_icectl )      CALL ice_prt( kt, iiceprt, jiceprt, 1, ' - ice thermodyn. - ' )   ! control print 
    269       IF( ln_ctl )         CALL ice_prt3D( 'icethd' )      ! Control print 
    270       ! 
    271       IF( nn_timing == 1 )  CALL timing_stop('icethd') 
     265      ! controls 
     266      IF( ln_icectl      )   CALL ice_prt    (kt, iiceprt, jiceprt, 1, ' - ice thermodyn. - ') ! prints 
     267      IF( ln_ctl         )   CALL ice_prt3D  ('icethd')                                        ! prints 
     268      IF( nn_timing == 1 )   CALL timing_stop('icethd')                                        ! timing 
    272269      ! 
    273270   END SUBROUTINE ice_thd  
     
    539536      INTEGER  ::   ios   ! Local integer output status for namelist read 
    540537      !! 
    541       NAMELIST/namice_thd/ ln_icethd, rn_kappa_i, ln_cndi_U64, ln_cndi_P07, ln_dqns_i, rn_cnd_s,   & 
     538      NAMELIST/namice_thd/ rn_kappa_i, ln_cndi_U64, ln_cndi_P07, ln_dqns_i, rn_cnd_s,   & 
    542539         &                 ln_icedH, rn_blow_s,                                                    & 
    543540         &                 ln_icedA, rn_beta, rn_dmin,                                             & 
     
    560557         WRITE(numout,*) '~~~~~~~~~~~~~' 
    561558         WRITE(numout,*) '   Namelist namice_thd' 
    562          WRITE(numout,*) '      Ice thermodynamics (T) or not (F)                       ln_icethd    = ', ln_icethd 
    563559         WRITE(numout,*) '   -- icethd_dif --' 
    564560         WRITE(numout,*) '      extinction radiation parameter in sea ice               rn_kappa_i   = ', rn_kappa_i 
Note: See TracChangeset for help on using the changeset viewer.