Serial.c

Go to the documentation of this file.

/* General notes on revisions to this source file. 
       Nov/Dec 2004 v2.3.2: documentation upgrade 
                - Unused functions removed 
                - Re-organized, clarified scopes, Doxygen tags added 
                - Functionality same as v2.1/v2.2 
        Oct 2006 v2.6.0: added model-experiment parameters (boundary condition experiments right now)
        June 2008 v2.8.2: expanded options for spatial map outputs
                - (re)implemented use of HDF (v. 4) for spatial output of uchar (4 byte unsigned) only
                - implemented 4 byte floating point output in generic binary maps
                - implemented use of NetCDF (v. 3) for floating point maps
        
*/

#include "serial.h"

byte readMap( char* filename, void* data)
{
  int gsize, nbytes; byte bsize;
  FILE* bFile;
  char binFileName[300];
  
  bsize = read_header(filename);

  sprintf(binFileName,"%s/%s/Data/Map_bin/%s",ModelPath,ProjName,filename);

  if((bFile = fopen(binFileName,"rb") ) == NULL ) { 
    fprintf(stderr,"Can't Open Mapfile: %s",binFileName); 
    exit(0); 
  }
  gsize = gbl_size[0] * gbl_size[1];
  nbytes = fread((char*)data, (int)bsize, gsize, bFile );
  if(debug) { 
    fprintf(dFile," %d of %d items Read for %s, size = %x bytes ",
            nbytes,gsize,filename,bsize); 
  }
  fclose(bFile);
  return bsize;
}


int read_header(char* filename)
{
  FILE* file;
  char headerFileName[300];
  
  int rows, cols, i, Dsize=1, test, offset;
  char ch, dirCh;
  long int start;
        
  sprintf(headerFileName,"%s/%s/Data/Map_head/%s",ModelPath,ProjName,filename);

  file = fopen(headerFileName,"r");
  if(file==NULL) { 
    fprintf(stderr,"Unable to open header file %s.\n",headerFileName); 
    fflush(stderr); 
    exit(0);
  }
  else { 
    fprintf(dFile,"Successfully opened header file %s.\n",headerFileName);  
    fflush(dFile); 
  } 
/* v2.8.2 - get the UTM coordinates of model, for use in self-documenting NetCDF map output */
  scan_forward(file,"north:");
  fscanf(file,"%f",&UTMnorth);
  fprintf(dFile,"UTMnorth = %f\n",UTMnorth); 
  scan_forward(file,"south:");
  fscanf(file,"%f",&UTMsouth);
  fprintf(dFile,"UTMsouth = %f\n",UTMsouth); 
  scan_forward(file,"east:");
  fscanf(file,"%f",&UTMeast);
  fprintf(dFile,"UTMeast = %f\n",UTMeast); 
  scan_forward(file,"west:");
  fscanf(file,"%f",&UTMwest);
  fprintf(dFile,"UTMwest = %f\n",UTMwest); 

  scan_forward(file,"cols:");
  fscanf(file,"%d",&cols);
  fprintf(dFile,"cols = %d\n",cols); 
  scan_forward(file,"rows:");
  fscanf(file,"%d",&rows);
  fprintf(dFile,"rows = %d\n",rows);
  fflush(dFile);
  if( !( rows == gbl_size[0] && cols == gbl_size[1] ) ) { 
    fprintf(stderr,"Error, Wrong map size: %s: (%d,%d) : (%d,%d)\n",
             headerFileName,rows,gbl_size[0],cols,gbl_size[1]); 
  }
  start = ftell(file);
  test = scan_forward(file,"format:");
  if(test < 1) { 
    clearerr(file); 
    fseek(file, start, 0); 
  }
  else {        
    fscanf(file,"%d",&Dsize);
    fprintf(dFile,"size = %d\n",Dsize);
    Dsize = Dsize + 1; /* using header generated by GRASS, the byte size is actually the value read from header plus 1 */
    if( Dsize < 1 || Dsize > 4 ) {  
      fprintf(stderr,"Error, Illegal map size: %s  (%d)\n",headerFileName,Dsize); 
      exit(0); 
    }
  }
  
/* v2.8.2 - convert the ELM UTM coordinates from NAD27 to NAD83, for use in the netCDF output */
  UTMnorthNAD83 = UTMnorth + conv_northing_NAD27toNAD83;
  UTMsouthNAD83 = UTMsouth + conv_northing_NAD27toNAD83;
  UTMeastNAD83  = UTMeast  + conv_easting_NAD27toNAD83;
  UTMwestNAD83  = UTMwest  + conv_easting_NAD27toNAD83;
  
  fclose(file);
  return Dsize;
}

void writeMap(char* filename, void* data, int bSize, unsigned char type, int index)
{
  int gsize, istat;
  FILE* bFile; 
  
if (index == 0 && type=='M') write_header(filename,bSize);
  
  if(type == 'H') {
#if HDF4
    if(index==0) { istat = DFR8putimage(filename,(char*)data,s1,s0,0); }
    else { istat = DFR8addimage(filename,(char*)data,s1,s0,0); }
    if(istat != 0) printf("\nerror writing HDF file: %s\n",filename);
#endif
  } 

  else if(type == 'C') {
    printf("\nError in configuring netCDF output, which can only be 4-byte floating point values: %s\n",filename);
      exit(0); 
  } 

  else {
    if((bFile = fopen(filename,"wb") ) == NULL ) { 
      fprintf(stderr,"Can't Open Mapfile: %s",filename); 
      exit(0); 
    }
    gsize = gbl_size[0]*gbl_size[1]*bSize;
    fwrite((char*)data, gsize, 1, bFile ); 
        
    fclose(bFile);
  }
}

void writeFloatMap(char* filename, char* VARNAME, float* data, int bSize, unsigned char type, int index, char* thisvarUnits, int outstep)
{
  int gsize;
  FILE* bFile; 
  
    
  if(type == 'H') { /* HDF */
    printf("\nError in configuring HDF output, which can only be 1-byte values (you requested 4): %s\n",filename);
      exit(0); 
  } 

  else if(type == 'C') { /* netCDF */
        writeFloatNetCDF(filename,VARNAME,data,bSize,type,index,thisvarUnits, outstep); 

  } 
  else if(type == 'M')  { /* generic binary Map */
    if (index == 0) write_header(filename,bSize);

        if((bFile = fopen(filename,"wb") ) == NULL ) { 
      fprintf(stderr,"Can't Open output Mapfile: %s",filename); 
      exit(0); 
    }
    gsize = (gbl_size[0]+2)*(gbl_size[1]+2)*bSize; /* note size is 2 rows and 2 cols larger than standard input/output arrays */
    fwrite(data, gsize, 1, bFile ); 
        
    fclose(bFile);
  }

  else { /* this error shouldn't be reached, but it is possible */
    printf("\nError in configuring 4-byte output request, which can only be netCDF ('C' argument in 'G' command) or generic map binary ('M' argument in 'G' command): output variable= %s\n",filename);
  } 

}

/* Handle errors by printing an error message and exiting with a
 * non-zero status. */
#define ERRCODE 2
#define ERR(e) {printf("\nError in writing netCDF output map: %s\n", nc_strerror(e)); exit(ERRCODE);}
#define NDIMS 3

void writeFloatNetCDF(char* filename, char* VARNAME, float* data, int bSize, unsigned char type, int index, char* thisvarUnits, int outstep)
{
  int retval;
  int ncid, dimids[NDIMS];
  size_t start[NDIMS], count[NDIMS];
  int time_dimid, x_dimid, y_dimid;
  int x_id, y_id, time_id, varid, date_id, TM_id;
  char s_mo[3], s_da[3], ELMdate[11], dateUnits[100], attribute_string[1000];
  float north[s0+2],east[s1+2];
  float northrev[s0+2];  /* for the reversal of the row ordering in having output w/ origin at lower left */

  time_t begin_sim, end_sim; /* Feb 2011 - just threw in here - poor practice!! */
  struct tm *local_begin_time; /* struct holding real-time calendar time from standard time.h */
  /* a couple of items that use the standard time.h to clock the simulation runtime*/
  begin_sim = time(NULL); 
  local_begin_time = localtime(&begin_sim);

  
  int i,j, irev;

/* this prevents code from being compiled when netCDF is not installed 
        (If you do not have netCDF installed on your system, you must set NCDF3 to false in globals.h) */
#if NCDF3
/*  append the date (yyyymmdd) to the root filename*/
  if (SimTime.mo[0]<10) sprintf(s_mo,"0%d", SimTime.mo[0]);
  else sprintf(s_mo,"%d", SimTime.mo[0]);
  if (SimTime.da[0]<10) sprintf(s_da,"0%d", SimTime.da[0]);
  else sprintf(s_da,"%d", SimTime.da[0]);
  sprintf(ELMdate,"%d-%s-%s",SimTime.yr[0],s_mo,s_da);
  
/*  Coordinates used for the netCDF cell attributes are the centroids:
        float maps are 1-cell larger on each side, so
        subtract 1 cell width to start at western edge, and 
        add 1 cell width to start at northern edge 
        (remember that origin of SME/ELM is upper left,
        while UTM goes in positive direction from lower left) 
        We then use the centroid (not upper and left edge) of each cell for coords,
        so (- celWid + 0.5*celWid) = (-0.5*celWid) */
for (i=0; i< (s0+2); i++) { 
    for (j=0; j< (s1+2); j++) { 
//               east[j] = UTMwest  - 0.5 * celWid + (j * celWid); 
//               north[i]= UTMnorth + 0.5 * celWid - (i * celWid); 
// Note that we've switched to having output in NAD83 horizontal datum - if you change this, be sure to alter the metadata strings below
         east[j] = UTMwestNAD83  - 0.5 * celWid + (j * celWid); 
         north[i]= UTMnorthNAD83 + 0.5 * celWid - (i * celWid); 
    }
  }

/* reverse row order (ELM origin is upper left, this outputs origin at lower left) */
for (i=0; i< (s0+2); i++) { 
    for (j=0; j< (s1+2); j++) { 
         irev = (s0+1) - i; 
         nc_dataRowRev[T(irev,j)] = data[T(i,j)]; 
         northrev[irev] = north[i]; 
    }
  }

   /* Create the file. The NC_CLOBBER parameter tells netCDF to
    * overwrite this file, if it already exists; NC_NOCLOBBER does not overwrite.  */
   if (index==0) {
                 if ( (retval = nc_create(filename, NC_CLOBBER, &ncid)))
                           ERR(retval);

   /* Define the dimensions: x, y, and time outstep-number (sequence). NetCDF will hand back an ID for each. */
           if ((retval = nc_def_dim(ncid, "time", NC_UNLIMITED, &time_dimid)))
              ERR(retval);
           if ((retval = nc_def_dim(ncid, "y", (s0+2), &y_dimid))) /* note the s0+2 array size */
              ERR(retval);
           if ((retval = nc_def_dim(ncid, "x", (s1+2), &x_dimid))) /* note the s1+2 array size */
              ERR(retval);

   /* The dimids array is used to pass the IDs of the dimensions of
    * the variable. */
           dimids[0] = time_dimid;
           dimids[1] = y_dimid;
           dimids[2] = x_dimid;

   /* Define the model output variable, time, x, y, and transverse_mercator (required for EverVIEW application) variables. 
   The type of the model output variable in this case is * NC_FLOAT (4-byte floating point). */
           if ((retval = nc_def_var(ncid, "transverse_mercator", NC_CHAR, 0, 
                            0, &TM_id)))
              {  printf("oops: %s",nc_strerror(retval) );  exit(-1); }

           if ((retval = nc_def_var(ncid, VARNAME, NC_FLOAT, NDIMS, 
                            dimids, &varid)))
              ERR(retval);
           if ((retval = nc_def_var(ncid, "time", NC_FLOAT, 1, 
                            &time_dimid, &time_id)))  
              ERR(retval);
           if ((retval = nc_def_var(ncid, "y", NC_FLOAT, 1, 
                            &y_dimid, &y_id)))
              ERR(retval);
           if ((retval = nc_def_var(ncid, "x", NC_FLOAT, 1, 
                            &x_dimid, &x_id)))
              ERR(retval);

   /* #### Attributes of the time variable. #### */
           if ((retval = nc_put_att_text(ncid, time_id, "long_name", 
                                         strlen("model time"), "model time")))
              ERR(retval);
   /* Assign units attributes to the julian time. In definition mode, the current date is initial day 0 of simulation */
           sprintf(dateUnits,"days since %sT12:00:00Z",ELMdate);
           if ((retval = nc_put_att_text(ncid, time_id, "units", 
                                         strlen(dateUnits), dateUnits)))
              ERR(retval);
           if ((retval = nc_put_att_text(ncid, time_id, "_CoordinateAxisType", 
                                         strlen("Time"), "Time")))
              ERR(retval);


   /* #### Attributes of the model output variable. #### */
   /* Assign "long_name" attribute to the variable - to do this right, need to pass the variable's description into here. */
           if ((retval = nc_put_att_text(ncid, varid, "long_name", 
                                         strlen(VARNAME), VARNAME)))
              ERR(retval);
           if ((retval = nc_put_att_text(ncid, varid, "units", 
                                         strlen(thisvarUnits), thisvarUnits)))
              ERR(retval);
           if ((retval = nc_put_att_float(ncid, varid, "missing_value",  /* note that "_FillValue" is updated attribute, but EverVIEW didn't respond to that, so using deprecated "missing_value" */
                                         NC_FLOAT, 1, &offMap_float)))
              ERR(retval);
           if ((retval = nc_put_att_int(ncid, varid, "OutputInterval_inDays", 
                                         NC_INT, 1, &outstep)))
              ERR(retval);
           if ((retval = nc_put_att_text(ncid, varid, "coordinates", 
                                         strlen("x y"), "x y")))
              ERR(retval);
           if ((retval = nc_put_att_text(ncid, varid, "grid_mapping", 
                                         strlen("transverse_mercator"), "transverse_mercator")))
              ERR(retval);

/* for compatability with other south Florida models, we switched (only the netCDF) output to NAD83 (in the above spatial loop) */
           sprintf(attribute_string,"PROJCS[\"NAD_1983_UTM_Zone_17N\",GEOGCS[\"GCS_North_American_1983\",DATUM[\"D_North_American_1983\",SPHEROID[\"GRS_1980\",6378137.0,298.257222101]],PRIMEM[\"Greenwich\",0.0],UNIT[\"Degree\",0.0174532925199433]],PROJECTION[\"Transverse_Mercator\"],PARAMETER[\"False_Easting\",500000.0],PARAMETER[\"False_Northing\",0.0],PARAMETER[\"Central_Meridian\",-81.0],PARAMETER[\"Scale_Factor\",0.9996],PARAMETER[\"Latitude_Of_Origin\",0.0],UNIT[\"Meter\",1.0]]");

           if ((retval = nc_put_att_text(ncid, varid, "esri_pe_string", 
                                         strlen(attribute_string), attribute_string)))
              ERR(retval);

   /* #### Attributes of the UTM coordinate variables. #### */
           if ((retval = nc_put_att_text(ncid, x_id, "long_name", 
                                         strlen("easting (x) coordinate"), "easting (x) coordinate")))
              ERR(retval);
           if ((retval = nc_put_att_text(ncid, y_id, "long_name", 
                                         strlen("northing (y) coordinate"), "northing (y) coordinate")))
              ERR(retval);
           if ((retval = nc_put_att_text(ncid, x_id, "standard_name", 
                                         strlen("projection_x_coordinate"), "projection_x_coordinate")))
              ERR(retval);
           if ((retval = nc_put_att_text(ncid, y_id, "standard_name", 
                                         strlen("projection_y_coordinate"), "projection_y_coordinate")))
              ERR(retval);
           if ((retval = nc_put_att_text(ncid, x_id, "units", 
                                         strlen("m"), "m")))
              ERR(retval);
           if ((retval = nc_put_att_text(ncid, y_id, "units", 
                                         strlen("m"), "m")))
              ERR(retval);
           if ((retval = nc_put_att_text(ncid, x_id, "_CoordinateAxisType", 
                                         strlen("GeoX"), "GeoX")))
              ERR(retval);
           if ((retval = nc_put_att_text(ncid, y_id, "_CoordinateAxisType", 
                                         strlen("GeoY"), "GeoY")))
              ERR(retval);
           if ((retval = nc_put_att_text(ncid, x_id, "CellRegistration", 
                                         strlen("coords are at cell centroids"), "coords are at cell centroids")))
              ERR(retval);
           if ((retval = nc_put_att_text(ncid, y_id, "CellRegistration", 
                                         strlen("coords are at cell centroids"), "coords are at cell centroids")))
              ERR(retval);

     /* #### Attributes of the transverse_mercator variable (done only for EverVIEW application). #### */
           if ((retval = nc_put_att_text(ncid, TM_id, "grid_mapping_name", 
                                         strlen("transverse_mercator"), "transverse_mercator")))
              ERR(retval);
           if ((retval = nc_put_att_text(ncid, TM_id, "longitude_of_central_meridian", 
                                         strlen(" -81."), " -81.")))
              ERR(retval);
           if ((retval = nc_put_att_text(ncid, TM_id, "latitude_of_projection_origin", 
                                         strlen("0."), "0.")))
              ERR(retval);
           if ((retval = nc_put_att_text(ncid, TM_id, "scale_factor_at_central_meridian", 
                                         strlen("0.9996"), "0.9996")))
              ERR(retval);
           if ((retval = nc_put_att_text(ncid, TM_id, "false_easting", 
                                         strlen("500000."), "500000.")))
              ERR(retval);
           if ((retval = nc_put_att_text(ncid, TM_id, "false_northing", 
                                         strlen("0."), "0.")))
              ERR(retval);
           if ((retval = nc_put_att_text(ncid, TM_id, "_CoordinateTransformType", 
                                         strlen("Projection"), "Projection")))
              ERR(retval);
           if ((retval = nc_put_att_text(ncid, TM_id, "_CoordinateAxisTypes", 
                                         strlen("GeoX GeoY"), "GeoX GeoY")))
              ERR(retval);
           if ((retval = nc_put_att_text(ncid, TM_id, "units", 
                                         strlen("m"), "m")))
              ERR(retval);


 /* #### Global attributes. #### */
           if ((retval = nc_put_att_text(ncid, NC_GLOBAL, "title", 
                                         strlen("Everglades Landscape Model output"), "Everglades Landscape Model output")))
              ERR(retval);
           if ((retval = nc_put_att_text(ncid, NC_GLOBAL, "ModelName", 
                                         strlen(modelName), modelName)))
              ERR(retval);
           if ((retval = nc_put_att_text(ncid, NC_GLOBAL, "ModelVersion", 
                                         strlen(modelVers), modelVers)))
              ERR(retval);
           if ((retval = nc_put_att_text(ncid, NC_GLOBAL, "ScenarioType", 
                                         strlen(SimAlt), SimAlt)))
              ERR(retval);
           if ((retval = nc_put_att_text(ncid, NC_GLOBAL, "Scenario", 
                                         strlen(SimModif), SimModif)))
              ERR(retval);
           if ((retval = nc_put_att_text(ncid, NC_GLOBAL, "Day0_SimulationPeriod", 
                                         strlen(initDateRead), initDateRead)))
              ERR(retval);
           if ((retval = nc_put_att_text(ncid, NC_GLOBAL, "Time_model_run_started", 
                                         strlen(asctime(local_begin_time)), asctime(local_begin_time))))
              ERR(retval);
           /* todo get rid of this hardcoding */
           sprintf(attribute_string,"H. Carl Fitz, carlfitz3@gmail.com, EcoLandMod, Inc, http://www.ecolandmod.com");

           if ((retval = nc_put_att_text(ncid, NC_GLOBAL, "modelDeveloper", 
                                         strlen(attribute_string), attribute_string)))
              ERR(retval);


   /* #### End define mode. This tells netCDF we are done defining metadata. #### */
           if ((retval = nc_enddef(ncid)))
              ERR(retval);



/* Still on the first outstep, write the coordinate variable data. This will put the northings
      and eastings of our data grid into the netCDF file. */
           if ((retval = nc_put_var_float(ncid, y_id, &northrev[0])))  /* note the use of reverse rows */
              ERR(retval);
           if ((retval = nc_put_var_float(ncid, x_id, &east[0])))
              ERR(retval);

           if ((retval = nc_close(ncid)))
              ERR(retval);


        } /* end of all defs. Done once, prior to first writing of data; have a closed file */

   /* These settings tell netcdf to write one timestep of data. (The
     value of start[0] below tells netCDF which (index)
     timestep to write.) */
   count[0] = 1;
   count[1] = s0+2; /* was s0+2 */
   count[2] = s1+2; /* was s1+2 */
   start[1] = 0;
   start[2] = 0;

        start[0] = index; /* tried to use 'SimTime.TIME' , but it filled in multiple blank days */
        /* now open the ncdf file, and write to it */
   if ((retval = nc_open(filename, NC_WRITE, &ncid)))
      ERR(retval);
   if ((retval = nc_inq_varid(ncid, VARNAME, &varid)))
      ERR(retval);
   if ((retval = nc_inq_varid(ncid, "time", &time_id)))  
      ERR(retval);

   /* Write the data to the file.  */
   if ((retval = nc_put_vara_float(ncid, varid, start, count, nc_dataRowRev)))
      ERR(retval);
   if ((retval = nc_put_var1_float(ncid, time_id,  start, &SimTime.TIME )))
      ERR(retval);

   /* Close the file. This frees up any internal netCDF resources
    * associated with the file, and flushes any buffers. */
   if ((retval = nc_close(ncid)))
      ERR(retval);

/* end of netCDF code */
#else

    printf("\nError in configuring 4-byte output request, you asked for netCDF ('C' argument in 'G' command) output, but your globals.h file has NCDF3 set to false: output variable= %s\n",filename);
    exit(0); 

#endif
}



void  write_header(char* mapFileName,int size)
{
  FILE* header;
  char headFIleName[320];
  int outrow,outcol;
  /* v2.8.2 modified - GRASS format header */
  sprintf(headFIleName,"%s_header.txt",mapFileName);
  fprintf(dFile,"Writing Header: filename = %s\n",headFIleName); 
  fflush(dFile);
  header = fopen(headFIleName,"w");
  if(header==NULL) { 
    fprintf(stderr,"Can't open header file: %s\n",headFIleName); 
    fflush(stderr); exit(-1); 
  }
  if(size == 4) { /* v2.8.2 - size of float arrays are 2 larger (both row and col) than the standard output */
        outrow = gbl_size[0]+2;
        outcol = gbl_size[1]+2;
  }
  else {
        outrow = gbl_size[0];
        outcol = gbl_size[1];
  }
  /* GRASS header - note that the "format" byte size is the size - 1 .  Also, the UTM zone is hardcoded to 17 here - lazy (unused anywhere in model) */
  fprintf(header,"proj:       1\nzone:       17\nnorth:      %.0f\nsouth:      %.0f\neast:       %.0f\nwest:       %.0f\n", 
                  UTMnorth, UTMsouth, UTMeast, UTMwest);
  fprintf(header,"cols:       %d\nrows:       %d\ne-w resol:  %.0f\nn-s resol:  %.0f\nformat:     %d\ncompressed: 0\n",
                  outcol,outrow,celWid,celWid,(size-1) ); 
  fclose(header);
}



float* readSeriesCol(char *filename, int format, int index, int* Npt, float* TStep, int col )
{
    int line = 0, cread=1, itest, j, sLen = 372 ;
    FILE *cfile;
    char  ctest, cmark[373], ret = '\n';
    unsigned char cmove[373], cmv=0;
    char tfilename[200], date_read[20],ss[200];
    int yyyy,mm, dd;
    double Jdate_read; 
    float  *tvar, *nullvar, first[10];

    sprintf(tfilename,"%s%d.ts",filename,index);
    if(debug>4) fprintf(dFile,"\nReading file %s, col %d, index = %d\n",tfilename,col,index); fflush(dFile);
    cfile = fopen(tfilename,"r"); 
    if(cfile==NULL) {  
        if( index > 0 ) {  fprintf(stdout,"\nWARNING: Unable to open timeseries file %s, using %s0.ts\n",tfilename,filename); return 0; } 
        else { fprintf(stdout,"\nERROR: Unable to open timeseries file %s\n",tfilename); exit(-1); } 
    }

    if (format != 2) { fprintf(stderr,"ERROR: only using floats in read_timeseries: %d\n",format); exit(0); }
    
           sLen = PORnumday; /* read/determined in genericDriver */
            tvar = (float*) nalloc( (sLen+2)*sizeof(float), "timeseries temp" );
            scan_forward(cfile,"{END COMMENT}\n");
            fgets(ss,200,cfile); /* skip  the header line */
            
            do {
                fscanf(cfile, "%d %d %d %f %f %f %f", &yyyy,&mm,&dd,&first[1],&first[2],&first[3],&first[4] ); 
        /* julian day, returns a double (args = month, day, year, hours, minutes, seconds) */
                Jdate_read = julday(mm, dd, yyyy, 0, 0, 0); 
                if (Jdate_read > Jdate_init) {printf (" \n***Error\nInit date in %s file > simulation start date.\n",tfilename); exit (-1);}
                if (feof(cfile)) {
        printf (" \n***Error\nNon-matching starting date for %s meteorological time series\n",tfilename); 
                        exit (-1);
                }
            } while (Jdate_read < Jdate_init);
            
            tvar[0] = first[col]; /* have now read the first data record, assign to array */
            if(debug>4) fprintf(dFile,"%d %d %d %f\n",yyyy,mm,dd,tvar[0]);
            for (line = 1; line<sLen; line++) {
            
            /* each case for reading a different column */
                switch (col) {
                    case 1:
                        itest = fscanf(cfile,"%d %d %d %f %f %f %f",&yyyy,&mm,&dd, &(tvar[ line ]),&nullvar, &nullvar, &nullvar);
                        break;
                    case 2:
                        itest = fscanf(cfile,"%d %d %d %f %f %f %f",&yyyy,&mm,&dd,&nullvar, &(tvar[ line ]), &nullvar, &nullvar);
                        break;
                    case 3:
                        itest = fscanf(cfile,"%d %d %d %f %f %f %f",&yyyy,&mm,&dd,&nullvar, &nullvar, &(tvar[ line ]), &nullvar);
                        break;
                    case 4:
                        itest = fscanf(cfile,"%d %d %d %f %f %f %f",&yyyy,&mm,&dd,&nullvar, &nullvar, &nullvar, &(tvar[ line ]));
                        break;
                    default:
                        printf ("\nError in interpolation time series data.\n"); 
                        exit(-1);
                        break;
                } /* end switch */

                if(debug>4) fprintf(dFile,"%d %d %d %f\n",yyyy,mm,dd,tvar[ line ]);
                if (feof(cfile)) {printf (" \n***Error\nExceeded number of records for %s meteorological time series\n",tfilename);  exit (-1); } 
            }
            tvar[ sLen ] = *TStep;
 
    if(debug>4) { fprintf(dFile,"\nDone Reading file %s\n",tfilename); fflush(dFile); }
    
    fflush(stdout);
    fclose(cfile);
    *Npt = sLen;
    return tvar;
}




float* get_hab_parm( char* s_parm_name, int s_parm_relval, char* parmName ) 
{
        int hIndex, end=0;
        FILE* parmFile;
        float* fTmp;
        char lline[2001], *line;
        int pLoc = -1;
        int foundParm= False;
        char parmNameHead[30];
        char modelFileName[300];
    char HabParmFile[300];
    char *fileAppend;
    extern ProgAttr *ProgExec;
    
        
        if ( (fTmp = (float*) malloc( MAX_NHAB * sizeof(float) ) ) == NULL) {
        sprintf(msgStr, "Failed to allocate memory for HabParm array.\n ") ;
        usrErr(msgStr);
        exit( -2 ) ;
        };
        fTmp[0] = 0.0;

    fileAppend = match_Sparm( s_parm_name, s_parm_relval, parmName);
        sprintf(HabParmFile,"HabParms_%s",fileAppend);
        sprintf(modelFileName,"%s/%s/Data/%s",ModelPath,ProjName,HabParmFile);
        parmFile = fopen(modelFileName,"r");
        if(parmFile==NULL) {fprintf(stdout,"\nERROR: Unable to open dBase file %s\n",modelFileName); exit(-1); }
        
     fgets(lline, 2000, parmFile); /* skip 1st header line  */
     fgets(lline, 2000, parmFile); /* read 2nd header line with column names */
     line=lline;

     while (!foundParm)  
     {
        sscanf( line,"%s",&parmNameHead);
        if (strcmp(parmName,parmNameHead) ==0) foundParm = True;
        pLoc++;
        if (*line == EOS || *line == EOL) {
           sprintf(msgStr,"ERROR: Could not find header tag of %s in HabParms dbase file.", parmName); 
           usrErr(msgStr); 
           exit(-1);
        }
        line = Scip( line, TAB );
     }
    
    sprintf(msgStr,"Parameter group: %s",parmName); 
    fflush(dFile); 
    WriteMsg(msgStr,1); 
        for( hIndex = 1; hIndex < MAX_NHAB; hIndex++) {
                goto_index (parmFile, gHabChar, hIndex); 
                fTmp[hIndex] = get_Nth_parm( parmFile, pLoc, &end, hIndex );
                if(end) break;
        }
        fclose(parmFile);
        habNumTot = hIndex-1;

        fflush(dFile); 
        /* TODO: this HARDCODES into the BIR (BIRavg) model output file ONLY, 
        the incorrect impression of running sensitivity on a single habitat ID 
        (hab 2 = sawgrass plain, ELMv2.4) */
    if (strcmp(s_parm_name,parmName) == 0) ProgExec->S_ParmVal = fTmp[2];
        return fTmp; 
}

char* match_Sparm( char* s_parm_name, int s_parm_relval, char* parmName) 
{
    char *fileAppend;
 
    if (strcmp(s_parm_name,parmName) == 0) /* match, it is a parameter being varied in this set of runs */
        { 
                switch (s_parm_relval) {
                   case 0: fileAppend = "NOM"; break;
                   case 1: fileAppend = "LO"; break;
                   case 2: fileAppend = "HI"; break;
                   default:  {
                      sprintf(msgStr,"ERROR: The relative parameter value (%d) is unknown for %s",s_parm_relval, s_parm_name); 
                      usrErr(msgStr);
                      exit(-1);
                   }
           }
        }
    else
        { /* no match, it is NOT a parameter being varied in this set of runs */
              fileAppend = "NOM"; 
            }
    return (fileAppend);
}


float get_global_parm( char* s_parm_name, int s_parm_relval, char* parmName) 
{
  int test;
  char modelFileName[300];
  FILE *parmFile;
  float parmVal;
  char parmNameHead[30];
  extern ProgAttr *ProgExec;


  char GlobalParmFile[50];
  char *fileAppend;
  
  fileAppend = match_Sparm( s_parm_name, s_parm_relval, parmName);
  sprintf(GlobalParmFile,"GlobalParms_%s",fileAppend);
  sprintf(modelFileName,"%s/%s/Data/%s",ModelPath,ProjName,GlobalParmFile);
  
  parmFile = fopen(modelFileName,"r");
  if(parmFile==NULL) { 
        sprintf(msgStr,"ERROR, can't open data file %s",modelFileName); 
    usrErr(msgStr);
        Exit(0); 
  }

  sprintf(parmNameHead,"%s=", parmName);
  scan_forward(parmFile,parmNameHead);
  if ( (test = fscanf(parmFile,"%f",&parmVal) ) <= 0 ) { 
    sprintf(msgStr,"ERROR in reading %s from %s; see Driver0.out debug file.", parmName, modelFileName); 
    usrErr(msgStr);
    Exit(0); 
  }
  sprintf(msgStr,"%s %f\n",parmNameHead, parmVal); 
  WriteMsg(msgStr,1); 

  fclose(parmFile);
  if (strcmp(s_parm_name,parmName) == 0) ProgExec->S_ParmVal = parmVal;
  return parmVal;
}

float get_modexperim_parm( char* s_parm_name, int s_parm_relval, char* parmName) 
{
  int test;
  char modelFileName[300];
  FILE *parmFile;
  float parmVal;
  char parmNameHead[30];
  extern ProgAttr *ProgExec;


  char ModExParmFile[50];
  char *fileAppend;
  
  fileAppend = match_Sparm( s_parm_name, s_parm_relval, parmName);
  sprintf(ModExParmFile,"ModExperimParms_%s",fileAppend);
  sprintf(modelFileName,"%s/%s/Data/%s",ModelPath,ProjName,ModExParmFile);
  
  parmFile = fopen(modelFileName,"r");
  if(parmFile==NULL) { 
        sprintf(msgStr,"ERROR, can't open data file %s",modelFileName); 
    usrErr(msgStr);
        Exit(0); 
  }

  sprintf(parmNameHead,"%s=", parmName);
  scan_forward(parmFile,parmNameHead);
  if ( (test = fscanf(parmFile,"%f",&parmVal) ) <= 0 ) { 
    sprintf(msgStr,"ERROR in reading %s from %s; see Driver0.out debug file.", parmName, modelFileName); 
    usrErr(msgStr);
    Exit(0); 
  }
  sprintf(msgStr,"%s %f\n",parmNameHead, parmVal); 
  WriteMsg(msgStr,1); 

  fclose(parmFile);
  if (strcmp(s_parm_name,parmName) == 0) ProgExec->S_ParmVal = parmVal;
  return parmVal;
}

void open_point_lists(SeriesParm *pSeries, int nSeries)
{
  FILE *oFile;
  int i, j, ix, iy, k= 0 ;

  for (i=0; i<nSeries; i++) {
     if( pSeries[i].data == NULL ) {
        if ( numPtFiles>0 ) fclose(oFile); /* in (rare?) case where user asks for no point output, no files stream to close */
        return;
      }
      if (strcmp(pSeries[i].name,pSeries[i-1].name) != 0) {
              /* this is new variable, thus new file.  numPtFiles== increment counter of number of variables (one file per variable) used for output */
          if (i>0) { 
            numPtFiles++; 
            fclose (oFile); 
          }
          
          sprintf(msgStr,"%s%s/Output/PtSer/%s.pts",OutputPath,ProjName,pSeries[i].name);
          if( (oFile = fopen(msgStr,"w") ) == NULL) { 
              fprintf(stderr,"\nERROR, unable to open %s point time series output file.",msgStr); 
              exit(-1); 
          }
          fprintf (oFile, "%s %s %s %s scenario: \nDate\t", &modelName, &modelVers, &SimAlt, &SimModif);
          
      }
/* populate the grid cell locations in the header for this variable/file */
      ix = pSeries[i].Loc.x;
      iy = pSeries[i].Loc.y;
      fprintf(oFile,"%s:(%d,%d)\t",pSeries[i].name,ix,iy);
      cell_pts[numPtFiles]++; /* count the number of cell point locs per variable (1 var per file) */
      
  }
}


void send_point_lists2(SeriesParm *pSeries, int nSeries)
{
  FILE *oFile;
  int i=0, ii, j, ix, iy, k= 0, last_pt;
  int yr_pt[2],mo_pt[2],da_pt[2],hr_pt[2],mi_pt[2];    
  double se_pt[2];           
  double Jdate_pt;
  
  if (numPtFiles==0) return;
  /* k = file/variable counter; i = indiv cellLoc time series counter */
  for (k = 0; k <= numPtFiles; k++) { /* loop over the number of point ts files */
      sprintf(msgStr,"%s%s/Output/PtSer/%s.pts",OutputPath,ProjName,pSeries[i].name);
      if( (oFile = fopen(msgStr,"a") ) == NULL) { 
          fprintf(stderr,"\nERROR, unable to open %s point time series output file.",msgStr); 
          exit(-1); 
      }

      for(j=pSeries[i].laststep; j<pSeries[i].Length; j++ ) { /*temporal loop */
          Jdate_pt = Jdate_init+j*pSeries[i].outstep;
          calcdate( Jdate_pt, mo_pt, da_pt, yr_pt, hr_pt, mi_pt, se_pt); /* get the calendar date info from the julian date */
          fprintf(oFile,"\n%d/%d/%d\t",yr_pt[0],mo_pt[0],da_pt[0] ); /* calendar date in col 1 */

          last_pt = i+cell_pts[k]; /* # of last cellLoc point of the current variable's file */
          for (ii=i; ii<last_pt; ii++) {/* loop over number of point locations per file */
              fprintf(oFile,"%f\t",pSeries[ii].data[j]);
          }
      }
      fclose (oFile); 
      pSeries[i].laststep = j; /* remember the last temporal outstep for this file */
      i += cell_pts[k];     /* increment the indiv cellLoc time series counter */ 
  }
  return;
}


void writeSeries(void* fValue, char* label, char* desc, int N0, int N1, byte Mtype, byte format )
{
        /* check on these! (who wrote this? Why? (HCF Dec'04) )*/
  int ix, iy; 
  static unsigned char first_write = 1;
  long int  ret, dimsizes[2];
  unsigned short int refnum;
  if(format == 'H') {
  
/* NOTE: June 2008 (v2.8.2): do not #define HDF in globals.h to true until until the code here is updated - 
   this code needs to be updated to HDF4 before it will compile and be used */
#if HDF   
          dimsizes[0] = N0;
          dimsizes[1] = N1;
          sprintf(cTemp,"%s%s:Output:Windows.hdf",OutputPath,ProjName);
          DFSDclear();
          
          switch(Mtype) {
                 case 'f': 
                    DFSDsetNT(DFNT_FLOAT32);  break;
                 case 'L': 
                    DFSDsetNT(DFNT_FLOAT32);  break;
                 case 'E': 
                    DFSDsetNT(DFNT_FLOAT32);  break;
                 case 'd': case 'i': 
                    DFSDsetNT(DFNT_INT32);  break;
                 case 'c': 
                    DFSDsetNT(DFNT_UINT8);  break;
          }
          if(first_write) {
                ret = DFSDputdata(cTemp,2,dimsizes,fValue);
                first_write = 0;
          } else {
                ret = DFSDadddata(cTemp,2,dimsizes,fValue);
          }
          HDF_VERIFY("DFSDadddata");
          refnum = DFSDlastref();
          ret = DFANputlabel(cTemp,DFTAG_SDG,refnum,label);
          HDF_VERIFY("DFANputlabel");
          ret = DFANputdesc(cTemp,DFTAG_SDG,refnum,desc,strlen(desc));
          HDF_VERIFY("DFANputdesc");
#endif
  }
  else {
    fprintf(Driver_outfile,"\n_________%s\n",label);    
    fprintf(Driver_outfile,"%s\n",desc);        
    for (ix=0; ix<N0; ix++) {
      fprintf(Driver_outfile,"\n");
      for (iy=0; iy<N1; iy++) {
        switch(Mtype) {
          case 'f': 
                if(format=='L')                 fprintf(Driver_outfile,"%f ",*(((float*)fValue)+iy+ix*N1)); 
                else if(format=='E')    fprintf(Driver_outfile,"%.3E ",*(((float*)fValue)+iy+ix*N1)); 
                        else                            fprintf(Driver_outfile,"%.3f ",*(((float*)fValue)+iy+ix*N1));
                        break;
        case 'd': case 'i': 
                        fprintf(Driver_outfile,"%d ",*(((int*)fValue)+iy+ix*N1)); 
                        break;
        case 'c': 
                        fprintf(Driver_outfile,"%x ",*(((UCHAR*)fValue)+iy+ix*N1)); 
                        break;
        }
      }
    }  
  }
}

void Combine(float* fValue, char* label, int nComp, int cType, int step)
{
        int print, cIndex, i, cum;
        static int type = 11;
        char tmpStr[50];
        if( ++type > 99) type = 11;

        cum = ((cType == kAVECUM) || (cType == kSUMCUM)) ? 1 : 0;       
        if(cum) { cIndex = getCombineIndex(label,step,cType,&print); }
        
        switch(cType) {
                case kMAXMIN: 
                        sprintf(msgStr,"\nMAXMIN(%d) for %s:",step,label); 
                        for(i=0; i<nComp; i++) { sprintf(tmpStr," %f",fValue[i]); strcat(msgStr,tmpStr); }
                        WriteMsg(msgStr,1);
                break;
                case kMAX:
                        sprintf(msgStr,"\nMAX(%d) for %s:",step,label); 
                        for(i=0; i<nComp; i++) { sprintf(tmpStr," %f",fValue[i]); strcat(msgStr,tmpStr); }
                        WriteMsg(msgStr,1);
                break;
                case kMIN:
                        sprintf(msgStr,"\nMIN(%d) for %s:",step,label); 
                        for(i=0; i<nComp; i++) { sprintf(tmpStr," %f",fValue[i]); strcat(msgStr,tmpStr); }
                        WriteMsg(msgStr,1);
                break;
                case kSUM:
                        sprintf(msgStr,"\nSUM(%d) for %s:",step,label); 
                        for(i=0; i<nComp; i++) { sprintf(tmpStr," %f",fValue[i]); strcat(msgStr,tmpStr); }
                        WriteMsg(msgStr,1);
                break;
                case kAVE:
                        sprintf(msgStr,"\nAVE(%d) for %s:",step,label); 
                        for(i=0; i<nComp; i+=2) { sprintf(tmpStr," %f",fValue[i]/fValue[i+1]); strcat(msgStr,tmpStr); }
                        WriteMsg(msgStr,1);
                break;
                case kSUMCUM:
                        if( print == -1) {
                                for(i=0; i<nComp; i++) ctable[cIndex].fvalue[i] = fValue[i]; 
                        
                        } else {
                                for(i=0; i<nComp; i++) ctable[cIndex].fvalue[i] += fValue[i]; 
                        }
                        if(print==1) {                          
                                sprintf(msgStr,"\nSUMCUM(%d) for %s:",step,label); 
                                for(i=0; i<nComp; i++) { sprintf(tmpStr," %f",ctable[cIndex].fvalue[i]); strcat(msgStr,tmpStr); }
                                WriteMsg(msgStr,1);
                        }
                break;
                case kAVECUM:
                        if( print == -1) {
                                for(i=0; i<nComp; i++) ctable[cIndex].fvalue[i] = fValue[i]; 
                        
                        } else {
                                for(i=0; i<nComp; i++) ctable[cIndex].fvalue[i] += fValue[i]; 
                        }
                        if(print==1) {                          
                                sprintf(msgStr,"\nAVECUM(%d) for %s:",step,label); 
                                for(i=0; i<nComp; i+=2) { sprintf(tmpStr," %f",ctable[cIndex].fvalue[i]/ctable[cIndex].fvalue[i+1]); strcat(msgStr,tmpStr); }
                                WriteMsg(msgStr,1);
                        }
                break;
        }
}       

int getCombineIndex( char* name, int step,int type,int *last)
{
        int count, index = -1;
        *last = 0;
        
        for(count = 0; count < max_combos_open; count++) {
        
                if( ctable[count].free )        {                       /* Determine if the slot is in use. */
                        if(index == -1) index = count;          /* Nope, keep track of first free global stream table slot. */ 
                        if(debug) { fprintf(dFile,"cnt = %d, index = %d, free combo stream\n",count,index); fflush(dFile); }
                }
                else if( strncmp(name,ctable[count].name,25) == 0 && ctable[count].type == type ) {
                    *last = 1; 
                    if(debug) { 
                      fprintf(dFile,"\n(%s)combo cnt = %d, step = %d, type = %d\n",
                              name,count,step,type); 
                        fflush(dFile); 
                    }
                    return count;
                }
        }
        if( index == -1 ) { 
                index = max_combos_open; 
                if( ++max_combos_open > MAXCOMBOS) { sprintf(msgStr,"Out of combo slots: %s, step=%d, type=%d\n",name,step,type);  usrErr(msgStr); return -1; }
        }
        if(debug) { fprintf(dFile,"\n(%s)index = %d, step=%d, set-up combo\n",name,index,step); fflush(dFile); }
        strncpy(ctable[index].name,name,25);
        ctable[index].step = step;
        ctable[index].type = type;
        ctable[index].free = 0;
        *last = -1;
        return index;
}
void open_debug_outFile(int index)
{
  char filename[120];
  
  sprintf(filename,"%s/%s/Output/Debug/Driver%d.out",OutputPath,ProjName,index);

  Driver_outfile = fopen(filename,"w");
  if(Driver_outfile == NULL) { 
        fprintf(stderr,"Error, unable to open %s file.\n",filename);
        fprintf(stderr,"OutputPath: %s\n",OutputPath);  
        fprintf(stderr,"Project: %s\n",ProjName);  
        exit(0); 
  }

}

int init_config_file(FILE *vpFile, char term1, char term2, char term3, char term4)
{
  char test;
  int format,size=-1;

  gTerm[0] = term1;
  gTerm[1] = term2;  
  gTerm[2] = term3;  
  gTerm[3] = term4;  
  skip_white(vpFile);
  test = fgetc(vpFile);
  if(test == gTerm[0])
      fscanf(vpFile,"%d",&format);
  else
      format = -1;
  return format;
}

int skip_cnfg_space(FILE *vpFile, char* tch)
{
        char ch; int rv = 1;
        ch = *tch;
        while( isspace(ch) ) { 
                if( (ch=fgetc(vpFile)) == EOF || ch == gTerm[3] || ch == '\0' ) {
                        fclose(vpFile); cnfgFile = NULL; return -3; 
                }
        }
        *tch = ch;
        return rv;
}

int parse_packet( FILE *vpFile, int* nArgs, char* test) 
{
    static      int     gVPindex = -1;
    static      char    gCnfg;
    
        char ch = ' ', eChar = ' '; 
        int btype=0, argc=2, go=1, i=0, j=0, new_name = 0;

        if(vpFile == NULL) return -3;
        while (1) {
                if ( skip_cnfg_space(vpFile, &ch) < 0) return -3; j=0;
                if( ch == gTerm[1] || ch == gTerm[2] ) {
                    gCnfg = ch; ch = ' ';
                        if ( skip_cnfg_space(vpFile, &ch) < 0 ) return -3;
                        while ( !isspace(ch) ) { cTemp[j++] = ch; ch=fgetc(vpFile); }
                        cTemp[j++] = '\0'; new_name = 1;
                        gVPindex++;
                }
                else break;
        }
        strcpy(gCArg[0],cTemp);
        *test = gCnfg;
        while ( isalnum( ch ) ) { gCArg[1][i++]=ch; ch=fgetc(vpFile); }
        gCArg[1][i]='\0'; i=0;

        while( 1 ) {
                if( ch == '(' ) { eChar = ')'; break;} 
                else if( ch == '[' ) { eChar = ']'; break;} 
                else if( ch == '{' ) { eChar = '}'; break;} 
                else if( !isspace(ch) ) { return -3;}
                ch=fgetc(vpFile);
        }
        while (go) {
                while( ch=fgetc(vpFile) ) {
                        if( ch == ',' ) { argc++; break; }
                        if( ch == ')' ||  ch == '}' ||  ch == ']')  { 
                                if( ch == eChar ) {  argc++; go=0; break;} 
                                else { printf( "\nWarning: Syntax error in parse_config, var: %s\n",gCArg[0]); return -2; } 
                        }
                        gCArg[argc][i++]=ch;
                        if( i==(kCArgWidth-1) ) break;
                }
                if(i==0) { printf( "\nWarning: Syntax error in parse_config, var: %s\n",gCArg[0]); return -2; } 
                else { gCArg[argc-1][i]='\0'; i=0; }
                if(argc == kCArgDepth) { 
                        while( ch != eChar ) ch=fgetc(vpFile);  
                        go = 0; 
                }
        }
        *nArgs = argc;
        if(argc==0) { printf( "\nWarning: Syntax error in parse_config, var: %s\n",gCArg[0]); return -2; }
        else return gVPindex;           
} 


int get_number( FILE *infile ) {
        char ch;
        int rv;
        
        ch = fgetc(infile);
        if( !isdigit( ch ) ) return(-1);
        rv = ch - '0';
        
        while( isdigit( ch = fgetc(infile) ) ) 
                rv = rv*10 + ( ch - '0' );
                
        return rv;
}

int goto_index ( FILE *infile, char tchar, int index) 
{
        int rv=1, current_index=-1, itest;

        while(current_index != index) { 
                itest = find_char (infile, tchar);
                if(itest <= 0) return 0;
                current_index = get_number(infile);
                if(current_index < 0) fatal("Bad number format in dBase");
        }
        return 1;
}

float get_Nth_parm( FILE *infile, int pIndex, int* end, int hIndex ) 
{
        int i=1, itest=1;
        float rv; char ch;
        
        while ( (i++ < pIndex) && itest > 0 ) { itest = find_char ( infile, '\t' ); }
        if( itest == 0 ) { *end = 1; return(0.0); }
        
        itest = fscanf(infile,"%f",&rv);
        if(itest==1) { 
        if(debug) { 
           sprintf(msgStr,"Habitat %d:\t%f",hIndex,rv); 
           WriteMsg(msgStr,1); 
        }
        }
        else if ( (ch = fgetc(infile)) == EOF ) { *end = 1; return(0.0); }
        else {
                sprintf(msgStr,"Read Error in dBase(%d)\n REad Dump:\n",itest); WriteMsg(msgStr,1); usrErr(msgStr);
                for(i=0; i<12; i++) { 
                        ch = fgetc(infile);
                        if(ch==EOF) { sprintf(msgStr,"\nAt EOF\n"); WriteMsg(msgStr,1); break; }
                        else fputc(ch,dFile);
                }
                exit(0);
        }
        *end = 0;
        return(rv);
}


float SMDRAND(float fminVal, float fmaxVal )
{ 
  float rv, tmp;
  tmp = rand();
  rv = fminVal + (fmaxVal-fminVal)*(tmp/RAND_MAX); 
  return(rv);
} 


void local_setup(int argc, char** argv)
{
  int i;
  char debugfileName[300];
  nprocs[0] = 1;
  nprocs[1] = 1;
  tramType =  0;
  procnum = 1;
  recpnum[0] = 0;
  recpnum[1] = 0;
  Lprocnum = 1;
  /* the combo_table (ctable) used in Combine operation for spatial array summaries (debug-related) */ 
  max_combos_open = 0;  
  for(i = 0; i < MAXCOMBOS; i++) {
    ctable[i].ocnt = 0; 
    ctable[i].free = 1;
  }

  /* this debug output file is not used extensively in current code - it is opened before we get
     environment vars, thus is written to user's current directory ($ModelPath/$ProjName/Load if
     ELM is executed from the "go" script) */ 
 dFile = fopen("ELM.debug","w");
  if(dFile == NULL) { 
    usrErr("Can't open ELM.debug file."); 
    exit(0); 
  } 
  /* there are no random processes in current ELM (v2.4) w/o fire */ 
  srand(seed);
  fprintf(dFile," RAND_MAX = %d\n",(int)RAND_MAX);
}

/*****/
/* begin parallel, effectively unused, code (ex== EXpress parallel compiler compatibility) */
/*****/

void exparam( struct nodenv* envInfo)
{       
        envInfo->procnum = 1;
        envInfo->nprocs = 1;
        envInfo->groupid = 0;
        envInfo->taskid = 0;
 
}

int  exgridinit(int dim, int* nprocs)   
{ return 0;}

void exgridsplit(int nprocs, int ndim,int nprocs2[2])
{
  nprocs2[0] = 1;
  nprocs2[1] = 1;
}

void exgridcoord(int pnum, int rnum[2])
{
    tramNum[0] = rnum[0] = 0;
    tramNum[1] = rnum[1] = 0;
}

void exgridsize( int pnum, int gsize[2],int lsize[2],int lstart[2])
{
        int rem[2], i, j;
        for( i=0; i<2; i++ ) {
                lsize[i] = gsize[i]/nprocs[i];
                rem[i] = gsize[i] - lsize[i]*nprocs[i];
                if(recpnum[i]<rem[i]) lsize[i]++;
                for(j=0; j<recpnum[i]; j++) if( j<rem[i] && recpnum[i] >= rem[i] ) lstart[i] += (lsize[i]+1); else lstart[i] += lsize[i];
        }
}

void set_async_mode(FILE* file) 
{ return;}

void fmulti(FILE* file) 
{ return;}

void fsingl(FILE* file) 
{ return;}

void fasync(FILE* file) 
{ return;}


void exchange_borders(UCHAR* map, int size) 
{return;}

int on_this_proc(int x,int y)
{ return 1; }

void Cplot(VOIDP Map, unsigned char Mtype, float max_value, float min_value) 
{}

void broadcastMsg( UCHAR* msgPtr)
{}

void broadcastInt(int* iValPtr )
{}

void broadcastChar( UCHAR* cPtr)
{}

void broadcastData( void* dataPtr, int *dataSize)
{}

void sync_processors() 
{return;}

void broadcastFloat(void* dataPtr)
{}

/*****/
/* end parallel, effectively unused, code (ex== EXpress parallel compiler compatibility) */
/*****/

/****************************************************/