Changes to directory structure in git repository, also added git tools directory

Author: Joe Hamman

readme.txt

@@ -0,0 +1,19 @@
+The following programs are found in the basin_prep.tgz tarball. 
+
+multiply_basin_param.c    
+
+	This program was written to read a runfile mask and multiply
+  	parameters in the activated grid cells by the given value.
+ 
+set_basin_param.c  
+
+  	This program was written to read a runfile mask and replaces
+  	parameters in the activated grid cells with the given value.
+
+sum_basin_param.c  
+
+  	This program was written to read a runfile mask and computes
+  	the sum and average parameter values for the activated grid
+  	cells.
+
+

global.param.sample renamed to samples/global.param.sample

@@ -0,0 +1,66 @@
+#include <stdio.h>
+#include <math.h>
+
+main(int argc, char *argv[]) {
+/****************************************************************
+  multiply_basin_param.c    Keith Cherkauer    August 14, 1998
+
+  This program was written to read a runfile mask and multiply
+  parameters in the activated grid cells by the given value.
+****************************************************************/
+
+  FILE *frun, *foldparam, *fnewparam;
+  int   i, j, ncol, nrow;
+  int   runflag;
+  char  tmpstr[512];
+  float multfact;
+  float param;
+
+  if(argc!=5) {
+    fprintf(stderr,"Usage: %s <run mask> <old param file> <new param file> <factor>\n",argv[0]);
+    fprintf(stderr,"\tThis program multiplies the parameters for the activated grid cells in the old parameter file by the given value, and writes the results to the new parameter file.  Activated cells have values other than 0 in the run mask file.\n");
+    exit(0);
+  }
+
+  if((frun=fopen(argv[1],"r"))==NULL) {
+    fprintf(stderr,"ERROR: Unable to open %s\n",argv[1]);
+    exit(0);
+  }
+
+  if((foldparam=fopen(argv[2],"r"))==NULL) {
+    fprintf(stderr,"ERROR: Unable to open %s\n",argv[2]);
+    exit(0);
+  }
+
+  if((fnewparam=fopen(argv[3],"w"))==NULL) {
+    fprintf(stderr,"ERROR: Unable to open %s\n",argv[3]);
+    exit(0);
+  }
+
+  multfact = atof(argv[4]);
+
+  /** Skip and copy header **/
+  for(i=0;i<6;i++) {
+    fgets(tmpstr,512,frun);
+    fgets(tmpstr,512,foldparam);
+    if(i==0) sscanf(tmpstr,"%*s %i",&ncol);
+    if(i==1) sscanf(tmpstr,"%*s %i",&nrow);
+    fprintf(fnewparam,"%s",tmpstr);
+  }
+
+  for(i=0;i<nrow;i++) {
+    for(j=0;j<ncol;j++) {
+      fscanf(frun,"%i",&runflag);
+      fscanf(foldparam,"%f",&param);
+      if(runflag!=0) {
+	fprintf(fnewparam,"%f",param*multfact);
+      }
+      else {
+	fprintf(fnewparam,"%f",param);
+      }
+      if(j==ncol-1) fprintf(fnewparam,"\n");
+      else fprintf(fnewparam," ");
+    }
+  }
+}
+

output.LDAS_OUTPUT.template renamed to samples/output.LDAS_OUTPUT.template

@@ -0,0 +1,66 @@
+#include <stdio.h>
+#include <math.h>
+
+main(int argc, char *argv[]) {
+/****************************************************************
+  set_basin_param.c        Keith Cherkauer       August 14, 1998
+
+  This program was written to read a runfile mask and replaces
+  parameters in the activated grid cells with the given value.
+****************************************************************/
+
+  FILE *frun, *foldparam, *fnewparam;
+  int   i, j, ncol, nrow;
+  int   runflag;
+  char  tmpstr[512];
+  float value;
+  float param;
+
+  if(argc!=5) {
+    fprintf(stderr,"Usage: %s <run mask> <old param file> <new param file> <value>\n",argv[0]);
+    fprintf(stderr,"\tThis program sets the activated grid cells in the old parameter file to the given value, and writes the results to the new parameter file.  Activated cells have values other than 0 in the run mask file.\n");
+    exit(0);
+  }
+
+  if((frun=fopen(argv[1],"r"))==NULL) {
+    fprintf(stderr,"ERROR: Unable to open %s\n",argv[1]);
+    exit(0);
+  }
+
+  if((foldparam=fopen(argv[2],"r"))==NULL) {
+    fprintf(stderr,"ERROR: Unable to open %s\n",argv[2]);
+    exit(0);
+  }
+
+  if((fnewparam=fopen(argv[3],"w"))==NULL) {
+    fprintf(stderr,"ERROR: Unable to open %s\n",argv[3]);
+    exit(0);
+  }
+
+  value = atof(argv[4]);
+
+  /** Skip and copy header **/
+  for(i=0;i<6;i++) {
+    fgets(tmpstr,512,frun);
+    fgets(tmpstr,512,foldparam);
+    if(i==0) sscanf(tmpstr,"%*s %i",&ncol);
+    if(i==1) sscanf(tmpstr,"%*s %i",&nrow);
+    fprintf(fnewparam,"%s",tmpstr);
+  }
+
+  for(i=0;i<nrow;i++) {
+    for(j=0;j<ncol;j++) {
+      fscanf(frun,"%i",&runflag);
+      fscanf(foldparam,"%f",&param);
+      if(runflag!=0) {
+	fprintf(fnewparam,"%f",value);
+      }
+      else {
+	fprintf(fnewparam,"%f",param);
+      }
+      if(j==ncol-1) fprintf(fnewparam,"\n");
+      else fprintf(fnewparam," ");
+    }
+  }
+}
+

output.OPTIMIZE.template renamed to samples/output.OPTIMIZE.template

@@ -0,0 +1,66 @@
+#include <stdio.h>
+#include <math.h>
+
+main(int argc, char *argv[]) {
+/****************************************************************
+  sum_basin_param.c        Keith Cherkauer       August 14, 1998
+
+  This program was written to read a runfile mask and computes
+  the sum and average parameter values for the activated grid 
+  cells.
+
+****************************************************************/
+
+  FILE *frun, *fparam;
+  int   i, j, ncol, nrow;
+  int   runflag;
+  int   N;
+  char  tmpstr[512];
+  float factor;
+  float param;
+  float sum;
+
+  if(argc!=4) {
+    fprintf(stderr,"Usage: %s <run mask> <param file> <factor>\n",argv[0]);
+    fprintf(stderr,"\tThis program sums the parameters divided by factor in the activated grid cells, and prints the results to stdout.  Activated cells have values other than 0 in the run mask file.\n");
+    exit(0);
+  }
+
+  if((frun=fopen(argv[1],"r"))==NULL) {
+    fprintf(stderr,"ERROR: Unable to open %s\n",argv[1]);
+    exit(0);
+  }
+
+  if((fparam=fopen(argv[2],"r"))==NULL) {
+    fprintf(stderr,"ERROR: Unable to open %s\n",argv[2]);
+    exit(0);
+  }
+
+  factor = atof(argv[3]);
+
+  /** Skip and copy header **/
+  for(i=0;i<6;i++) {
+    fgets(tmpstr,512,frun);
+    fgets(tmpstr,512,fparam);
+    if(i==0) sscanf(tmpstr,"%*s %i",&ncol);
+    if(i==1) sscanf(tmpstr,"%*s %i",&nrow);
+  }
+
+  sum = 0.;
+  N   = 0;
+  for(i=0;i<nrow;i++) {
+    for(j=0;j<ncol;j++) {
+      fscanf(frun,"%i",&runflag);
+      fscanf(fparam,"%f",&param);
+      if(runflag!=0) {
+	sum += param / factor;
+	N ++;
+      }
+    }
+  }
+
+  fprintf(stdout,"Parameter Sum = %f\nParameter Average = %f\n",
+    sum,sum/(float)N);
+
+}
+

output.OUTPUT_FORCE.ALMA.template renamed to samples/output.OUTPUT_FORCE.ALMA.template

@@ -0,0 +1,32 @@
+#!/bin/csh
+
+# This script is designed to change all of the files associated with
+# running the model for the current subbasin.  The parameters calibrated
+# are b_infilt, Ws, Ds, depth2, and depth3.
+
+set BASINNAME = wis
+set CALIBDIR = ARCINFO/
+set CALIBPREFIX = WIS_OPTI/
+
+if ( $5 != "" ) then
+
+	cd $CALIBDIR
+
+	set_basin_param.linux $BASINNAME"_run.asc" b_infilt.upm $CALIBPREFIX/b_infilt.upm $1
+
+	set_basin_param.linux $BASINNAME"_run.asc" Ds.upm $CALIBPREFIX/Ds.upm $2
+
+	set_basin_param.linux $BASINNAME"_run.asc" Ws.upm $CALIBPREFIX/Ws.upm $3
+
+	set_basin_param.linux $BASINNAME"_run.asc" depth2.upm $CALIBPREFIX/depth2.upm $4
+
+	set_basin_param.linux $BASINNAME"_run.asc" depth3.upm $CALIBPREFIX/depth3.upm $5
+
+	cd ..
+
+else
+
+	echo "Usage: $0 <b_infilt> <Ds> <Ws> <depth2> <depth3>"
+
+endif
+

output.OUTPUT_FORCE.template renamed to samples/output.OUTPUT_FORCE.template

@@ -0,0 +1,113 @@
+#include <stdio.h>
+#include <stdlib.h>
+#include <math.h>
+
+/** Define the units conversion factor between the simulated 
+    discharge (cfs), and the observed dischare (???) **/
+#define sim_to_obs 1
+
+/** Define the conversion factor from observed simulation rate to total
+  volume **/
+
+
+main(int argc, char *argv[]) {
+/************************************************************************
+  compute_R2.c              Keith Cherkauer            January 10, 1999
+
+  This program was written to comute the R^2 value for the current 
+  simulated daily discharge.
+
+************************************************************************/
+
+  FILE *fvic, *fobs, *fout;
+  int startrec, endrec;
+  int Nrecs, i;
+  float basinfactor;
+  float meanobs;
+  float meanvic;
+  float Rsqr;
+  float tmpdata;
+  float *vicdata;
+  float *obsdata;
+  float vicflow, obsflow;
+  float numsum, densum;
+
+  if(argc!=7) {
+    fprintf(stderr,"Usage: %s <VIC discharge> <observed discharge> <start rec> <end rec> <basin size factor> <R2 outfile>\n",argv[0]);
+    exit(0);
+  }
+
+  if((fvic=fopen(argv[1],"r"))==NULL) {
+    fprintf(stderr,"ERROR: Unable to open VIC simulation file %s.\n",argv[1]);
+    exit(1);
+  }
+
+  if((fobs=fopen(argv[2],"r"))==NULL) {
+    fprintf(stderr,"ERROR: Unable to open observed flow file %s.\n",argv[2]);
+    exit(2);
+  }
+
+  if((fout=fopen(argv[6],"w"))==NULL) {
+    fprintf(stderr,"ERROR: Unable to open output file %s.\n","R2.out");
+    exit(1);
+  }
+
+  startrec = atoi(argv[3]);
+  endrec = atoi(argv[4]);
+  basinfactor = atof(argv[5]);
+
+  /** Process Data **/
+  Nrecs   = endrec - startrec + 1;
+  vicdata = (float*)calloc(Nrecs,sizeof(float));
+  obsdata = (float*)calloc(Nrecs,sizeof(float));
+
+  for(i=0;i<startrec;i++) {
+    fscanf(fvic,"%*s %f\n",&tmpdata);
+    fscanf(fobs,"%*s %f\n",&tmpdata);
+  }
+
+  meanvic = meanobs = 0;
+  vicflow = obsflow = 0;
+  for(i=0;i<Nrecs;i++) {
+    fscanf(fvic,"%*s %f\n",&vicdata[i]);
+    fscanf(fobs,"%*s %f\n",&obsdata[i]);
+
+    /* Convert Simulated Flow to Match Observed */
+    vicdata[i] *= sim_to_obs;
+
+    vicdata[i] *= basinfactor;
+    meanvic += vicdata[i];
+    meanobs += obsdata[i];
+    /* Convert discharge from cubic feet per sec to cubic meters 
+       (per time step) to look at total flow volumes */
+    vicflow += vicdata[i] * 0.02831685 * 24. * 60. * 60.;
+    obsflow += obsdata[i] * 0.02831685 * 24. * 60. * 60.;
+  }
+  meanvic /= (float)Nrecs;
+  meanobs /= (float)Nrecs;
+  vicflow /= 1.e9;
+  obsflow /= 1.e9;
+
+  numsum = densum = 0;
+  for(i=0;i<Nrecs;i++) {
+    numsum += pow( obsdata[i] - vicdata[i], 2. );
+    densum += pow( obsdata[i] - meanobs, 2. );
+  }
+  Rsqr = -1. * ( 1. - ( numsum / densum) );
+
+  fprintf(fout,"Total observed flow  = %f km^3.\n",obsflow);
+  fprintf(fout,"Total simulated flow = %f km^3.\n",vicflow);
+  fprintf(fout,"Mean observed flow   = %f cfs.\n",meanobs);
+  fprintf(fout,"Mean simulated flow  = %f cfs.\n",meanvic);
+  fprintf(fout,"R^2 model error      = %f.\n",Rsqr);
+
+  fprintf(stdout,"%f\n", Rsqr);
+
+  free((char*)vicdata);
+  free((char*)obsdata);
+
+  fclose(fvic);
+  fclose(fobs);
+  fclose(fout);
+
+}