Ignore:
Timestamp:
2020-04-09T16:30:48+02:00 (8 months ago)
Author:
clem
Message:

clean useless variables

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2020/r4.0-HEAD_r12713_clem_dan_fixcpl/src/ICE/icethd_pnd.F90

    r12725 r12726  
    8888         ! 
    8989         IF( a_i_1d(ji) > 0._wp .AND. t_su_1d(ji) >= rt0 ) THEN 
    90             a_ip_frac_1d(ji) = rn_apnd 
    9190            h_ip_1d(ji)      = rn_hpnd     
    92             a_ip_1d(ji)      = a_ip_frac_1d(ji) * a_i_1d(ji) 
     91            a_ip_1d(ji)      = rn_apnd * a_i_1d(ji) 
    9392            h_il_1d(ji)      = 0._wp    ! no pond lids whatsoever 
    9493         ELSE 
    95             a_ip_frac_1d(ji) = 0._wp 
    9694            h_ip_1d(ji)      = 0._wp     
    9795            a_ip_1d(ji)      = 0._wp 
     
    184182            !--- Remove ponds on thin ice or tiny ice fractions 
    185183            a_ip_1d(ji)      = 0._wp 
    186             a_ip_frac_1d(ji) = 0._wp 
    187184            h_ip_1d(ji)      = 0._wp 
    188185            h_il_1d(ji)      = 0._wp 
     
    261258            ! a_ip/a_i = a_ip_frac = h_ip / zaspect (cf Holland 2012, fitting SHEBA so that knowing v_ip we can distribute it to a_ip and h_ip) 
    262259            a_ip_1d(ji)      = MIN( a_i_1d(ji), SQRT( v_ip_1d(ji) * z1_aspect * a_i_1d(ji) ) ) ! make sure a_ip < a_i 
    263             a_ip_frac_1d(ji) = a_ip_1d(ji) / a_i_1d(ji) 
    264             h_ip_1d(ji)      = zaspect * a_ip_frac_1d(ji) 
     260            h_ip_1d(ji)      = zaspect * a_ip_1d(ji) / a_i_1d(ji) 
    265261 
    266262            !---------------!             
     
    269265            IF( ln_pnd_flush ) THEN 
    270266               ! height of top of the pond above sea-level 
    271                zhp = ( h_i_1d(ji) * ( rau0 - rhoi ) + h_ip_1d(ji) * ( rau0 - rhow * a_ip_frac_1d(ji) ) ) * r1_rau0 
     267               zhp = ( h_i_1d(ji) * ( rau0 - rhoi ) + h_ip_1d(ji) * ( rau0 - rhow * a_ip_1d(ji) / a_i_1d(ji) ) ) * r1_rau0 
    272268 
    273269               ! Calculate the permeability of the ice (Assur 1958) 
     
    288284               !--- Set new pond area and depth ---! assuming linear relation between h_ip and a_ip_frac 
    289285               a_ip_1d(ji)      = MIN( a_i_1d(ji), SQRT( v_ip_1d(ji) * z1_aspect * a_i_1d(ji) ) ) ! make sure a_ip < a_i 
    290                a_ip_frac_1d(ji) = a_ip_1d(ji) / a_i_1d(ji) 
    291                h_ip_1d(ji)      = zaspect * a_ip_frac_1d(ji) 
     286               h_ip_1d(ji)      = zaspect * a_ip_1d(ji) / a_i_1d(ji) 
    292287 
    293288            ENDIF 
     
    302297               IF ( h_il_1d(ji) > h_ip_1d(ji) * 10._wp ) THEN 
    303298                  a_ip_1d(ji)      = 0._wp 
    304                   a_ip_frac_1d(ji) = 0._wp 
    305299                  h_ip_1d(ji)      = 0._wp 
    306300                  h_il_1d(ji)      = 0._wp 
Note: See TracChangeset for help on using the changeset viewer.