/[lmdze]/trunk/phylmd/fisrtilp.f
ViewVC logotype

Diff of /trunk/phylmd/fisrtilp.f

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 22 by guez, Fri Jul 31 15:18:47 2009 UTC revision 52 by guez, Fri Sep 23 12:28:01 2011 UTC
# Line 11  SUBROUTINE fisrtilp(dtime,paprs,pplay,t, Line 11  SUBROUTINE fisrtilp(dtime,paprs,pplay,t,
11    USE dimens_m    USE dimens_m
12    USE dimphy    USE dimphy
13    USE tracstoke    USE tracstoke
14    USE yomcst    USE suphec_m
15    USE yoethf    USE yoethf_m
16    USE fcttre    USE fcttre
17    USE comfisrtilp    USE comfisrtilp
18    use numer_rec, only: erf    use numer_rec, only: nr_erf
19    
20    IMPLICIT NONE    IMPLICIT NONE
21    
# Line 24  SUBROUTINE fisrtilp(dtime,paprs,pplay,t, Line 24  SUBROUTINE fisrtilp(dtime,paprs,pplay,t,
24    REAL, INTENT (IN) :: & ! intervalle du temps (s)                    REAL, INTENT (IN) :: & ! intervalle du temps (s)                
25         dtime         dtime
26    REAL, INTENT (IN) :: paprs(klon,klev+1) ! pression a inter-couche      REAL, INTENT (IN) :: paprs(klon,klev+1) ! pression a inter-couche  
27    REAL, INTENT (IN) :: pplay(klon,klev) ! pression au milieu de couche    REAL, INTENT (IN):: pplay(klon,klev) ! pression au milieu de couche
28    REAL t(klon,klev) ! temperature (K)                                  REAL, INTENT (IN):: t(klon,klev) ! temperature (K)
29    REAL q(klon,klev) ! humidite specifique (kg/kg)                      REAL q(klon,klev) ! humidite specifique (kg/kg)                  
30    REAL d_t(klon,klev) ! incrementation de la temperature (K)            REAL d_t(klon,klev) ! incrementation de la temperature (K)        
31    REAL d_q(klon,klev) ! incrementation de la vapeur d'eau              REAL d_q(klon,klev) ! incrementation de la vapeur d'eau          
# Line 338  SUBROUTINE fisrtilp(dtime,paprs,pplay,t, Line 338  SUBROUTINE fisrtilp(dtime,paprs,pplay,t,
338                zpdf_b(i) = zpdf_k(i)/(2.*sqrt(2.))                zpdf_b(i) = zpdf_k(i)/(2.*sqrt(2.))
339                zpdf_e1(i) = zpdf_a(i) - zpdf_b(i)                zpdf_e1(i) = zpdf_a(i) - zpdf_b(i)
340                zpdf_e1(i) = sign(min(abs(zpdf_e1(i)),5.),zpdf_e1(i))                zpdf_e1(i) = sign(min(abs(zpdf_e1(i)),5.),zpdf_e1(i))
341                zpdf_e1(i) = 1. - erf(zpdf_e1(i))                zpdf_e1(i) = 1. - nr_erf(zpdf_e1(i))
342                zpdf_e2(i) = zpdf_a(i) + zpdf_b(i)                zpdf_e2(i) = zpdf_a(i) + zpdf_b(i)
343                zpdf_e2(i) = sign(min(abs(zpdf_e2(i)),5.),zpdf_e2(i))                zpdf_e2(i) = sign(min(abs(zpdf_e2(i)),5.),zpdf_e2(i))
344                zpdf_e2(i) = 1. - erf(zpdf_e2(i))                zpdf_e2(i) = 1. - nr_erf(zpdf_e2(i))
345                IF (zpdf_e1(i)<1.E-10) THEN                IF (zpdf_e1(i)<1.E-10) THEN
346                   rneb(i,k) = 0.                   rneb(i,k) = 0.
347                   zqn(i) = zqs(i)                   zqn(i) = zqs(i)
Legend:
Removed from v.22  
changed lines
  Added in v.52
  ViewVC Help
Powered by ViewVC 1.1.21