#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include <string.h>
#include "mucal.h"
#include "panel.h"
/* ********************************************** */

int init_facet_REFL(
   reflect REFL[],
   material MATER[],
   facet  FACET[],
   int    LEN[],
   int    L_order[numLayers][N_CryPos],
   int    numOptions)
{
   int ifac, iifac, nl, iOpt;
   static int ill[3] = {1, 1, 1};
/*
   printf("init_facet_REFL: LEN: %6d %6d %6d, numOptions: %3d, First facet options, layer 0/1/2: %2d %2d %2d\n", 
      LEN[0], LEN[1], LEN[2], numOptions, FACET[L_order[0][0]].nOption, 
      FACET[(L_order[1][0]+N_CryPos)].nOption, FACET[(L_order[2][0]+(N_CryPos*2))].nOption);
*/
   for (nl=0; nl<numLayers; nl++) {
      for (iifac=0; iifac<LEN[nl]; iifac++) {
         ifac = L_order[nl][iifac] + nl * N_CryPos;
         if (FACET[ifac].nOption < 0) continue;
         iOpt = FACET[ifac].nOption;
         if ((REFL[iOpt]).Z <= 0) continue;
         reflect2(&(REFL[iOpt]), &(FACET[ifac]), iOpt);  
         reflect2(&(REFL[2]), &(FACET[ifac]), 2);  
         reflect2(&(REFL[5]), &(FACET[ifac]), 5);  
/*         if (ill[nl]-- > 0) printf("init_facet_REFL: Layer %1d NumOpt: %2d, Opt %2d, d_spacing %8.6lf (AA)\n", 
                                    nl, numOptions, iOpt, REFL[iOpt].dist_geo*1.0e+8);
*/
      }
   }
   return(0);
}