Changeset 1071 for trunk

Timestamp:

Oct 14, 2013, 7:15:25 PM (11 years ago)

Author:

slebonnois

Message:

SL: still some corrections to version 1056

Location:

trunk

Files:

• LMDZ.COMMON/libf/dyn3d/calfis.F (modified) (3 diffs)
• LMDZ.COMMON/libf/dyn3dpar/calfis_p.F (modified) (4 diffs)
• LMDZ.TITAN/libf/phytitan/physiq.F (modified) (1 diff)
• LMDZ.TITAN/libf/phytitan/phytrac.F (modified) (1 diff)

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

@@ -195 +195 @@
 ! to allow for 2D computation of microphys and chemistry
       LOGICAL,save :: flag_moyzon
-      REAL,dimension(iip1,llm) :: tmpvar
-      REAL,dimension(iip1,llmp1) :: tmpvarp1
-      REAL,dimension(llm) :: tmpvarbar
-      REAL,dimension(llmp1) :: tmpvarbarp1
+      REAL,allocatable :: tmpvar(:,:)
+      REAL,allocatable :: tmpvarp1(:,:)
+      REAL,allocatable :: tmpvarbar(:)
+      REAL,allocatable :: tmpvarbarp1(:)
 
 c-----------------------------------------------------------------------
@@ -222 +222 @@
         if(moyzon_ch.or.moyzon_mu) then
          flag_moyzon = .true.
+         allocate(tmpvar(iip1,llm))
+         allocate(tmpvarp1(iip1,llmp1))
+         allocate(tmpvarbar(llm))
+         allocate(tmpvarbarp1(llmp1))
         endif
+
         if (flag_moyzon) call moyzon_init
 
 c----------------------------------------------
 c moyennes globales pour le profil de pression
+      if(planet_type.eq."titan") then
         ALLOCATE(plevmoy(llm+1))
         ALLOCATE(playmoy(llm))
@@ -252 +258 @@
         playmoy = preff * (pkmoy/cpp) ** unskap
         call tpot2t(llm,tetamoy,tmoy,pkmoy)
-c----------------------------------------------
+c-------------------
 c + lat index
-      allocate(klat(ngridmx))
-      klat=0
-      klat(1)  = 1
-      ig0  = 2
-      DO j = 2,jjm
-         do i=0,iim-1
-          klat(ig0+i) = j
-         enddo
-         ig0 = ig0+iim
-      ENDDO
-      klat(ngridmx)  = jjp1
+        allocate(klat(ngridmx))
+        klat=0
+        klat(1)  = 1
+        ig0  = 2
+        DO j = 2,jjm
+           do i=0,iim-1
+            klat(ig0+i) = j
+           enddo
+           ig0 = ig0+iim
+        ENDDO
+        klat(ngridmx)  = jjp1
+      endif   ! planet_type=titan
 c----------------------------------------------
       ELSE

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

@@ -256 +256 @@
 ! to allow for 2D computation of microphys and chemistry
       LOGICAL,save :: flag_moyzon
-      REAL,dimension(iip1,llm) :: tmpvar
-      REAL,dimension(iip1,llmp1) :: tmpvarp1
-      REAL,dimension(llm) :: tmpvarbar
-      REAL,dimension(llmp1) :: tmpvarbarp1
+      REAL,allocatable :: tmpvar(:,:)
+      REAL,allocatable :: tmpvarp1(:,:)
+      REAL,allocatable :: tmpvarbar(:)
+      REAL,allocatable :: tmpvarbarp1(:)
 
 c-----------------------------------------------------------------------
@@ -287 +287 @@
         if(moyzon_ch.or.moyzon_mu) then
          flag_moyzon = .true.
+         allocate(tmpvar(iip1,llm))
+         allocate(tmpvarp1(iip1,llmp1))
+         allocate(tmpvarbar(llm))
+         allocate(tmpvarbarp1(llmp1))
         endif
 
@@ -310 +314 @@
 c------------------------------------------------------------------
 c moyennes globales pour les profils de pression et de temperature
+      if(planet_type.eq."titan") then
         call AllGather_Field(pp,iip1*jjp1,llmp1)
         call AllGather_Field(pteta,iip1*jjp1,llm)
@@ -338 +343 @@
         playmoy = preff * (pkmoy/cpp) ** unskap
         call tpot2t_p(1,llm,tetamoy,tmoy,pkmoy)
-c------------------------------------------------------------------
+c-------------------
 c + lat index
-      allocate(klat(klon))
-      do ig0=1,klon
-        j=index_j(ig0)
-        klat(ig0)=j
-      enddo
+        allocate(klat(klon))
+        do ig0=1,klon
+           j=index_j(ig0)
+           klat(ig0)=j
+        enddo
+      endif   ! planet_type=titan
 c------------------------------------------------------------------
 

trunk/LMDZ.TITAN/libf/phytitan/physiq.F

@@ -408 +408 @@
 c appel a la lecture du physiq.def 
 c
-         call conf_phys(ok_journe, ok_mensuel,
+         call conf_phys(ok_mensuel,ok_journe,
      .                  ok_instan,
      .                  if_ebil)

trunk/LMDZ.TITAN/libf/phytitan/phytrac.F

@@ -471 +471 @@
 c qaer est ce qui entre dans muphy, donc la moyenne zonale
 c qaer0 est la valeur initiale du champ
-c qaer0 est la moyenne zonale initiale
+c qaer0bar est la moyenne zonale initiale
 c la variation relative pour une bande de latitude est donc (qaer/qaer0bar)
 c la nouvelle valeur en un point (3D) est donc qaer0*(qaer/qaer0bar)