7 extern char forecast_basin_file[100];
9 extern struct map map[1500];
13 int ib,ier,l,numpts,ip,x,y,ip2,ip3,ip4,zones,
i,maxib,dum;
18 char hb5[10],parm[10],wfo[10];
27 if(basin_file[0]==0 ) {
29 printf(
"no basin_file specified\n");
36 printf(
"no basin_file specified\n");
41 fr=fopen(basin_file,
"r");
45 printf(
"could not open basin_file %s\n",basin_file);
79 ier=sscanf(ibuf,
"%d %d %s",&numpts,&dum,
map[
ib].hb5);
87 for(l=0;l<numpts;l++) {
99 ier=sscanf(ibuf,
"%lf %lf",&
lon,&
lat);
129 map[
ib].
gz=calloc(200,
sizeof(
float));
130 map[
ib].
uz=calloc(200,
sizeof(
float));
131 map[
ib].
mz=calloc(200,
sizeof(
float));
132 map[
ib].
lz=calloc(200,
sizeof(
float));
133 map[
ib].
gzc=calloc(24,
sizeof(
float));
134 map[
ib].
uzc=calloc(24,
sizeof(
float));
135 map[
ib].
mzc=calloc(24,
sizeof(
float));
136 map[
ib].
lzc=calloc(24,
sizeof(
float));
137 map[
ib].
zgz=calloc(200,
sizeof(
float));
138 map[
ib].
zuz=calloc(200,
sizeof(
float));
139 map[
ib].
zmz=calloc(200,
sizeof(
float));
140 map[
ib].
zlz=calloc(200,
sizeof(
float));
141 map[
ib].
tgz=calloc(200,
sizeof(
float));
142 map[
ib].
tuz=calloc(200,
sizeof(
float));
143 map[
ib].
tmz=calloc(200,
sizeof(
float));
144 map[
ib].
tlz=calloc(200,
sizeof(
float));
149 for(m=0;m < 200;m++) {
183 ier=sscanf(ibuf,
"%d %d %*s %s",&numpts,&zones,bchar);
185 strcpy(
map[
ib].bchar,bchar);
192 map[
ib].zones[mm]=-1;
194 for(l=0;l<numpts;l++) {
199 p=fgets(ibuf,100,
fp);
204 ier=sscanf(ibuf,
"%d %d %d %d %d %d\n",&x,&y,&ip,&ip2,&ip3,&ip4);
206 if(ip < 0 || ip > 4) {
220 if(ip2 < 0 || ip2 > 4) {
235 if(ip3 < 0 || ip3 > 4) {
249 if(ip4 < 0 || ip4 > 4) {
274 for(
i=0;
i<maxib;
i++) {
287 if(climo_file[0] != 0) {
289 fp=fopen(climo_file,
"r");
293 printf(
"could not open %s\n",climo_file);
301 p=fgets(ibuf,200,
fp);
306 ier=sscanf(ibuf,
"%s %s %f %f %f %f %f %f %f %f %f %f %f %f",
307 hb5,parm,&f[0],&f[1],&f[2],&f[3],&f[4],&f[5],&f[6],&f[7],&f[8],
308 &f[9],&f[10],&f[11]);
314 for(
i=0;
i<maxib;
i++) {
316 len=strlen(
map[
i].hb5);
318 if(strncmp(
map[
i].hb5,hb5,len)==0) {
324 if(
map[
i].zones[l] < 0)
353 amapvalue=amapvalue +
map[
i].
lzc[m];
383 if(forecast_basin_file[0] == 0 )
386 fr=fopen(forecast_basin_file,
"r");
401 ier=sscanf(ibuf,
"%s %s",hb5,wfo);
405 if(strcasecmp(
map[
ib].hb5,hb5)==0) {
409 if(strcmp(wfo,
tag[m].wfo)==0) {
void get_basin_data(char *basin_file, char *climo_file, char *hrap_file)
printf("fbuf is %s\n", fbuf)
struct hrap_data * hrap_data