Ignore:
Timestamp:
2017-09-22T16:55:24+02:00 (3 years ago)
Author:
clem
Message:

almost useless commits

File:
1 edited

Legend:

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

    r8534 r8559  
    5959      REAL(wp) ::   iflush, igravdr              ! local scalars 
    6060      REAL(wp) ::   zs_sni, zsm_i_gd, zsm_i_fl, zsm_i_si, zsm_i_bg   ! local scalars 
     61      REAL(wp) ::   z1_time_gd, z1_time_fl 
    6162      !!--------------------------------------------------------------------- 
    6263 
     
    6667      CASE( 2 )       !  time varying salinity with linear profile  ! 
    6768      !               !---------------------------------------------! 
     69         z1_time_gd = 1._wp / rn_time_gd * rdt_ice 
     70         z1_time_fl = 1._wp / rn_time_fl * rdt_ice 
     71         ! 
    6872         DO ji = 1, nidx 
    6973 
     
    7175            !  Update ice salinity from snow-ice and bottom growth 
    7276            !--------------------------------------------------------- 
    73             zs_sni   = sss_1d(ji) * ( rhoic - rhosn ) * r1_rhoic   ! Salinity of snow ice 
    74             rswitch  = MAX( 0._wp , SIGN( 1._wp , ht_i_1d(ji) - epsi20 ) ) 
    75             zsm_i_si = ( zs_sni      - sm_i_1d(ji) ) *             dh_snowice(ji)  / MAX( ht_i_1d(ji), epsi20 ) * rswitch ! snow-ice     
    76             zsm_i_bg = ( s_i_new(ji) - sm_i_1d(ji) ) * MAX( 0._wp, dh_i_bott(ji) ) / MAX( ht_i_1d(ji), epsi20 ) * rswitch ! bottom growth 
    77  
    78             ! Update salinity (nb: salt flux already included in icethd_dh) 
    79             sm_i_1d(ji) = sm_i_1d(ji) + zsm_i_bg + zsm_i_si 
     77            IF( ht_i_1d(ji) > 0._wp ) THEN 
     78               zs_sni   = sss_1d(ji) * ( rhoic - rhosn ) * r1_rhoic                                 ! Salinity of snow ice 
     79               zsm_i_si = ( zs_sni      - sm_i_1d(ji) ) *             dh_snowice(ji)  / ht_i_1d(ji) ! snow-ice     
     80               zsm_i_bg = ( s_i_new(ji) - sm_i_1d(ji) ) * MAX( 0._wp, dh_i_bott(ji) ) / ht_i_1d(ji) ! bottom growth 
     81               ! Update salinity (nb: salt flux already included in icethd_dh) 
     82               sm_i_1d(ji) = sm_i_1d(ji) + zsm_i_bg + zsm_i_si 
     83            ENDIF 
    8084 
    8185            IF( ld_sal ) THEN 
     
    8589               iflush   = MAX( 0._wp , SIGN( 1._wp , t_su_1d(ji) - rt0         ) )  ! =1 if summer  
    8690               igravdr  = MAX( 0._wp , SIGN( 1._wp , t_bo_1d(ji) - t_su_1d(ji) ) )  ! =1 if t_su < t_bo 
    87                zsm_i_gd = - igravdr * MAX( sm_i_1d(ji) - rn_sal_gd , 0._wp ) / rn_time_gd * rdt_ice  ! gravity drainage  
    88                zsm_i_fl = - iflush  * MAX( sm_i_1d(ji) - rn_sal_fl , 0._wp ) / rn_time_fl * rdt_ice  ! flushing 
     91 
     92               zsm_i_gd = - igravdr * MAX( sm_i_1d(ji) - rn_sal_gd , 0._wp ) * z1_time_gd  ! gravity drainage  
     93               zsm_i_fl = - iflush  * MAX( sm_i_1d(ji) - rn_sal_fl , 0._wp ) * z1_time_fl  ! flushing 
    8994                
    9095               ! Update salinity    
Note: See TracChangeset for help on using the changeset viewer.