Ignore:
Timestamp:
2017-09-14T17:52:02+02:00 (3 years ago)
Author:
clem
Message:

changes in style - part6 - ice diffusion (hope the scheme still converges)

File:
1 edited

Legend:

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

    r8514 r8522  
    7878      REAL(wp) ::   zgrr         ! bottom growth rate 
    7979      REAL(wp) ::   zt_i_new     ! bottom formation temperature 
     80      REAL(wp) ::   z1_rho       ! 1/(rhosn+rau0-rhoic) 
    8081 
    8182      REAL(wp) ::   zQm          ! enthalpy exchanged with the ocean (J/m2), >0 towards the ocean 
     
    596597      ! When snow load excesses Archimede's limit, snow-ice interface goes down under sea-level,  
    597598      ! flooding of seawater transforms snow into ice dh_snowice is positive for the ice 
     599      z1_rho = 1._wp / ( rhosn+rau0-rhoic ) 
    598600      DO ji = 1, nidx 
    599601         ! 
    600          dh_snowice(ji) = MAX(  0._wp , ( rhosn * ht_s_1d(ji) + (rhoic-rau0) * ht_i_1d(ji) ) / ( rhosn+rau0-rhoic ) ) 
     602         dh_snowice(ji) = MAX(  0._wp , ( rhosn * ht_s_1d(ji) + (rhoic-rau0) * ht_i_1d(ji) ) * z1_rho ) 
    601603 
    602604         ht_i_1d(ji)    = ht_i_1d(ji) + dh_snowice(ji) 
     
    646648            e_s_1d(ji,jk) = rswitch * e_s_1d(ji,jk) 
    647649            ! recalculate t_s_1d from e_s_1d 
    648             t_s_1d(ji,jk) = rt0 + rswitch * ( - e_s_1d(ji,jk) / ( rhosn * cpic ) + lfus / cpic ) 
     650            t_s_1d(ji,jk) = rt0 + rswitch * ( - e_s_1d(ji,jk) * r1_rhosn * r1_cpic + lfus * r1_cpic ) 
    649651         END DO 
    650652      END DO 
Note: See TracChangeset for help on using the changeset viewer.