Changeset 1454 for LMDZ5/trunk/libf/dyn3d/addfi.F

Timestamp:

Nov 18, 2010, 1:01:24 PM (14 years ago)

Author:

Laurent Fairhead

Message:

Merge of LMDZ5V1.0-dev branch r1453 into LMDZ5 trunk r1434

---

Fusion entre la version r1453 de la branche de développement LMDZ5V1.0-dev
et le tronc LMDZ5 (r1434)

Location:

LMDZ5/trunk

Files:

• . (modified) (1 prop)
• libf/dyn3d/addfi.F (modified) (5 diffs)

LMDZ5/trunk/libf/dyn3d/addfi.F

@@ -1 +1 @@
 !
-! $Header$
+! $Id$
 !
       SUBROUTINE addfi(pdt, leapf, forward,
@@ -7 +7 @@
 
       USE infotrac, ONLY : nqtot
+      USE control_mod, ONLY : planet_type
       IMPLICIT NONE
 c
@@ -116 +117 @@
       ENDDO
  
-      DO iq = 1, 2
+      if (planet_type=="earth") then
+      ! earth case, special treatment for first 2 tracers (water)
+       DO iq = 1, 2
          DO k = 1,llm
             DO j = 1,ip1jmp1
@@ -123 +126 @@
             ENDDO
          ENDDO
-      ENDDO
+       ENDDO
 
-      DO iq = 3, nqtot
+       DO iq = 3, nqtot
          DO k = 1,llm
             DO j = 1,ip1jmp1
@@ -132 +135 @@
             ENDDO
          ENDDO
-      ENDDO
+       ENDDO
+      else
+      ! general case, treat all tracers equally)
+       DO iq = 1, nqtot
+         DO k = 1,llm
+            DO j = 1,ip1jmp1
+               pq(j,k,iq)= pq(j,k,iq) + pdqfi(j,k,iq) * pdt
+               pq(j,k,iq)= AMAX1( pq(j,k,iq), qtestt )
+            ENDDO
+         ENDDO
+       ENDDO
+      endif ! of if (planet_type=="earth")