#######################################################################
# VIC Model Parameters - 4.1.x
#######################################################################
# $Id: global.param.sample,v 5.7.2.12 2009/06/16 06:27:02 vicadmin Exp $
#######################################################################
# Simulation Parameters
#######################################################################
NLAYER          3       # number of soil layers
NODES          	5       # number of soil thermal nodes  #not known the function of number of nodes. Used 5 as used in previous version.
TIME_STEP       3       # model time step in hours (set to 24 if FULL_ENERGY = FALSE, set to < 24 if FULL_ENERGY = TRUE)
SNOW_STEP       3       # time step in hours for which to solve the snow model (should = TIME_STEP if TIME_STEP < 24)
STARTYEAR       1980    # year model simulation starts
STARTMONTH      01      # month model simulation starts
STARTDAY        01      # day model simulation starts
STARTHOUR       00      # hour model simulation starts
ENDYEAR         2012    # year model simulation ends
ENDMONTH        12      # month model simulation ends
ENDDAY          31      # day model simulation ends
FULL_ENERGY     TRUE    # TRUE = calculate full energy balance; FALSE = compute water balance only
FROZEN_SOIL     TRUE    # TRUE = calculate frozen soils #used true
QUICK_FLUX      FALSE   # TRUE = use simplified ground heat flux method of Liang et al (1999); FALSE = use finite element method of Cherkauer et al (1999) #used false as in older version
QUICK_SOLVE     FALSE   # TRUE = Use Liang et al., 1999 formulation for iteration, but explicit finite difference method for final step.
NO_FLUX         FALSE   # TRUE = use no flux lower boundary for ground heat flux computation; FALSE = use constant flux #lower boundary condition.  If NO_FLUX = TRUE, QUICK_FLUX MUST = FALSE
IMPLICIT        FALSE   # TRUE = use implicit solution for soil heat flux equation of Cherkauer et al (1999), otherwise uses		#original explicit solution.
EXP_TRANS       FALSE   # TRUE = exponentially distributes the thermal nodes in the Cherkauer et al. (1999) finite difference 
				#algorithm, otherwise uses linear distribution
SNOW_ALBEDO     USACE   # USACE = use traditional VIC algorithm based on US Army Corps of Engineers empirical snow albedo 
				#decay curves, using hard-coded dates for transitions from snow accumulation to melting; 
				#SUN1999 = use algorithm of Sun et al 1999, in which albedo decay depends on snow cold content 
				#(more appropriate for simulations outside the US).
SNOW_DENSITY    DENS_BRAS       # DENS_BRAS = use traditional VIC algorithm taken from Bras, 1990; DENS_SNTHRM = use 
				#algorithm taken from SNTHRM model.
BLOWING         FALSE   # TRUE = compute evaporative fluxes due to blowing snow
DIST_PRCP       FALSE   # TRUE = use distributed precipitation
PREC_EXPT       0.6     # exponent for use in distributed precipitation eqn (only used if DIST_PRCP is TRUE)
CORRPREC        FALSE   # TRUE = correct precipitation for gauge undercatch
MIN_WIND_SPEED  0.1     # minimum allowable wind speed (m/s)
MAX_SNOW_TEMP   1.5     # maximum temperature (C) at which snow can fall
MIN_RAIN_TEMP   -0.5    # minimum temperature (C) at which rain can fall
CONTINUEONERROR FALSE   # TRUE = continue simulation on cell error; TFALLBACK = when energy balance fails to converge, use previous T value
COMPUTE_TREELINE        FALSE   # Can be either FALSE or the id number of an understory veg class; FALSE = turn 
				#treeline computation off; VEG_CLASS_ID= replace any overstory veg types with the this 
				#understory veg type in all snow bands for which the average July Temperature <= 10 C 
				#(e.g. "COMPUTE_TREELINE 10" replaces any overstory veg cover with class 10)
EQUAL_AREA      FALSE   # TRUE = grid cells are from an equal-area projection; FALSE = grid cells are on a regular lat-lon grid
RESOLUTION      0.041666666666668   # Grid cell resolution (degrees if EQUAL_AREA is FALSE, km^2 if EQUAL_AREA is TRUE); ignored if LAKES is FALSE
#AERO_RESIST_CANSNOW     AR_410  # Options for aerodynamic resistance in snow-filled canopy:
                                # AR_406 = multiply by 10 for latent heat but do NOT multiply by 10 for sensible heat
                                #          and do NOT apply stability correction (as in VIC 4.0.6).
                                #          When no snow in canopy, use surface aero_resist for ET.
                                # AR_406_LS = multiply by 10 for latent heat AND sensible heat
                                #             and do NOT apply stability correction.
                                #             When no snow in canopy, use surface aero_resist for ET.
                                # AR_406_FULL = multiply by 10 for latent heat AND sensible heat
                                #               and do NOT apply stability correction;
                                #               Additionally, always use overstory aero_resist for
                                #               ET (as in 4.1.0).
                                # AR_410 = apply stability correction but do NOT multiply by 10 (as in VIC 4.1.0)
                                #          Additionally, always use overstory aero_resist for
                                #          ET (as in 4.1.0).
                                # AR_COMBO = multiply by 10 AND apply stability correction
                                #            Additionally, always use overstory aero_resist for
                                #            ET (as in 4.1.0).
                                # default = AR_410
# GRND_FLUX_TYPE  GF_FULL # Options for ground flux:
                        # GF_406  = use (flawed) formulas for ground flux, deltaH, and fusion
                        #           from VIC 4.0.6 and earlier
                        # GF_410  = use formulas from VIC 4.1.0 (ground flux is correct,
                        #           but deltaH and fusion ignore surf_atten)
                        # GF_FULL = use correct ground flux formula from VIC 4.1.0 and
                        #           also take surf_atten into account in deltaH and fusion
#MIN_LIQ         FALSE   # TRUE = use temperature-dependent minimum for liquid soil moisture instead of residual moisture
#PLAPSE          FALSE   # This controls how VIC computes air pressure when air pressure is not supplied as an input forcing:
                        # TRUE = set air pressure to sea level pressure, lapsed to grid cell average elevation;
                        # FALSE = set air pressure to constant 95.5 kPa (as in all versions of VIC pre-4.1.1)

#######################################################################
# State Files and  Parameters
#######################################################################
#INIT_STATE     (put the initial state path/filename here)      # Initial state path/file
#STATENAME      (put the path/prefix of output state file here) # Output state file path/prefix.  The date 
			#(STATEYEAR,STATEMONTH,STATEDAY) will be appended to the prefix automatically in the format yyyymmdd.
#STATEYEAR      2000    # year to save model state
#STATEMONTH     12      # month to save model state
#STATEDAY       31      # day to save model state
#BINARY_STATE_FILE       FALSE  # TRUE if state file should be binary format; FALSE if ascii

#######################################################################
# Forcing Files and Parameters
#
#       All FORCING filenames are actually the pathname, and prefix
#               for gridded data types: ex. DATA/forcing_
#               Latitude and longitude index suffix is added by VIC
#
#       There must be 1 FORCE_TYPE entry for each variable (column) in the forcing file
#
#       If FORCE_TYPE is BINARY, each FORCE_TYPE must be followed by:
#                       SIGNED/UNSIGNED SCALE_FACTOR
#               For example (BINARY):
#                       FORCE_TYPE      PREC    UNSIGNED        40
#               or (ASCII):
#                       FORCE_TYPE      PREC
#######################################################################
FORCING1        ../../basins/5020004/forcings/Obs2012C/data_ # Forcing file path and prefix, ending in "_"
FORCE_FORMAT    ASCII  # BINARY or ASCII
#FORCE_ENDIAN    LITTLE  # LITTLE (PC/Linux) or BIG (SUN)
N_TYPES         4       # Number of variables (columns)
FORCE_TYPE      PREC    #UNSIGNED        40
FORCE_TYPE      TMAX    #SIGNED  100
FORCE_TYPE      TMIN    #SIGNED  100
FORCE_TYPE      WIND    #SIGNED  100
FORCE_DT        24      # Forcing time step length (hours)
FORCEYEAR       1980    # Year of first forcing record
FORCEMONTH      01      # Month of first forcing record
FORCEDAY        01      # Day of first forcing record
FORCEHOUR       00      # Hour of first forcing record
GRID_DECIMAL    4       # Number of digits after decimal point in forcing file names
WIND_H          10.0    # height of wind speed measurement (m)
MEASURE_H       2.0     # height of humidity measurement (m)
ALMA_INPUT      FALSE   # TRUE = ALMA-compliant input variable units; FALSE = standard VIC units

#######################################################################
# Land Surface Files and Parameters
#######################################################################
SOIL            ../../basins/5020004/params/vic/soil170406/soil_5020004_cal7835.asc # Soil parameter path/file
#ARC_SOIL        FALSE   # TRUE = read soil parameters from ARC/INFO ASCII grids
#SOIL_DIR        (soil param directory)   # Directory containing ARC/INFO ASCII grids of soil parameters - only valid if ARC_SOIL is TRUE
BASEFLOW        ARNO    # ARNO = columns 5-8 are the standard VIC baseflow parameters; NIJSSEN2001 = columns 5-8 
				#of soil file are baseflow parameters from Nijssen et al (2001)
#JULY_TAVG_SUPPLIED      FALSE   # TRUE = final column of the soil parameter file will contain average July 
				#air temperature, for computing treeline; this will be ignored if COMPUTE_TREELINE 
				#is FALSE; FALSE = compute the treeline based on the average July air temperature of 
				#the forcings over the simulation period
VEGPARAM        ../../basins/5020004/params/vic/veg_20160308/veg_5020004_cal7835.asc # Veg parameter path/file
GLOBAL_LAI      TRUE    # TRUE = read LAI from veg param file; FALSE = read LAI from veg library file
ROOT_ZONES      3       # Number of root zones (must match format of veg param file)
VEGLIB          ../../basins/5020004/params/vic/veg_lib.asc # Veg library path/file
SNOW_BAND       5       ../../basins/5020004/params/vic/elev5/elev_5020004_cal7835.asc    # number of snow bands and snow band filename # Number of snow bands; if number of snow bands > 1, you must insert the snow band 
				#path/file after the number of bands (e.g. SNOW_BAND 5 my_path/my_snow_band_file)
#LAKES           (put lake parameter path/file here)     # Lake parameter path/file
#LAKE_PROFILE    FALSE   # TRUE = User-specified depth-area parameters in lake parameter file; FALSE = VIC 
				#computes a parabolic depth-area profile

#######################################################################
# Output Files and Parameters
#######################################################################
RESULT_DIR      ../../basins/5020004/results/vic/output140704      # Results directory path (with trailing )
OUT_STEP        24       # Output interval (hours); if 0, OUT_STEP = TIME_STEP
SKIPYEAR        0       # Number of years of output to omit from the output files
#COMPRESS        FALSE   # TRUE = compress input and output files when done
BINARY_OUTPUT   FALSE   # TRUE = binary output files
ALMA_OUTPUT     FALSE   # TRUE = ALMA-format output files; FALSE = standard VIC units
#MOISTFRACT      FALSE   # TRUE = output soil moisture as volumetric fraction; FALSE = standard VIC units
#PRT_HEADER      FALSE   # TRUE = insert a header at the beginning of each output file; FALSE = no header
PRT_SNOW_BAND   TRUE   # TRUE = write a "snowband" output file, containing band-specific values of snow variables
                        # NOTE: this is ignored if N_OUTFILES is specified below.
						
# If FROZEN_SOIL is TRUE, increase N_OUTFILES by 1
# If PRT_SNOW_BAND is TRUE, increase N_OUTFILES by 1

N_OUTFILES	2

OUTFILE         fluxes  24
OUTVAR          OUT_PREC        %.4f    OUT_TYPE_FLOAT  1
OUTVAR          OUT_EVAP        %.4f    OUT_TYPE_FLOAT  1
OUTVAR          OUT_RUNOFF      %.4f    OUT_TYPE_FLOAT  1
OUTVAR          OUT_BASEFLOW    %.4f    OUT_TYPE_FLOAT  1
OUTVAR          OUT_AIR_TEMP    %.4f    OUT_TYPE_FLOAT  1
OUTVAR          OUT_SURF_TEMP   %.4f    OUT_TYPE_FLOAT  1
OUTVAR          OUT_EVAP_BARE   %.4f    OUT_TYPE_FLOAT  1
OUTVAR          OUT_PET_SATSOIL %.4f    OUT_TYPE_FLOAT  1
OUTVAR          OUT_PET_H2OSURF %.4f    OUT_TYPE_FLOAT  1
OUTVAR          OUT_LONGWAVE    %.4f    OUT_TYPE_FLOAT  1
OUTVAR          OUT_SHORTWAVE   %.4f    OUT_TYPE_FLOAT  1
OUTVAR          OUT_SENSIBLE    %.4f    OUT_TYPE_FLOAT  1
OUTVAR          OUT_LATENT      %.4f    OUT_TYPE_FLOAT  1
OUTVAR          OUT_REL_HUMID   %.4f    OUT_TYPE_FLOAT  1
OUTVAR          OUT_VP          %.4f    OUT_TYPE_FLOAT  1
OUTVAR          OUT_VPD         %.4f    OUT_TYPE_FLOAT  1
OUTVAR          OUT_WIND        %.4f    OUT_TYPE_FLOAT  1
OUTVAR          OUT_PRESSURE    %.4f    OUT_TYPE_FLOAT  1
OUTVAR          OUT_DENSITY     %.4f    OUT_TYPE_FLOAT  1
OUTVAR          OUT_QAIR        %.4f    OUT_TYPE_FLOAT  1
OUTVAR          OUT_ALBEDO      %.4f    OUT_TYPE_FLOAT  1
OUTVAR          OUT_TRANSP_VEG  %.4f    OUT_TYPE_FLOAT  1
OUTVAR          OUT_INFLOW      %.4f    OUT_TYPE_FLOAT  1
OUTVAR          OUT_SOIL_MOIST  %.4f    OUT_TYPE_FLOAT  1

OUTFILE         snow  5
OUTVAR          OUT_SWE         %.4f    OUT_TYPE_FLOAT  1
OUTVAR          OUT_SNOW_CANOPY %.4f    OUT_TYPE_FLOAT  1
OUTVAR          OUT_SNOW_COVER  %.4f    OUT_TYPE_FLOAT  1
OUTVAR          OUT_SNOW_DEPTH  %.4f    OUT_TYPE_FLOAT  1
OUTVAR          OUT_SWE_BAND    %.4f    OUT_TYPE_FLOAT  1
