Changeset 688

Timestamp:

03/05/01 16:41:16 (24 years ago)

Author:

harald

Message:

There was an error if one starts simulation from zenithangle = 0.
Ciro send me the solution.

File:

• trunk/MagicSoft/Simulation/Corsika/Mmcs/main.f (modified) (1 diff)

trunk/MagicSoft/Simulation/Corsika/Mmcs/main.f

@@ -739 +739 @@
 c And then, RE-calculate the GLOBAL direction in CORSIKA
 c We use formulae for spherical triangles
-        theprim = acos( cos(THETAP)*cos(spinthe)+
+
+        if ( THETAP .EQ. 0.0 ) then                   
+           theprim = spinthe 
+           phiprim = spinphi
+        elseif ( spinthe .eq. 0.0 ) then 
+           theprim = THETAP 
+           phiprim = 0.0
+        else 
+           theprim = acos( cos(THETAP)*cos(spinthe)+
      $                  sin(THETAP)*sin(spinthe)*cos(spinphi) )
-        phiprim = asin( sin(spinthe)*sin(spinphi)/sin(theprim) )
+           phiprim = asin( sin(spinthe)*sin(spinphi)/sin(theprim) )
+        endif 
         THETAP = theprim
         EVTH(140) = spinthe