/*    
 * Vis5D system for visualizing five dimensional gridded data sets.
 * Copyright (C) 1990 - 2000 Bill Hibbard, Johan Kellum, Brian Paul,
 * Dave Santek, and Andre Battaiola.
 *
 * This program is free software; you can redistribute it and/or modify
 * it under the terms of the GNU General Public License as published by
 * the Free Software Foundation; either version 2 of the License, or
 * (at your option) any later version.
 *
 * As a special exception to the terms of the GNU General Public
 * License, you are permitted to link Vis5D with (and distribute the
 * resulting source and executables) the LUI library (copyright by
 * Stellar Computer Inc. and licensed for distribution with Vis5D),
 * the McIDAS library, and/or the NetCDF library, where those
 * libraries are governed by the terms of their own licenses.
 *
 * This program is distributed in the hope that it will be useful,
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 * GNU General Public License for more details.
 *
 * You should have received a copy of the GNU General Public License
 * along with this program; if not, write to the Free Software
 * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
 *
 */

#ifndef V5D_H
#define V5D_H

/* SGJ: use extern "C" if included from a C++ file: */
#ifdef __cplusplus
extern "C" {
#endif                /* __cplusplus */

/*
 * A numeric version number which we can test for in utility programs which
 * use the v5d functions.  For example, we can do tests like this:
 * #if V5D_VERSION > 42
 *      do something
 * #else
 *      do something else
 * #endif
 *
 * If V5D_VERSION is not defined, then its value is considered to be zero.
 */

#define V5D_VERSION 42


/*
 * Define our own 1 and 2-byte data types.  We use these names to avoid
 * collisions with types defined by the OS include files.
 */
typedef unsigned char V5Dubyte;     /* Must be 1 byte, except for cray */
typedef unsigned short V5Dushort;   /* Must be 2 byte, except for cray */



#define MISSING 1.0e35
#define IS_MISSING(X)  ( (X) >= 1.0e30 )


/* Limits on 5-D grid size:  (must match those in v5df.h!!!) */
#define MAXVARS     200
#define MAXTIMES    400
#define MAXROWS     400
#define MAXCOLUMNS  400 
#define MAXLEVELS   400

#define MAXRECS     10000

/************************************************************************/
/***                                                                  ***/
/*** Functions for writing v5d files.  See README file for details.   ***/
/*** These are the functions user's will want for writing file        ***/
/*** converters, etc.                                                 ***/
/***                                                                  ***/
/************************************************************************/
extern int v5dCreateSimple( const char *name,
                            int numtimes, int numvars,
                            int nr, int nc, int nl,
                            const char varname[MAXVARS][10],
                            const int timestamp[],
                            const int datestamp[],
                            float northlat, float latinc,
                            float westlon, float loninc,
                            float bottomhgt, float hgtinc );


extern int v5dCreate( const char *name,
                      int numtimes, int numvars,
                      int nr, int nc, const int nl[],
                      const char varname[MAXVARS][10],
                      const int timestamp[],
                      const int datestamp[],
                      int compressmode,
                      int projection,
                      const float proj_args[],
                      int vertical,
                      const float vert_args[] );


extern int v5dWrite( int time, int var, const float data[] );


extern int v5dClose( void );


extern int v5dSetLowLev( int lowlev[] );

extern int v5dSetUnits( int var, const char *units );



/************************************************************************/
/***                                                                  ***/
/*** Definition of v5d struct and function prototypes.                ***/
/*** These functions are used by vis5d and advanced v5d utilities.    ***/
/***                                                                  ***/
/************************************************************************/

/* ZLB 02-09-2000 */
#define MAXPROJARGS (MAXROWS+MAXCOLUMNS+1)
#define MAXVERTARGS (MAXLEVELS+1)

/*
 * This struct describes the structure of a .v5d file.
 */
typedef struct {
    /* PUBLIC (user can freely read, sometimes write, these fields) */
        int NumTimes;                   /* Number of time steps */
        int NumVars;                    /* Number of variables */
        int Nr;                         /* Number of rows */
        int Nc;                         /* Number of columns */
        int Nl[MAXVARS];                /* Number of levels per variable */
        int LowLev[MAXVARS];            /* Lowest level per variable */
        char VarName[MAXVARS][10];      /* 9-character variable names */
        char Units[MAXVARS][20];        /* 19-character units for variables */
        int TimeStamp[MAXTIMES];        /* Time in HHMMSS format */
        int DateStamp[MAXTIMES];        /* Date in YYDDD format */
        float MinVal[MAXVARS];          /* Minimum variable data values */
        float MaxVal[MAXVARS];          /* Maximum variable data values */

        /* This info is used for external function computation */
        short McFile[MAXTIMES][MAXVARS];/* McIDAS file number in 1..9999 */
        short McGrid[MAXTIMES][MAXVARS];/* McIDAS grid number in 1..? */

        int VerticalSystem;             /* Which vertical coordinate system */
        float VertArgs[MAXVERTARGS];    /* Vert. Coord. Sys. arguments... */

        /*
        IF VerticalSystem==0 THEN
                -- Linear scale, equally-spaced levels in generic units
                VertArgs[0] = Height of bottom-most grid level in generic units
                VertArgs[1] = Increment between levels in generic units
        ELSE IF VerticalSystem==1 THEN
                -- Linear scale, equally-spaced levels in km
                VertArgs[0] = Height of bottom grid level in km
                VertArgs[1] = Increment between levels in km
        ELSE IF VerticalSystem==2 THEN
                -- Linear scale, Unequally spaced levels in km
                VertArgs[0] = Height of grid level 0 (bottom) in km
                ...                ...
                VertArgs[n] = Height of grid level n in km
        ELSE IF VerticalSystem==3 THEN
                -- Linear scale, Unequally spaced levels in mb
                VertArgs[0] = Pressure of grid level 0 (bottom) in mb
                ...             ...
                VertArgs[n] = Pressure of grid level n in mb
        ENDIF
        */

        int Projection;                     /* Which map projection */
        float ProjArgs[MAXPROJARGS];        /* Map projection arguments... */

        /*
        IF Projection==0 THEN
                -- Rectilinear grid, generic units
                ProjArgs[0] = North bound, Y coordinate of grid row 0
                ProjArgs[1] = West bound, X coordiante of grid column 0
                ProjArgs[2] = Increment between rows
                ProjArgs[3] = Increment between colums
                NOTES: X coordinates increase to the right, Y increase upward.
                NOTES: Coordinate system is right-handed.
        ELSE IF Projection==1 THEN
                -- Cylindrical equidistant (Old Vis5d)
                -- Rectilinear grid in lat/lon
                ProjArgs[0] = Latitude of grid row 0, north bound, in degrees
                ProjArgs[1] = Longitude of grid column 0, west bound, in deg.
                ProjArgs[2] = Increment between rows in degrees
                ProjArgs[3] = Increment between rows in degrees
                NOTES: Coordinates (degrees) increase to the left and upward.
        ELSE IF Projection==2 THEN
                -- Lambert conformal
                ProjArgs[0] = Standared Latitude 1 of conic projection
                ProjArgs[1] = Standared Latitude 2 of conic projection
                ProjArgs[2] = Row of North/South pole
                ProjArgs[3] = Column of North/South pole
                ProjArgs[4] = Longitude which is parallel to columns
                ProjArgs[5] = Increment between grid columns in km
        ELSE IF Projection==3 THEN
                -- Polar Stereographic
                ProjArgs[0] = Latitude of center of projection
                ProjArgs[1] = Longitude of center of projection
                ProjArgs[2] = Grid row of center of projection
                ProjArgs[3] = Grid column of center of projection
                ProjArgs[4] = Increment between grid columns at center in km
        ELSE IF Projection==4 THEN
                -- Rotated
                ProjArgs[0] = Latitude on rotated globe of grid row 0
                ProjArgs[1] = Longitude on rotated globe of grid column 0
                ProjArgs[2] = Degrees of latitude on rotated globe between
                                grid rows
                ProjArgs[3] = Degrees of longitude on rotated globe between
                                grid columns
                ProjArgs[4] = Earth latitude of (0, 0) on rotated globe
                ProjArgs[5] = Earth longitude of (0, 0) on rotated globe
                ProjArgs[6] = Clockwise rotation of rotated globe in degrees
        ELSE IF Projection == 5 THEN
                -- Mercator
                ProjArgs[0] = Latitude of center of projection
                ProjArgs[1] = Longitude of center of projection
                ProjArgs[2] = Row Increment in Kilometers
                ProjArgs[3] = Column Increment in Kilometers
        ELSE IF Projection==-1 THEN
            ProjArgs[0]    = coordinate of latitude[0] (row[0])
            ...
            ProjArgs[n]    = coordinate of latitude[n] (row[n])

            ProjArgs[NR]   = coordinate of latitude[0] (column[0])
            ...
            ProjArgs[NR+n] = coordinate of latitude[n] (column[n])
        ENDIF
        */

        int CompressMode;        /* 1, 2 or 4 = # bytes per grid point */
        char FileVersion[10];    /* 9-character version number */

    /* PRIVATE (not to be touched by user code) */
        unsigned int FileFormat; /* COMP5D file version or 0 if .v5d */
        int FileDesc;            /* Unix file descriptor */
        char Mode;               /* 'r' = read, 'w' = write */
        int CurPos;              /* current position of file pointer */
        int FirstGridPos;        /* position of first grid in file */
        int GridSize[MAXVARS];   /* size of each grid */
        int SumGridSizes;        /* sum of GridSize[0..NumVars-1] */
} v5dstruct;



extern float pressure_to_height( float pressure);

extern float height_to_pressure( float height );




extern int v5dYYDDDtoDays( int yyddd );

extern int v5dHHMMSStoSeconds( int hhmmss );

extern int v5dDaysToYYDDD( int days );

extern int v5dSecondsToHHMMSS( int seconds );


extern void v5dPrintStruct( const v5dstruct *v );


extern v5dstruct *v5dNewStruct( void );

extern void v5dFreeStruct( v5dstruct* v );

extern void v5dInitStruct( v5dstruct *v );

extern int v5dVerifyStruct( const v5dstruct *v );


extern void v5dCompressGrid( int nr, int nc, int nl, int compressmode,
                             const float data[], void *compdata,
                             float ga[], float gb[],
                             float *minval, float *maxval );


extern void v5dDecompressGrid( int nr, int nc, int nl, int compressmode,
                               void *compdata,
                               float ga[], float gb[],
                               float data[] );


extern int v5dSizeofGrid( const v5dstruct *v, int time, int var );


extern v5dstruct *v5dOpenFile( const char *filename, v5dstruct *v );


extern int v5dCreateFile( const char *filename, v5dstruct *v );


extern v5dstruct *v5dUpdateFile( const char *filename, v5dstruct *v );


extern int v5dCloseFile( v5dstruct *v );


extern int v5dReadCompressedGrid( v5dstruct *v,
                                  int time, int var,
                                  float *ga, float *gb,
                                  void *compdata );


extern int v5dReadGrid( v5dstruct *v, int time, int var, float data[] );


extern int v5dWriteCompressedGrid( const v5dstruct *v,
                                   int time, int var,
                                   const float *ga, const float *gb,
                                   const void *compdata );


extern int v5dWriteGrid( v5dstruct *v, int time, int var, const float data[] );

  /* JPE added 09-19-2000  */
extern int v5dCreateStruct( v5dstruct *v, int numtimes, int numvars,
                                        int nr, int nc, const int nl[],
                                        const char varname[MAXVARS][10],
                                        const int timestamp[], const int datestamp[],
                                        int compressmode,
                                        int projection,
                                        const float proj_args[],
                                        int vertical,
                                        const float vert_args[] );



#ifdef __cplusplus
}                               /* extern "C" */
#endif                /* __cplusplus */

#endif
