Changeset 953 for LMDZ4/trunk

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

Location:

LMDZ4/trunk/libf/phylmd

Files:

• calwake.F (modified) (5 diffs)
• physiq.F (modified) (1 diff)
• wake.F (modified) (10 diffs)

LMDZ4/trunk/libf/phylmd/calwake.F

@@ -1 +1 @@
       SUBROUTINE CALWAKE(paprs,pplay,dtime
-     :             ,t,q,omgb,ibas_con
+     :             ,t,q,omgb
      :             ,dt_dwn,dq_dwn,M_dwn,M_up
      :             ,dt_a,dq_a,sigd
@@ -43 +43 @@
       REAL paprs(klon,klev+1),pplay(klon,klev)
       REAL t(klon,klev), q(klon,klev), omgb(klon,klev)
-      INTEGER ibas_con(klon)
       REAL dt_dwn(klon,klev), dq_dwn(klon,klev),M_dwn(klon,klev)
       REAL M_up(klon,klev)
@@ -73 +72 @@
       REAL p(klev),ph(klev+1),pi(klev)
       REAL te(klev),qe(klev),omgbe(klev),dtdwn(klev),dqdwn(klev)
-      INTEGER  ibas
       REAL dta(klev),dqa(klev)
       REAL wdtPBL(klev),wdqPBL(klev)
@@ -121 +119 @@
 c      print*, 'sigd0,sigd', sigd0, sigd(i)
       ph(klev+1)=0.
-cIM bug indice ibas_con: ibas = ibas_con(1)
-      ibas = ibas_con(i)
 
       ktopw = wake_k(i)
@@ -156 +152 @@
 
       call WAKE(p,ph,pi,dtime,sigd0
-     $                ,te,qe,omgbe,ibas
+     $                ,te,qe,omgbe
      $                ,dtdwn,dqdwn,amdwn,amup,dta,dqa
      $                ,wdtPBL,wdqPBL,udtPBL,udqPBL

LMDZ4/trunk/libf/phylmd/physiq.F

@@ -2207 +2207 @@
 ccalcul caracteristiques de la poche froide
       call calWAKE (paprs,pplay,dtime
-     :               ,t_seri,q_seri,omega,ibas_con
+     :               ,t_seri,q_seri,omega
      :               ,dt_dwn,dq_dwn,M_dwn,M_up
      :               ,dt_a,dq_a,sigd

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)