#include <stdio.h>
#include <string.h>
double atof();
/* #define DEBUG */
#define TOKEN 40
#define NVAR 50
#define COMMENTSIZE 2430
#define SEPARATOR " ,=\t\n"

char names[NVAR][TOKEN];
char values[NVAR][TOKEN];
char attr[NVAR][TOKEN];
int map[NVAR];
int nheader, ndata;
FILE *fl;
int maxlev,newlev;
int eofflag=0;
char str[1025],*strp;
char tabcomments[COMMENTSIZE];

err(s,t)
char *s,*t;
{
char tmp[80];
printf("&x error - %s%t\r\n",s,t);
ioclose_();
exit(1);
}

int fg(str)
char *str;
{
  *str = 0;
  while(*str == 0){
    if(fgets(str,1024,fl) == NULL) return 1;
    if(strspn(str,SEPARATOR)==strlen(str))*str = 0;
  };
  return 0;
}

int readdata(flag)
int flag;
{
int state,nn,i;
char *pstr;

state=1;
if(*str == 0)
  if(fg(str))return -1;

  if (str[0] == '#'){
    *tabcomments=0;
    while (str[0] == '#') {
      strcat(tabcomments,str+1);
      if (fg(str)) return -1;
      };
    };

if(flag<=1) while(strchr(str,'=')==NULL) if(fg(str))return -1;

while(strchr(str,'=')){
  state=2;
  strp=strtok(str,SEPARATOR);
  while(strp){
    if(flag==0){
      if(pstr=strchr(strp,'[')){
	*(pstr+strlen(pstr)-1)=0;
	strcpy(attr[nheader],pstr+1);
	*pstr=0;
      } else attr[nheader][0]=0;
      strcpy(names[nheader],strp);
      nn=nheader++;
    } else {
      if(pstr=strchr(strp,'['))*pstr=0;
      for(nn=0;nn<nheader;nn++)if(strcmp(names[nn],strp)==0)break;
      if(nn==nheader)err("unknown header variable name found - ",strp);
    };
    strp=strtok(NULL,SEPARATOR);
    if(strp==NULL)err("no value assigned - ",names[nn]);
    strcpy(values[nn],strp);
    strp=strtok(NULL,SEPARATOR);
  };
  if(fg(str))return -1;
};

if (state==2 || flag<=1){
  strp=strtok(str,SEPARATOR);
  i=0;
  while(strp){
    if(flag==0){
      if(pstr=strchr(strp,'[')){
	*(pstr+strlen(pstr)-1)=0;
	strcpy(attr[nheader+ndata],pstr+1);
	*pstr=0;
      } else attr[nheader+ndata][0]=0;
      strcpy(names[nheader+ndata],strp);
      nn=nheader+ndata;
      ndata++;
    } else {
      if(pstr=strchr(strp,'['))*pstr=0;
      for(nn=nheader;nn<nheader+ndata;nn++)
	if(strcmp(names[nn],strp)==0)break;
      if(nn==(nheader+ndata))err("unknown data variable - ",strp);
    };
    map[i++]=nn;
    strp=strtok(NULL,SEPARATOR);
  };
  *str=0;
  return 1;
};

strp=strtok(str,SEPARATOR);
for(nn=0;nn<ndata;nn++){
  if(strp==NULL){
    if(fg(str)) return -1;
    strp=strtok(str,SEPARATOR);
  };
  strcpy(values[map[nn]],strp);
  strp=strtok(NULL,SEPARATOR);
};
*str = 0;
return 2;
}

int ioopen_(s,nparams,ntotal)
char *s[];
int *nparams,*ntotal;
{
char tmp[258],*tok;
int ncnt,i,j;
if((fl=fopen(s[0],"r")) == NULL)err("file name ",s);
s[0][0] = 0;
str[0]=0;
strp=str;
nheader=0;
ndata=0;
*tabcomments=0;
if(readdata(0)<=0)err("file structure","");
*ntotal=nheader+ndata;
if (nheader>0){newlev=3;return 1;};
newlev=2;return 0;
}

ioname_(vn,s)
int *vn;
char *s;
{
strcpy(s,names[*vn]);
}

int iovarlevel_(vn)
int *vn;
{
if(*vn<nheader) return 0;
return 1;
}

ioattrout_(vn,str)
int *vn;
char *str;
{
char *at;
int j;
j= *vn;
if (attr[j][0]){
  at=strchr(attr[j],';');
  if(at){
    *at = '\0';
    strcpy(str,attr[j]);
    strcpy(attr[j],at+1);
  } else {
    strcpy(str,attr[j]);
    attr[j][0]='\0';
  };
  return 1;
} else return 0;
}

int iocommout_(str)
char *str;
{
char *at;
if (tabcomments[0]){
  at=strchr(tabcomments,'\n');
  if(at){
    *at = '\0';
    strcpy(str,tabcomments);
    strcpy(tabcomments,at+1);
  } else {
    strcpy(str,tabcomments);
    tabcomments[0]='\0';
  };
  return 1;
} else return 0;
}
  

iovalreal_(vn,f)
int *vn;
float *f;
{
int i;
i= *vn;
if(i<0){*f= -9999.0;return;};
if(strspn(values[i],"0123456789.+-"))
   *f=atof(values[i]);
else
   *f = -9999.0;
#ifdef DEBUG
printf("iovalreal %d %f\n",*vn,f);
#endif
}

iovalstr_(vn,tmp)
int *vn;
char *tmp;
{
static char *s;
int i;
i= *vn;
if(i<0){strcpy(tmp,"nd");return;};
s=values[i];
s=s+strspn(s," ");
strcpy(tmp,s);
#ifdef DEBUG
printf("iovalstr %d %s\n",*vn,tmp);
#endif
}

int ioreadrec_(level)
int *level;
{
if(eofflag)return 0;
switch(*level){
 case 0:
  if(nheader){
    if (newlev & 1) {newlev=0;return 1;};
    newlev=readdata(1);
    if(newlev<0){eofflag=1;return 0;};
    newlev=newlev-1;
    return 1;
  } else {
    if (newlev & 2){newlev=0;return 1;};
    newlev=readdata(2);
    if(newlev<0){eofflag=1;return 0;};
    newlev=newlev-1;
    return 1;
  };
 case 1:
  if(newlev & 2){newlev=newlev-2; return 1;};
  newlev=readdata(2);
  if(newlev<0){eofflag=1;return 0;};
  if(newlev & 1) return 0;
  newlev=newlev-2;
  return 1;
};
}

ioclose_(maxlevel)
int *maxlevel;
{
fclose(fl);
}