ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/radiance/ray/src/gen/gendaylit.c
(Generate patch)

Comparing ray/src/gen/gendaylit.c (file contents):
Revision 2.1 by greg, Sat Jun 6 20:18:32 2009 UTC vs.
Revision 2.4 by greg, Wed Aug 10 22:28:45 2011 UTC

# Line 21 | Line 21 | static const char RCSid[] = "$Id$";
21   #include  <string.h>
22   #include  <math.h>
23   #include  <stdlib.h>
24 + #include  <ctype.h>
25  
26   #include  "rtio.h"
27   #include  "fvect.h"
# Line 71 | Line 72 | void   theta_phi_to_dzeta_gamma(double theta,double phi,
72   double  integ_lv(float *lv,float *theta);
73   float   *theta_ordered(char *filename);
74   float   *phi_ordered(char *filename);
75 + void    skip_comments(FILE *fp);
76  
77  
78  
# Line 271 | Line 273 | char  *argv[];
273   }
274  
275  
276 + void
277   computesky()                    /* compute sky parameters */
278   {
279  
# Line 470 | Line 473 | computesky()                   /* compute sky parameters */
473   /* Compute the ground radiance */
474   zenithbr=calc_rel_lum_perez(0.0,radians(sunzenith),radians(sunzenith),skyclearness,skybrightness,coeff_perez);
475   zenithbr*=diffnormalization;
476 + /*
477   fprintf(stderr, "gendaylit : the actual zenith radiance(W/m^2/sr) or luminance(cd/m^2) is : %.0lf\n", zenithbr);
478 <        
478 > */
479 >
480   if (skyclearness==1)
481          normfactor = 0.777778;
482                  
# Line 608 | Line 613 | register char  **av;
613  
614  
615  
616 + void
617 + skip_comments(FILE *fp)         /* skip comments in file */
618 + {
619 +        int     c;
620 +        
621 +        while ((c = getc(fp)) != EOF)
622 +                if (c == '#') {
623 +                        while ((c = getc(fp)) != EOF)
624 +                                if (c == '\n')
625 +                                        break;
626 +                } else if (!isspace(c)) {
627 +                        ungetc(c, fp);
628 +                        break;
629 +                }
630 + }
631  
632  
633  
614
615
616
617
618
634   /* Perez models */
635  
636   /* Perez global horizontal luminous efficacy model */
# Line 998 | Line 1013 | int lect_coeff_perez(char *filename,float **coeff_pere
1013          {
1014                  /*printf("file %s  open\n", filename);*/
1015          }
1016 +        
1017 +        skip_comments(fcoeff_perez);
1018  
1002        rewind(fcoeff_perez); /* on se place en debut de fichier */
1003
1019          for (i=0;i<8;i++)
1020                  for (j=0;j<20;j++)
1021                  {
# Line 1183 | Line 1198 | float *theta_ordered(char *filename)
1198                  fprintf(stderr,"Cannot open file %s in function theta_ordered\n",filename);
1199                  exit(1);
1200          }
1201 +        
1202 +        skip_comments(file_in);
1203  
1187        rewind(file_in);
1188
1204          if ( (ptr = malloc(145*sizeof(float))) == NULL )
1205          {
1206                  fprintf(stderr,"Out of memory in function theta_ordered\n");
# Line 1229 | Line 1244 | float *phi_ordered(char *filename)
1244                  fprintf(stderr,"Cannot open file %s in function phi_ordered\n",filename);
1245                  exit(1);
1246          }
1247 <
1248 <        rewind(file_in);
1247 >        
1248 >        skip_comments(file_in);
1249  
1250          if ( (ptr = malloc(145*sizeof(float))) == NULL )
1251          {

Diff Legend

Removed lines
+ Added lines
< Changed lines
> Changed lines