Changeset 1909

Timestamp:

Nov 27, 2013, 11:13:29 AM (11 years ago)

Author:

acozic

Message:

Fix a bug in calcul of tr_seri for the injection of source in the first model layer

File:

• LMDZ5/trunk/libf/phylmd/phytrac_mod.F90 (modified) (1 diff)

LMDZ5/trunk/libf/phylmd/phytrac_mod.F90

@@ -590 +590 @@
     DO it=1,nbtr
        d_tr_cl(:,1,it)=source(:,it)*rg/delp(:,1)*pdtphys
+       tr_seri(:,1,it)=tr_seri(:,1,it)+d_tr_cl(:,1,it)
     ENDDO
-    tr_seri(:,1,it)=tr_seri(:,1,it)+d_tr_cl(:,1,it)
     d_tr_cl(:,2:klev,1:nbtr)=0.
   ELSE IF (iflag_vdf_trac==-1) THEN