Changeset 953 for LMDZ4/trunk/libf/phylmd/wake.F

Timestamp:

Apr 21, 2008, 6:46:08 PM (17 years ago)

Author:

lmdzadmin

Message:

On supprime variables pbase, ibas qui servent pas ; print sous clef JYG/FH
IM

File:

• LMDZ4/trunk/libf/phylmd/wake.F (modified) (10 diffs)

LMDZ4/trunk/libf/phylmd/wake.F

@@ -1 +1 @@
       Subroutine WAKE (p,ph,ppi,dtime,sigd_con
-     :                ,te0,qe0,omgb,ibas
+     :                ,te0,qe0,omgb
      :                ,dtdwn,dqdwn,amdwn,amup,dta,dqa
      :                ,wdtPBL,wdqPBL,udtPBL,udqPBL
@@ -61 +61 @@
 C                        qe0  : humidite dans l'environnement     (kg/kg)
 C                        omgb : vitesse verticale moyenne sur la maille (Pa/s)
-C                        ibas : cloud base level number
 C                        dtdwn: source de chaleur due aux descentes (K/s)
 C                        dqdwn: source d'humidite due aux descentes (kg/kg/s)
@@ -116 +115 @@
 #include "YOMCST.h"
 #include "cvthermo.h"
+#include "iniprint.h"
 
 c Arguments en entree
@@ -124 +124 @@
       REAL te0(klev),qe0(klev)
       REAL omgb(klev+1)
-      INTEGER ibas
       REAL dtdwn(klev), dqdwn(klev)
       REAL wdtPBL(klev),wdqPBL(klev)
@@ -198 +197 @@
       REAL epaisseur1(klev), epaisseur2(klev)
 
-      REAL pbase
-      
       REAL the(klev), thu(klev)
 
@@ -264 +261 @@
 
       delta_t_min = 0.2
-
-
-C Cloud base
-c-------------------------------------------------------------------------
-
-      Pbase = P(ibas)
 
 
@@ -544 +535 @@
 
       IF ( wape .LT. 0.) THEN
-        print*,'wape<0'
+        if(prt_level.ge.10) print*,'wape<0'
         wape = 0.
         hw = hwmin
@@ -555 +546 @@
         ENDDO
       ELSE
-        print*,'wape>0'
+        if(prt_level.ge.10) print*,'wape>0'
         Cstar = stark*sqrt(2.*wape)
       ENDIF
@@ -1012 +1003 @@
 
       IF ( wape2 .LT. 0.) THEN
-        print*,'wape2<0'
+        if(prt_level.ge.10) print*,'wape2<0'
         wape2 = 0.
         hw = hwmin
@@ -1023 +1014 @@
         ENDDO
       ELSE
-        print*,'wape2>0'
+        if(prt_level.ge.10) print*,'wape2>0'
         Cstar2 = stark*sqrt(2.*wape2)