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

Comparing ray/src/cv/epw2wea.c (file contents):
Revision 2.1 by greg, Wed Sep 25 16:52:28 2013 UTC vs.
Revision 2.5 by greg, Fri Aug 2 21:18:40 2024 UTC

# Line 6 | Line 6 | static const char RCSid[] = "$Id$";
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 < */
9 > /* epw2wea: daylight analysis subprogram of DAYSIM */
10 > /* Program converts EnergyPlus weather format (*.ppw) into DAYSIM format (*.wea) */
11  
13 – #include <stdlib.h>
12   #include <stdio.h>
13 + #include <stdlib.h>
14   #include <string.h>
16 – #include <errno.h>
15  
16 < int
19 < main(int argc, char *argv[])
16 > int main( int argc, char  *argv[])
17   {
18 <        FILE *EPW_FILE;
19 <        FILE* WEA_FILE;
20 <        int year, month,day, hour_in,minute,i;
18 >        FILE *EPW_FILE = NULL;
19 >        FILE* WEA_FILE = stdout;
20 >        int year, month,day, hour_in,minute=60,i;
21 >        int minute_message=1;
22 >        int get_ac = 0;
23          float dir_norm_rad, dif_or_rad,dummy_float;
24 <        char keyword[2000]="";
25 <        char epw_file[200]="";
26 <        char wea_file[200] ="";
24 >        float aod;
25 >        int cc;
26 >    char keyword[2000]="";
27 >        char minute_string[2000]="";
28 >        char *epw_file = NULL;
29 >        char *wea_file = NULL;
30          char city[200] ="";
31          char country[200] ="";
32          char latitude[200] ="",longitude[200] ="",time_zone[200] ="",elevation[200] ="";
33  
34 <        if (argc < 3)
34 >        for ( ; argc > 2 && argv[1][0] == '-'; argc--, argv++)
35 >                if (argv[1][1] == 'a') {
36 >                        get_ac = !get_ac;
37 >                } else {
38 >                        fprintf(stderr, "epw2wea: unknown option: %s\n", argv[1]);
39 >                        exit(1);
40 >                }
41 >
42 >        if ((argc < 2) | (argc > 3))
43          {
44                  fprintf(stderr,"epw2wea: FATAL ERROR - wrong number of arguments\n");
45 <                fprintf(stderr,"start program with: epw2wea file-name.epw file-name.wea\n");
46 <                return(1);
45 >                fprintf(stderr,"Usage: epw2wea [-a] file-name.epw [file-name.wea]\n");
46 >                exit(1);
47          }
48 <        if (argc >= 3)
49 <        {
50 <                strcpy(epw_file, argv[1]);
51 <                strcpy(wea_file, argv[2]);
52 <                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
48 >        epw_file = argv[1];
49 >        EPW_FILE=fopen(epw_file, "r");
50 >        if (!EPW_FILE) {
51 >                fprintf(stderr, "%s: cannot open for reading\n", epw_file);
52 >                exit(1);
53          }
54 <
55 <        EPW_FILE=fopen(epw_file, "r");;
56 <        WEA_FILE=fopen(wea_file, "w");;
57 <        fscanf(EPW_FILE,"%[^,]s",keyword);
58 <        if(strcmp(keyword,"LOCATION") ){
59 <                fprintf(stderr,"epw2wea: FATAL ERROR - this does not seem to be an epw file \n");
60 <                return(1);
54 >        if (argc == 3) {
55 >                wea_file = argv[2];
56 >                WEA_FILE=fopen(wea_file, "w");
57 >                if (!WEA_FILE) {
58 >                        fprintf(stderr, "%s: cannot open for writing\n", wea_file);
59 >                        exit(1);
60 >                }
61          }
62 +        fscanf(EPW_FILE,"%[^,]s",keyword);
63 +        if( !strcmp(keyword,"LOCATION") ){
64          fscanf(EPW_FILE,",%[^,]s",city);
65          fscanf(EPW_FILE,",%[^,]s",country);
66          fscanf(EPW_FILE,",%[^,]s",country);
# Line 68 | Line 69 | main(int argc, char *argv[])
69          fprintf(WEA_FILE,"%s",keyword);
70  
71          fscanf(EPW_FILE,",%[^,]s",country);
72 <        fscanf(EPW_FILE,",%[^,]s",country);//printf("%s\n",country);
72 >        fscanf(EPW_FILE,",%[^,]s",country);
73          fscanf(EPW_FILE,",%[^,]s",latitude);
74          printf("latitude %s\n",latitude);
75          fprintf(WEA_FILE,"latitude %s\n",latitude);
# Line 77 | Line 78 | main(int argc, char *argv[])
78          printf("longitude %.2f\n",-1.0*atof(longitude));
79          fprintf(WEA_FILE,"longitude %.2f\n",-1.0*atof(longitude));
80          fscanf(EPW_FILE,",%[^,]s",time_zone);
81 <        printf("time_zone %.0f\n",-15.0*atoi(time_zone));
81 >        printf("time_zone %.2f\n",-15.0*atof(time_zone));
82          fprintf(WEA_FILE,"time_zone %.0f\n",-15.0*atoi(time_zone));
83          fscanf(EPW_FILE,",%s[^\n]",elevation);
84          printf("site_elevation %s\nweather_data_file_units 1\n",elevation);
# Line 90 | Line 91 | main(int argc, char *argv[])
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 +        /* read in time step interval */
96 +        fscanf(EPW_FILE,"%[^,]s",keyword);
97 +        fscanf(EPW_FILE,",%[^,]s",keyword);
98 +        fscanf(EPW_FILE,",%[^,]s",minute_string);
99 +        minute=atoi(minute_string);
100 +        if(minute==1)   /* one measurement per hour equals a 60 minute time step */
101 +                minute=60;
102          fscanf(EPW_FILE,"%*[^\n]");fscanf(EPW_FILE,"%*[\n\r]");
103 <
104 <
105 <        while( EOF != fscanf(EPW_FILE,"%d,%d,%d,%d,%d",&year,&month,&day, &hour_in,&minute)){
103 >    
104 >        while( EOF != fscanf(EPW_FILE,"%d,%d,%d,%d",&year,&month,&day, &hour_in)){
105 >                
106                  fprintf(WEA_FILE,"%d %d %.3f ",month,day,hour_in*1.0-minute*(0.5/60));
107 <
107 >                
108 >                fscanf(EPW_FILE,",%f",&dummy_float);
109                  fscanf(EPW_FILE,",%[^,]s",city);
110                  fscanf(EPW_FILE,",%f",&dummy_float);
111                  fscanf(EPW_FILE,",%f",&dummy_float);
# Line 108 | Line 118 | main(int argc, char *argv[])
118  
119                  fscanf(EPW_FILE,",%f,%f",&dir_norm_rad, &dif_or_rad);
120                  fprintf(WEA_FILE,"%.0f %.0f",dir_norm_rad, dif_or_rad);
121 +                if (get_ac){
122 +                        fscanf(EPW_FILE,",%f",&dummy_float);
123 +                        fscanf(EPW_FILE,",%f",&dummy_float);
124 +                        fscanf(EPW_FILE,",%f",&dummy_float);
125 +                        fscanf(EPW_FILE,",%f",&dummy_float);
126 +                        fscanf(EPW_FILE,",%f",&dummy_float);
127 +                        fscanf(EPW_FILE,",%f",&dummy_float);
128 +                        fscanf(EPW_FILE,",%d",&cc);
129 +                        fscanf(EPW_FILE,",%f",&dummy_float);
130 +                        fscanf(EPW_FILE,",%f",&dummy_float);
131 +                        fscanf(EPW_FILE,",%f",&dummy_float);
132 +                        fscanf(EPW_FILE,",%f",&dummy_float);
133 +                        fscanf(EPW_FILE,",%f",&dummy_float);
134 +                        fscanf(EPW_FILE,",%f",&dummy_float);
135 +                        fscanf(EPW_FILE,",%f",&aod);
136  
137 +                        fprintf(WEA_FILE," %.3f %.1f",aod, cc/10.);
138 +                        
139 +                }
140 +
141                  fscanf(EPW_FILE,"%*[^\n]");fscanf(EPW_FILE,"%*[\n\r]");
142                  fprintf(WEA_FILE,"\n");
143          }
144  
145 <        fclose(EPW_FILE);
146 <        fclose(WEA_FILE);
145 > }else{
146 >                fprintf(stderr,"epw2wea: FATAL ERROR - this does not seem to be an epw file \n");exit(1);
147 > }
148  
149 <        return(0);
149 > fclose(EPW_FILE);
150 > fclose(WEA_FILE);
151 > return 0;
152   }
153  
154  

Diff Legend

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