Main Subroutines:
| aero_flux_pilp_gra.c | aero_flux_pilp_gra | |||
| twolayer_vic_roff.c | twolayer_vic_roff | Computes runoff from VIC grid cells. | ||
| melt_2.c | melt_2 | Calculates snow accumulation and melt using an energy balance. | ||
| snow_met.c | snow_met | air_temperature, ppt, cold_content, snow_table | snow_albedo, isnow, last_snow | Computes meteorological variables and snow albedo for snow accumulation and melt computations |
| snow_table_albedo.c | snow_table_albedo | none | snow_table | Creates table for snow albedos with snow age |
| long_c.c | long_c | rhm, temp, cloud_cover | long_wave | Computes long wave radiation using relative humidity, temperature, and cloud cover. |
| rad_c.c | rad_c | lat, lon, elev, ntime, julian, cloud_cover | short_wave | Computes short wave radiation. |
| met_inline.c | met_inline | Calculates meteorological variables from base weather station data. | ||
| flux_rtsafe_pilp.c | flux_rtsafe_pilp | x1 ,x2, i, layer, Tsur0, k | rts | ??? |
| flux_Tsur_fun_pilp.c | flux_Tsur_fun_pilp | Tsur, i, layer, Tsur0, k | f, df | ??? |
| evap_coef_pilp_snow.c | evap_coef_pilp_snow | Compute evaporation for snow?? | ||
| evap_coef_pilp_gra_Ts.c | evap_coef_pilp_gra_Ts | Compute surface temperature and evaporation from grass?? | ||
| evap_coef_pilp_gra.c | evap_coef_pilp_gra | Compute evaporation from grass |
Utility Subroutines:
| end_of_month.c | end_of_month | day, month, year | TRUE or FALSE | Determines whether the end of the month has been reached. |
Input and Output Subroutines:
| read_pilp_data.c | read_pilp_data | Reads in PILPs forcing data | ||
Snow Melt Subroutine:
Variables:
| Variable Type | Variable Name | Description |
| double | Rs | Incoming shortwave radiation at snow surface |
| double | Ls | Downward longwave radiation from the sky |
| double | aero_cond | Aerodynamic conductance over snow |
| double | air_temperature | Air temperature |
| double | air_dens | Density of air |
| double | gamma | Psychometric constant |
| double | ppt | Precipitation depth after interception |
| double | pack_temp_min | Minimum allowable snow pack temperature |
| double | vel_2m | Wind velocity at 2 m above the surface |
| double | air_pressure | Air pressure |
| double | air_vp | Air vapor pressure |
| double | snow_coeff | Coefficient used in determining snow_temp (value 0.0 - 1,0) |
| double * | swq | Snow water equivalent |
| double * | snow_vapor_flux | Depth of water evaporation, sublimation, or condensation |
| double * | pack_water | Liquid water content of the snow pack |
| double * | cold_content | Cold content of the snow pack (<= 0C) |
| double * | snow_bal | Sum of mass balance errors |
| double * | Qs | Sensible heat transfer |
| double * | Qe | Latent heat transfer |
| double * | snow_melt | Amount of snow melted in depth |
| double * | snow_temp | Temperature of the snow surface |
| double * | net_rad | Net radiation exchange at the snow surface |
Called Subroutines: NONE
Met Inline Subroutine:
Variables:
| Variable Type | Variable Name | Description |
| double | PPT | Base station precipitation |
| double | dew_temp | Base station dew point temperature |
| double | qh | Specific humidity |
| double | air_temperature | Air temperature |
| double | air_pressure | Air pressure |
| double * | gamma | Psychometric constant |
| double | hv | Latent heat of vaporization |
| double * | air_dens | Air density |
| double * | air_vp | Air vapor density |
Called Subroutines: NONE
Read PILP Subroutine:
Variables:
| Variable Type | Variable Name | Description |
| int | num_recs | Number of records to be read |
| double * | prec1hr | 1 hour precipitation in mm/hr |
| double * | soldn | Solar radiation??? |
| double * | lwdn | Long wave radiation ??? |
| double * | wind | Wind |
| double * | pressure | Air pressure |
| double * | specific_humid | Specific humidity |
| double * | air_temp | Air temperature |
| double * | penman_temp | Temperature for evaporation ??? (penmen eqn) |
Called Subroutines: NONE
Snow Met Subroutine:
Variables:
| Variable Type | Variable Name | Description |
| double | air_temperature | Air temperature above snow pack |
| double | ppt | Depth of precipitation reaching snow pack |
| double | cold_content | Cold content of snow pack |
| double ** | snow_table | Look-up table of snow albedo based on age and season |
| int * | isnow | Snow cover flag ( 1 = snow cover, 0 = no snow) |
| double * | snow_albedo | albedo of the snow surface |
| int * | last_snow | Number of days since last snow |
Called Subroutines: NONE
Evaporation Coefficients for PILP Snow Subroutine:
Variables:
| Variable Type | Variable Name | Description |
| int | num_recs | Number of records |
| double * | soldn | ?? |
| double * | lwdn | |
| double * | pan_evap1hr | 1 hour pan evaporation rate?, amount? |
| double * | aero_resist | Aerodynamic resistance |
| double * | wind | Wind |
| double * | specific_humid | Specific humidity |
| double * | pressure | Air pressure |
| double * | air_temp | Air temperature |
| double * | penman_temp | Evaporation temperature |
| double | hour_coef[25] |
Called Subroutines: NONE
Evapoartion Coefficients for PILP grass Ts Subroutine:
Variables:
| Variable Type | Variable Name | Description |
| int | num_recs | Number of records |
| int | iglobal | ?? |
| int | iyear | |
| int | vegetat_type_num | |
| int * | month | |
| double * | rhum | Relative humidity |
| double * | soldn | Solar downwelling radiation?? |
| double * | lwdn | Long wave downwelling radiation?? |
| double * | pan_evap1hr | 1 hour pan evaporation |
| double * | wind | Wind |
| double * | specific_humid | Specific humidity |
| double * | pressure | Air pressure |
| double * | penman_temp | Temperature |
| double * | vege_par | Vegitation parameter? |
| double ** | ratio_archi_to_aero | Ratio of architectural resistance to aerodynamic resistance |
| double ** | ratio_canop_to_aero | Ratio of canopy resistance to aerodynamic resistance |
| double ** | groot_percent1 | Percent of roots in layer 1? |
| double ** | groot_percent2 | Percent of roots in layer 2? |
| double ** | Cv | Specific heat |
| double ** | LAI | Leaf Area Index |
| double | Fraction | Fraction of vegetiation type coverage? |
| double ** | z0 | |
| double ** | d0 | |
| double ** | raph | Average albedo? |
| double ** | z2 | |
| double ** | Cv_aero_resist | |
| double ** | Cv_mod_Tsur | |
| double * | g_Cv_aero_resist | |
| double * | g_wind_z2_vege | |
| double * | g_wind_z2_soil |
Called Subroutines: NONE
Evapoartion Coefficients for PILP grass Subroutine:
Seems to be variation of previous subroutine.
Variables:
| Variable Type | Variable Name | Description |
| int | num_recs | Number of records |
| int | iglobal | ?? |
| int | iyear | |
| int | vegetat_type_num | |
| int * | month | |
| double * | rhum | Relative humidity |
| double * | soldn | Solar downwelling radiation?? |
| double * | lwdn | Long wave downwelling radiation?? |
| double * | pan_evap1hr | 1 hour pan evaporation |
| double ** | xrnet | Writes net radiation |
| double * | wind | Wind |
| double * | specific_humid | Specific humidity |
| double * | pressure | Air pressure |
| double * | penman_temp | Temperature |
| double * | vege_par | Vegitation parameter? |
| double ** | ratio_archi_to_aero | Ratio of architectural resistance to aerodynamic resistance |
| double ** | ratio_canop_to_aero | Ratio of canopy resistance to aerodynamic resistance |
| double ** | groot_percent1 | Percent of roots in layer 1? |
| double ** | groot_percent2 | Percent of roots in layer 2? |
| double ** | Cv | Specific heat |
| double ** | LAI | Leaf Area Index |
| double | Fraction | Fraction of vegetiation type coverage? |
| double ** | z0 | |
| double ** | d0 | |
| double ** | raph | Average albedo? |
| double ** | z2 | |
| double ** | Cv_aero_resist | |
| double ** | Cv_mod_Tsur | |
| double * | g_Cv_aero_resist | |
| double * | g_wind_z2_vege | |
| double * | g_wind_z2_soil |
Called Subroutines: NONE
Aerodynamic flux for PILP grass Subroutine:
Variables:
| Variable Type | Variable Name | Description |
| int | itry | |
| int | num_recs | Number of records |
| int | iglobal | ?? |
| int | iyear | |
| int | vegetat_type_num | |
| int * | day | |
| int * | month | |
| int * | snow_index | |
| double * | soldn | Solar downwelling radiation?? |
| double * | lwdn | Long wave downwelling radiation?? |
| double | prec1hr | Precipitation per time step |
| double * | evap1hr | Evaporation per time step |
| double | C | Heat capacity of air (1204.8 Jm-3K-1) |
| double * | K | Saturated soil conductivity? (0.866 for trees, 1.03 for grass) |
| double | D | Depth of lower soil layer |
| double | Za | Height of the air layer in the energy balance equation |
| double | DT | Time step? (set to 3600*TS) |
| double | avg_soil_temp | Average soil temperature at depth D |
| double | mixalbedo | |
| double * | air_temp | Air temperature |
| double ** | Cv_aero_resist | |
| double ** | Cv_mod_Tsur | |
| double ** | Cv_mod_sensible | |
| double ** | Cv_mod_grnd | |
| double * | mod_sensible | |
| double * | mod_grnd | |
| double * | mod_Tsur | |
| double ** | raph | Average albedo? |
| double ** | Cv | |
| double ** | const_k | const_k[k]=Cs/2/DT |
| double * | pilp_rnet | |
| FILE * | out_flux | |
| double | sum_year_rnet | |
| double | sum_year_evap | |
| double | sum_year_sensible | |
| double | sum_year_grnd | |
| double | sum_year_rnet1 | |
| double | sum_year_evap1 | |
| double | sum_year_sensible1 | |
| double | sum_year_grnd1 |
Called Subroutines: flux_rtsafe_pilp
Flux RTSAFE PILP Subroutine:
Variables:
| Variable Type | Variable Name | Description |
| double | x1 | 270.0 |
| double | x2 | 315.0 |
| int | i | ?? |
| int | layer | layer = 1: assume there is an air-layer in the energy balance
layer = 0: assume an ideal energy balance equation applies |
| double | Tsur0 | Surface temperature for vegetation type |
| int | k | vegetation type |
Called Subroutines: flux_Tsur_fun_pilp, nrerror
Flux Tsur fun PILP Subroutine:
Variables:
| Variable Type | Variable Name | Description |
| double | Tsur | Surface temperature?? |
| double * | f | ?? |
| double * | df | ?? |
| int | i | ?? |
| int | layer | ?? |
| double | Tsur0 | ?? |
| int | k | ?? |
| int | month | ?? |
| double | C | Heat capacity of air (1204.8 Jm-3K-1) |
| double | K[8] | |
| double | D | Depth of the lower layer |
| double | Za | |
| double | DT | |
| double | avg_soil_temp | Average soil temperature at depth D |
| double ** | Cv_evap1hr | |
| double ** | Cv_mod_Tsur | |
| double * | soldn | |
| double * | lwdn | |
| double * | air_temp | |
| double ** | Cv_aero_resist | |
| double | raph[13][8] | Average albedo |
| double | const_k[8] |
Called Subroutines: NONE
Read PET Data Subroutine:
Variables:
| Variable Type | Variable Name | Description |
| int | num_recs | Number of records |
| double * | soldn | Short wave radiation |
| double * | lwdn | Long wave radiation |
| double * | prec1hr | Precipitation on bare ground (plus melt)? |
| double * | wind | Wind |
| double * | air_temp | Air temperature |
| double * | penman_temp | Penman temperature |
| double * | melt | Snow melt |
| double * | rhum | Relative humidity |
| double * | tskc | Total cloud cover percentage |
| char ** | date | |
| int * | julian_day | Julian day of year |
| char * | infile | Input data file |
| char * | infile_sa | Surface airways input data file |
Called Subroutines: NONE
Snow Melt Subroutine:
Variables:
| Variable Type | Variable Name | Description |
| double | Rs | Incoming shortwave radiation at snow surface |
| double | Ls | Downward longwave radiation from the sky |
Called Subroutines: NONE