Changeset 1431 for trunk/MagicSoft/Simulation/Detector/ReflectorII/attenu.f

Timestamp:

07/24/02 15:36:10 (22 years ago)

Author:

bigongia

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Simulation/Detector/ReflectorII/attenu.f (modified) (3 diffs)

trunk/MagicSoft/Simulation/Detector/ReflectorII/attenu.f

@@ -360 +360 @@
       x4 = sqrt((2. * h      + rt)    / (2. * hscale))
 
+c--   AM Dec 2001, to avoid crash! A few photons seem to be "corrupted" 
+c--    (have absurd value) in a cer file...
+
+      if (abs(x3-x4) .lt. 1.e-10) then
+         tr_atmos = -1.
+         RETURN
+      endif
+
+
       e1 = derfc(x1)
       e2 = derfc(x2)
@@ -366 +375 @@
 
       m = exp(-Rsin2 / (2. * hscale)) * ((e1 - e2) / (e3 - e4))
+
 
 **********************************************************************   
@@ -442 +452 @@
 
       CON_OZ = 2
+
  2001 IF (LONG(CON_OZ).GT.wavelength) THEN
-         A = (OZ_ABI(ROW,CON_OZ)-OZ_ABI(ROW,CON_OZ-1))/
+             A = (OZ_ABI(ROW,CON_OZ)-OZ_ABI(ROW,CON_OZ-1))/
      +        (LONG(CON_OZ) - LONG(CON_OZ-1))
          B = OZ_ABI(ROW,CON_OZ) - (A*LONG(CON_OZ))