#From glenn@lake.mit.edu Fri Jul  7 11:37 EDT 1995
#Date: Fri, 7 Jul 95 11:38:02 EDT
#From: Glenn Flierl <glenn@lake.mit.edu>
#To: jmanning@whsun1.wh.whoi.edu, rgroman@whoi.edu
#Subject: Matlab method

#include <stdio.h>
#include <string.h>
#include <math.h>
#include "/usr/local/matlab/extern/include/mat.h"
double atof();
#define MAXVAR 20

MATFile *fl;
char **names;
static int nn;
static int lev[MAXVAR],ptr[MAXVAR];
static int maxlev=0;
static float oldval[20];
double *mp[MAXVAR];
static int index= -1;
static int maxindex=0;

/* flags for controlling reading sequence */
#define OK0 -999
#define ERR1 -9999

/* error reporting */
err(s,t)
char *s,*t;
{
printf("&x error - &s &s\n",s,t);
exit(1);
}

/* 
     open subroutine 
*/
int ioopen_(s,nparams,ntotal)
char *s[];
int *nparams,*ntotal;
{
char tmp[80];
FILE *fh;
int i,j,k;
double *mpi;

/* open file named in first parameter */
strcpy(tmp,s[0]);
strcat(tmp,".mat");
fl=matOpen(tmp,"r");
if(fl==NULL)err("MATLAB file name ",s[0]);

names=matGetDir(fl,&nn);
if(names==NULL)err("No Variables in MATLAB FILE","");
*ntotal=nn;

for(i=0;i<nn;i++){lev[i]=0;oldval[i]= -9999.0;ptr[i]=i;};

strcpy(tmp,s[0]);
strcat(tmp,".hdr");
if((fh=fopen(tmp,"r")) != NULL){
  while(!feof(fh)){
    fscanf(fh,"%d %s",&i,tmp);
    for(j=0;j<nn;j++)if(strcmp(names[j],tmp)==0)break;
    if(j==nn)err("Header structure ",tmp);
    lev[j]=i;
    if(i>maxlev)maxlev=i;
  };
  fclose(fh);
};

s[0][0] = 0;

if(maxlev){
  k=0;
  for(i=maxlev;i>=0;i--){
    for(j=0;j<nn;j++)if(lev[j]==i)ptr[k++]=j;
  };
};

for(i=0;i<nn;i++){
  if(matGetFull(fl,names[ptr[i]],&j,&k,&(mp[i]),&mpi))
         err("Missing matrix ",names[ptr[i]]);
  if(k != 1)err("Can't do matrices yet","");
  if(maxindex==0)maxindex=j;
  if(j != maxindex)err("Not uniform size","");
};

matClose(fl);

return maxlev;
}

/* 
     return name of variable number *vn 
*/
ioname_(vn,s)
int *vn;
char *s;
{
strcpy(s,names[ptr[*vn]]);
}

/* 
      return level of variable number *vn 
*/
int iovarlevel_(vn)
int *vn;
{
return lev[ptr[*vn]];
}

/* 
      output attributes of variable number *vn 
*/
int ioattrout_(vn)
int *vn;
{
return 0;
}

/* 
    return value of variable *vn as a real number 
*/
iovalreal_(vn,f)
int *vn;
float *f;
{
*f=mp[*vn][index];
}

/* 
      return value of variable *vn as a string 
*/
iovalstr_(vn,tmp)
int *vn;
char *tmp;
{
char s[80];
float f;
f=mp[*vn][index];
if(f == -9999.0){
  strcpy(tmp,"nd");
  return;
};
sprintf(s,"%g",f);
strcpy(tmp,s+strspn(s," "));
}

/* 
      get next record at requested level 
*/
int ioreadrec_(level)
int *level;
{
if(index>=maxindex-1)return 0;
index++;
return 1;
}

/* 
    close file 
*/
ioclose_()
{
}

/*
   comments
*/
int iocommout_()
{
return 0;
}