本文档详细介绍了SWAT模型中控制天气输入的子程序 clicon，涉及降水、温度数据处理、辐射、风速和湿度生成，以及气候变化研究中的参数调整。它解释了输入变量如站点高度、降雨量、温度等，以及输出变量如每日辐射、降水和气温数据。
由于这个子程序不能一下全部粘贴上去，所以我只要分开附上：
subroutine clicon
!! ~ ~ ~ PURPOSE ~ ~ ~
!! this subroutine controls weather inputs to SWAT. Precipitation and
!! temperature data is read in and the weather generator is called to
!! fill in radiation, wind speed and relative humidity as well as
!! missing precipitation and temperatures. Adjustments for climate
!! changes studies are also made in this subroutine.
!! ~ ~ ~ INCOMING VARIABLES ~ ~ ~
!! name |units |definition
!! ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~
!! elevb(:,:) |m |elevation at center of band
!! elevb_fr(:,:)|none |fraction of subbasin area within elevation
!! |band
!! elevp(:) |m |elevation of precipitation gage station
!! elevt(:) |m |elevation of temperature gage station
!! hru_sub(:) |none |subbasin in which HRU is located
!! huminc(:,:) |none |monthly humidity adjustment. Daily values
!! |for relative humidity within the month are
!! |raised or lowered by the specified amount.
!! |(used in climate change studies)
!! id1 |julian date |first day of simulation in year
!! ifirstpet |none |potential ET data search code
!! |0 first day of potential ET data located in
!! | file
!! |1 first day of potential ET data not located
!! | in file
!! ipet |none |code for potential ET method
!! |0 Priestley-Taylor method
!! |1 Penman/Monteith method
!! |2 Hargreaves method
!! |3 read in daily potential ET values
!! irgage(:) |none |HRU rain gage data code (gage # for rainfall
!! |data used in HRU)
!! itgage(:) |none |HRU temperature gage data code (gage # for
!! |temperature data used in HRU)
!! iyr |year |year currently being simulated (eg 1980)
!! i_mo |none |current month of simulation
!! nhru |none |number of HRUs in watershed
!! nstep |none |number of lines of rainfall data for each
!! |day
!! pcpsim |none |rainfall input code
!! |1 gage read for each subbasin
!! |2 gage simulated for each subbasin
!! plaps(:) |mm H2O/km |precipitation lapse rate: precipitation
!! |increase due to increase in elevation
!! radinc(:,:) |MJ/m^2 |monthly solar radiation adjustment. Daily
!! |radiation within the month is raised or
!! |lowered by the specified amount. (used in
!! |climate change studies)
!! rfinc(:,:) |% |monthly rainfall adjustment. Daily rainfall
!! |within the month is adjusted to the specified
!! |percentage of the original value (used in
!! |climate change studies)
!! rhsim |none |relative humidity input code
!! |1 measured data read for each subbasin
!! |2 data simulated for each subbasin
!! slrsim |none |solar radiation input code
!! |1 measured data read for each subbasin
!! |2 data simulated for each subbasin
!! tlaps(:) |deg C/km |temperature lapse rate: temperature increase
!! |due to increase in elevation
!! tmpinc(:,:) |deg C |monthly temperature adjustment. Daily maximum
!! |and minimum temperatures within the month are
!! |raised or lowered by the specified amount
!! |(used in climate change studies)
!! tmpsim |none |temperature input code
!! |1 daily max/min read for each subbasin
!! |2 daily max/min simulated for each subbasin
!! welev(:) |m |elevation of weather station used to compile
!! |weather generator data
!! wndsim |none |wind speed input code
!! |1 measured data read for each subbasin
!! |2 data simulated for each subbasin
!! ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~
!! ~ ~ ~ OUTGOING VARIABLES ~ ~ ~
!! name |units |definition
!! ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~
!! frad(:,:) |none |fraction of solar radiation occuring during
!! |hour in day in HRU
!! hhsubp(:,:) |mm H2O |precipitation falling during hour in day in
!! |HRU
!! hru_ra(:) |MJ/m^2 |solar radiation for the day in HRU
!! hru_rmx(:) |MJ/m^2 |maximum solar radiation for the day in HRU
!! ifirstpet |none |potential ET data search code
!! |0 first day of potential ET data located in
!! | file
!! |1 first day of potential ET data not located
!! | in file
!! pcpband(:,:)|mm H2O |precipitation for the day in band in HRU
!! petmeas |mm H2O |potential ET value read in for day
!! rainsub(:,:)|mm H2O |precipitation for the time step during the
!! |day in HRU
!! rhd(:) |none |relative humidity for the day in HRU
!! subp(:) |mm H2O |precipitation for the day in HRU
!! tavband(:,:)|deg C |average temperature for the day in band in HRU
!! tmn(:) |deg C |minimum temperature for the day in HRU
!! tmnband(:,:)|deg C |minimum temperature for the day in band in HRU
!! tmpav(:) |deg C |average temperature for the day in HRU
!! tmx(:) |deg C |maximum temperature for the day in HRU
!! tmxband(:,:)|deg C |maximum temperature for the day in band in HRU
!! u10(:) |m/s |wind speed for the day in HRU
!! ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~
!! ~ ~ ~ LOCAL DEFINITIONS ~ ~ ~
!! name |units |definition
!! ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~
!! fradbsb(:) |none |hourly solar radiation fractions for subbasin
!! ib |none |counter
!! idap |julain date |day currently being simulated
!! ii |none |counter
!! inum3sprev |none |subbasin number of previous HRU
!! iyp |none |year currently being simulated
!! k |none |counter
!! pdif |mm H2O |difference in precipitation for station and
!! |precipitation for elevation band
!! rabsb |MJ/m^2 |generated solar radiation for subbasin
!! ratio |none |fraction change in precipitation due to
!! |elevation changes
!! rbsb |mm H2O |generated precipitation for subbasin
!! rhdbsb |none |generated relative humidity for subbasin
!! rmxbsb |MJ/m^2 |generated maximum solar radiation for subbasin
!! tdif |deg C |difference in temperature for station and
!! |temperature for elevation band
!! tmnbsb |deg C |generated minimum temperature for subbasin
!! tmxbsb |deg C |generated maximum temperature for subbasin
!! u10bsb |m/s |generated wind speed for subbasin
扫一扫
8969
被折叠的 条评论
为什么被折叠?
到【灌水乐园】发言
