Changeset 2262 for LMDZ5/trunk/libf/phylmd/physiq.F90

Timestamp:

Apr 17, 2015, 4:34:43 PM (9 years ago)

Author:

jyg

Message:

Correction of a bug concerning the number 'nqo' of
water phases transported by the dynamic : the
default value (= 2, corresponding to vapour and
liquid phases) was still explicitely present in
various places.

Modified files:

infotrac.F90,
physiq.F90

File:

• LMDZ5/trunk/libf/phylmd/physiq.F90 (modified) (2 diffs)

LMDZ5/trunk/libf/phylmd/physiq.F90

@@ -1557 +1557 @@
      ENDDO
 !!! RomP >>>   td dyn traceur
-     IF (nqtot.GE.3) THEN
-        DO iq = 3, nqtot
+!!     IF (nqtot.GE.3) THEN       ! jyg
+!!        DO iq = 3, nqtot        ! jyg
+     IF (nqtot.GE.nqo+1) THEN     ! jyg
+        DO iq = nqo+1, nqtot      ! jyg
            DO k = 1, klev
               DO i = 1, klon
-                 d_tr_dyn(i,k,iq-2)= &
-                      (tr_seri(i,k,iq-2)-tr_ancien(i,k,iq-2))/dtime
+!!                 d_tr_dyn(i,k,iq-2)= &                                 ! jyg
+!!                      (tr_seri(i,k,iq-2)-tr_ancien(i,k,iq-2))/dtime    ! jyg
+                 d_tr_dyn(i,k,iq-nqo)= &                                 ! jyg
+                      (tr_seri(i,k,iq-nqo)-tr_ancien(i,k,iq-nqo))/dtime  ! jyg
                  !         iiq=niadv(iq)
-                 !         print*,i,k," d_tr_dyn",d_tr_dyn(i,k,iq-2),"tra:",iq,tname(iiq)
+                 !         print*,i,k," d_tr_dyn",d_tr_dyn(i,k,iq-nqo),"tra:",iq,tname(iiq)
               ENDDO
            ENDDO
@@ -1580 +1584 @@
      ENDDO
 !!! RomP >>>   td dyn traceur
-     IF (nqtot.GE.3) THEN
-        DO iq = 3, nqtot
+!!     IF (nqtot.GE.3) THEN                                            ! jyg
+!!        DO iq = 3, nqtot                                             ! jyg
+     IF (nqtot.GE.nqo+1) THEN                                          ! jyg
+        DO iq = nqo+1, nqtot                                           ! jyg
            DO k = 1, klev
               DO i = 1, klon
-                 d_tr_dyn(i,k,iq-2)= 0.0
+!!                 d_tr_dyn(i,k,iq-2)= 0.0                             ! jyg
+                 d_tr_dyn(i,k,iq-nqo)= 0.0                             ! jyg
               ENDDO
            ENDDO