Changeset 3781 for LMDZ6/trunk/libf/phylmd/dyn1d/scm.F90

Timestamp:

Oct 26, 2020, 10:33:06 AM (4 years ago)

Author:

evignon

Message:

Initialisation de la TKE pour les cas 1D (important pour GABLS1), Etienne

File:

• LMDZ6/trunk/libf/phylmd/dyn1d/scm.F90 (modified) (3 diffs)

LMDZ6/trunk/libf/phylmd/dyn1d/scm.F90

@@ -412 +412 @@
 !          (phys_state_var_init is called again in physiq)
       read_climoz = 0
-!
+      nsw=6
+
       call phys_state_var_init(read_climoz)
 
@@ -624 +625 @@
         zpic = zpicinp
         ftsol=tsurf
-        nsw=6 ! on met le nb de bandes SW=6, pour initialiser
-              ! 6 albedo, mais on peut quand meme tourner avec
-              ! moins. Seules les 2 ou 4 premiers seront lus
-              print*, 'les albedos sont', albedo
         falb_dir=albedo
         falb_dif=albedo
@@ -639 +636 @@
         prw_ancien = 0.
 !jyg<
+! Etienne: comment those lines since now the TKE is inialized in 1D_read_forc_cases
 !!      pbl_tke(:,:,:)=1.e-8
-        pbl_tke(:,:,:)=0.
-! EV: pourquoi????
+!        pbl_tke(:,:,:)=0.
 !        pbl_tke(:,2,:)=1.e-2
-        PRINT *, ' pbl_tke dans lmdz1d '
-        if (prt_level .ge. 5) then 
-         DO nsrf = 1,4
-           PRINT *,'pbl_tke(1,:,',nsrf,') ',pbl_tke(1,:,nsrf)
-         ENDDO
-        end if
-
 !>jyg
         rain_fall=0.