Changeset 1126 for trunk/LMDZ.COMMON/libf

Timestamp:

Dec 17, 2013, 1:02:44 PM (11 years ago)

Author:

slebonnois

Message:

SL: update of Titan photochemical module to include computation of chemistry up to 1300 km

Location:

trunk/LMDZ.COMMON/libf

Files:

• dyn3d/calfis.F (modified) (6 diffs)
• dyn3d/moyzon_mod.F90 (modified) (1 diff)
• dyn3dpar/calfis_p.F (modified) (5 diffs)
• dyn3dpar/moyzon_mod.F90 (modified) (1 diff)

trunk/LMDZ.COMMON/libf/dyn3d/calfis.F

@@ -199 +199 @@
       REAL,allocatable,save :: tmpvarbar(:)
       REAL,allocatable,save :: tmpvarbarp1(:)
+      real :: zz1,zz2
 
 c-----------------------------------------------------------------------
@@ -232 +233 @@
 c----------------------------------------------
 c moyennes globales pour le profil de pression
-      if(planet_type.eq."titan") then
+       if(planet_type.eq."titan".or.planet_type.eq."venus") then
         ALLOCATE(plevmoy(llm+1))
         ALLOCATE(playmoy(llm))
@@ -238 +239 @@
         ALLOCATE(tetamoy(llm))
         ALLOCATE(pkmoy(llm))
+        ALLOCATE(phimoy(0:llm))
+        ALLOCATE(zlevmoy(llm+1))
+        ALLOCATE(zlaymoy(llm))
         plevmoy=0.
         do l=1,llmp1
@@ -248 +252 @@
         tetamoy=0.
         pkmoy=0.
+        phimoy=0.
+        do i=1,iip1
+         do j=1,jjp1
+            phimoy(0)=phimoy(0)+pphis(i,j)/(iip1*jjp1)
+         enddo
+        enddo
         do l=1,llm
          do i=1,iip1
@@ -253 +263 @@
             tetamoy(l)=tetamoy(l)+pteta(i,j,l)/(iip1*jjp1)
             pkmoy(l)=pkmoy(l)+ppk(i,j,l)/(iip1*jjp1)
+            phimoy(l)=phimoy(l)+pphi(i,j,l)/(iip1*jjp1)
           enddo
          enddo
         enddo
-        playmoy = preff * (pkmoy/cpp) ** unskap
+        playmoy(:) = preff * (pkmoy(:)/cpp) ** unskap
         call tpot2t(llm,tetamoy,tmoy,pkmoy)
+c SI ON TIENT COMPTE DE LA VARIATION DE G AVEC L'ALTITUDE:
+        zlaymoy(:) = g*rad*rad/(g*rad-phimoy(:))-rad
+        zlevmoy(1) = phimoy(0)/g
+        DO l=2,llm
+            zz1=(playmoy(l-1)+plevmoy(l))/(playmoy(l-1)-plevmoy(l))
+            zz2=(plevmoy(l)  +playmoy(l))/(plevmoy(l)  -playmoy(l))
+            zlevmoy(l)=(zz1*zlaymoy(l-1)+zz2*zlaymoy(l))/(zz1+zz2)
+        ENDDO
+        zlevmoy(llmp1)=zlaymoy(llm)+(zlaymoy(llm)-zlevmoy(llm))
 c-------------------
 c + lat index
@@ -265 +285 @@
         ig0  = 2
         DO j = 2,jjm
-           do i=0,iim-1
-            klat(ig0+i) = j
-           enddo
-           ig0 = ig0+iim
+         do i=0,iim-1
+          klat(ig0+i) = j
+         enddo
+         ig0 = ig0+iim
         ENDDO
         klat(ngridmx)  = jjp1
-      endif   ! planet_type=titan
+       endif   ! planet_type=titan
 c----------------------------------------------
       ELSE

trunk/LMDZ.COMMON/libf/dyn3d/moyzon_mod.F90

@@ -18 +18 @@
 ! to be changed...
       REAL,ALLOCATABLE,SAVE :: plevmoy(:),playmoy(:)
+      REAL,ALLOCATABLE,SAVE :: zlevmoy(:),zlaymoy(:),phimoy(:)
       REAL,ALLOCATABLE,SAVE :: tmoy(:),tetamoy(:),pkmoy(:)
       INTEGER,ALLOCATABLE,SAVE :: klat(:)

trunk/LMDZ.COMMON/libf/dyn3dpar/calfis_p.F

@@ -260 +260 @@
       REAL,allocatable,save :: tmpvarbar(:)
       REAL,allocatable,save :: tmpvarbarp1(:)
+      real :: zz1,zz2
 
 c-----------------------------------------------------------------------
@@ -314 +315 @@
 c------------------------------------------------------------------
 c moyennes globales pour les profils de pression et de temperature
-      if(planet_type.eq."titan") then
+      if(planet_type.eq."titan".or.planet_type.eq."venus") then
         call AllGather_Field(pp,iip1*jjp1,llmp1)
         call AllGather_Field(pteta,iip1*jjp1,llm)
         call AllGather_Field(ppk,iip1*jjp1,llm)
+        call AllGather_Field(pphi,iip1*jjp1,llm)
+        call AllGather_Field(pphis,iip1*jjp1,1)
         ALLOCATE(plevmoy(llm+1))
         ALLOCATE(playmoy(llm))
@@ -323 +326 @@
         ALLOCATE(tetamoy(llm))
         ALLOCATE(pkmoy(llm))
+        ALLOCATE(phimoy(0:llm))
+        ALLOCATE(zlevmoy(llm+1))
+        ALLOCATE(zlaymoy(llm))
         plevmoy=0.
         do l=1,llmp1
@@ -333 +339 @@
         tetamoy=0.
         pkmoy=0.
+        phimoy=0.
+        do i=1,iip1
+         do j=1,jjp1
+            phimoy(0)=phimoy(0)+pphis(i,j)/(iip1*jjp1)
+         enddo
+        enddo
         do l=1,llm
          do i=1,iip1
@@ -338 +350 @@
             tetamoy(l)=tetamoy(l)+pteta(i,j,l)/(iip1*jjp1)
             pkmoy(l)=pkmoy(l)+ppk(i,j,l)/(iip1*jjp1)
+            phimoy(l)=phimoy(l)+pphi(i,j,l)/(iip1*jjp1)
           enddo
          enddo
         enddo
-        playmoy = preff * (pkmoy/cpp) ** unskap
+        playmoy(:) = preff * (pkmoy(:)/cpp) ** unskap
         call tpot2t_p(1,llm,tetamoy,tmoy,pkmoy)
+c SI ON TIENT COMPTE DE LA VARIATION DE G AVEC L'ALTITUDE:
+        zlaymoy(:) = g*rad*rad/(g*rad-phimoy(:))-rad
+        zlevmoy(1) = phimoy(0)/g
+        DO l=2,llm
+            zz1=(playmoy(l-1)+plevmoy(l))/(playmoy(l-1)-plevmoy(l))
+            zz2=(plevmoy(l)  +playmoy(l))/(plevmoy(l)  -playmoy(l))
+            zlevmoy(l)=(zz1*zlaymoy(l-1)+zz2*zlaymoy(l))/(zz1+zz2)
+        ENDDO
+        zlevmoy(llmp1)=zlaymoy(llm)+(zlaymoy(llm)-zlevmoy(llm))
 c-------------------
 c + lat index
         allocate(klat(klon))
         do ig0=1,klon
-           j=index_j(ig0)
-           klat(ig0)=j
+          j=index_j(ig0)
+          klat(ig0)=j
         enddo
       endif   ! planet_type=titan

trunk/LMDZ.COMMON/libf/dyn3dpar/moyzon_mod.F90

@@ -24 +24 @@
 ! to be changed...
       REAL,ALLOCATABLE,SAVE :: plevmoy(:),playmoy(:)
+      REAL,ALLOCATABLE,SAVE :: zlevmoy(:),zlaymoy(:),phimoy(:)
       REAL,ALLOCATABLE,SAVE :: tmoy(:),tetamoy(:),pkmoy(:)
       INTEGER,ALLOCATABLE,SAVE :: klat(:)