Changeset 1026 for LMDZ4/trunk/libf/phylmd/thermcell_init.F90

Timestamp:

Oct 24, 2008, 5:09:39 PM (16 years ago)

Author:

lmdzadmin

Message:

Modifs thermiques
FH/IM

File:

• LMDZ4/trunk/libf/phylmd/thermcell_init.F90 (modified) (3 diffs)

LMDZ4/trunk/libf/phylmd/thermcell_init.F90

@@ -7 +7 @@
       IMPLICIT NONE
 #include "iniprint.h"
+#include "thermcell.h"
 
       INTEGER l,ig
@@ -34 +35 @@
          lalim(ig)=1
       enddo
+
+      if (iflag_thermals_ed.ge.1) then
+!si la première couche est instable, on declenche un thermique
+         do ig=1,ngrid
+            if (ztv(ig,1).gt.ztv(ig,2)) then
+               lmin(ig)=1
+               lalim(ig)=2
+               alim_star(ig,1)=1.
+               alim_star_tot(ig)=alim_star(ig,1)
+               print*,'init',alim_star(ig,1),alim_star_tot(ig)
+            else
+                lmin(ig)=1
+                lalim(ig)=1
+                alim_star(ig,1)=0.
+                alim_star_tot(ig)=0. 
+            endif
+         enddo
+ 
+         else
+!else iflag_thermals_ed=0 ancienne def de l alim 
 
 !on ne considere que les premieres couches instables
@@ -131 +152 @@
       enddo
        
+!on remet alim_star_tot a 1
+      do ig=1,ngrid 
+         alim_star_tot(ig)=1.
+      enddo
+
+      endif
+!endif iflag_thermals_ed
       return 
       end