Changeset 420 for trunk/LMDZ.MARS/libf/phymars/physiq.F

Timestamp:

Nov 24, 2011, 9:54:55 AM (13 years ago)

Author:

tnavarro

Message:

24/11/11 == TN

corrected minor bug in updatereffrad.F : reffdust was not saved

ccn_factor as not correctly used in sedimentation.

It is now initialized in inifis.F, declared in tracer.h and
used in both simpleclouds.F & callsedim.F to update ice radius.

Commented diagfi outputs in aeropacity.F & improvedclouds.F for non scavenging users.

File:

• trunk/LMDZ.MARS/libf/phymars/physiq.F (modified) (2 diffs)

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

@@ -318 +318 @@
 c Test 1d/3d scavenging
       real h2o_tot
-      real ccndust_mass(nlayermx)
-      real ccndust_number(nlayermx)
-      real rescale                ! to rescale GCM dust quantities
 
       REAL time_phys
@@ -1994 +1991 @@
            do l=1,nlayer
              h2o_tot = h2o_tot + 
-     &           (zq(ig,l,igcm_h2o_vap) + zq(ig,l,igcm_h2o_ice))
-     &                     * (pplev(ig,l) - pplev(ig,l+1)) / g
-             ccndust_mass(l) = 
-     &       pq(1,l,igcm_dust_mass)+pq(1,l,igcm_ccn_mass)
-             ccndust_number(l) = 
-     &       pq(1,l,igcm_dust_number)+pq(1,l,igcm_ccn_number)
-     
+     &           (zq(1,l,igcm_h2o_vap) + zq(1,l,igcm_h2o_ice))
+     &                     * (pplev(1,l) - pplev(1,l+1)) / g
            end do