Changeset 1650 for LMDZ5/trunk/libf/phylmd/cv3a_uncompress.F

Timestamp:

Sep 6, 2012, 3:25:45 PM (12 years ago)

Author:

Laurent Fairhead

Message:

Inclusion de modifications pour régler le problème convection/traceurs dans la nouvelle
physique

1. Cozic

---

Modifications needed to correct the convection/tracers problem with the new physics

1. Cozic

File:

• LMDZ5/trunk/libf/phylmd/cv3a_uncompress.F (modified) (6 diffs)

LMDZ5/trunk/libf/phylmd/cv3a_uncompress.F

@@ -9 +9 @@
      :         ,Plim1,Plim2,asupmax,supmax0
      :         ,asupmaxmin
+!AC!
+     :         ,da,phi
+!AC!
      o         ,iflag1,kbas1,ktop1
      :         ,precip1,cbmf1,plcl1,plfc1,wbeff1,sig1,w01,ptop21
@@ -17 +20 @@
      :         ,ftd1,fqd1
      :         ,Plim11,Plim21,asupmax1,supmax01
-     :         ,asupmaxmin1     )
+     :         ,asupmaxmin1     
+!AC!
+     :         ,da1,phi1  )
+!AC!
 ***************************************************************
 *                                                             *
@@ -50 +56 @@
       real asupmax(nloc,nd),supmax0(nloc)
       real asupmaxmin(nloc)
-
+!AC!
+      real da(nloc,nd),phi(nloc,nd,nd)
+!AC!
 c outputs:
       integer iflag1(len),kbas1(len),ktop1(len)
@@ -68 +76 @@
       real asupmax1(len,nd),supmax01(len)
       real asupmaxmin1(len)
+!AC!
+      real da1(nloc,nd),phi1(nloc,nd,nd)
+!AC!
 c
 c local variables:
@@ -111 +122 @@
             fqd1(idcum(i),k)=fqd(i,k)
             asupmax1(idcum(i),k)=asupmax(i,k)
- 2010     continue
+!AC!
+            da1(idcum(i),k)=da(i,k)
+!AC!
+ 2010    continue
  2020   continue
 
@@ -119 +133 @@
 
 
-        do 2100 j=1,ntra
-c oct3         do 2110 k=1,nl
-         do 2110 k=1,nd ! oct3
-          do 2120 i=1,ncum
-            ftra1(idcum(i),k,j)=ftra(i,k,j)
- 2120     continue
- 2110    continue
- 2100   continue
+!AC!        do 2100 j=1,ntra
+!AC!c oct3         do 2110 k=1,nl
+!AC!         do 2110 k=1,nd ! oct3
+!AC!          do 2120 i=1,ncum
+!AC!            ftra1(idcum(i),k,j)=ftra(i,k,j)
+!AC! 2120     continue
+!AC! 2110    continue
+!AC! 2100   continue
+
+!AC!
+       do j=1,nd
+         do k=1,nd 
+          do i=1,ncum
+            phi1(idcum(i),k,j)=phi(i,k,j)
+          end do
+         end do
+        end do
+!AC!
+
 c
 c        do 2220 k2=1,nd