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

Comparing src/gen/gensky.c (file contents):
Revision 2.2 by greg, Thu Dec 19 15:09:28 1991 UTC vs.
Revision 2.4 by greg, Fri Mar 20 12:23:38 1992 UTC

# Line 35 | Line 35 | extern double  s_latitude;
35   extern double  s_longitude;
36   extern double  s_meridian;
37                                          /* required values */
38 < int  month, day;
39 < double  hour;
38 > int  month, day;                                /* date */
39 > double  hour;                                   /* standard time */
40 > double  altitude, azimuth;                      /* or solar angles */
41                                          /* default values */
42 < int  cloudy = 0;
42 > int  cloudy = 0;                                /* 1=standard, 2=uniform */
43   int  dosun = 1;
44   double  zenithbr = -1.0;
45   double  turbidity = 2.75;
# Line 67 | Line 68 | char  *argv[];
68          }
69          if (argc < 4)
70                  userror("arg count");
71 <        month = atoi(argv[1]);
72 <        day = atoi(argv[2]);
73 <        hour = atof(argv[3]);
71 >        if (!strcmp(argv[1], "-ang")) {
72 >                altitude = atof(argv[2]) * (PI/180);
73 >                azimuth = atof(argv[3]) * (PI/180);
74 >                month = 0;
75 >        } else {
76 >                month = atoi(argv[1]);
77 >                day = atoi(argv[2]);
78 >                hour = atof(argv[3]);
79 >        }
80          for (i = 4; i < argc; i++)
81                  if (argv[i][0] == '-' || argv[i][0] == '+')
82                          switch (argv[i][1]) {
# Line 78 | Line 85 | char  *argv[];
85                                  dosun = argv[i][0] == '+';
86                                  break;
87                          case 'c':
88 <                                cloudy = 1;
88 >                                cloudy = argv[i][0] == '+' ? 2 : 1;
89                                  dosun = 0;
90                                  break;
91                          case 't':
# Line 120 | Line 127 | char  *argv[];
127  
128   computesky()                    /* compute sky parameters */
129   {
123 –        int  jd;
124 –        double  sd, st;
125 –        double  altitude, azimuth;
130                                          /* compute solar direction */
131 <        jd = jdate(month, day);                 /* Julian date */
132 <        sd = sdec(jd);                          /* solar declination */
133 <        st = hour + stadj(jd);                  /* solar time */
134 <        altitude = salt(sd, st);
135 <        azimuth = sazi(sd, st);
131 >        if (month) {                    /* from date and time */
132 >                int  jd;
133 >                double  sd, st;
134 >
135 >                jd = jdate(month, day);         /* Julian date */
136 >                sd = sdec(jd);                  /* solar declination */
137 >                st = hour + stadj(jd);          /* solar time */
138 >                altitude = salt(sd, st);
139 >                azimuth = sazi(sd, st);
140 >        }
141          sundir[0] = -sin(azimuth)*cos(altitude);
142          sundir[1] = -cos(azimuth)*cos(altitude);
143          sundir[2] = sin(altitude);
# Line 146 | Line 155 | computesky()                   /* compute sky parameters */
155                  zenithbr = 0.0;
156                                          /* Compute horizontal radiance */
157          if (cloudy) {
158 <                groundbr = zenithbr*0.777778;
158 >                if (cloudy == 2)
159 >                        groundbr = zenithbr;
160 >                else
161 >                        groundbr = zenithbr*0.777778;
162                  printf("# Ground ambient level: %f\n", groundbr);
163          } else {
164                  F2 = 0.274*(0.91 + 10.0*exp(-3.0*(PI/2.0-altitude)) +
# Line 182 | Line 194 | printsky()                     /* print out sky */
194          printf("2 skybright skybright.cal\n");
195          printf("0\n");
196          if (cloudy)
197 <                printf("3 1 %.2e %.2e\n", zenithbr, groundbr);
197 >                printf("3 %d %.2e %.2e\n", cloudy, zenithbr, groundbr);
198          else
199                  printf("7 -1 %.2e %.2e %.2e %f %f %f\n", zenithbr, groundbr,
200                                  F2, sundir[0], sundir[1], sundir[2]);
# Line 191 | Line 203 | printsky()                     /* print out sky */
203  
204   printdefaults()                 /* print default values */
205   {
206 <        if (cloudy)
206 >        if (cloudy == 1)
207                  printf("-c\t\t\t\t# Cloudy sky\n");
208 +        else if (cloudy == 2)
209 +                printf("+c\t\t\t\t# Uniform cloudy sky\n");
210          else if (dosun)
211                  printf("+s\t\t\t\t# Sunny sky with sun\n");
212          else
# Line 214 | Line 228 | char  *msg;
228          if (msg != NULL)
229                  fprintf(stderr, "%s: Use error - %s\n", progname, msg);
230          fprintf(stderr, "Usage: %s month day hour [options]\n", progname);
231 +        fprintf(stderr, "   Or: %s -ang altitude azimuth [options]\n", progname);
232          fprintf(stderr, "   Or: %s -defaults\n", progname);
233          exit(1);
234   }

Diff Legend

– Removed lines
+ Added lines
< Changed lines
> Changed lines