Changeset 2961 for LMDZ5/trunk


Ignore:
Timestamp:
Jul 20, 2017, 8:55:20 PM (7 years ago)
Author:
jbmadeleine
Message:

Added iflag_t_glace=3 which uses a better function of the ice fraction:
iflag_t_glace.LT.3 uses icefrac = (1.-(T-Tmin)/(Tmax-Tmin))nx
iflag_t_glace.GE.3 uses icefrac = 1. - ((T-Tmin)/(Tmax-Tmin))
nx

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ5/trunk/libf/phylmd/icefrac_lsc_mod.F90

    r2507 r2961  
    2828  REAL, DIMENSION(np), INTENT(OUT) :: icefrac
    2929
    30   REAL :: sig0,www,tmin_tmp,icefrac_tmp
     30  REAL :: sig0,www,tmin_tmp,liqfrac_tmp
    3131  INTEGER :: ip
    3232
     
    4343     ENDIF
    4444     tmin_tmp=www*t_glace_max+(1.-www)*t_glace_min
    45      icefrac_tmp= 1.0 - (temp(ip)-tmin_tmp) / (t_glace_max-tmin_tmp)
    46      icefrac_tmp = MIN(MAX(icefrac_tmp,0.0),1.0)
    47      icefrac(ip) = icefrac_tmp**exposant_glace
     45     liqfrac_tmp=  (temp(ip)-tmin_tmp) / (t_glace_max-tmin_tmp)
     46     liqfrac_tmp = MIN(MAX(liqfrac_tmp,0.0),1.0)
     47     IF (iflag_t_glace.GE.3) THEN
     48       icefrac(ip) = 1.0-liqfrac_tmp**exposant_glace
     49     ELSE
     50       icefrac(ip) = (1.0-liqfrac_tmp)**exposant_glace
     51     ENDIF
    4852  ENDDO
    4953
Note: See TracChangeset for help on using the changeset viewer.