Changeset 2565 for LMDZ5/trunk/libf/phylmd/physiq_mod.F90

Timestamp:

Jun 10, 2016, 4:01:28 PM (8 years ago)

Author:

jyg

Message:

Some small improvements to the 1D model:
(1) initialization of pbl_tke (-> 1+1=2);
(2) bug fix concerning nudge_tsoil=y;
(3) possibility of a constant moisture convergence
when forcing_type=0.

File:

• LMDZ5/trunk/libf/phylmd/physiq_mod.F90 (modified) (2 diffs)

LMDZ5/trunk/libf/phylmd/physiq_mod.F90

@@ -1165 +1165 @@
 
        CALL phyetat0 ("startphy.nc",clesphy0,tabcntr0)
+!jyg<
        IF (klon_glo==1) THEN
-          coefh=0. ; coefm=0. ; pbl_tke=0.
-          coefh(:,2,:)=1.e-2 ; coefm(:,2,:)=1.e-2 ; pbl_tke(:,2,:)=1.e-2
-          PRINT*,'FH WARNING : lignes a supprimer'
+          pbl_tke(:,:,is_ave) = 0.
+          DO nsrf=1,nbsrf
+            DO k = 1,klev+1
+                 pbl_tke(:,k,is_ave) = pbl_tke(:,k,is_ave) &
+                     +pctsrf(:,nsrf)*pbl_tke(:,k,nsrf)
+            ENDDO
+          ENDDO
+!>jyg
        ENDIF
        !IM begin
@@ -2776 +2782 @@
 !jyg<
 !
-          CALL alpale_th( dtime, lmax_th, t_seri,  &
+          CALL alpale_th( dtime, lmax_th, t_seri, cell_area,  &
                           cin, s2, n2,  &
                           ale_bl_trig, ale_bl_stat, ale_bl,  &