static void*
pty_slave(void *arg)
{
char *cmd;
char *av[20];
int ac = 0;
cmdStub *ep = NULL;
int matched = 0;
while (!cmd_quit) {
cmd = readline("HV > " , rls);
fprintf(fptys, "\r\n");
fflush(fptys);
usleep(1);
if (!cmd)
continue;
ttrim(cmd);
ac = mkargv(cmd, (char **)av, 20);
if (ac == 0)
continue;
for (ep = cmd_entry, matched = 0; ep->name != NULL; ep++) {
if(!strncmp(av[0], ep->name, strlen(av[0]))) {
matched = 1;
break;
}
}
if (matched) {
ep->f(ac, av);
} else
ERR(fptys, "Invalid or incomplete command\r\n");
}
return NULL;
}
veg_con_struct *read_vegparam(FILE *vegparam,
int gridcel,
int Nveg_type)
/**********************************************************************
read_vegparam.c Keith Cherkauer and Dag Lohmann 1997
This routine reads in vegetation parameters for the current grid cell.
It also relates each vegetation class in the cell to the appropriate
parameters in the vegetation library.
Modifications:
09-24-98 Modified to read root zone distribution information so
that soil layer root fractions can be computed for new
soil layer depths - see calc_root_fractions.c KAC
07-15-99 Modified to read LAI values from a new line in the vegetation
parameter file. Added specifically to work with the new
global LAI files.
11-18-02 Added code to read in blowing snow parameters. LCB
03-27-03 Modified code to update Wdmax based on LAI values read in
for the current grid cell. If LAI is not obtained from this
function, then the values cacluated in read_veglib.c are
left unchanged. DP & KAC
2006-Nov-07 Allocates MaxVeg+1 veg tiles. TJB
2007-May-11 Changed some 'fscanf' statements to 'fgets' and 'sscanf'
to count rootzone and BLOWING fields. Also tests for
fetch < 1. GCT
2007-Oct-31 Added missing brackets in if(options.GLOBAL_LAI) block. TJB
2008-Oct-23 Added blocks to free vegarr[]. LCB via TJB
2009-Jan-16 Added logic for COMPUTE_TREELINE option. TJB
2009-Jun-09 Modified to use extension of veg_lib structure to contain
bare soil information. TJB
2009-Jun-17 Modified to understand both tabs and spaces as delimiters.TJB
2009-Jun-17 Fixed incorrect placement of free vegarr[] for case of
GLOBAL_LAI==FALSE. TJB
2009-Jul-26 Allocate extra veg tile for COMPUTE_TREELINE case. TJB
2009-Jul-31 Removed extra veg tile for lake/wetland case. TJB
2009-Sep-14 Made error messages clearer. TJB
2009-Oct-01 Added error message for case of LAI==0 and overstory==1. TJB
2010-Apr-28 Replaced GLOBAL_LAI with VEGPARAM_LAI and LAI_SRC. TJB
2012-Jan-16 Removed LINK_DEBUG code BN
2013-Jul-25 Added photosynthesis parameters. TJB
2013-Dec-28 Removed NO_REWIND option. TJB
2014-Apr-25 Added optional albedo values; added VEGPARAM_ALB and
ALB_SRC. TJB
2014-Apr-25 Added optional vegcover values; added VEGPARAM_VEGCOVER
and VEGCOVER_SRC. TJB
**********************************************************************/
{
void ttrim( char *string );
extern veg_lib_struct *veg_lib;
extern option_struct options;
veg_con_struct *temp;
int vegcel, i, j, k, vegetat_type_num, skip, veg_class;
int MaxVeg;
int Nfields, NfieldsMax;
int NoOverstory;
float depth_sum;
float sum;
char str[500];
char ErrStr[MAXSTRING];
char line[MAXSTRING];
char tmpline[MAXSTRING];
const char delimiters[] = " \t";
char *token;
char *vegarr[500];
size_t length;
int cidx;
double tmp;
skip = 1;
if(options.VEGPARAM_LAI) skip++;
if(options.VEGPARAM_VEGCOVER) skip++;
if(options.VEGPARAM_ALB) skip++;
NoOverstory = 0;
while ( ( fscanf(vegparam, "%d %d", &vegcel, &vegetat_type_num) == 2 ) && vegcel != gridcel ){
if (vegetat_type_num < 0) {
sprintf(ErrStr,"ERROR number of vegetation tiles (%i) given for cell %i is < 0.\n",vegetat_type_num,vegcel);
nrerror(ErrStr);
}
for (i = 0; i <= vegetat_type_num * skip; i++){
if ( fgets(str, 500, vegparam) == NULL ){
sprintf(ErrStr,"ERROR unexpected EOF for cell %i while reading root zones and LAI\n",vegcel);
nrerror(ErrStr);
}
}
}
fgets(str, 500, vegparam); // read newline at end of veg class line to advance to next line
if (vegcel != gridcel) {
fprintf(stderr, "Error in vegetation file. Grid cell %d not found\n",
gridcel);
exit(99);
}
if(vegetat_type_num >= MAX_VEG) {
sprintf(ErrStr,"Vegetation parameter file wants more vegetation tiles in grid cell %i (%i) than are allowed by MAX_VEG (%i) [NOTE: bare soil class is assumed]. Edit vicNl_def.h and recompile.",gridcel,vegetat_type_num+1,MAX_VEG);
nrerror(ErrStr);
}
// Make sure to allocate extra memory for bare soil tile
// and optionally an above-treeline veg tile
MaxVeg = vegetat_type_num+1;
if ( options.AboveTreelineVeg >= 0 )
MaxVeg++;
/** Allocate memory for vegetation grid cell parameters **/
temp = (veg_con_struct*) calloc( MaxVeg, sizeof(veg_con_struct));
temp[0].Cv_sum = 0.0;
for (i = 0; i < vegetat_type_num; i++) {
temp[i].zone_depth = calloc(options.ROOT_ZONES,sizeof(float));
temp[i].zone_fract = calloc(options.ROOT_ZONES,sizeof(float));
temp[i].vegetat_type_num = vegetat_type_num;
/* Upper boundaries of canopy layers, expressed in terms of fraction of total LAI */
if (options.CARBON) {
temp[i].CanopLayerBnd = calloc(options.Ncanopy,sizeof(double));
for (cidx=0; cidx<options.Ncanopy; cidx++) {
/* apportion LAI equally among layers */
temp[i].CanopLayerBnd[cidx] = (double)((cidx+1))/(double)(options.Ncanopy);
}
}
// Read the root zones line
if ( fgets( line, MAXSTRING, vegparam ) == NULL ){
sprintf(ErrStr,"ERROR unexpected EOF for cell %i while reading vegetat_type_num %d\n",vegcel,vegetat_type_num);
nrerror(ErrStr);
}
strcpy(tmpline, line);
ttrim( tmpline );
token = strtok (tmpline, delimiters); /* token => veg_class, move 'line' pointer to next field */
Nfields = 0;
vegarr[Nfields] = calloc( 500, sizeof(char));
strcpy(vegarr[Nfields],token);
Nfields++;
token = strtok (NULL, delimiters);
while (token != NULL && (length=strlen(token))==0) token = strtok (NULL, delimiters);
while ( token != NULL ) {
vegarr[Nfields] = calloc( 500, sizeof(char));
strcpy(vegarr[Nfields],token);
Nfields++;
token = strtok (NULL, delimiters);
while (token != NULL && (length=strlen(token))==0) token = strtok (NULL, delimiters);
}
NfieldsMax = 2 + 2 * options.ROOT_ZONES; /* Number of expected fields this line */
if( options.BLOWING ){
NfieldsMax += 3;
}
if ( Nfields != NfieldsMax ) {
sprintf(ErrStr,"ERROR - cell %d - expecting %d fields but found %d in veg line %s\n",gridcel,NfieldsMax, Nfields, line);
nrerror(ErrStr);
}
temp[i].LAKE = 0;
temp[i].veg_class = atoi( vegarr[0] );
temp[i].Cv = atof( vegarr[1] );
depth_sum = 0;
sum = 0.;
for(j=0;j<options.ROOT_ZONES;j++) {
temp[i].zone_depth[j] = atof( vegarr[2 + j*2] );
temp[i].zone_fract[j] = atof( vegarr[3 + j*2] );
depth_sum += temp[i].zone_depth[j];
sum += temp[i].zone_fract[j];
}
if(depth_sum <= 0) {
sprintf(str,"Root zone depths must sum to a value greater than 0.");
nrerror(str);
}
if(sum != 1.) {
fprintf(stderr,"WARNING: Root zone fractions sum to more than 1 ( = %f), normalizing fractions. If the sum is large, check that your vegetation parameter file is in the form - <zone 1 depth> <zone 1 fract> <zone 2 depth> <zone 2 fract> ...\n", sum);
for(j=0;j<options.ROOT_ZONES;j++) {
temp[i].zone_fract[j] /= sum;
}
}
if(options.BLOWING) {
j = 2 * options.ROOT_ZONES;
temp[i].sigma_slope = atof( vegarr[2 + j] );
temp[i].lag_one = atof( vegarr[3 + j] );
temp[i].fetch = atof( vegarr[4 + j]) ;
if( temp[i].sigma_slope <= 0. || temp[i].lag_one <= 0.) {
sprintf(str,"Deviation of terrain slope must be greater than 0.");
nrerror(str);
}
if( temp[i].fetch < 1.0 ) {
sprintf(str,"ERROR - BLOWING parameter fetch should be >> 1 but cell %i has fetch = %.2f\n", gridcel, temp[i].fetch );
nrerror(str);
}
}
veg_class = MISSING;
for(j=0;j<Nveg_type;j++)
if(temp[i].veg_class == veg_lib[j].veg_class)
veg_class = j;
if(veg_class == MISSING) {
sprintf(ErrStr,"The vegetation class id %i in vegetation tile %i from cell %i is not defined in the vegetation library file.", temp[i].veg_class, i, gridcel);
nrerror(ErrStr);
}
else
temp[i].veg_class = veg_class;
temp[0].Cv_sum += temp[i].Cv;
for(k=0; k<Nfields; k++)
free(vegarr[k]);
if ( options.VEGPARAM_LAI ) {
// Read the LAI line
if ( fgets( line, MAXSTRING, vegparam ) == NULL ){
sprintf(ErrStr,"ERROR unexpected EOF for cell %i while reading LAI for vegetat_type_num %d\n",vegcel,vegetat_type_num);
nrerror(ErrStr);
}
Nfields = 0;
vegarr[Nfields] = calloc( 500, sizeof(char));
strcpy(tmpline, line);
ttrim( tmpline );
token = strtok (tmpline, delimiters);
strcpy(vegarr[Nfields],token);
Nfields++;
while( ( token = strtok (NULL, delimiters)) != NULL ){
vegarr[Nfields] = calloc( 500, sizeof(char));
strcpy(vegarr[Nfields],token);
Nfields++;
}
NfieldsMax = 12; /* For LAI */
if ( Nfields != NfieldsMax ) {
sprintf(ErrStr,"ERROR - cell %d - expecting %d LAI values but found %d in line %s\n",gridcel, NfieldsMax, Nfields, line);
nrerror(ErrStr);
}
for ( j = 0; j < MONTHSPERYEAR; j++ ) {
temp[i].albedo[j] = veg_lib[temp[i].veg_class].albedo[j];
temp[i].LAI[j] = veg_lib[temp[i].veg_class].LAI[j];
temp[i].vegcover[j] = veg_lib[temp[i].veg_class].vegcover[j];
temp[i].Wdmax[j] = veg_lib[temp[i].veg_class].Wdmax[j];
}
if (options.LAI_SRC == FROM_VEGPARAM) {
for ( j = 0; j < 12; j++ ) {
tmp = atof( vegarr[j] );
if (tmp != NODATA_VH)
temp[i].LAI[j] = tmp;
if (veg_lib[temp[i].veg_class].overstory && temp[i].LAI[j] == 0) {
sprintf(ErrStr,"ERROR: cell %d, veg tile %d: the specified veg class (%d) is listed as an overstory class in the veg LIBRARY, but the LAI given in the veg PARAM FILE for this tile for month %d is 0.\n",gridcel, i+1, temp[i].veg_class+1, j+1);
nrerror(ErrStr);
}
temp[i].Wdmax[j] = LAI_WATER_FACTOR * temp[i].LAI[j];
}
}
for(k=0; k<Nfields; k++)
free(vegarr[k]);
}
if ( options.VEGPARAM_VEGCOVER ) {
// Read the vegcover line
if ( fgets( line, MAXSTRING, vegparam ) == NULL ){
sprintf(ErrStr,"ERROR unexpected EOF for cell %i while reading vegcover for vegetat_type_num %d\n",vegcel,vegetat_type_num);
nrerror(ErrStr);
}
Nfields = 0;
vegarr[Nfields] = calloc( 500, sizeof(char));
strcpy(tmpline, line);
ttrim( tmpline );
token = strtok (tmpline, delimiters);
strcpy(vegarr[Nfields],token);
Nfields++;
while( ( token = strtok (NULL, delimiters)) != NULL ){
vegarr[Nfields] = calloc( 500, sizeof(char));
strcpy(vegarr[Nfields],token);
Nfields++;
}
NfieldsMax = 12; /* For vegcover */
if ( Nfields != NfieldsMax ) {
sprintf(ErrStr,"ERROR - cell %d - expecting %d vegcover values but found %d in line %s\n",gridcel, NfieldsMax, Nfields, line);
nrerror(ErrStr);
}
if (options.VEGCOVER_SRC == FROM_VEGPARAM) {
for ( j = 0; j < 12; j++ ) {
tmp = atof( vegarr[j] );
if (tmp != NODATA_VH)
temp[i].vegcover[j] = tmp;
}
}
for(k=0; k<Nfields; k++)
free(vegarr[k]);
}
if ( options.VEGPARAM_ALB ) {
// Read the albedo line
if ( fgets( line, MAXSTRING, vegparam ) == NULL ){
sprintf(ErrStr,"ERROR unexpected EOF for cell %i while reading albedo for vegetat_type_num %d\n",vegcel,vegetat_type_num);
nrerror(ErrStr);
}
Nfields = 0;
vegarr[Nfields] = calloc( 500, sizeof(char));
strcpy(tmpline, line);
ttrim( tmpline );
token = strtok (tmpline, delimiters);
strcpy(vegarr[Nfields],token);
Nfields++;
while( ( token = strtok (NULL, delimiters)) != NULL ){
vegarr[Nfields] = calloc( 500, sizeof(char));
strcpy(vegarr[Nfields],token);
Nfields++;
}
NfieldsMax = 12; /* For albedo */
if ( Nfields != NfieldsMax ) {
sprintf(ErrStr,"ERROR - cell %d - expecting %d albedo values but found %d in line %s\n",gridcel, NfieldsMax, Nfields, line);
nrerror(ErrStr);
}
if (options.ALB_SRC == FROM_VEGPARAM) {
for ( j = 0; j < 12; j++ ) {
tmp = atof( vegarr[j] );
if (tmp != NODATA_VH)
temp[i].albedo[j] = tmp;
}
}
for(k=0; k<Nfields; k++)
free(vegarr[k]);
}
// Determine if cell contains non-overstory vegetation
if ( options.COMPUTE_TREELINE && !veg_lib[temp[i].veg_class].overstory )
NoOverstory++;
}
// Determine if we have bare soil
if(temp[0].Cv_sum>1.0){
fprintf(stderr,"WARNING: Cv exceeds 1.0 at grid cell %d, fractions being adjusted to equal 1\n", gridcel);
for(j=0;j<vegetat_type_num;j++)
temp[j].Cv = temp[j].Cv / temp[0].Cv_sum;
temp[0].Cv_sum = 1.;
}
else if(temp[0].Cv_sum>0.99 && temp[0].Cv_sum<1.0){
fprintf(stderr,"WARNING: Cv > 0.99 and Cv < 1.0 at grid cell %d, model assuming that bare soil is not to be run - fractions being adjusted to equal 1\n", gridcel);
for(j=0;j<vegetat_type_num;j++)
temp[j].Cv = temp[j].Cv / temp[0].Cv_sum;
temp[0].Cv_sum = 1.;
}
// Handle veg above the treeline
if ( options.SNOW_BAND > 1 && options.COMPUTE_TREELINE
&& ( !NoOverstory && temp[0].Cv_sum == 1. ) ) {
// All vegetation in the current cell is defined with overstory.
// Add default non-overstory vegetation so that snow bands above treeline
// can be sucessfully simulated.
if ( options.AboveTreelineVeg < 0 ) {
// Above treeline snowband should be treated as bare soil
for ( j = 0; j < vegetat_type_num; j++ )
temp[j].Cv -= ( 0.001 / (float)vegetat_type_num );
temp[0].Cv_sum -= 0.001;
}
else {
// Above treeline snowband should use the defined vegetation
// add vegetation to typenum
// check that veg type exists in library and does not have overstory
if(vegetat_type_num > 0) {
for ( j = 0; j < vegetat_type_num; j++ ) {
temp[j].Cv -= ( 0.001 / (float)vegetat_type_num );
temp[j].vegetat_type_num++;
}
temp[vegetat_type_num].Cv = 0.001;
temp[vegetat_type_num].veg_class = options.AboveTreelineVeg;
temp[vegetat_type_num].Cv_sum = temp[vegetat_type_num-1].Cv_sum;
temp[vegetat_type_num].zone_depth = calloc( options.ROOT_ZONES,
sizeof(float));
temp[vegetat_type_num].zone_fract = calloc( options.ROOT_ZONES,
sizeof(float));
temp[vegetat_type_num].vegetat_type_num = vegetat_type_num+1;
// Since root zones are not defined they are copied from the last
// vegetation type.
for ( j = 0; j < options.ROOT_ZONES; j++ ) {
temp[vegetat_type_num].zone_depth[j]
= temp[vegetat_type_num-1].zone_depth[j];
temp[vegetat_type_num].zone_fract[j]
= temp[vegetat_type_num-1].zone_fract[j];
}
}
// Identify current vegetation class
veg_class = MISSING;
for ( j = 0; j < Nveg_type; j++ ) {
if(temp[vegetat_type_num].veg_class == veg_lib[j].veg_class) {
veg_class = j;
break;
}
}
if ( veg_class == MISSING ) {
sprintf(ErrStr,"The vegetation class id %i defined for above-treeline from cell %i is not defined in the vegetation library file.", temp[vegetat_type_num].veg_class, gridcel);
nrerror(ErrStr);
}
else {
temp[vegetat_type_num].veg_class = veg_class;
}
if ( veg_lib[veg_class].overstory ) {
sprintf(ErrStr,"Vegetation class %i is defined to have overstory, so it cannot be used as the default vegetation type for above canopy snow bands.", veg_lib[veg_class].veg_class );
nrerror(ErrStr);
}
}
}
// Bare soil tile
if (temp[0].Cv_sum < 1.) {
j = vegetat_type_num;
temp[j].veg_class = Nveg_type; // Create a veg_class ID for bare soil, which is not mentioned in the veg library
temp[j].Cv = 1.0 - temp[0].Cv_sum;
// Don't allocate any root-zone-related arrays
if(options.BLOWING) {
if (vegetat_type_num > 0) {
temp[j].sigma_slope = temp[0].sigma_slope;
temp[j].lag_one = temp[0].lag_one;
temp[j].fetch = temp[0].fetch;
}
else {
temp[j].sigma_slope = 0.005;
temp[j].lag_one = 0.95;
temp[j].fetch = 2000;
}
}
}
return temp;
}
