struct zone_object *construct_zone(
struct command_line_object *command_line,
FILE *world_file,
int *num_world_base_stations,
struct base_station_object **world_base_stations,
struct default_object *defaults,
struct base_station_ncheader_object *base_station_ncheader,
struct world_object *world)
{
/*--------------------------------------------------------------*/
/* Local function definition. */
/*--------------------------------------------------------------*/
struct base_station_object *assign_base_station(
int ,
int ,
struct base_station_object **);
struct base_station_object *assign_base_station_xy(
float ,
float ,
int ,
int *,
struct base_station_object **);
struct patch_object *construct_patch(
struct command_line_object *,
FILE *,
int num_world_base_stations,
struct base_station_object **world_base_stations,
struct default_object *defaults);
struct base_station_object *construct_netcdf_grid(
struct base_station_ncheader_object *,
int *,
float,
float,
float,
struct date,
struct date);
int get_netcdf_xy(char *, char *, char *, float, float, float, float *, float *);
void *alloc(size_t, char *, char *);
double atm_pres( double );
/*--------------------------------------------------------------*/
/* Local variable definition. */
/*--------------------------------------------------------------*/
int base_stationID;
int i, k, j;
int default_object_ID;
int notfound;
int *dum;
float base_x, base_y;
char record[MAXSTR];
struct zone_object *zone;
notfound = 0;
base_x = 0.0;
base_y = 0.0;
j = 0;
k = 0;
/*--------------------------------------------------------------*/
/* Allocate a zone object. */
/*--------------------------------------------------------------*/
zone = (struct zone_object *) alloc( 1 *
sizeof( struct zone_object ),"zone","construct_zone" );
/*--------------------------------------------------------------*/
/* Read in the next zone record for this hillslope. */
/*--------------------------------------------------------------*/
fscanf(world_file,"%d",&(zone[0].ID));
read_record(world_file, record);
fscanf(world_file,"%lf",&(zone[0].x));
read_record(world_file, record);
fscanf(world_file,"%lf",&(zone[0].y));
read_record(world_file, record);
fscanf(world_file,"%lf",&(zone[0].z));
read_record(world_file, record);
fscanf(world_file,"%d",&(default_object_ID));
read_record(world_file, record);
fscanf(world_file,"%lf",&(zone[0].area));
read_record(world_file, record);
fscanf(world_file,"%lf",&(zone[0].slope));
read_record(world_file, record);
fscanf(world_file,"%lf",&(zone[0].aspect));
read_record(world_file, record);
fscanf(world_file,"%lf",&(zone[0].precip_lapse_rate));
read_record(world_file, record);
fscanf(world_file,"%lf",&(zone[0].e_horizon));
read_record(world_file, record);
fscanf(world_file,"%lf",&(zone[0].w_horizon));
read_record(world_file, record);
fscanf(world_file,"%d",&(zone[0].num_base_stations));
read_record(world_file, record);
/*--------------------------------------------------------------*/
/* convert from degrees to radians for slope and aspects */
/*--------------------------------------------------------------*/
zone[0].aspect = zone[0].aspect * DtoR;
zone[0].slope = zone[0].slope * DtoR;
/*--------------------------------------------------------------*/
/* Define cos and sine of slopes and aspects. */
/*--------------------------------------------------------------*/
zone[0].cos_aspect = cos(zone[0].aspect);
zone[0].cos_slope = cos(zone[0].slope);
zone[0].sin_aspect = sin(zone[0].aspect);
zone[0].sin_slope = sin(zone[0].slope);
/*--------------------------------------------------------------*/
/* Initialize accumulator variables */
/*--------------------------------------------------------------*/
zone[0].acc_month.K_direct = 0.0;
zone[0].acc_month.K_diffuse = 0.0;
zone[0].acc_month.tmax = 0.0;
zone[0].acc_month.tmin = 0.0;
zone[0].acc_month.precip = 0.0;
zone[0].acc_month.length = 0;
/*--------------------------------------------------------------*/
/* Define hourly array zone */
/*--------------------------------------------------------------*/
zone[0].hourly = (struct zone_hourly_object *) alloc(
sizeof( struct zone_hourly_object ),
"hourly","zone_hourly");
/*--------------------------------------------------------------*/
/* Assign defaults for this zone */
/*--------------------------------------------------------------*/
zone[0].defaults = (struct zone_default **)
alloc( sizeof(struct zone_default *),"defaults",
"construct_zone" );
i = 0;
while (defaults[0].zone[i].ID != default_object_ID) {
i++;
/*--------------------------------------------------------------*/
/* Report an error if no match was found. Otherwise assign */
/* the default to point to this zone. */
/*--------------------------------------------------------------*/
if ( i>= defaults[0].num_zone_default_files ){
fprintf(stderr,
"\nFATAL ERROR: in construct_zone, zone default ID %d not found.\n",
default_object_ID);
exit(EXIT_FAILURE);
}
} /* end-while */
zone[0].defaults[0] = &defaults[0].zone[i];
if ( command_line[0].verbose_flag == -3 ){
printf("\n>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>\n");
printf("Constructing zone base stations %d \n", zone[0].ID);
}
/*--------------------------------------------------------------*/
/* Allocate a list of base stations for this zone. */
/*--------------------------------------------------------------*/
zone[0].base_stations = (struct base_station_object **)
alloc(zone[0].num_base_stations *
sizeof(struct base_station_object *),
"base_stations","construct_zone" );
/*--------------------------------------------------------------*/
/* NON NETCDF BASE STATIONS */
/*--------------------------------------------------------------*/
/* Read each base_station ID and then point to that base station */
/*--------------------------------------------------------------*/
if (command_line[0].gridded_netcdf_flag == 0){
for (i = 0; i < zone[0].num_base_stations; i++ ){
fscanf(world_file, "%d", &(base_stationID));
read_record(world_file, record);
/*--------------------------------------------------------------*/
/* Point to the appropriate base station in the base */
/* station list for this world. */
/* */
/*--------------------------------------------------------------*/
zone[0].base_stations[i] = assign_base_station(
base_stationID,
*num_world_base_stations,
world_base_stations);
} /*end for*/
}
else {
fscanf(world_file, "%d", &(dum));
read_record(world_file, record);
}
/*--------------------------------------------------------------*/
/* NETCDF BASE STATIONS */
/*--------------------------------------------------------------*/
/* Construct the list of base stations. */
/*--------------------------------------------------------------*/
if (command_line[0].dclim_flag == 0) {
if(command_line[0].gridded_netcdf_flag == 1){
/*--------------------------------------------------------------*/
/* Construct the base_stations. */
/*--------------------------------------------------------------*/
if ( command_line[0].verbose_flag == -3 ){
printf("Constructing base stations from %s: zoney=%lf zonex=%lf num=%d sdist=%lf base_y=%lf base_x=%lf \n",
base_station_ncheader[0].netcdf_tmax_filename,
zone[0].y,
zone[0].x,
*num_world_base_stations,
base_station_ncheader[0].sdist,
base_y,
base_x);
printf("STARTING CLOSEST CELL: y=%lf x=%lf \n",base_y,base_x);
}
/* Identify centerpoint coords for closest netcdf cell to zone x, y */
k = get_netcdf_xy(base_station_ncheader[0].netcdf_tmax_filename,
base_station_ncheader[0].netcdf_y_varname,
base_station_ncheader[0].netcdf_x_varname,
zone[0].y,
zone[0].x,
base_station_ncheader[0].sdist,
&(base_y),
&(base_x));
if ( command_line[0].verbose_flag == -3 ){
printf("CLOSEST CELL: y=%lf x=%lf num=%d \n",base_y,base_x,*num_world_base_stations);
}
if(k == -1){
fprintf(stderr,"Can't locate station data in netcdf for var tmax\n");
exit(0);
}
/* Assign base station based on closest found coordinates */
if (*num_world_base_stations > 0) {
zone[0].base_stations[0] = assign_base_station_xy(
base_x,
base_y,
*num_world_base_stations,
&(notfound),
world_base_stations);
}
else notfound = 1;
/* If station is not already in list, add it */
if (notfound == 1) {
if ( command_line[0].verbose_flag == -3 ){
printf("Starting construct_netcdf_grid: file=%s num=%d \n",
world[0].base_station_files[0],
*num_world_base_stations);
}
j = *num_world_base_stations;
world_base_stations[j] = construct_netcdf_grid(
base_station_ncheader,
num_world_base_stations,
base_x,
base_y,
zone[0].z,
world[0].start_date,
world[0].duration);
world[0].num_base_stations = *num_world_base_stations;
if ( command_line[0].verbose_flag == -3 ){
printf("Assigning base station: file=%s num=%d lai=%lf notfound=%d lastid=%d numworld=%d ID=%d LAI=%lf \n",
world[0].base_station_files[0],
*num_world_base_stations,
(*(world_base_stations[j])).effective_lai,
notfound,
base_station_ncheader[0].lastID,
world[0].num_base_stations,
(*(world_base_stations[j])).ID,
(*(world_base_stations[j])).effective_lai);
}
/* Now link the zone to the newly added base station */
zone[0].base_stations[0] = world_base_stations[j];
if ( command_line[0].verbose_flag == -3 ){
printf("Zone base station: lai=%lf tmin=%lf tmax=%lf rain=%lf wind=%lf \n",
(*(zone[0].base_stations[0])).effective_lai,
(*(zone[0].base_stations[0])).daily_clim[0].tmin[0],
(*(zone[0].base_stations[0])).daily_clim[0].tmax[0],
(*(zone[0].base_stations[0])).daily_clim[0].rain[0]
//(*(zone[0].base_stations[0])).daily_clim[0].wind[0]
); // T.N
}
}
}
}
if ( command_line[0].verbose_flag == -3 ){
printf("Ending zone base station: num=%d worldnum=%d \n", *num_world_base_stations, world[0].num_base_stations);
printf("------------------------------\n");
}
/*--------------------------------------------------------------*/
/* Read in number of patches in this zone. */
/*--------------------------------------------------------------*/
fscanf(world_file,"%d",&(zone[0].num_patches));
read_record(world_file, record);
/*--------------------------------------------------------------*/
/* Allocate list of pointers to patch objects . */
/*--------------------------------------------------------------*/
zone[0].patches = ( struct patch_object ** )
alloc( zone[0].num_patches * sizeof( struct patch_object *),
"patches","construct_zone");
/*--------------------------------------------------------------*/
/* Initialize any variables that should be initialized at */
/* the start of a simulation run for the zone. */
/* for temperaturature we initialize to 999 so that they will be set on the first day based */
/* on air temperatures on that day - which we don't know at this point */
/*--------------------------------------------------------------*/
zone[0].metv.pa = atm_pres( zone[0].z );
zone[0].metv.tsoil_sum = 0.0;
zone[0].metv.tsoil = 0.0;
zone[0].metv.tmin_ravg = 3.0;
zone[0].metv.vpd_ravg = 900;
zone[0].metv.dayl_ravg = 38000;
/*--------------------------------------------------------------*/
/* Construct the intervals in this zone. */
/*--------------------------------------------------------------*/
for ( i=0 ; i<zone[0].num_patches ; i++ ){
zone[0].patches[i] = construct_patch(
command_line,
world_file,
*num_world_base_stations,
world_base_stations,
defaults);
zone[0].patches[i][0].zone = zone;
} /*end for*/
return(zone);
} /*end construct_zone.c*/
struct zone_object *construct_zone(
struct command_line_object *command_line,
FILE *world_file,
int num_world_base_stations,
struct base_station_object **world_base_stations,
struct default_object *defaults)
{
/*--------------------------------------------------------------*/
/* Local function definition. */
/*--------------------------------------------------------------*/
struct base_station_object *assign_base_station(
int ,
int ,
struct base_station_object **);
struct patch_object *construct_patch(
struct command_line_object *,
FILE *,
int num_world_base_stations,
struct base_station_object **world_base_stations,
struct default_object *defaults);
void *alloc(size_t, char *, char *);
double atm_pres( double );
/*--------------------------------------------------------------*/
/* Local variable definition. */
/*--------------------------------------------------------------*/
int base_stationID;
int i;
int default_object_ID;
char record[MAXSTR];
struct zone_object *zone;
/*--------------------------------------------------------------*/
/* Allocate a zone object. */
/*--------------------------------------------------------------*/
zone = (struct zone_object *) alloc( 1 *
sizeof( struct zone_object ),"zone","construct_zone" );
/*--------------------------------------------------------------*/
/* Read in the next zone record for this hillslope. */
/*--------------------------------------------------------------*/
fscanf(world_file,"%d",&(zone[0].ID));
read_record(world_file, record);
fscanf(world_file,"%lf",&(zone[0].x));
read_record(world_file, record);
fscanf(world_file,"%lf",&(zone[0].y));
read_record(world_file, record);
fscanf(world_file,"%lf",&(zone[0].z));
read_record(world_file, record);
fscanf(world_file,"%d",&(default_object_ID));
read_record(world_file, record);
fscanf(world_file,"%lf",&(zone[0].area));
read_record(world_file, record);
fscanf(world_file,"%lf",&(zone[0].slope));
read_record(world_file, record);
fscanf(world_file,"%lf",&(zone[0].aspect));
read_record(world_file, record);
fscanf(world_file,"%lf",&(zone[0].precip_lapse_rate));
read_record(world_file, record);
fscanf(world_file,"%lf",&(zone[0].e_horizon));
read_record(world_file, record);
fscanf(world_file,"%lf",&(zone[0].w_horizon));
read_record(world_file, record);
fscanf(world_file,"%d",&(zone[0].num_base_stations));
read_record(world_file, record);
/*--------------------------------------------------------------*/
/* convert from degrees to radians for slope and aspects */
/*--------------------------------------------------------------*/
zone[0].aspect = zone[0].aspect * DtoR;
zone[0].slope = zone[0].slope * DtoR;
/*--------------------------------------------------------------*/
/* Define cos and sine of slopes and aspects. */
/*--------------------------------------------------------------*/
zone[0].cos_aspect = cos(zone[0].aspect);
zone[0].cos_slope = cos(zone[0].slope);
zone[0].sin_aspect = sin(zone[0].aspect);
zone[0].sin_slope = sin(zone[0].slope);
/*--------------------------------------------------------------*/
/* Initialize accumulator variables */
/*--------------------------------------------------------------*/
zone[0].acc_month.K_direct = 0.0;
zone[0].acc_month.K_diffuse = 0.0;
zone[0].acc_month.tmax = 0.0;
zone[0].acc_month.tmin = 0.0;
zone[0].acc_month.precip = 0.0;
zone[0].acc_month.length = 0;
/*--------------------------------------------------------------*/
/* Define hourly array zone */
/*--------------------------------------------------------------*/
zone[0].hourly = (struct zone_hourly_object *) alloc(
sizeof( struct zone_hourly_object ),
"hourly","zone_hourly");
/*--------------------------------------------------------------*/
/* Assign defaults for this zone */
/*--------------------------------------------------------------*/
zone[0].defaults = (struct zone_default **)
alloc( sizeof(struct zone_default *),"defaults",
"construct_zone" );
i = 0;
while (defaults[0].zone[i].ID != default_object_ID) {
i++;
/*--------------------------------------------------------------*/
/* Report an error if no match was found. Otherwise assign */
/* the default to point to this zone. */
/*--------------------------------------------------------------*/
if ( i>= defaults[0].num_zone_default_files ){
fprintf(stderr,
"\nFATAL ERROR: in construct_zone, zone default ID %d not found.\n",
default_object_ID);
exit(EXIT_FAILURE);
}
} /* end-while */
zone[0].defaults[0] = &defaults[0].zone[i];
/*--------------------------------------------------------------*/
/* Allocate a list of base stations for this zone. */
/*--------------------------------------------------------------*/
zone[0].base_stations = (struct base_station_object **)
alloc(zone[0].num_base_stations *
sizeof(struct base_station_object *),
"base_stations","construct_zone" );
/*--------------------------------------------------------------*/
/* Read each base_station ID and then point to that base station */
/*--------------------------------------------------------------*/
for (i=0 ; i<zone[0].num_base_stations ; i++ ){
fscanf(world_file,"%d",&(base_stationID));
read_record(world_file, record);
/*--------------------------------------------------------------*/
/* Point to the appropriate base station in the base */
/* station list for this world. */
/* */
/*--------------------------------------------------------------*/
zone[0].base_stations[i] = assign_base_station(
base_stationID,
num_world_base_stations,
world_base_stations);
} /*end for*/
/*--------------------------------------------------------------*/
/* Read in number of patches in this zone. */
/*--------------------------------------------------------------*/
fscanf(world_file,"%d",&(zone[0].num_patches));
read_record(world_file, record);
/*--------------------------------------------------------------*/
/* Allocate list of pointers to patch objects . */
/*--------------------------------------------------------------*/
zone[0].patches = ( struct patch_object ** )
alloc( zone[0].num_patches * sizeof( struct patch_object *),
"patches","construct_zone");
/*--------------------------------------------------------------*/
/* Initialize any variables that should be initialized at */
/* the start of a simulation run for the zone. */
/* for temperaturature we initialize to 999 so that they will be set on the first day based */
/* on air temperatures on that day - which we don't know at this point */
/*--------------------------------------------------------------*/
zone[0].metv.pa = atm_pres( zone[0].z );
zone[0].metv.tsoil_sum = 0.0;
zone[0].metv.tsoil = 0.0;
zone[0].metv.tmin_ravg = 3.0;
zone[0].metv.vpd_ravg = 900;
zone[0].metv.dayl_ravg = 38000;
/*--------------------------------------------------------------*/
/* Construct the intervals in this zone. */
/*--------------------------------------------------------------*/
for ( i=0 ; i<zone[0].num_patches ; i++ ){
zone[0].patches[i] = construct_patch(
command_line,
world_file,
num_world_base_stations,
world_base_stations,
defaults);
zone[0].patches[i][0].zone = zone;
} /*end for*/
return(zone);
} /*end construct_zone.c*/
