Changeset 1501 for LMDZ5/trunk/libf/phylmd

Timestamp:

Mar 15, 2011, 5:12:29 PM (14 years ago)

Author:

idelkadi

Message:

Rajout de la nouvelle formulation de l'efficacite des precipitations fonction de la temperature et contenu en eau liquide

File:

• LMDZ5/trunk/libf/phylmd/cv3_routines.F (modified) (2 diffs)

LMDZ5/trunk/libf/phylmd/cv3_routines.F

@@ -936 +936 @@
       logical lcape(nloc)
       integer iposit(nloc)
+      real elcrit, tlcrit, sigs
+      parameter(sigs=0.15)
+      integer flag_epKEorig
 
 !=====================================================================
@@ -1044 +1047 @@
 ! --- THESE MAY BE FUNCTIONS OF TP(I), P(I) AND CLW(I)
 !=====================================================================
-c
-c ori      do 320 k=minorig+1,nl
-      do 320 k=1,nl ! convect3
+cIM beg: ajout fis. reglage ep
+      flag_epKEorig=0
+      elcrit=0.0011
+      tlcrit=-55.0
+cIM Lecture du fichier ep_param.data
+      OPEN(99,file='ep_param.data',status='old',
+     $          form='formatted',err=9999)
+      READ(99,*,end=9998) flag_epKEorig
+      READ(99,*,end=9998) elcrit
+      READ(99,*,end=9998) tlcrit
+9998  Continue
+      CLOSE(99)
+9999  Continue
+      WRITE(*,*)'flag_epKEorig',flag_epKEorig
+      WRITE(*,*)'elcrit=',elcrit
+      WRITE(*,*)'tlcrit=',tlcrit
+cIM end: ajout fis. reglage ep
+
+      if(flag_epKEorig.ne.1) THEN
+        do 320 k=1,nl ! convect3
         do 310 i=1,ncum
            pden=ptcrit-pbcrit
            ep(i,k)=(plcl(i)-p(i,k)-pbcrit)/pden*epmax
-           ep(i,k)=max(ep(i,k),0.0)
+           ep(i,k)=amax1(ep(i,k),0.0)
            ep(i,k)=amin1(ep(i,k),epmax)
            sigp(i,k)=spfac
-c ori          if(k.ge.(nk(i)+1))then
-c ori            tca=tp(i,k)-t0
-c ori            if(tca.ge.0.0)then
-c ori              elacrit=elcrit
-c ori            else
-c ori              elacrit=elcrit*(1.0-tca/tlcrit)
-c ori            endif
-c ori            elacrit=max(elacrit,0.0)
-c ori            ep(i,k)=1.0-elacrit/max(clw(i,k),1.0e-8)
-c ori            ep(i,k)=max(ep(i,k),0.0 )
-c ori            ep(i,k)=min(ep(i,k),1.0 )
-c ori            sigp(i,k)=sigs
-c ori          endif
  310    continue
- 320  continue
-c
+ 320    continue
+      else
+        do 325 k=1,nl
+        do 315 i=1,ncum
+          if(k.ge.(nk(i)+1))then
+              tca=tp(i,k)-t0
+              if(tca.ge.0.0)then
+               elacrit=elcrit
+              else
+                elacrit=elcrit*(1.0-tca/tlcrit)
+              endif
+              elacrit=max(elacrit,0.0)
+              ep(i,k)=1.0-elacrit/max(clw(i,k),1.0e-8)
+              ep(i,k)=max(ep(i,k),0.0 )
+              ep(i,k)=min(ep(i,k),epmax )
+              sigp(i,k)=sigs
+          endif
+ 315    continue
+ 325    continue
+      endif
 !=====================================================================
 ! --- CALCULATE VIRTUAL TEMPERATURE AND LIFTED PARCEL