/* *******************************************************************
 *                                                                   *
 * Copyright (c) L-DGO/MIT/JGOFS                                     *
 *                                                                   *
 *                                                                   *
 * File             : tabp.c                                         *
 *                                                                   *
 * Purpose          :                                                *         
 *                                                                   *
 * Version Number   : 1.1                                            *
 *                                                                   *
 * Revision History :                                                *
 *                                                                   *
 * Date                              Developer                       *
 * ----                              ---------                       *
 *                                                                   *
 * Sat Oct 17  1992 10.00           Glenn Flierl                     *
 * Oct 1993                         tabattributes                       *
 *                                                                   *   
 *                                                                   *
 ******************************************************************* */

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

char tabnames[NVAR][TOKEN];
char tabvalues[NVAR][TOKEN];
char tabattributes[NVAR][TOKEN];
char tabcomments[COMMENTSIZE];

int nvarlevel[6];
int fileopen[5] = {0,0,0,0,0};
FILE *fl[5];
int maxlev;
int pointers[NVAR];
char dirstring[80];

err(s,t)
char *s,*t;
{
printf("&x %s%s\n",s,t);
exit(1);
}

int scanheader(lev)
int lev;
{
int ncnt,i,j,k,m;
char tmp[1258],*tok,attrsav[TOKEN],*attr;
ncnt= 0;
i=0;j=1;
if(lev>0){
  for(i=nvarlevel[lev-1];i<nvarlevel[lev];i++)pointers[i]= -1;
  m=nvarlevel[lev-1];
  };
while(j){
  if (fgets(tmp,1257,fl[lev]) == NULL)err("Bad file structure","");
  if (tmp[0] == '#'){
    while (tmp[0] == '#') {
      strcat(tabcomments,tmp+1);
      if (fgets(tmp,1257,fl[lev]) == NULL)err("Bad file structure","");
      };
    };
  tok=strtok(tmp,SEPARATOR);
  while (tok != NULL){
    if(tok[0] == '"'){tok++,tok[strlen(tok)-1]=0;};
    if(attr=strchr(tok,'[')){
      strcpy(attrsav,attr+1);
      attrsav[strlen(attrsav)-1]='\0';
      *attr = '\0';
    } else attrsav[0]='\0';
    if(lev==0){
      strcpy(tabnames[ncnt],tok);
      strcpy(tabattributes[ncnt],attrsav);
      pointers[ncnt]=ncnt;ncnt++;}
    else if(strcmp(tok,">")){
      for(k=nvarlevel[lev-1];k<nvarlevel[lev];k++)
	if(strcmp(tabnames[k],tok)==0)break;
      if(k==nvarlevel[lev])err("New variable name in file ",tok);
      pointers[k]=m++;
      };
    tok=strtok(NULL,SEPARATOR);
    };
  if(lev==0){
    if(tabnames[ncnt-1][0]=='>') ncnt--; else j=0;
    nvarlevel[i]=ncnt;
    if(j) i++;
    } else j=0;
  };
if(lev==0)maxlev=i; else
for(i=lev;i<maxlev;i++)
  if(fgets(tmp,1257,fl[lev])==NULL)err("Bad file structure","");
return ncnt;
}

int ioopen_(s,nparams,ntotal)
char *s[];
int *nparams,*ntotal;
{
char tmp[1258],*tok;
int ncnt,i,j;
strcpy(dirstring,s[0]);
if((tok=strrchr(dirstring,DIRSEP))!=NULL)*(tok+1)='\0';
  else dirstring[0]='\0';
if((fl[0]=fopen(s[0],"r")) == NULL)err("Bad header file name ",s[0]);
s[0][0] = 0;
*tabcomments = '\0';
*ntotal=scanheader(0);
fileopen[0]=1;
return maxlev;
}

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

int iovarlevel_(vn)
int *vn;
{
int i;
for(i=maxlev-1;i>=0;i--)if (*vn >= nvarlevel[i]) return i+1;
return 0;
}

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

iovalreal_(vn,f)
int *vn;
float *f;
{
int i;
i=pointers[*vn];
if(i<0){*f= -9999.0;return;};
if(strspn(tabvalues[i],"0123456789.+-"))
   *f=atof(tabvalues[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=pointers[*vn];
if(i<0){strcpy(tmp,"nd");return;};
s=tabvalues[i];
s=s+strspn(s," ");
strcpy(tmp,s);
#ifdef DEBUG
printf("iovalstr %d %s\n",*vn,tmp);
#endif
}

static char subfilename[80];

int ioreadrec_(level)
int *level;
{
char tmp[1258],*tok;
int i,j,k;
#ifdef DEBUG
printf("ioreadrec %d\n",*level);
#endif

k= *level;

if (! fileopen[k]) {
#ifdef DEBUG
printf("open %s\n",subfilename);
#endif
  if((fl[k]=fopen(subfilename,"r"))  == NULL)err("Missing file ",subfilename);
  scanheader(k);
  fileopen[k]=1;
  };

j=1;
while(j){
  if(fgets(tmp,1257,fl[k]) == NULL){
#ifdef DEBUG
printf("eof at level %d\n",k);
#endif
    fclose(fl[k]);
    fileopen[k]=0;
    return 0;
    };
  i=strspn(tmp,SEPARATOR);
  if(i<strlen(tmp))j=0;
};

if(k==0)j=0; else j=nvarlevel[k-1];

tok=strtok(tmp,SEPARATOR);
for(i=j;i<nvarlevel[k];i++){
#ifdef DEBUG
printf("assign %d %s\n",i,tok);
#endif
    strcpy(tabvalues[i],tok);
    tok=strtok(NULL,SEPARATOR);
    };

if(k<maxlev){strcpy(subfilename,dirstring);strcat(subfilename,tok);};
return 1;
}

ioclose_()
{
int i;
for(i=0;i<=maxlev;i++)if (fileopen[i])fclose(fl[i]);
}

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;
}