vicnl.h

超强的大尺度水文模拟工具

double linear_interp(double,double,double,double,double);

cell_data_struct **make_cell_data(int, int);
dist_prcp_struct make_dist_prcp(int, int *);
dmy_struct *make_dmy(global_param_struct *);
energy_bal_struct **make_energy_bal(int, int *);
filenames_struct make_in_and_outfiles(infiles_struct *, filenames_struct *, 
				      soil_con_struct *, outfiles_struct *);
out_data_struct *make_out_data(int);
snow_data_struct **make_snow_data(int);
veg_var_struct **make_veg_var(int);
void   MassRelease(double *,double *,double *,double *);
double maximum_unfrozen_water(double, double, double, double);
#if QUICK_FS
double maximum_unfrozen_water_quick(double, double, double **);
#endif
double modify_Ksat(double);
void mtclim42_wrapper(int, int, double, double, double, double, 
		      global_param_struct *, dmy_struct *, double *, 
		      double *, double *, double *, double *, double *);

double new_snow_density(double);
void   nrerror(char *);

void   open_debug();
FILE  *open_file(char string[], char type[]);
#if SAVE_STATE
FILE  *open_state_file(global_param_struct *, int, int);
#endif

double penman(double, double, double, double, double, double, double, 
	      double, double, float, float);
void   prepare_full_energy(int, int, int, dist_prcp_struct *,
			   soil_con_struct *, double *, double *);
double priestley(double, double);
void   put_data(dist_prcp_struct *, atmos_data_struct *, veg_con_struct *, 
		outfiles_struct *, double *, double *, double, double *, 
		char *, dmy_struct *, int, int, int, int); 

double quick_penman(double, double, double, double, double, double, 
		    double, double);
double read_arcinfo_value(char *, double, double);
int    read_arcinfo_info(char *, double **, double **, int **);
void   read_atmos_data(FILE *, global_param_struct, int, int, double **);
double **read_forcing_data(FILE **, global_param_struct);
void   read_initial_model_state(FILE *, dist_prcp_struct *, 
				global_param_struct *, int, int, int, 
				soil_con_struct *, int, char *, int *);
void   read_PILPS2c(atmos_data_struct *, FILE *, int *, int, int, int);
void   read_rosemount(atmos_data_struct *, FILE *, int *, int, int, int, int);
void   read_sawd(atmos_data_struct *, FILE *, int *, int, int, int);
void   read_sawd_binary(atmos_data_struct *, FILE *, int *, int, int, int);
void   read_snowband(FILE *, int, double, double **, double **, double **, 
                     char **);
void   read_snowmodel(atmos_data_struct *, FILE *, int, int, int, int);
soil_con_struct read_soilparam(FILE *, int);
soil_con_struct read_soilparam_arc(FILE *, char *, int *, int *, int);
veg_lib_struct *read_veglib(FILE *, int *);
veg_con_struct *read_vegparam(FILE *, int, int);
void   redistribute_during_storm(cell_data_struct ***, veg_var_struct ***,
				 int, int, int, double, double, double, 
				 double *);
void   redistribute_moisture(layer_data_struct *, double *, double *,
			     double *, double *, double *, int);
unsigned char redistribute_moisture_for_storm(double *, double *, double, 
					      double, double);
double root_brent(double, double, char *, double (*Function)(double, va_list), ...);
void   runoff(layer_data_struct *, layer_data_struct *, energy_bal_struct *, 
	      soil_con_struct *, double *, double *, double *, double *, 
	      double *, double, int, int, int, int, int);

void set_max_min_hour(double *, int, int *, int *);
void set_node_parameters(double *, double *, double *, double *, double *,
			 double *, double *, double *, double *, double *,
			 double *, double *, float **,
#if QUICK_FS
			 double ***,
#endif
			 int, int, char);
void   setup_frozen_soil(soil_con_struct *, layer_data_struct  *,
			 layer_data_struct *, layer_data_struct *,
			 energy_bal_struct, int, int, int, double,
			 double *, double *, double *);
double shrad(double,double,double,double,double,int,double);
double snow_albedo(double, double, double, double, int, char);
double snow_density(int, double, double, double, double, double, double, 
		    double);
void   snow_intercept(double, double, double, double, double, double, double,
		      double, double, double, double, double, double, double,
		      double *, double *, double *, double *, double *, 
		      double *, double *, double *, int, int, int);
void   snow_melt(soil_con_struct *, int, int, double, double, double,
		 snow_data_struct *, double, double, double, double, double, 
		 double, double, double, double, double, double, double, 
		 double, double, double, double *, double *, double *,
		 double *, double *, double *, double *, double *);
double SnowPackEnergyBalance(double, va_list);
double soil_conductivity(double, double, double, double, double);
void   soil_thermal_calc(soil_con_struct *, layer_data_struct *,
			 energy_bal_struct, double *, double *, double *,
			 int, int);
double soil_thermal_eqn(double, va_list);
double solve_snow(snow_data_struct *, layer_data_struct *, 
		  layer_data_struct *, veg_var_struct *, veg_var_struct *, 
		  int, int, energy_bal_struct *, soil_con_struct *, 
		  char, int, int, int, int, int, int, int, int, double, 
		  double, double, double, double, double, double, double, 
		  double, double, double, double, double, double, double, 
		  double, double, double, double, double, double *, double *,
		  double *, double *, double *, double *, double *, double *,
		  double *, double *, double *, double *, double *, double *, 
		  float *);
double solve_snow_ground_flux(double Tsurf, ...);
double solve_surf_energy_bal(double Tsurf, ...);
#if QUICK_FS
void   solve_T_profile(double *, double *, double *, double *,double *,
		       double *, double, double *, double *, double *,
		       double *, double *, double *, double *, double ***, 
		       int, char *, char, int);
#else
void   solve_T_profile(double *, double *, double *, double *,double *,
		       double *, double, double *, double *, double *,
		       double *, double *, double *, double *, int, char *, 
		       char, int);
#endif
double StabilityCorrection(double, double, double, double, double, double);
void   store_moisture_for_debug(int,int,double *,cell_data_struct ***,
				veg_var_struct ***,snow_data_struct **,
				soil_con_struct *);
void   surface_fluxes(char, int, int, int, int, int, int, int, int, int, 
		      double, double, double, double, double, double, 
		      double, double, double, double *, double *, double *, 
		      double *, double *, double *, double *, double *, 
		      double *, double *, double *, double *, double *, 
		      double *, float *, atmos_data_struct *, 
		      soil_con_struct *, dmy_struct *, global_param_struct *,
		      energy_bal_struct *, snow_data_struct *, 
		      layer_data_struct *, layer_data_struct *, 
		      veg_var_struct *, veg_var_struct *);
double svp(double);
double svp_slope(double);

void transpiration(layer_data_struct *, int, int, double, double, double, 
		   double, double, double, double, double, double, double, 
		   double *, double *, double *, double *, double *, 
		   double *, float *);
void tridag(double *,double *,double *,double *,double *,int);

void usage(char *);

void   vicerror(char *);
double volumetric_heat_capacity(double,double,double);

void write_atmosdata(atmos_data_struct *, int);
void write_data(out_data_struct *, outfiles_struct *, dmy_struct *, int);
void write_debug(atmos_data_struct *, soil_con_struct *, cell_data_struct *,
                 energy_bal_struct *, snow_data_struct *, veg_var_struct *,
                 dmy_struct *, global_param_struct *,
                 double, double, int, int, int, int, int, char);
void write_dist_prcp(dist_prcp_struct *);
#if OUTPUT_FORCE
void write_forcing_file(atmos_data_struct *, int nrecs, outfiles_struct *);
#endif
void write_layer(layer_data_struct *, int, int, double *);
#if SAVE_STATE
void write_model_state(dist_prcp_struct *, global_param_struct *, int, 
		       int, outfiles_struct *, soil_con_struct *, char, int);
#endif
void write_soilparam(soil_con_struct *);
void write_vegparam(veg_con_struct *);
void write_vegvar(veg_var_struct *, int);