source: LMDZ5/trunk/libf/phylmd/icefrac_lsc_mod.F90 @ 2451

Last change on this file since 2451 was 2109, checked in by jbmadeleine, 10 years ago

nettoyage de icefrac_lsc: les paramètres t_glace_min, t_glace_max et
exposant_glace sont à présents chargés à partir de nuage.h et non
donnés en argument;

cleaning of icefrac_lsc: parameters t_glace_min, t_glace_max and
exposant_glace are loaded from the nuage.h file instead of being
given as arguments;

  • Property svn:executable set to *
File size: 1.2 KB
Line 
1!
2! $Header$
3!
4MODULE icefrac_lsc_mod
5
6IMPLICIT NONE
7
8CONTAINS
9!*******************************************************************
10
11SUBROUTINE icefrac_lsc(np,temp, sig, icefrac)
12  !
13  ! Compute the ice fraction 1-xliq (see e.g.
14  ! Doutriaux-Boucher & Quaas 2004, section 2.2.)
15  !
16  ! (JBM 3/14 8/14)
17 
18  INCLUDE "nuage.h"
19
20  ! nuage.h contains:
21  ! t_glace_min: if T < Tmin, the cloud is only made of water ice
22  ! t_glace_max: if T > Tmax, the cloud is only made of liquid water
23  ! exposant_glace: controls the sharpness of the transition
24  INTEGER :: np
25  REAL, DIMENSION(np), INTENT(IN) :: temp ! temperature
26  REAL, DIMENSION(np), INTENT(IN) :: sig
27  REAL, DIMENSION(np), INTENT(OUT) :: icefrac
28
29  REAL :: sig0,www,tmin_tmp,icefrac_tmp
30  INTEGER :: ip
31
32  sig0=0.8
33
34  DO ip=1,np
35     www=(max(sig(ip)-sig0,0.))/(1.-sig0) ! w=1 at the surface and 0 for sig < sig0
36     tmin_tmp=www*t_glace_max+(1.-www)*t_glace_min
37     icefrac_tmp= 1.0 - (temp(ip)-tmin_tmp) / (t_glace_max-tmin_tmp)
38     icefrac_tmp = MIN(MAX(icefrac_tmp,0.0),1.0)
39     icefrac(ip) = icefrac_tmp**exposant_glace
40  ENDDO
41
42  RETURN
43END SUBROUTINE icefrac_lsc
44
45!*******************************************************************
46!
47END MODULE icefrac_lsc_mod
Note: See TracBrowser for help on using the repository browser.