Changeset 1496 for LMDZ5/trunk/libf/phylmd/physiq.F

Timestamp:

Mar 11, 2011, 10:44:05 AM (14 years ago)

Author:

Laurent Fairhead

Message:

Modified closure for the new physics package, new values for the iflag_coupl parameter
FH

---

Fermeture modifiée pour la nouvelle physique, nouvelles valeurs définies pour
le paramètre iflag_coupl
FH

File:

• LMDZ5/trunk/libf/phylmd/physiq.F (modified) (4 diffs)

LMDZ5/trunk/libf/phylmd/physiq.F

@@ -1234 +1234 @@
      .     flag_aerosol, new_aod,
      .     bl95_b0, bl95_b1,
-     .     iflag_thermals,nsplit_thermals,tau_thermals,
-     .     iflag_thermals_ed,iflag_thermals_optflux,
 c     nv flags pour la convection et les poches froides
-     .     iflag_coupl,iflag_clos,iflag_wake, read_climoz,
+     .     read_climoz,
      &     alp_offset)
       call phys_state_var_init(read_climoz)
@@ -2846 +2844 @@
       if (iflag_cldcon>=5) then
 
+         ptconvth(:,:)=fm_therm(:,:)>0.
+
+       if (iflag_coupl==4) then
+
+! Dans le cas iflag_coupl==4, on prend la somme des convertures
+! convectives et lsc dans la partie des thermiques
+! Le controle par iflag_coupl est peut etre provisoire.
+         do k=1,klev
+            do i=1,klon
+               if (ptconv(i,k).and.ptconvth(i,k)) then
+                   cldliq(i,k)=cldliq(i,k)+rnebcon(i,k)*clwcon(i,k)
+                   cldfra(i,k)=min(cldfra(i,k)+rnebcon(i,k),1.)
+               else if (ptconv(i,k)) then
+                   cldfra(i,k)=rnebcon(i,k)
+                   cldliq(i,k)=rnebcon(i,k)*clwcon(i,k)
+               endif
+            enddo
+         enddo
+
+         else
 ! Si on est sur un point touche par la convection profonde et pas
 ! par les thermiques, on prend la couverture nuageuse et l'eau nuageuse
@@ -2852 +2870 @@
 !IM/FH: 2011/02/23 
 ! definition des points sur lesquels ls thermiques sont actifs
-         if (prt_level>9)write(*,*)'TEST SCHEMA DE NUAGES '
-         ptconvth(:,:)=fm_therm(:,:)>0.
+
          do k=1,klev
             do i=1,klon
@@ -2862 +2879 @@
             enddo
          enddo
+
+        endif
+
       else
+
 ! Ancienne version
       cldfra(:,:)=min(max(cldfra(:,:),rnebcon(:,:)),1.)