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

Timestamp:

Jun 19, 2008, 12:25:57 PM (16 years ago)

Author:

lmdzadmin

Message:

Initialisations : concvl, cv3_routines, cva_driver, physiq
Correction bug i0 + ajout tests : cv3p1_closure
Ajout sorties : ale, alp, cin, wape
Ajout variables wake : phyetat0, phyredem
IM

File:

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

LMDZ4/trunk/libf/phylmd/physiq.F

@@ -117 +117 @@
       PARAMETER (ok_stratus=.FALSE.)
 c======================================================================
-      LOGICAL :: rnpb=.TRUE.
+      LOGICAL, SAVE :: rnpb=.TRUE.
 cIM "slab" ocean
       REAL tslab(klon)    !Temperature du slab-ocean
@@ -1135 +1135 @@
       call phys_state_var_init
       print*, '================================================='
-                  
-        paire_ter(:)=0.    
-        clwcon(:,:)=0.
-        rnebcon(:,:)=0.
-        ratqs(:,:)=0.
-        sollw(:)=0.
+
+cIM beg
+          dnwd0=0.0
+          ftd=0.0
+          fqd=0.0
+          cin=0.
 cym Attention pbase pas initialise dans concvl !!!!
-        pbase(:)=0
-        
+          pbase=0
+          paire_ter(:)=0.    
+cIM 180608
+c         pmflxr=0.
+c         pmflxs=0.
         first=.false.
 
@@ -1229 +1232 @@
      .                  ok_ade, ok_aie, aerosol_couple, 
      .                  bl95_b0, bl95_b1,
-     .                  iflag_thermals,nsplit_thermals,
+     .                  iflag_thermals,nsplit_thermals,tau_thermals,
 cnv flags pour la convection et les poches froides
      .                   iflag_coupl,iflag_clos,iflag_wake)
@@ -1262 +1265 @@
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
 
-
          CALL phyetat0 ("startphy.nc",ocean, ok_veget,clesphy0,tabcntr0)
-
+cIM begin
+          print*,'physiq: clwcon rnebcon ratqs',clwcon(1,1),rnebcon(1,1)
+     $,ratqs(1,1)
+cIM end
 
 
@@ -1358 +1363 @@
 
         do i = 1,klon
-         wake_s(i) = 0.
-         wake_fip(i) = 0.
-         wake_cstar(i) = 0.
-         DO k=1,klev
-          wake_deltat(i,k)=0.
-          wake_deltaq(i,k)=0.
-         ENDDO
+         Ale_bl(i)=0.
+         Alp_bl(i)=0.
         enddo
+
 c================================================================================
 
@@ -2032 +2033 @@
      .        ftd,fqd,lalim_conv,wght_th)
 
+cIM begin
+        print*,'physiq: cin pbase dnwd0 ftd fqd ',cin(1),pbase(1),
+     .dnwd0(1,1),ftd(1,1),fqd(1,1)
+cIM end
 cIM cf. FH
               clwcon0=qcondc
@@ -2184 +2189 @@
         DO i=1,klon
           dt_dwn(i,k)  = ftd(i,k) 
-          wdt_PBL(i,k) = 0.
+          wdt_PBL(i,k) = 0.
           dq_dwn(i,k)  = fqd(i,k) 
-          wdq_PBL(i,k) = 0.
+          wdq_PBL(i,k) = 0.
           M_dwn(i,k)   = dnwd0(i,k)
           M_up(i,k)    = upwd(i,k)
           dt_a(i,k)    = d_t_con(i,k)/dtime - ftd(i,k) 
-          udt_PBL(i,k) = 0.
+          udt_PBL(i,k) = 0.
           dq_a(i,k)    = d_q_con(i,k)/dtime - fqd(i,k)
-          udq_PBL(i,k) = 0.
+          udq_PBL(i,k) = 0.
         ENDDO
       ENDDO
@@ -2244 +2249 @@
       clwcon0th(:,:)=0.
 c
+      fm_therm(:,:)=0.
+      entr_therm(:,:)=0.
+      detr_therm(:,:)=0.
+c
       IF(prt_level>9)WRITE(lunout,*)
      .    'AVANT LA CONVECTION SECHE , iflag_thermals='
@@ -2266 +2275 @@
      s      ,u_seri,v_seri,t_seri,q_seri,zqsat,debut
      s      ,d_u_ajs,d_v_ajs,d_t_ajs,d_q_ajs
-     s      ,fm_therm,entr_therm,zqasc,clwcon0th,lmax_th,ratqscth
+     s      ,fm_therm,entr_therm,detr_therm
+     s      ,zqasc,clwcon0th,lmax_th,ratqscth
      s      ,ratqsdiff,zqsatth
 con rajoute ale et alp, et les caracteristiques de la couche alim