stdaa.f

水文模型的原始代码

!!                               |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(1)  |mm H2O        |average amount of precipitation in watershed
!!                               |for the simulation
!!    wshdaao(3)  |mm H2O        |surface runoff in watershed for simulation
!!    wshdaao(4)  |mm H2O        |lateral flow contribution to streamflow in
!!                               |watershed for simulation
!!    wshdaao(5)  |mm H2O        |water percolation past bottom of soil profile
!!                               |in watershed for simulation
!!    wshdaao(6)  |mm H2O        |water yield to streamflow from HRUs in
!!                               |watershed for simulation
!!    wshdaao(7)  |mm H2O        |actual evapotranspiration in watershed
!!                               |for simulation
!!    wshdaao(11) |metric tons/ha|net change in sediment of reservoirs in
!!                               |watershed during simulation
!!    wshdaao(12) |metric tons/ha|sediment yield from HRUs in watershed for
!!                               |the simulation
!!    wshdaao(13) |metric tons/ha|sediment loading to ponds in watershed 
!!                               |during simulation
!!    wshdaao(14) |metric tons/ha|sediment loading from ponds in watershed
!!                               |during simulation
!!    wshdaao(15) |metric tons/ha|net change in sediment level in ponds in
!!                               |watershed during simulation
!!    wshdaao(19) |mm H2O        |evaporation from ponds in watershed during
!!                               |simulation
!!    wshdaao(20) |mm H2O        |seepage from ponds in watershed during
!!                               |simulation
!!    wshdaao(21) |mm H2O        |precipitation on ponds in watershed during
!!                               |simulation
!!    wshdaao(22) |mm H2O        |volume of water entering ponds in watershed
!!                               |during simulation
!!    wshdaao(23) |mm H2O        |volume of water leaving ponds in watershed
!!                               |during simulation
!!    wshdaao(33) |mm H2O        |net change in water volume of ponds in
!!                               |watershed during simulation
!!    wshdaao(34) |mm H2O        |net change in water volume of reservoirs in
!!                               |watershed during simulation
!!    wshdaao(36) |mm H2O        |snow melt in watershed for simulation
!!    wshdaao(38) |mm H2O        |average amount of tributary channel
!!                               |transmission losses in watershed during
!!                               |simulation
!!    wshdaao(39) |mm H2O        |freezing rain/snow fall in watershed for 
!!                               |the simulation
!!    wshdaao(40) |kg N/ha       |organic N loading to stream in watershed for
!!                               |the simulation
!!    wshdaao(41) |kg P/ha       |organic P loading to stream in watershed for
!!                               |the simulation
!!    wshdaao(42) |kg N/ha       |nitrate loading to stream in surface runoff
!!                               |in watershed for the simulation
!!    wshdaao(43) |kg P/ha       |soluble P loading to stream in watershed for
!!                               |the simulation
!!    wshdaao(44) |kg N/ha       |plant uptake of N in watershed for the 
!!                               |simulation
!!    wshdaao(45) |kg N/ha       |nitrate loading to stream in lateral flow
!!                               |in watershed for the simulation
!!    wshdaao(46) |kg N/ha       |nitrate percolation past bottom of soil
!!                               |profile in watershed for the simulation
!!    wshdaao(104)|mm H2O        |groundwater contribution to stream in
!!                               |watershed for the simulation
!!    wshdaao(105)|mm H2O        |amount of water moving from shallow aquifer
!!                               |to plants/soil profile in watershed during
!!                               |simulation
!!    wshdaao(106)|mm H2O        |deep aquifer recharge in watershed during
!!                               |simulation
!!    wshdaao(107)|mm H2O        |total amount of water entering both aquifers
!!                               |in watershed during simulation
!!    wshdaao(108)|mm H2O        |potential evapotranspiration in watershed
!!                               |for the simulation
!!    wshdaao(109)|mm H2O        |drainage tile flow contribution to stream
!!                               |in watershed for the simulation
!!    yldaa(:)    |metric tons/ha|average annual yield (dry weight) in HRU
!!    yldn(:,:,:) |kg/ha         |average value for yield of crop
!!    ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~

!!    ~ ~ ~ LOCAL DEFINITIONS ~ ~ ~
!!    name        |units         |definition
!!    ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~
!!    j           |none          |counter
!!    k           |none          |counter
!!    nicr        |none          |sequence number for crop in year
!!    nnro        |none          |sequence number for year in rotation
!!    sumpady     |none          |number of HRUs with potholes
!!    xirr        |mm H2O        |average annual amount of irrigation water
!!                               |applied to watershed
!!    ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~

!!    ~ ~ ~ ~ ~ ~ END SPECIFICATIONS ~ ~ ~ ~ ~ ~

      use parm
 
      real :: xirr
      integer :: j, nnro, nicr, k, sumpady, ncrp

      if (iscen == 1) then
      write (2,1000) prog
      write (2,1100) title
      write (2,1200)
      else if (isproj == 1) then
      write (19,1000) prog
      write (19,1100) title
      write (19,1200)
      endif
        xirr = 0.
        xirr = Sum(aairr)
        !! print irrigation data
        if (xirr > 0.) then
          if (iscen == 1) then
          write (2,1300)
          else if (isproj == 1) then
          write (19,1300)
          endif
          do j = 1, nhru
            if (iscen == 1) then
            if (aairr(j) > 1.e-4) write (2,1400) j, irn(j), aairr(j)
            else if (isproj == 1) then
            if (aairr(j) > 1.e-4) write (19,1400) j, irn(j), aairr(j)
            endif
          end do
        end if
     
!! write average crop information to std output file
      ncrp = 0
      if (mcr < 3) then
        ncrp = mcr
      else
        ncrp = 3
      end if
      if (iscen == 1) then
      write (2,1000) prog
      write (2,1100) title
      write (2,1500)
      else if (isproj == 1) then
      write (19,1000) prog
      write (19,1100) title
      write (19,1500)
      endif
      do j = 1, nhru
        do nnro = 1, nrot(j)
          if (iscen == 1) then
          write (2,1600) cpnm(idplt(nro(j),icr(j),j)), j, nnro,         &
     &        (yldn(nnro,nicr,j), bio_aahv(nnro,nicr,j),nicr = 1, ncrp)
          else if (isproj == 1) then
          write (19,1600) cpnm(idplt(nro(j),icr(j),j)), j, nnro,        &
     &        (yldn(nnro,nicr,j), bio_aahv(nnro,nicr,j),nicr = 1, ncrp)
          endif
        end do
      end do

!! write average annual HRU data

      if (iscen == 1) then
      write (2,1700)
      write (2,1800)
      do j = 1, nhru
        write (2,1900) j, hru_sub(j), cpnm(idplt(nro(j),icr(j),j)),     &
     &     snam(j), hru_km(j), cn2(j), sol_sumfc(j), usle_ls(j),        &
     &     hruaao(22,j), hruaao(28,j), hruaao(29,j), sumix(j),          &
     &     hruaao(1,j), hruaao(4,j), hruaao(5,j) + hruaao(6,j),         &
     &     hruaao(12,j), hruaao(14,j), hruaao(37,j) + hruaao(38,j),     &
     &     hruaao(35,j), bio_aams(j), yldaa(j)
      end do
      else if (isproj == 1) then
      write (19,1700)
      write (19,1800)
      do j = 1, nhru
        write (19,1900) j, hru_sub(j), cpnm(idplt(nro(j),icr(j),j)),    &
     &     snam(j), hru_km(j), cn2(j), sol_sumfc(j), usle_ls(j),        &
     &     hruaao(22,j), hruaao(28,j), hruaao(29,j), sumix(j),          &
     &     hruaao(1,j), hruaao(4,j), hruaao(5,j) + hruaao(6,j),         &
     &     hruaao(12,j), hruaao(14,j), hruaao(37,j) + hruaao(38,j),     &
     &     hruaao(35,j), bio_aams(j), yldaa(j)
      end do
      endif

!! write average annual watershed monthly values
      if (iscen == 1) then
      write (2,2000)
      do j = 1, 12
        write (2,2100) j, (wshd_aamon(j,k), k = 1, 8)
      end do
      else if (isproj == 1) then
      write (19,2000)
      do j = 1, 12
        write (19,2100) j, (wshd_aamon(j,k), k = 1, 8)
      end do
      endif
	
!! write average annual stress values
      if (iscen == 1) then
      write (2,2200) wshd_wstrs, wshd_tstrs, wshd_nstrs, wshd_pstrs
      else if (isproj == 1) then
      write (19,2200) wshd_wstrs, wshd_tstrs, wshd_nstrs, wshd_pstrs
      endif

!! watershed summary water balance table