Changeset 102 for trunk/libf/phytitan/ajsec.F

Timestamp:

Mar 24, 2011, 11:46:27 AM (14 years ago)

Author:

slebonnois

Message:

SL : corrections et modifications dans phytitan correspondant a celles
faites apres compilation Venus. Titan pas encore compile.

File:

• trunk/libf/phytitan/ajsec.F (modified) (10 diffs)

trunk/libf/phytitan/ajsec.F

@@ -3 +3 @@
 !
 ! ADAPTATION GCM POUR CP(T)
-      SUBROUTINE ajsec(paprs, pplay, ppk, tfi, ufi, vfi, nq, qfi,
+      SUBROUTINE ajsec(paprs, pplay, ppk, tfi, ufi, vfi, nq, qfi, 
      .                             d_tfi, d_ufi, d_vfi, d_qfi)
+
+      use dimphy
       IMPLICIT none
 c======================================================================
@@ -23 +25 @@
 c d_vfi-----output-R-Incrementation du vent meridien
 c d_qfi-----output-R-Incrementation des traceurs
-
 c======================================================================
 #include "dimensions.h"
-#include "dimphy.h"
 #include "YOMCST.h"
       REAL paprs(klon,klev+1), pplay(klon,klev)
@@ -36 +36 @@
       REAL qfi(klon,klev,nq), d_qfi(klon,klev,nq)
 c
-      INTEGER limbas, limhau ! les couches a ajuster
-ccc      PARAMETER (limbas=klev-3, limhau=klev)
-      PARAMETER (limbas=1, limhau=klev)
+      INTEGER,save :: limbas, limhau ! les couches a ajuster
 c
       REAL zh(klon,klev)
@@ -51 +49 @@
 c Initialisation:
 c
+      limbas=1
+      limhau=klev
+
       DO k = 1, klev
       DO i = 1, klon
@@ -56 +57 @@
          d_ufi(i,k) = 0.0
          d_vfi(i,k) = 0.0
-         d_qfi(i,k,:) = 0.0
+         d_qfi(i,k,:) = 0.0
          zu(i,k)    = ufi(i,k)
          zv(i,k)    = vfi(i,k)
-         zq(i,k,:)  = qfi(i,k,:)
+         zq(i,k,:)  = qfi(i,k,:)
       ENDDO
       ENDDO
@@ -119 +120 @@
               zum=0.
               zvm=0.
-              zqm=0.
+              zqm=0.
               DO k = k1, k2
                 zalpha=zalpha+ABS(zh(i,k)-hm)*zdp(i,k)
@@ -125 +126 @@
                 zum=zum+ufi(i,k)*zdp(i,k)
                 zvm=zvm+vfi(i,k)*zdp(i,k)
-                do iq=1,nq
-                   zqm(iq)=zqm(iq)+qfi(i,k,iq)*zdp(i,k)
-                enddo
+                do iq=1,nq
+                  zqm(iq)=zqm(iq)+qfi(i,k,iq)*zdp(i,k)
+                enddo
               ENDDO
               zalpha=zalpha/(hm*(paprs(i,k1)-paprs(i,k2+1)))
               zum=zum/(paprs(i,k1)-paprs(i,k2+1))
               zvm=zvm/(paprs(i,k1)-paprs(i,k2+1))
-              do iq=1,nq
-                 zqm(iq)=zqm(iq)/(paprs(i,k1)-paprs(i,k2+1))
-              enddo
+              do iq=1,nq
+                zqm(iq)=zqm(iq)/(paprs(i,k1)-paprs(i,k2+1))
+              enddo
 
               IF(zalpha.GT.1.) THEN
@@ -153 +154 @@
                  zv(i,k)=vfi(i,k)+zalpha*(zvm-vfi(i,k))
                  do iq=1,nq
-                   zq(i,k,iq)=qfi(i,k,iq)+zalpha*(zqm(iq)-qfi(i,k,iq))
+                   zq(i,k,iq)=qfi(i,k,iq)+zalpha*(zqm(iq)-qfi(i,k,iq))
                  enddo
               ENDDO
@@ -187 +188 @@
          do iq=1,nq
            d_qfi(i,k,iq) = 0.0
-         enddo 
+         enddo
       ENDDO
       ENDDO
@@ -200 +201 @@
          do iq=1,nq
            d_qfi(i,k,iq) = 0.0
-         enddo 
+         enddo
       ENDDO
       ENDDO