Changeset 1627 for trunk/LMDZ.MARS/libf/phymars

Timestamp:

Nov 16, 2016, 11:24:01 AM (8 years ago)

Author:

emillour

Message:

Mars GCM:
Bug fix in 1D in physiq, missing "if (co2clouds)" around some operations to
be done only when that option is set.
EM

File:

• trunk/LMDZ.MARS/libf/phymars/physiq_mod.F (modified) (3 diffs)

trunk/LMDZ.MARS/libf/phymars/physiq_mod.F

@@ -2812 +2812 @@
            endif ! of if (scavenging)
 
-          if (scavenging) then
+          if (co2clouds) then
+           if (scavenging) then
              Nccntot(:)= 0
              Mccntot(:)= 0
@@ -2842 +2843 @@
                 if (icetotco2(ig)*1e3.lt.0.01) raveco2(ig)=0.
              enddo
-          else                  ! of if (scavenging)
+           else                  ! of if (scavenging)
              raveco2(:)=0
              do ig=1,ngrid 
@@ -2854 +2855 @@
      &               max(icetotco2(ig),1.e-30),1.e-30) ! mass weight
              enddo
-          endif                 ! of if (scavenging)
+           endif                 ! of if (scavenging)
+              
+           call WRITEdiagfi(ngrid,"raveco2","ice eff radius","m",0
+     &        ,raveco2)
           
-             
-         call WRITEdiagfi(ngrid,"raveco2","ice eff radius","m",0
-     &        ,raveco2)
+          endif ! of if (co2clouds)
 
            CALL WRITEDIAGFI(ngrid,'reffice',