Changeset 557 for LMDZ4/trunk/libf/phylmd/physiq.F

Timestamp:

Oct 27, 2004, 12:14:46 PM (20 years ago)

Author:

lmdzadmin

Message:

Initialisations diverses heritees de LMDZ.3.3
LF

File:

• LMDZ4/trunk/libf/phylmd/physiq.F (modified) (5 diffs)

LMDZ4/trunk/libf/phylmd/physiq.F

@@ -34 +34 @@
 c   CLEFS CPP POUR LES IO
 c   =====================
-c #define histhf
+#define histhf
 #define histday
 #define histmth
 #define histins
-c #define histISCCP
-c #define histREGDYN
-c #define histmthNMC
+c#define histISCCP
+#define histREGDYN
+#define histmthNMC
 c======================================================================
 c    modif   ( P. Le Van ,  12/10/98 )
@@ -982 +982 @@
 #include "FCTTRE.h"
 c======================================================================
+!rv
+      u10m(:,:)=0.
+      v10m(:,:)=0.
+      t2m(:,:)=0.
+      q2m(:,:)=0.
+      ffonte(:,:)=0.
+      fqcalving(:,:)=0.
+      piz_ae(:,:,:)=0.
+      tau_ae(:,:,:)=0.
+      cg_ae(:,:,:)=0.
+      rain_con(:)=0.
+      snow_con(:)=0.
+      bl95_b0=0.
+      bl95_b1=0.
+!rv
       modname = 'physiq'
       IF (if_ebil.ge.1) THEN
@@ -1004 +1019 @@
 c======================================================================
 cIM BEG
-        DO k=1, nlev
+        DO k=1, nlevENS
           DO l=1, nlevSTD
 c
@@ -1870 +1885 @@
       entr_therm(:,:)=0.
 c
-      print*,'AVANT LA CONVECTION SECHE , iflag_thermals='
+      IF(prt_level>9)WRITE(lunout,*)
+     .    'AVANT LA CONVECTION SECHE , iflag_thermals='
      s   ,iflag_thermals,'   nsplit_thermals=',nsplit_thermals
       if(iflag_thermals.lt.0) then
 c  Rien
 c  ====
-         print*,'pas de convection'
+         IF(prt_level>9)WRITE(lunout,*)'pas de convection'
       else if(iflag_thermals.eq.0) then
 
 c  Ajustement sec
 c  ==============
-         print*,'ajsec'
+         IF(prt_level>9)WRITE(lunout,*)'ajsec'
          CALL ajsec(paprs, pplay, t_seri,q_seri, d_t_ajs, d_q_ajs)
          t_seri(:,:) = t_seri(:,:) + d_t_ajs(:,:)
@@ -1887 +1903 @@
 c  Thermiques
 c  ==========
-         print*,'JUSTE AVANT , iflag_thermals='
+         IF(prt_level>9)WRITE(lunout,*)'JUSTE AVANT , iflag_thermals='
      s   ,iflag_thermals,'   nsplit_thermals=',nsplit_thermals
          call calltherm(pdtphys