Changeset 1726 for trunk/LMDZ.VENUS/libf/phyvenus/phytrac_chimie.F

Timestamp:

Jul 22, 2017, 11:34:22 PM (7 years ago)

Author:

slebonnois

Message:

SL: adjustments after revision 1723, + some debug for cloud microphysics config and rcm1d

File:

• trunk/LMDZ.VENUS/libf/phyvenus/phytrac_chimie.F (modified) (2 diffs)

trunk/LMDZ.VENUS/libf/phyvenus/phytrac_chimie.F

@@ -273 +273 @@
          DO ilev=1, n_lev
 
+         if (temp(ilon,ilev).lt.500.) then
            CALL MAD_MUPHY(pdtphys,                               ! Timestep
      &  temp(ilon,ilev),pplay(ilon,ilev),                        ! Temperature and pressure
@@ -283 +284 @@
      &  trac(ilon,ilev,i_m3_mode2sa),trac(ilon,ilev,i_m3_mode2w),      ! Moments of mode 2
      &  trac(ilon,ilev,i_m3_mode2ccn))                                 ! Moments of mode 2
-
+          else
+           trac(ilon,ilev,i_m0_aer)=0. 
+           trac(ilon,ilev,i_m3_aer)=0.
+           trac(ilon,ilev,i_m0_mode1drop)=0.
+           trac(ilon,ilev,i_m0_mode1ccn)=0.
+           trac(ilon,ilev,i_m3_mode1sa)=0.
+           trac(ilon,ilev,i_m3_mode1w)=0.
+           trac(ilon,ilev,i_m3_mode1ccn)=0.
+           trac(ilon,ilev,i_m0_mode2drop)=0.
+           trac(ilon,ilev,i_m0_mode2ccn)=0.
+           trac(ilon,ilev,i_m3_mode2sa)=0.
+           trac(ilon,ilev,i_m3_mode2w)=0.
+           trac(ilon,ilev,i_m3_mode2ccn)=0.
+          endif
          ENDDO
         ENDDO