writeaa.f

水文模型的原始代码

      subroutine writeaa

!!    ~ ~ ~ PURPOSE ~ ~ ~
!!    this subroutine writes average annual output

!!    ~ ~ ~ INCOMING VARIABLES ~ ~ ~
!!    name         |units         |definition
!!    ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~
!!    aairr(:)     |mm H2O        |average annual amount of irrigation water
!!                                |applied to HRU
!!    bio_aams(:)  |metric tons/ha|average annual biomass (dry weight) in HRU
!!    bio_hv(:,:,:)|kg/ha         |harvested biomass (dry weight)
!!    da_ha        |ha            |area of watershed in hectares
!!    hru_dafr(:)  |none          |fraction of watershed area in HRU
!!    hruaao       |varies        |HRU average annual output array
!!    hrupest(:)   |none          |pesticide use flag:
!!                                | 0: no pesticides used in HRU
!!                                | 1: pesticides used in HRU
!!    hrupsta(:,1,:)|mg pst       |amount of pesticide type in surface runoff
!!                                |contribution to stream from HRU during 
!!                                |simulation (in solution)
!!    hrupsta(:,2,:)|mg pst       |amount of pesticide type in surface runoff
!!                                |contribution to stream from HRU during
!!                                |simulation (sorbed to sediment)
!!    icodes(:)     |none         |routing command code:
!!                                |0 = finish       9 = save
!!                                |1 = subbasin    10 = recday
!!                                |2 = route       11 = reccnst
!!                                |3 = routres     12 = structure
!!                                |4 = transfer    13 = 
!!                                |5 = add         14 = saveconc
!!                                |6 = rechour     15 = 
!!                                |7 = recmon      16 = autocal
!!                                |8 = recyear
!!    idaf         |julian date   |beginning day of simulation
!!    idal         |julian date   |ending day of simulation
!!    ihouts(:)    |none          |For ICODES equal to
!!                                |0: not used
!!                                |1,2,3,5,7,8,10,11: hydrograph storage
!!                                |                     location number
!!                                |4: water source type
!!                                |   (1=reach)
!!                                |   (2=reservoir)
!!                                |9: hydrograph storage location of data to
!!                                |   be printed to event file
!!                                |14:hydrograph storage location of data to
!!                                |   be printed to saveconc file
!!    inum1s(:)    |none          |For ICODES equal to
!!                                |0: not used
!!                                |1: HRU number
!!                                |2: reach number
!!                                |3: reservoir number
!!                                |4: reach or res # flow is diverted from
!!                                |5: hydrograph storage location of 1st
!!                                |   dataset to be added
!!                                |7,8,9,10,11,14: file number
!!    inum2s(:)    |none          |For ICODES equal to
!!                                |0,1,7,8,10,11: not used
!!                                |2,3: inflow hydrograph storage location
!!                                |4: destination type
!!                                |   (1=reach)
!!                                |   (2=reservoir)
!!                                |5: hydrograph storage location of 2nd
!!                                |   dataset to be added
!!                                |9,14:print frequency
!!                                |   (0=daily)
!!                                |   (1=hourly)
!!    inum3s(:)    |none          |For ICODES equal to
!!                                |0,1,2,3,5,7,8,10,11: not used
!!                                |4: destination number. Reach or
!!                                |   reservoir receiving water
!!                                |9: print format
!!                                |   (0=normal, fixed format)
!!                                |   (1=txt format for AV interface,recday)
!!    iprp         |none          |print code for output.pst file
!!                                |0 do not print pesticide output
!!                                |1 print pesticide output
!!    irn(:)       |none          |average annual number of irrigation
!!                                |applications in HRU
!!    lai_aamx(:)  |none          |average annual maximum leaf area index in
!!                                |HRU
!!    leapyr       |none          |leap year flag:
!!                                |0  leap year
!!                                |1  regular year
!!    mcr          |none          |max number of crops grown per year
!!    nbyr         |none          |number of calendar years simulated
!!    ncrops(:,:,:)|
!!    ndays(:)     |julian date   |julian date for last day of preceding
!!                                |month (where the array location is the
!!                                |number of the month). The dates are for
!!                                |leap years
!!    ndmo(:)      |days          |cumulative number of days accrued in the
!!                                |month since the simulation began where the
!!                                |array location number is the number of the
!!                                |month
!!    nhru         |none          |number of HRUs in watershed
!!    npmx         |none          |number of different pesticides used in
!!                                |the simulation
!!    nres         |none          |number of reservoirs in watershed
!!    nrot(:)      |none          |number of years of rotation
!!    nyskip       |none          |number of years to skip output
!!                                |summarization and printing
!!    rchaao       |varies        |reach average annual output array
!!    resouta(3,:) |metric tons   |sediment entering reservoir during simulation
!!    resouta(4,:) |metric tons   |sediment leaving reservoir during simulation
!!    resouta(17,:)|m^3 H2O       |evaporation from reservoir during simulation
!!    resouta(18,:)|m^3 H2O       |seepage from reservoir during simulation
!!    resouta(19,:)|m^3 H2O       |precipitation on reservoir during simulation
!!    resouta(20,:)|m^3 H2O       |water flowing into reservoir during 
!!                                |simulation
!!    resouta(21,:)|m^3 H2O       |water flowing out of reservoir during
!!                                |simulation
!!    sbactlchlp   |# colonies/ha |average annual number of less persistent
!!                                |bacteria lost from soil surface layer by
!!                                |percolation
!!    sbactlchp    |# colonies/ha |average annual number of persistent bacteria
!!                                |lost from soil surface layer by percolation
!!    sbactrolp    |# colonies/ha |average annual number of less persistent
!!                                |bacteria transported to main channel
!!                                |with surface runoff in solution
!!    sbactrop     |# colonies/ha |average annual number of persistent bacteria
!!                                |transported to main channel with surface
!!                                |runoff in solution
!!    sbactsedlp   |# colonies/ha |average annual number of less persistent
!!                                |bacteria transported with sediment in
!!                                |surface runoff
!!    sbactsedp    |# colonies/ha |average annual number of persistent bacteria
!!                                |transported with sediment in surface runoff
!!    sdiegrolpq   |# colonies/ha |average annual change in the number of
!!                                |less persistent bacteria colonies in soil
!!                                |solution in watershed
!!    sdiegrolps   |# colonies/ha |average annual change in the number of
!!                                |less persistent bacteria colonies on soil
!!                                |particles in watershed
!!    sdiegropq    |# colonies/ha |average annual change in the number of
!!                                |persistent bacteria colonies in soil solution
!!                                |in watershed
!!    sdiegrops    |# colonies/ha |average annual change in the number of
!!                                |persistent bacteria colonies on soil particles
!!                                |in watershed
!!    shyd(1,:)    |m^3 H2O       |water
!!    shyd(2,:)    |metric tons   |sediment or suspended solid load
!!    shyd(3,:)    |kg N          |organic nitrogen
!!    shyd(4,:)    |kg P          |organic phosphorus
!!    shyd(5,:)    |kg N          |nitrate
!!    shyd(6,:)    |kg P          |soluble phosphorus
!!    shyd(7,:)    |mg pst        |soluble pesticide
!!    sno3up       |kg N/ha       |amount of nitrate moving upward in the soil
!!                                |profile in watershed
!!    sol_actp(:,:)|kg P/ha       |amount of phosphorus stored in the
!!                                |active mineral phosphorus pool
!!    sol_aorgn(:,:)|kg N/ha      |amount of nitrogen stored in the active
!!                                |organic (humic) nitrogen pool
!!    sol_fon(:,:) |kg N/ha       |amount of nitrogen stored in the fresh
!!                                |organic (residue) pool
!!    sol_fop(:,:) |kg P/ha       |amount of phosphorus stored in the fresh
!!                                |organic (residue) pool
!!    sol_nly(:)   |none          |number of layers in soil profile
!!    sol_no3(:,:) |kg N/ha       |amount of nitrogen stored in the
!!                                |nitrate pool.
!!    sol_orgn(:,:)|kg N/ha       |amount of nitrogen stored in the stable
!!                                |organic N pool
!!    sol_orgp(:,:)|kg P/ha       |amount of phosphorus stored in the organic
!!                                |P pool
!!    sol_solp(:,:)|kg P/ha       |amount of phosohorus stored in solution
!!    sol_stap(:,:)|kg P/ha       |amount of phosphorus in the soil layer
!!                                |stored in the stable mineral phosphorus pool
!!    spadyev      |mm H2O        |average annual amount of water removed
!!                                |from potholes by evaporation in watershed
!!    spadyo       |mm H2O        |average annual amount of water released to
!!                                |main channel from potholes in watershed
!!    spadyrfv     |mm H2O        |average annual amount of precipitation on
!!                                |potholes in watershed
!!    spadysp      |mm H2O        |average annual amount of water removed
!!                                |from potholes by seepage in watershed
!!    subaao       |varies        |subbasin average annual output array
!!    sumix(:)     |none          |sum of mixing efficiencies in HRU
!!    wshd_aamon(:,:)|varies        |array of watershed monthly average values
!!    wshd_dnit    |kg N/ha       |average annual amount of nitrogen lost from
!!                                |nitrate pool due to denitrification in
!!                                |watershed
!!    wshd_fixn    |kg N/ha       |average annual amount of nitrogen added to
!!                                |plant biomass via fixation
!!    wshd_fminp   |kg P/ha       |average annual amount of mineral P applied
!!                                |in watershed
!!    wshd_fnh3    |kg N/ha       |average annual amount of NH3-N applied in
!!                                |watershed
!!    wshd_fno3    |kg N/ha       |average annual amount of NO3-N applied in
!!                                |watershed
!!    wshd_forgn   |kg N/ha       |average annual amount of organic N applied
!!                                |in watershed
!!    wshd_forgp   |kg P/ha       |average annual amount of organic P applied
!!                                |in watershed
!!    wshd_ftotn   |kg N/ha       |average annual amount of N (mineral &
!!                                |organic) applied in watershed
!!    wshd_ftotp   |kg P/ha       |average annual amount of P (mineral &
!!                                |organic) applied in watershed
!!    wshd_hmn     |kg N/ha       |average annual amount of nitrogen moving
!!                                |from active organic to nitrate pool in
!!                                |watershed
!!    wshd_hmp     |kg P/ha       |average annual amount of phosphorus moving
!!                                |from organic to labile pool in watershed
!!    wshd_nitn    |kg N/ha       |average annual amount of nitrogen moving
!!                                |from the NH3 to the NO3 pool by
!!                                |nitrification in the watershed
!!    wshd_nstrs   |stress units  |average annual number of nitrogen stress
!!                                |units in watershed
!!    wshd_pal     |kg P/ha       |average annual amount of phosphorus moving
!!                                |from labile mineral to active mineral pool
!!                                |in watershed
!!    wshd_pas     |kg P/ha       |average annual amount of phosphorus moving
!!                                |from active mineral to stable mineral pool
!!                                |in watershed
!!    wshd_plch    |kg P/ha       |average annual amount of phosphorus leached
!!                                |into second soil layer
!!    wshd_pstap(:)|kg pst/ha     |total amount of pesticide type applied in
!!                                |watershed during simulation
!!    wshd_pstdg(:)|kg pst/ha     |amount of pesticide lost through degradation
!!                                |in watershed
!!    wshd_pstrs   |stress units  |average annual number of phosphorus stress
!!                                |units in watershed
!!    wshd_pup     |kg P/ha       |average annual amount of plant uptake of
!!                                |phosphorus
!!    wshd_raino3  |kg N/ha       |average annual amount of NO3 added to soil
!!                                |by rainfall in watershed
!!    wshd_rmn     |kg N/ha       |average annual amount of nitrogen moving
!!                                |from fresh organic (residue) to nitrate
!!                                |and active organic pools in watershed
!!    wshd_rmp     |kg P/ha       |average annual amount of phosphorus moving
!!                                |from fresh organic (residue) to labile
!!                                |and organic pools in watershed
!!    wshd_rwn     |kg N/ha       |average annual amount of nitrogen moving
!!                                |from active organic to stable organic pool
!!                                |in watershed
!!    wshd_tstrs   |stress units  |average annual number of temperature stress
!!                                |units in watershed
!!    wshd_voln    |kg N/ha       |average annual amount if nitrogen lost by
!!                                |ammonia volatilization in watershed
!!    wshd_wstrs   |stress units  |average annual number of water stress units
!!                                |in watershed
!!    wshd_yldn    |kg N/ha       |amount of nitrogen removed from soil in
!!                                |watershed in the yield
!!    wshd_yldp    |kg P/ha       |amount of phosphorus removed from soil in
!!                                |watershed in the yield
!!    wshdaao      |varies        |watershed average annual output array
!!    wtraa(:,:)   |varies        |HRU impoundment average annual output array
!!    yldaa(:)     |metric tons/ha|average annual yield (dry weight) in HRU
!!    yldkg(:,:,:) |kg/ha         |yield (dry weight) by crop type in the HRU
!!    ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~

!!    ~ ~ ~ OUTGOING VARIABLES ~ ~ ~
!!    name        |units         |definition
!!    ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~
!!    basminpf    |kg P/ha       |final average amount of phosphorus in
!!                               |the mineral P pool in watershed soil
!!    basno3f     |kg N/ha       |final average amount of nitrogen in the
!!                               |nitrate pool in watershed soil
!!    basorgnf    |kg N/ha       |final average amount of nitrogen in the
!!                               |organic N pool in watershed soil