#include <stdio.h>#include <stdlib.h>#include <vicNl.h>static char vcid[] = "$Id: read_vegparam.c,v 4.2.2.4 2004/05/06 19:57:34 tbohn Exp $";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.  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  09-02-2003 Moved COMPUTE_TREELINE flag from user_def.h to the              options structure.  Now when not set to FALSE, the              value indicates the default above treeline vegetation             if no usable vegetation types are in the grid cell              (i.e. everything has a canopy).  A negative value               will cause the model to use bare soil.  Make sure that              positive index value refer to a non-canopied vegetation             type in the vegetation library.                   KAC  08-Dec-03 Applied Alan Hamlet's fix for COMPUTE_TREELINE option,	    which fixed a segmentation fault when COMPUTE_TREELINE=TRUE.	    This consisted of removing the call to realloc and	    instead allocating an extra veg class to begin with,	    as well as assigning this extra veg class a very small	    fraction of the grid cell's area to avoid changing the	    results for areas below the treeline.		TJB**********************************************************************/{  extern veg_lib_struct *veg_lib;  extern option_struct   options;#if LINK_DEBUG  extern debug_struct    debug;#endif  veg_con_struct *temp;  int             vegcel, i, j, vegetat_type_num, skip, veg_class;  int             NoOverstory;  float           depth_sum;  float           sum;  char            str[500];  char            ErrStr[MAXSTRING];  if(options.GLOBAL_LAI) skip=2;  else skip=1;  NoOverstory = 0;#if !NO_REWIND  rewind(vegparam);#endif        while ((fscanf(vegparam, "%d %d", &vegcel, &vegetat_type_num)) == 2 &&          vegcel != gridcel) {    for (i = 0; i <= vegetat_type_num * skip; i++)      fgets(str, 500, vegparam);  }  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 types in grid cell %i (%i) than are defined 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);  }  /** Allocate memory for vegetation grid cell parameters **/  if(vegetat_type_num>0)    temp = (veg_con_struct*) calloc(vegetat_type_num+1,                                     sizeof(veg_con_struct));  else    temp = (veg_con_struct*) calloc(1, 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;    fscanf(vegparam, "%d",  &temp[i].veg_class);    fscanf(vegparam, "%lf", &temp[i].Cv);        depth_sum = 0;    sum = 0.;    for(j=0;j<options.ROOT_ZONES;j++) {      fscanf(vegparam,"%f %f",&temp[i].zone_depth[j], &temp[i].zone_fract[j]);      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;      }    }    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,"Vegetation class %i from cell %i is not defined in the vegetation library file.", temp[i].veg_class, gridcel);      nrerror(ErrStr);    }    else      temp[i].veg_class = veg_class;    temp[0].Cv_sum += temp[i].Cv;    if ( options.GLOBAL_LAI ) {      for ( j = 0; j < 12; j++ ) {	fscanf(vegparam,"%lf",&veg_lib[temp[i].veg_class].LAI[j]);	veg_lib[temp[i].veg_class].Wdmax[j] = 	  LAI_WATER_FACTOR * veg_lib[temp[i].veg_class].LAI[j];      }    }    if ( options.COMPUTE_TREELINE && !veg_lib[temp[i].veg_class].overstory )       // Determine if cell contains non-overstory vegetation      NoOverstory++;      }  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.;  }  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.;  }  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;        for ( j = 0; j < options.ROOT_ZONES; j++ ) {	  // Since root zones are not defined they are copied from the last	  // vegetation type.	  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];        }      }      veg_class = MISSING;      for ( j = 0; j < Nveg_type; j++ ) {	// Identify current vegetation class	if(temp[vegetat_type_num].veg_class == veg_lib[j].veg_class) {	  veg_class = j;	  break;	}      }      if ( veg_class == MISSING ) {	sprintf(ErrStr,"Vegetation class %i from cell %i is not defined in the vegetation library file.", temp[i].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);      }    }  }  return temp;}