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

#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>
/* #define DEBUG */
#define TOKEN 40
#define NVAR 150
#define COMMSIZE 2048

int nlevels;
int firstvar[10];
char names[NVAR][TOKEN];
int namesize=TOKEN;
char values[NVAR][TOKEN];
int valuesize=TOKEN;
int lev[NVAR],vpntr[NVAR];
int handle;
int minlevelread;

char *nm[] = {"sal","temp","press","theta","sigth","svan","dynhgt"};
int nmpt[] = {-1,-1,-1,-1,-1,-1,-1};
float nmvl[7];
int newvarpt[] = {-1,-1,-1,-1,-1,-1,-1};
int newvarmin = 10000;
int newvarmax = 10000;
float sigma,th;
float thref = 0.0;
float dyn,p0,del0;

float atg();
float theta();
float svan();

char pstrings[10][255];
int npstrings;
int ncrit;
char comments[COMMSIZE];

docalc()
{
int i;
if(minlevelread<ncrit){dyn=0.0;p0=0.0;};
for (i = 0;i<3;i++)nmvl[i] = atof(values[nmpt[i]]);
    nmvl[5] = svan(nmvl[0],nmvl[1],nmvl[2],&sigma);
    th = theta(nmvl[0],nmvl[1],nmvl[2],thref);
    nmvl[3] = th;
    th = svan(nmvl[0],th,thref,&sigma);
    nmvl[4] = sigma;
    dyn = dyn+(del0+nmvl[5])*(nmvl[2]-p0)*5.0e-6;
    p0 = nmvl[2];
    del0 = nmvl[5];
    nmvl[6] = dyn;
for (i = 3;i<7;i++)
  sprintf(values[nmpt[i]],"%7.4f",nmvl[i]);
}

doinit(nparams,s)
char *s[];
int *nparams;
{
int i,j,k;

dyn = 0.0;
p0 = 0.0;

ncrit= 0;
for(k=0;k<=nlevels;k++)
for(i=firstvar[k];i<firstvar[k+1];i++){
  for(j=0;j<=6;j++) if(!strcmp(nm[j],names[i])){
    nmpt[j]=i;
    if(j<=2)ncrit=k;};
};
if(nmpt[0]<0 || nmpt[1]<0 || nmpt[2]<0){
  printf("&x error bad variable names for dyn"); 
  exit(1);};

j = firstvar[ncrit+1];
k=0;
for(i = 3;i<7;i++) if (nmpt[i]<0){
  strcpy(names[j],nm[i]);
  nmpt[i]=j;
  lev[j++]=ncrit;
  k++;  
};
for(i = ncrit+1;i<=nlevels+1;i++)firstvar[i] += k;
}

int ioreadrec_(level)
int *level;
{
  char tmp[80];
  int i,clevel;
  
  /*
     printf("** level %d  minlevelread %d \n",*level,minlevelread);
     */
  
  if(*level == minlevelread)
    {
      minlevelread= *level+1;
      return 1;
    }
  else if(*level >minlevelread) return 0;
  do {
    minlevelread=jdbreada_(&handle,values,&valuesize); 
  } while (minlevelread > *level);
  if(minlevelread<=ncrit) docalc();
  if (minlevelread < *level) return 0;
  else {
    minlevelread++;
    return 1;
  };
}

ioclose_()
{
  jdbclose_(&handle);
}

int ioopen_(s,nparams,ntotal)
char *s[];
int *nparams;
int *ntotal;
{
  char tmp[255],*sp;
  int i,j,k,m,maxclev,nread;
  
  for (i=0;i<NVAR;i++)lev[i]= -1;
  strcpy(tmp,s[0]);
  nread = -NVAR;
  putenv("QUERY_STRING=");
  maxclev=jdbopen_(&handle,tmp,names,&namesize,&nread);
  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");};
  nlevels=0;
  firstvar[nlevels]=0;
  for(i=0;i<nread;i++){
    lev[i]=jdblevel_(&handle,&i);
    if(lev[i]>firstvar[nlevels])firstvar[++nlevels]=lev[i];
  };
  firstvar[nlevels+1]=nread;

  doinit(nparams,s);

  k=0;
  *ntotal=firstvar[nlevels+1];
  for(i=0;i<=nlevels;i++)
    for(j=0;j< *ntotal;j++)if(lev[j]==i)vpntr[k++]=j;
  minlevelread = nlevels+1;
  return nlevels;
}

int iovarlevel_(vn)
int *vn;
{
  return lev[vpntr[*vn]];
}

ioattrout_(vn,str)
int *vn;
char *str;
{
  return 0;
}

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

char *iovalstr_(vn,tmp)
int *vn;
char *tmp;
{
  char *s;
  s=values[vpntr[*vn]];
  s=s+strspn(s," ");
  strcpy(tmp,s);
}

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

int iocommout_(s)
char *s;
{
  return 0;
}