Changeset 616 for LMDZ4/trunk/libf/dyn3d

Timestamp:

Apr 13, 2005, 10:58:34 AM (20 years ago)

Author:

lmdzadmin

Message:

Mise a jour pour INCA.2.0 Anne C
MAFi+LF

Location:

LMDZ4/trunk/libf/dyn3d

Files:

• advtrac.F (modified) (6 diffs)
• caladvtrac.F (modified) (3 diffs)
• calfis.F (modified) (4 diffs)
• iniadvtrac.F (modified) (2 diffs)
• iniinterp_horiz.F (added)
• interp_horiz.F (added)
• leapfrog.F (modified) (3 diffs)

LMDZ4/trunk/libf/dyn3d/advtrac.F

@@ -4 +4 @@
 c
 c
-#ifdef INCA_CH4
+#ifdef INCA
       SUBROUTINE advtrac(pbaru,pbarv ,
      *                   p,  masse,q,iapptrac,teta,
@@ -16 +16 @@
      *                  pk)
 #endif
-
 c     Auteur :  F. Hourdin
 c
@@ -51 +50 @@
       REAL p( ip1jmp1,llmp1 ),teta(ip1jmp1,llm)
       REAL pk(ip1jmp1,llm)
-#ifdef INCA_CH4
+#ifdef INCA
       INTEGER            :: hadv_flg(nqmx)
-      REAL               :: mmt_adj(ip1jmp1,llm)
+      REAL               :: mmt_adj(ip1jmp1,llm,1)
       REAL               :: flxw(ip1jmp1,llm)
 #endif
@@ -129 +128 @@
         CALL groupe( massem, pbaruc,pbarvc, pbarug,pbarvg,wg )
 
-#ifdef INCA_CH4
+#ifdef INCA
       ! ... Flux de masse diaganostiques traceurs
       flxw = wg / FLOAT(iapp_tracvl)
@@ -214 +213 @@
         else if (iadv(iq).eq.20) then
             call pentes_ini (q(1,1,iq),wg,massem,pbarug,pbarvg,0)
-#ifdef INCA_CH4
+#ifdef INCA
        do iiq = iq+1, iq+3
          q(:,:,iiq)=q(:,:,iiq)*mmt_adj(:,:,1)
@@ -232 +231 @@
            call  prather(q(1,1,iq),wg,massem,pbarug,pbarvg,
      s                     n,dtbon)
-#ifdef INCA_CH4
+#ifdef INCA
        do iiq = iq+1, iq+9
          q(:,:,iiq)=q(:,:,iiq)*mmt_adj(:,:,1)

LMDZ4/trunk/libf/dyn3d/caladvtrac.F

@@ -4 +4 @@
 c
 c
-#ifdef INCA_CH4
+#ifdef INCA
             SUBROUTINE caladvtrac(q,pbaru,pbarv ,
      *                   p ,masse, dq ,  teta,
@@ -41 +41 @@
       REAL p( ip1jmp1,llmp1),q( ip1jmp1,llm,nqmx),dq( ip1jmp1,llm,2 )
       REAL teta( ip1jmp1,llm),pk( ip1jmp1,llm)
-#ifdef INCA_CH4
+#ifdef INCA
       INTEGER            :: hadv_flg(nqmx)
       REAL               :: mmt_adj(iip1,jjp1,llm,1)
@@ -73 +73 @@
 c   advection
 
-#ifdef INCA_CH4
+#ifdef INCA
       CALL advtrac( pbaru,pbarv, 
      *             p,  masse,q,iapptrac, teta,

LMDZ4/trunk/libf/dyn3d/calfis.F

@@ -23 +23 @@
      $                  pdq,
      $                  pw,
-#ifdef INCA_CH4
+#ifdef INCA
      $                  flxw,
 #endif
@@ -163 +163 @@
       REAL unskap, pksurcp
 
-#ifdef INCA_CH4
+#ifdef INCA
       REAL flxw(iip1,jjp1,llm)
       REAL flxwfi(ngridmx,llm)
@@ -432 +432 @@
 
 
-#ifdef INCA_CH4
+#ifdef INCA
       CALL gr_dyn_fi(llm,iip1,jjp1,ngridmx,flxw,flxwfi)
 #endif
@@ -461 +461 @@
      .             zqfi,
      .             pvervel,
-#ifdef INCA_CH4
+#ifdef INCA
      .             flxwfi,
 #endif

LMDZ4/trunk/libf/dyn3d/iniadvtrac.F

@@ -8 +8 @@
 #ifdef INCA
       USE transport_controls, only : hadv_flg, vadv_flg
+      USE species_names
       USE chemshut
 #endif
@@ -192 +193 @@
          str1=tnom(iq)
          tname(iiq)=tnom(iq)
-         ttext(iiq)=str1(1:lnblnk(str1))//descrq(iadv(iiq))
+         IF (iadv(iiq).eq.0) THEN
+           ttext(iiq)=str1(1:lnblnk(str1))
+         ELSE
+           ttext(iiq)=str1(1:lnblnk(str1))//descrq(iadv(iiq))
+         endif
          str2=ttext(iiq)
 c   schemas tenant compte des moments d'ordre superieur.

LMDZ4/trunk/libf/dyn3d/leapfrog.F

@@ -132 +132 @@
 
       REAL alpha(ip1jmp1,llm),beta(ip1jmp1,llm)
-#ifdef INCA_CH4
+#ifdef INCA
       REAL :: flxw(ip1jmp1,llm)
 #endif
@@ -274 +274 @@
 
 c
-#ifdef INCA_CH4
+#ifdef INCA
              CALL caladvtrac(q,pbaru,pbarv,
      *                      p, masse, dq,  teta,
@@ -358 +358 @@
      $               ucov,vcov,teta,q,masse,ps,p,pk,phis,phi ,
      $               du,dv,dteta,dq,w,
-#ifdef INCA_CH4
+#ifdef INCA
      $               flxw,
 #endif