--- ray/src/gen/gensky.c	1994/01/31 12:53:49	2.14
+++ ray/src/gen/gensky.c	1996/07/16 14:14:50	2.17
@@ -17,12 +17,16 @@ static char SCCSid[] = "$SunId$ LBL";
 
 #include  <math.h>
 
+#include  <ctype.h>
+
 #include  "color.h"
 
 extern char  *strcpy(), *strcat(), *malloc();
-extern double  stadj(), sdec(), sazi(), salt();
+extern double  stadj(), sdec(), sazi(), salt(), tz2mer();
 
-#define  PI		3.141592654
+#ifndef  PI
+#define  PI		3.14159265358979323846
+#endif
 
 #define  DOT(v1,v2)	(v1[0]*v2[0]+v1[1]*v2[1]+v1[2]*v2[2])
 
@@ -38,6 +42,28 @@ double  normsc();
 extern double  s_latitude;
 extern double  s_longitude;
 extern double  s_meridian;
+
+#undef  toupper
+#define  toupper(c)	((c) & ~0x20)	/* ASCII trick to convert case */
+
+					/* European and North American zones */
+struct {
+	char	zname[8];	/* time zone name (all caps) */
+	float	zmer;		/* standard meridian */
+} tzone[] = {
+	"YST", 135, "YDT", 120,
+	"PST", 120, "PDT", 105,
+	"MST", 105, "MDT", 90,
+	"CST", 90, "CDT", 75,
+	"EST", 75, "EDT", 60,
+	"AST", 60, "ADT", 45,
+	"NST", 52.5, "NDT", 37.5,
+	"GMT", 0, "BST", -15,
+	"WET", -15, "WETDST", -30,
+	"MET", -30, "METDST", -45,
+	"MEZ", -30, "MESZ", -45,
+	"", 0
+};
 					/* required values */
 int  month, day;				/* date */
 double  hour;					/* time */
@@ -85,10 +111,7 @@ char  *argv[];
 		day = atoi(argv[2]);
 		if (day < 1 || day > 31)
 			userror("bad day");
-		hour = atof(argv[3]);
-		if (hour < 0 || hour >= 24)
-			userror("bad hour");
-		tsolar = argv[3][0] == '+';
+		cvthour(argv[3]);
 	}
 	for (i = 4; i < argc; i++)
 		if (argv[i][0] == '-' || argv[i][0] == '+')
@@ -148,6 +171,8 @@ char  *argv[];
 
 	computesky();
 	printsky();
+
+	exit(0);
 }
 
 
@@ -167,6 +192,7 @@ computesky()			/* compute sky parameters */
 			st = hour + stadj(jd);
 		altitude = salt(sd, st);
 		azimuth = sazi(sd, st);
+		printf("# Local solar time: %.2f\n", st);
 		printf("# Solar altitude and azimuth: %.1f %.1f\n",
 				180./PI*altitude, 180./PI*azimuth);
 	}
@@ -323,6 +349,47 @@ normsc()			/* compute normalization factor (E0*F2/L0) 
 		nsc = nsc*x + nf[i];
 
 	return(nsc);
+}
+
+
+cvthour(hs)			/* convert hour string */
+char  *hs;
+{
+	register char  *cp = hs;
+	register int	i, j;
+
+	if (tsolar = *cp == '+') cp++;		/* solar time? */
+	while (isdigit(*cp)) cp++;
+	if (*cp == ':')
+		hour = atoi(hs) + atoi(++cp)/60.0;
+	else {
+		hour = atof(hs);
+		if (*cp == '.') cp++;
+	}
+	while (isdigit(*cp)) cp++;
+	if (!*cp)
+		return;
+	if (tsolar || !isalpha(*cp)) {
+		fprintf(stderr, "%s: bad time format: %s\n", progname, hs);
+		exit(1);
+	}
+	i = 0;
+	do {
+		for (j = 0; cp[j]; j++)
+			if (toupper(cp[j]) != tzone[i].zname[j])
+				break;
+		if (!cp[j] && !tzone[i].zname[j]) {
+			s_meridian = tzone[i].zmer * (PI/180);
+			return;
+		}
+	} while (tzone[i++].zname[0]);
+
+	fprintf(stderr, "%s: unknown time zone: %s\n", progname, cp);
+	fprintf(stderr, "Known time zones:\n\t%s", tzone[0].zname);
+	for (i = 1; tzone[i].zname[0]; i++)
+		fprintf(stderr, " %s", tzone[i].zname);
+	putc('\n', stderr);
+	exit(1);
 }