vicnl.h

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

/* RCS Id String
 * $Id: vicNl.h,v 4.4.2.4 2004/09/22 00:40:34 vicadmin Exp $
 */
#include <math.h>
#include <vicNl_def.h>

/*** SubRoutine Prototypes ***/

void alloc_atmos(int, atmos_data_struct **);
double arno_evap(layer_data_struct *, layer_data_struct *, double, double, 
		 double, double, double, double, double, double, double, 
		 double, double, double, double, double, double, double);
unsigned char average_moisture_for_storm(double *, double *, double, double);

void CalcAerodynamic(char, int, int, double, double, double, double,
		     double *, double *, double *, double, double *,
		     double *);
void   calc_cloud_cover_fraction(atmos_data_struct *, dmy_struct *, int,
				 int, int, double *);
void   calc_energy_balance_error(int, double, double, double, double, double);
void   calc_longwave(double *, double, double, double);
void   calc_netlongwave(double *, double, double, double);
double calc_netshort(double, int, double, double *);
double calc_rainonly(double,double,double,double,double);
void   calc_root_fractions(veg_con_struct *, soil_con_struct *);
double calc_snow_ground_flux(int, int, int, int, double, double, double, 
			     double, double, double *, double *, double *, 
			     double *, energy_bal_struct *, 
			     snow_data_struct *, layer_data_struct *,
                             layer_data_struct *, soil_con_struct *, char *);
#if QUICK_FS
int    calc_soil_thermal_fluxes(int, double *, double *, double *, double *, 
				double *, double *, double *, double *, 
				double *, double *, double *, double *, 
				double *, double *, double ***, char);
#else
int    calc_soil_thermal_fluxes(int, double *, double *, double *, double *, 
				double *, double *, double *, double *, 
				double *, double *, double *, double *, 
				double *, double *, char);
#endif
double CalcSnowPackEnergyBalance(double Tsurf, ...);
double calc_surf_energy_bal(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 *, float *,
			    atmos_data_struct *, veg_var_struct *, 
			    veg_var_struct *, energy_bal_struct *, 
			    snow_data_struct *,
			    layer_data_struct *, layer_data_struct *, 
			    soil_con_struct *, dmy_struct *);
double calc_trans(double, double);
double calc_veg_displacement(double);
double calc_veg_height(double);
double calc_veg_roughness(double);
void   calc_water_balance_error(int, double, double, double);
double canopy_evap(layer_data_struct *, layer_data_struct *,
		   veg_var_struct *, veg_var_struct *, char, int, int, 
		   double, double *, double, double, double, double, 
		   double, double,  double, double, double, double, 
		   double *, double *, double *, double *, float *);
void   check_files(infiles_struct *, filenames_struct *);
FILE  *check_state_file(char *, dmy_struct *, global_param_struct *, int, int, 
                        int *);
void   close_files(infiles_struct *, outfiles_struct *, filenames_struct *);
filenames_struct cmd_proc(int argc, char *argv[]);
void   compress_files(char string[]);
void   compute_dz(double *, double *, int, double);
void   compute_penman_constants(double, double, double, double, double, 
				double, double, float, float, double *, 
				double *, double *, double *, double *);

void   correct_precip(double *, double, double, double, double);
void   compute_soil_layer_thermal_properties(layer_data_struct *, double *,
					     double *, double *, double *, 
					     int);
void   compute_treeline(atmos_data_struct *, dmy_struct *, double, double *, char *);
void   display_current_settings(int, filenames_struct *, global_param_struct *);
void   dist_prec(atmos_data_struct *,dist_prcp_struct *,soil_con_struct *,
		 veg_con_struct *,dmy_struct *,global_param_struct *,
		 outfiles_struct *,int,int,char,char, char, int);
#if QUICK_FS
void distribute_node_moisture_properties(double *, double *, double *,
					 double *, double *, double *,
					 double *, double ***, 
					 double *, double *, double *,
					 double *, double *, int, int, char);
#else
void distribute_node_moisture_properties(double *, double *, double *,
					 double *, double *, double *,
					 double *, double *, double *,
					 double *, double *, double *,
					 double *, double *, int, int, char);
#endif
void   distribute_soil_property(double *,double,double,
				double **l_param,
				int, int, double *, double *);

double error_calc_snow_ground_flux(double Tsurf, ...);
double error_calc_surf_energy_bal(double Tsurf, ...);
double ErrorSnowPackEnergyBalance(double Tsurf, ...);
double error_print_snow_ground_flux(double, va_list);
double ErrorPrintSnowPackEnergyBalance(double, va_list);
double error_print_solve_T_profile(double, va_list);
double error_print_surf_energy_bal(double, va_list);
double error_solve_T_profile(double Tsurf, ...);
double estimate_dew_point(double, double, double, double, double);
#if QUICK_FS
void estimate_layer_ice_content(layer_data_struct *, double *, double *,
				double *, double ***, double *,
				double *, double ***, double *,
				double *, double *, float **, int, int, char);
#else
void estimate_layer_ice_content(layer_data_struct *, double *, double *,
				double *, double *, double *, double *,
				double *, double *, double *, double *,
				double *, double *, float **, int, int, char);
#endif
double estimate_T1(double, double, double, double, double, double, double, 
		   double, double, double, double);
double exp_interp(double,double,double,double,double);

double f(double, double, double, double, double, double, double, double,
         double, double, int, double *, double, double, double, double *,
         double *, double *, double *, double *, double *);
void   find_0_degree_fronts(energy_bal_struct *,double *, double *, int);
layer_data_struct find_average_layer(layer_data_struct *, layer_data_struct *,
				     double, double);
void   find_sublayer_temperatures(layer_data_struct *, double *, double *,
				  double *, double, double, int, int);
void   finish_frozen_soil_calcs(energy_bal_struct *, layer_data_struct *,
				layer_data_struct *, layer_data_struct *,
				soil_con_struct *, int, int, double, 
				double *, double *, double *, double *);
void   free_atmos(int nrecs, atmos_data_struct **atmos);
void   free_dist_prcp(dist_prcp_struct *, int);
void   free_vegcon(veg_con_struct **);
void   full_energy(int, atmos_data_struct *, soil_con_struct *,
		   veg_con_struct *, dist_prcp_struct *,
		   dmy_struct *,global_param_struct *,int,char);
double func_aero_resist(double,double,double,double,double);
double func_snow_ground_flux(double, va_list);
double func_surf_energy_bal(double, va_list);

double get_avg_temp(double, double, double *, double *, int);
void   get_force_type(char *, int, int *);
global_param_struct get_global_param(filenames_struct *, FILE *);
void   get_next_time_step(int *, int *, int *, int *, int *, int);

double hermint(double, int, double *, double *, double *, double *, double *);
void   hermite(int, double *, double *, double *, double *, double *);
void   HourlyT(int, int, int *, double *, int *, double *, double *);

void   initialize_atmos(atmos_data_struct *, dmy_struct *, FILE **, double, 
			double, double, double, double, double, double, 
			double, double *, 
#if OUTPUT_FORCE
			char *, outfiles_struct *);
#else
			char *);
#endif

void   initialize_global();
void   initialize_model_state(dist_prcp_struct *, dmy_struct, double,
			      global_param_struct *, infiles_struct, int,
			      int, int, int, soil_con_struct *, 
			      veg_con_struct *, char *, int *);
void   initialize_new_storm(cell_data_struct ***, veg_var_struct ***,
			    int, int, int, double, double);
void   initialize_snow(snow_data_struct **,int,FILE *,int);
void   initialize_soil(cell_data_struct **, soil_con_struct *, int);
void   initialize_veg( veg_var_struct **, veg_con_struct *,
		       global_param_struct *);