Changeset 1722 for trunk/MagicSoft/Simulation/Detector/ReflectorII/reflector.c

Timestamp:

01/21/03 15:02:35 (22 years ago)

Author:

moralejo

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Simulation/Detector/ReflectorII/reflector.c (modified) (3 diffs)

trunk/MagicSoft/Simulation/Detector/ReflectorII/reflector.c

@@ -47 +47 @@
 extern CerEventHeader *cheadp;  /*  var inited in header.c      */
 extern RflEventHeader *rheadp;  /*  var inited in header.c      */
+extern void SetAtmModel(int model, float ol); /*  from atm.c    */
 
 /*  Prototypes  */
@@ -470 +471 @@
   RflRunHeader RflRunHead;
 
-  extern int atmModel;      /* current atm. model */
+  const char *AtmModelNames[]={"ATM_NOATMOSPHERE","ATM_90PERCENT","ATM_CORSIKA"};
+  extern char atmosphere[256];  /* current atm. model */
+  int model=0;
 
   do
@@ -490 +493 @@
           if (newFile)
             {
+              while (strcmp(atmosphere, AtmModelNames[model]))
+                if (++model == sizeof(AtmModelNames)/sizeof(AtmModelNames[0]))
+                  {   
+                    model = 0;
+                    Error(" *** Atm model \"%s\" not found.\n", atmosphere);
+                    break;      
+                  }
+
               /* Write Reflector "run header" (one per cer file!): */
               memcpy(&RflRunHead, cheadp, sizeof(RflRunHead));
               RflRunHead.wobble_mode = wobble_position;
-              RflRunHead.atmospheric_model = atmModel;
+              RflRunHead.atmospheric_model = model;
               fwrite(&RflRunHead, sizeof(RflRunHead), 1, rflfile);
+
+              /* Set up atmosphere (do some necessary calculations) once we 
+               * read in the observation level from the run header: 
+               */
+
+              Log("Setting atm model \"%s\" and observation level %.0f cm.\n", AtmModelNames[model], RflRunHead.HeightLev[0]);
+              SetAtmModel(model,RflRunHead.HeightLev[0]);
+
             }