ViewVC Help
View File | Revision Log | Show Annotations | Download File | Root Listing
root/radiance/ray/src/cv/epw2wea.c
Revision: 2.1
Committed: Wed Sep 25 16:52:28 2013 UTC (10 years, 7 months ago) by greg
Content type: text/plain
Branch: MAIN
Log Message:
Added epw2wea by Christoph Reinhart

File Contents

# Content
1 #ifndef lint
2 static const char RCSid[] = "$Id$";
3 #endif
4 /* Copyright (c) 2003
5 * National Research Council Canada
6 * written by Christoph Reinhart
7 */
8
9 /* epw2wea: daylight analysis subprogram of DAYSIM
10 Program converts EnergyPlus weather format (*.ppw) into DAYSIM format (*.wea)
11 */
12
13 #include <stdlib.h>
14 #include <stdio.h>
15 #include <string.h>
16 #include <errno.h>
17
18 int
19 main(int argc, char *argv[])
20 {
21 FILE *EPW_FILE;
22 FILE* WEA_FILE;
23 int year, month,day, hour_in,minute,i;
24 float dir_norm_rad, dif_or_rad,dummy_float;
25 char keyword[2000]="";
26 char epw_file[200]="";
27 char wea_file[200] ="";
28 char city[200] ="";
29 char country[200] ="";
30 char latitude[200] ="",longitude[200] ="",time_zone[200] ="",elevation[200] ="";
31
32 if (argc < 3)
33 {
34 fprintf(stderr,"epw2wea: FATAL ERROR - wrong number of arguments\n");
35 fprintf(stderr,"start program with: epw2wea file-name.epw file-name.wea\n");
36 return(1);
37 }
38 if (argc >= 3)
39 {
40 strcpy(epw_file, argv[1]);
41 strcpy(wea_file, argv[2]);
42 for (i = 3; i < argc; i++)
43 if (argv[i][0] == '-' )
44 switch (argv[i][1])
45 {
46 case 'h': /* scaling factor */
47 break;
48 }//end switch (argv[i][1])
49 else
50 {
51 fprintf(stdout,"epw2wea: fatal error - %s bad option for input arguments\n", argv[i]);
52 return(1);
53 }//end else
54 }
55
56 EPW_FILE=fopen(epw_file, "r");;
57 WEA_FILE=fopen(wea_file, "w");;
58 fscanf(EPW_FILE,"%[^,]s",keyword);
59 if(strcmp(keyword,"LOCATION") ){
60 fprintf(stderr,"epw2wea: FATAL ERROR - this does not seem to be an epw file \n");
61 return(1);
62 }
63 fscanf(EPW_FILE,",%[^,]s",city);
64 fscanf(EPW_FILE,",%[^,]s",country);
65 fscanf(EPW_FILE,",%[^,]s",country);
66 sprintf(keyword,"place %s_%s\n",city,country);
67 printf("%s",keyword);
68 fprintf(WEA_FILE,"%s",keyword);
69
70 fscanf(EPW_FILE,",%[^,]s",country);
71 fscanf(EPW_FILE,",%[^,]s",country);//printf("%s\n",country);
72 fscanf(EPW_FILE,",%[^,]s",latitude);
73 printf("latitude %s\n",latitude);
74 fprintf(WEA_FILE,"latitude %s\n",latitude);
75 fscanf(EPW_FILE,",%[^,]s",longitude);
76
77 printf("longitude %.2f\n",-1.0*atof(longitude));
78 fprintf(WEA_FILE,"longitude %.2f\n",-1.0*atof(longitude));
79 fscanf(EPW_FILE,",%[^,]s",time_zone);
80 printf("time_zone %.0f\n",-15.0*atoi(time_zone));
81 fprintf(WEA_FILE,"time_zone %.0f\n",-15.0*atoi(time_zone));
82 fscanf(EPW_FILE,",%s[^\n]",elevation);
83 printf("site_elevation %s\nweather_data_file_units 1\n",elevation);
84 fprintf(WEA_FILE,"site_elevation %s\nweather_data_file_units 1\n",elevation);
85
86 fscanf(EPW_FILE,"%*[^\n]");fscanf(EPW_FILE,"%*[\n\r]");
87 fscanf(EPW_FILE,"%*[^\n]");fscanf(EPW_FILE,"%*[\n\r]");
88 fscanf(EPW_FILE,"%*[^\n]");fscanf(EPW_FILE,"%*[\n\r]");
89 fscanf(EPW_FILE,"%*[^\n]");fscanf(EPW_FILE,"%*[\n\r]");
90 fscanf(EPW_FILE,"%*[^\n]");fscanf(EPW_FILE,"%*[\n\r]");
91 fscanf(EPW_FILE,"%*[^\n]");fscanf(EPW_FILE,"%*[\n\r]");
92 fscanf(EPW_FILE,"%*[^\n]");fscanf(EPW_FILE,"%*[\n\r]");
93 fscanf(EPW_FILE,"%*[^\n]");fscanf(EPW_FILE,"%*[\n\r]");
94
95
96 while( EOF != fscanf(EPW_FILE,"%d,%d,%d,%d,%d",&year,&month,&day, &hour_in,&minute)){
97 fprintf(WEA_FILE,"%d %d %.3f ",month,day,hour_in*1.0-minute*(0.5/60));
98
99 fscanf(EPW_FILE,",%[^,]s",city);
100 fscanf(EPW_FILE,",%f",&dummy_float);
101 fscanf(EPW_FILE,",%f",&dummy_float);
102 fscanf(EPW_FILE,",%f",&dummy_float);
103 fscanf(EPW_FILE,",%f",&dummy_float);
104 fscanf(EPW_FILE,",%f",&dummy_float);
105 fscanf(EPW_FILE,",%f",&dummy_float);
106 fscanf(EPW_FILE,",%f",&dummy_float);
107 fscanf(EPW_FILE,",%f",&dummy_float);
108
109 fscanf(EPW_FILE,",%f,%f",&dir_norm_rad, &dif_or_rad);
110 fprintf(WEA_FILE,"%.0f %.0f",dir_norm_rad, dif_or_rad);
111
112 fscanf(EPW_FILE,"%*[^\n]");fscanf(EPW_FILE,"%*[\n\r]");
113 fprintf(WEA_FILE,"\n");
114 }
115
116 fclose(EPW_FILE);
117 fclose(WEA_FILE);
118
119 return(0);
120 }
121
122