Index: trunk/MagicSoft/Mars/Changelog
===================================================================
--- trunk/MagicSoft/Mars/Changelog	(revision 9428)
+++ trunk/MagicSoft/Mars/Changelog	(revision 9429)
@@ -100,4 +100,8 @@
    * mbase/MParSpline.[h,cc]:
      - added
+
+   * msim/MSimAtmosphere.cc:
+     - fixed a bug which caused the aerosol coefficient to be
+       incorrectly initialized
 
 
Index: trunk/MagicSoft/Mars/NEWS
===================================================================
--- trunk/MagicSoft/Mars/NEWS	(revision 9428)
+++ trunk/MagicSoft/Mars/NEWS	(revision 9429)
@@ -35,4 +35,7 @@
    * There is no default pulse shape anymoe it must be given in 
      the resource file
+
+   * Fixed a bug which caused the aerosol coefficient to be incorrectly
+     initialized. The effect might be random.
 
  ;callisto:
Index: trunk/MagicSoft/Mars/msim/MSimAtmosphere.cc
===================================================================
--- trunk/MagicSoft/Mars/msim/MSimAtmosphere.cc	(revision 9428)
+++ trunk/MagicSoft/Mars/msim/MSimAtmosphere.cc	(revision 9429)
@@ -554,5 +554,5 @@
 
         if (IsValid())
-            PreCalcOzone();
+            PreCalcAerosol();
 
         return kTRUE;
@@ -621,5 +621,4 @@
             h = 30.e5;
 
-        // Mie absorption coefficient at sea level through interpolation
         //const float aero_betap[15] = {0.27, 0.26, 0.25, 0.24, 0.24, 0.23, 0.20, 0.180, 0.167, 0.158, 0.150, 0.142, 0.135, 0.127, 0.120};
         //const Double_t beta0 = getbeta(wavelength, aero_betap);
@@ -630,4 +629,5 @@
         const UInt_t H = TMath::Nint(h/1e4);
         const UInt_t T = TMath::Min(89, TMath::Nint(theta/STEPTHETA));
+
 
         const Double_t path = aerosol_path[H][T];