C语言站点数据转换为Grads数据格式代码

daveryjl

int i;
int j;
int ci;
float cf;
char stcode[8];
char endstn[8];
FILE *infile;
if( BuildValue==0 )
{
  if( ( infile=fopen(filename,"wb") )==NULL )
  {
   printf("the %s file don't be build.\n",filename);
   return -1;
  }
}
else
{
  if( ( infile=fopen(filename,"ab") )==NULL )
  {
   printf("the %s file don't be opened.\n",filename);
   return -11;
  }
}

//for(i=0;i<hn;i++)
for(i=0;i<HstNum;i++)
{
  j=i;//stinf.No;
  //****************每站一头*********************
  sprintf(stcode,"%d",st[j].stcode);
  fwrite(stcode,sizeof(char)*8,1,infile);  
  fwrite(&st[j].lat,sizeof(float),1,infile);
  fwrite(&st[j].lon,sizeof(float),1,infile);
  cf=0.0;
  fwrite(&cf,sizeof(float),1,infile);//时间偏差
  ci=1;
  fwrite(&ci,sizeof(int),1,infile);//总层数
  ci=1;
  fwrite(&ci,sizeof(int),1,infile);//是否有地面报
  //************地面报**************************
  fwrite(&st[j].height,sizeof(float),1,infile);
  fwrite(&st[j].wv,sizeof(float),1,infile);
  fwrite(&st[j].wd,sizeof(float),1,infile);
  fwrite(&st[j].temp,sizeof(float),1,infile);
  fwrite(&st[j].td,sizeof(float),1,infile);
  fwrite(&st[j].SeaPressue,sizeof(float),1,infile);
  //************end 地面报**************************
}
///*************每一时的结尾*********************
sprintf(endstn,"end_time");
fwrite(endstn,sizeof(char)*8,1,infile);
cf=0.0;
fwrite(&cf,sizeof(float),1,infile);
fwrite(&cf,sizeof(float),1,infile);
fwrite(&cf,sizeof(float),1,infile);
ci=0;
fwrite(&ci,sizeof(int),1,infile);
fwrite(&ci,sizeof(int),1,infile);
///*************end  每一时的结尾*********************
fclose(infile);
return 0;