/* 
rs - regrouping method
version 1.0 grf  used to establish levels in a flat file
version 1.1 clh  add subroutine to return column widths - iowidth_()
version 1.11 clh fixed rsw to respect the re-ordering of variables
version 2.0 Mar 18, 1996 clh to prepare for the major update distribution
        use minbufsize.h to determine size of buffers, from wjs.
*/
#if HP || IBM
#define jdbopen_ jdbopen
#define jdbreada_ jdbreada
#define jdbread_ jdbread
#define jdbclose_ jdbclose
#define jdbcomments_ jdbcomments
#define jdblevel_ jdblevel
#define jdbattributes_ jdbattributes
#endif
#include <stdio.h>
#include <string.h>
double atof();
#include <math.h>
#include "minbufsize.h"

/* #define DEBUG */
#define TOKEN 40

#ifndef NVAR
#define NVAR 250
#endif

#ifndef COMMENTSIZE
#define COMMENTSIZE 30*MINCOMMENTLINE
#endif

#ifndef VARNAMESIZE
#define VARNAMESIZE MINVARNAMESIZE
#endif

#ifndef ATTRSIZE
#define ATTRSIZE MINATTRSIZE
#endif

#ifndef DATUMSIZE
#define DATUMSIZE MINDATUMSIZE
#endif

#ifndef INBUFSIZE
#define INBUFSIZE 255
#endif

int nlevels,nnames;
int firstvar[6];
char names[NVAR][VARNAMESIZE];
int fldwidths[NVAR];
int namesize=VARNAMESIZE;
char values[NVAR][DATUMSIZE];
char comp[NVAR][DATUMSIZE];
int valuesize=DATUMSIZE;
int handle;
int minlevelread;
char comments[COMMENTSIZE];

/* #define TOK2 42  */
char ss[NVAR][DATUMSIZE];
int ns=0;
int newlev[NVAR];
int newpntr[NVAR];
int newlpntr[11];
int maxoutlev;

int ioreadrec_(level)
int *level;
{
  int clevel,i;

  if(*level == minlevelread)
    {
      minlevelread= *level+1;
      return 1;
    }
  else if(*level >minlevelread) return 0;

  while(1){
    clevel=jdbreada_(&handle,values,&valuesize);
    /* printf("clevel %d\n",clevel); */
    if(clevel<0){
      minlevelread= -1;
      return 0;
    };
    
    for(i=firstvar[clevel];i<nnames;i++){
      if(newlev[i]<minlevelread)
	if(strcmp(values[i],comp[i]))minlevelread=newlev[i];
      strcpy(comp[i],values[i]);
    };
    
    if (minlevelread < *level) return 0;
    else if(minlevelread == *level){
      minlevelread++;
      return 1;
    };
  };
}

void ioclose_()
{
  jdbclose_(&handle);
}

int ioopen_(s,nparams,ntotal)
char *s[];
int *nparams;
int *ntotal;
{
  char tmp[INBUFSIZE],*sp;
  int i,j,k,m,maxclev;

  strcpy(tmp,s[0]);
  nnames = -NVAR;
  putenv("PATH_INFO=");

  maxclev=jdbopen_(&handle,tmp,names,&namesize,&nnames);

  if (maxclev<0){printf("&x Bad Object:%s\n",s[0]);exit(1);};
  s[0][0]=0;
  *comments=0;
  while(jdbcomments_(&handle,tmp))
    {strcat(comments,tmp);strcat(comments,"\n");};
  for(i=1;i<*nparams;i++)if(s[i][0]>='0' && s[i][0]<='9' && s[i][1]==':'){
    strcpy(ss[ns++],s[i]);

    s[i][0]='\0';
  };
  
  j=0;
  firstvar[j]=0;
  for(i=0;i<nnames;i++){
    newlev[i]= -1;
    comp[i][0]=0;
    if((k=iovarlevel_(&handle,&i))>j)
      firstvar[++j]=i;
  };

  maxoutlev= -1;
  m=0;
  for(k=0;k<ns;k++){
    strcpy(tmp,ss[k]);
    i=tmp[0]-'0';
    if(i>maxoutlev){maxoutlev=i;newlpntr[maxoutlev]=m;};
    if(tmp[2]=='*'){
      for(j=0;j<nnames;j++)if(newlev[j]== -1){
	newlev[j]=maxoutlev;
	newpntr[m++]=j;
      };
    } else {
      for(j=0;j<nnames;j++)
	if(strcmp(names[j],tmp+2)==NULL)break;
      if(j<nnames){
	newlev[j]=maxoutlev;
	newpntr[m++]=j;
      };
    };
  };
  
  if(maxoutlev<0)maxoutlev=0;
  
  for(j=0;j<nnames;j++)
    if(newlev[j]<0){
      newlev[j]=maxoutlev;
      newpntr[m++]=j;
    };
  
  *ntotal=nnames;

  minlevelread = maxoutlev+1;
  return maxoutlev;
}

int iovarlevel_(vn)
int *vn;
{
  return newlev[newpntr[*vn]];
}

int ioattrout_(vn,str)
int *vn;
char *str;
{
   int j;

   j = newpntr[*vn];

   if ((jdbattributes_(&handle,&j,str)) == NULL)
       return 0;
   else {
       if (strncmp(str,"width=",6) == 0) {
         getwidth(&j,str);
         str = '\0';        /* don't show the width attr */
       }      
       return 1;
   }
}

getwidth(vn, ptr)
/* read the width from the string 'ptr' (an attribute that has
        been tested prior to coming here for containing the 
        string 'width=') and write it to the fldwidths array
        element '*vn'
*/
int *vn;
char *ptr;
{
  char *p;
  int len;

  p = strchr(ptr,'=');

  sscanf(p+1,"%d", &len);

  fldwidths[*vn] = len;
  return;
}

void iovalreal_(vn,f)
int *vn;
float *f;
{
  if(strspn(values[*vn],"0123456789.+-"))
    *f=atof(values[newpntr[*vn]]);
  else
    *f= -9999.0;
}

void iovalstr_(vn,tmp)
int *vn;
char *tmp;
{
  char *s;

  s=values[newpntr[*vn]];
  s=s+strspn(s," ");
  strcpy(tmp,s);
}

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

int iocommout_(str)
char *str;
{
char *at;

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

int iowidth_(vn)
int *vn;
{
   int i;

   i = fldwidths[newpntr[*vn]];
return i;
}