Changeset 2801 for LMDZ5/trunk/libf/phylmd/physiq_mod.F90

Timestamp:

Feb 28, 2017, 12:48:45 PM (7 years ago)

Author:

jyg

Message:

Bug fix in a call to add_phys_tend.

File:

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

LMDZ5/trunk/libf/phylmd/physiq_mod.F90

@@ -821 +821 @@
     !
     ! tendance nulles
-    REAL, dimension(klon,klev):: du0, dv0, dt0, dq0, dql0, dqi0
-    REAL, dimension(klon)     :: dsig0, ddens0
-    INTEGER, dimension(klon)  :: wkoccur1
+    REAL, dimension(klon,klev)  :: du0, dv0, dt0, dq0, dql0, dqi0
+    REAL, dimension(klon)       :: dsig0, ddens0
+    INTEGER, dimension(klon)    :: wkoccur1
+    ! tendance buffer pour appel de add_phys_tend
+    REAL, DIMENSION(klon,klev)  :: d_q_ch4_dtime
     !
     ! Flag pour pouvoir ne pas ajouter les tendances. 
@@ -3988 +3990 @@
        CALL METHOX(1,klon,klon,klev,q_seri,d_q_ch4,pplay)
        ! ajout de la tendance d'humidite due au methane
-       CALL add_phys_tend(du0, dv0, dt0, d_q_ch4*dtime, dql0, dqi0, paprs, &
+       d_q_ch4_dtime(:,:) = d_q_ch4(:,:)*dtime
+       CALL add_phys_tend(du0, dv0, dt0, d_q_ch4_dtime, dql0, dqi0, paprs, &
             'q_ch4', abortphy,flag_inhib_tend,itap)
+       d_q_ch4(:,:) = d_q_ch4_dtime(:,:)/dtime
     ENDIF
     !