⭐ 欢迎来到虫虫下载站! | 📦 资源下载 📁 资源专辑 ℹ️ 关于我们
⭐ 虫虫下载站

📄 stdaa.f

📁 水文模型的原始代码
💻 F
📖 第 1 页 / 共 3 页
字号:
      if (iscen == 1) then
      write (2,1000) prog
      write (2,1100) title
      write (2,2300) wshdaao(1), wshdaao(39), wshdaao(36), wshdaao(37), &
     &    wshdaao(3), wshdaao(4), wshdaao(109), wshdaao(104),           &
     &    wshdaao(105), wshdaao(106), wshdaao(107), wshdaao(6),         &
     &    wshdaao(5), wshdaao(7), wshdaao(108), wshdaao(38),            &
     &    wshdaao(12), (wshdaao(k),k = 19,22), wshdaao(13),             &
     &    wshdaao(23), wshdaao(14), (resdata(k),k = 1,7)
      write (2,2400) wshdaao(33), wshdaao(15), wshdaao(34), wshdaao(11)
      else if (isproj == 1) then
      write (19,1000) prog
      write (19,1100) title
      write (19,2300) wshdaao(1), wshdaao(39), wshdaao(36), wshdaao(37),&
     &    wshdaao(3), wshdaao(4), wshdaao(109), wshdaao(104),           &
     &    wshdaao(105), wshdaao(106), wshdaao(107), wshdaao(6),         &
     &    wshdaao(5), wshdaao(7), wshdaao(108), wshdaao(38),            &
     &    wshdaao(12), (wshdaao(k),k = 19,22), wshdaao(13),             &
     &    wshdaao(23), wshdaao(14), (resdata(k),k = 1,7)
      write (19,2400) wshdaao(33), wshdaao(15), wshdaao(34), wshdaao(11)
      endif

!! watershed pothole summary values
      sumpady = 0
      sumpady = Sum(ipot)
      if (sumpady > 0) then
        if (iscen == 1) then
        write (2,2500) spadyo, spadyev, spadysp, spadyrfv
        else if (isproj == 1) then
        write (19,2500) spadyo, spadyev, spadysp, spadyrfv
        endif
      end if

!! watershed summary nutrient table
      if (iscen == 1) then
      write (2,1000) prog
      write (2,1100) title
      write (2,2600)
      write (2,2700) wshdaao(40), wshdaao(41), wshdaao(42), wshdaao(45),&
     &    wshdaao(43), wshdaao(46), wshd_plch, wshdaao(44), wshd_pup,   &
     &    wshdaao(110)
      write (2,2800) wshd_pal, wshd_pas, wshd_ftotn, wshd_ftotp,        &
     &    wshd_fixn, wshd_dnit
      write (2,2900) wshd_hmn, wshd_rwn, wshd_hmp, wshd_rmn, wshd_rmp
      write (2,3000) wshd_raino3, basno3i, basno3f, basorgni, basorgnf, &
     &    basminpi, basminpf, basorgpi, basorgpf, wshd_fno3,            &
     &    wshd_fnh3, wshd_forgn, wshd_fminp, wshd_forgp, wshd_yldn,     &
     &    wshd_yldp, wshd_voln, wshd_nitn, sno3up
      else if (isproj == 1) then
      write (19,1000) prog
      write (19,1100) title
      write (19,2600)
      write (19,2700) wshdaao(40), wshdaao(41), wshdaao(42),            &
     &    wshdaao(45), wshdaao(43), wshdaao(46), wshd_plch, wshdaao(44),&
     &    wshd_pup, wshdaao(110)
      write (19,2800) wshd_pal, wshd_pas, wshd_ftotn, wshd_ftotp,       &
     &    wshd_fixn, wshd_dnit
      write (19,2900) wshd_hmn, wshd_rwn, wshd_hmp, wshd_rmn, wshd_rmp
      write (19,3000) wshd_raino3, basno3i, basno3f, basorgni, basorgnf,&
     &    basminpi, basminpf, basorgpi, basorgpf, wshd_fno3,            &
     &    wshd_fnh3, wshd_forgn, wshd_fminp, wshd_forgp, wshd_yldn,     &
     &    wshd_yldp, wshd_voln, wshd_nitn, sno3up
      endif

!! watershed bacteria summary table
      if (iscen == 1) then
      write (2,3100) sdiegropq,sdiegrolpq,sdiegrops,sdiegrolps,         &
     &     sbactrop,sbactrolp,sbactsedp,sbactsedlp,sbactlchp,sbactlchlp
      else if (isproj == 1) then
      write (19,3100) sdiegropq,sdiegrolpq,sdiegrops,sdiegrolps,        &
     &     sbactrop,sbactrolp,sbactsedp,sbactsedlp,sbactlchp,sbactlchlp
      endif

      return

 1000 format ('1',/t5,a80,t105,2(a2,'/'),a2,5x,2(i2,':'),i2)
 1100 format (/(t5,20a4))
 1200 format (t5,'FINAL VALUES'/)
 1300 format (//t21,'IRRIGATION - AVE. ANNUAL',/,t14,'HRU',t27,         &
     &    'NO.OF ',t44,'VOLUME',/,t14,'NO.',t24,'APPLICATIONS',t42,     &
     &    'APPLIED(MM)',/)
 1400 format (11x,i6,9x,i4,13x,f8.3)
 1500 format (41x,'Average Crop Values',/,26x,'Crop 1',15x,'Crop 2',15x,&
     &    'Crop 3',/,20x,'Yld',t30,'Biomass',t41,'Yld',t51,'Biomass',   &
     &    t62,'Yld',t72,'Biomass',/,t19,'(kg/ha)',t30,'(kg/ha)',t40,    &
     &    '(kg/ha)',t51,'(kg/ha)',t61,'(kg/ha)',t72,'(kg/ha)',/)
 1600 format (1x,a4,2x,'HRU ',i6,' Rot ',i3,6(f8.1,2x))
 1700 format (/t5,'HRU STATISTICS'//t17,'AVE ANNUAL VALUES'/)
 1800 format (1x,'HRU',t6,' SUB',t11,'CPMN',t16,'SOIL',t24,'AREAkm2',   &
     & t36,'CN',                                                        &
     & t42,'AWCmm',t48,'USLE_LS',t58,'IRRmm',t64,'AUTONkh ',t72,        &
     & 'AUTOPkh ',t82,'MIXEF',t89,'PRECmm',t97,'SURQmm',t106,           &
     & 'GWQmm',t115,'ETmm',t122,'SEDth ',t129,'NO3kgh ',t136,           &
     & 'ORGNkgh ',t145,'BIOMth',t154,'YLDth')
 1900 format (i6,i4,1x,a4,a8,e8.3,16f8.2)
 2000 format (///,t17,'AVE MONTHLY BASIN VALUES',/t20,'SNOW',t46,       &
     &   'WATER',t66,'SED',/t3,'MON',t11,'RAIN',t20,'FALL',t27,'SURF Q',&
     &    t37,'LAT Q',t46,'YIELD',t58,'ET',t64,'YIELD',t75,'PET',/t11,  &
     &    '(MM)',t20,'(MM)',t29,'(MM)',t38,'(MM)',t47,'(MM)',t56,'(MM)',&
     &    t63,'(T/HA)',t74,'(MM)')
 2100 format (i5,14f9.2)
 2200 format (/,t5,' AVE ANNUAL BASIN STRESS DAYS',/,t15,               &
     &    ' WATER STRESS DAYS = ',f8.2,/,t15,                           &
     &    ' TEMPERATURE STRESS DAYS = ',f8.2,/,t15,                     &
     &    ' NITROGEN STRESS DAYS = ',f8.2,/,t15,                        &
     &    ' PHOSPHORUS STRESS DAYS = ',f8.2)
 2300 format (t10,'AVE ANNUAL BASIN VALUES'//                           &
     &        t15,'PRECIP = ',f8.1,' MM'/                               &
     &        t15,'SNOW FALL =',f8.2,' MM'/                             &
     &        t15,'SNOW MELT = ',f8.2,' MM'/                            &
     &        t15,'SUBLIMATION = ',f8.2,' MM'/                          &
     &        t15,'SURFACE RUNOFF Q = ',f8.2,' MM'/                     &
     &        t15,'LATERAL SOIL Q =',f8.2,' MM'/                        &
     &        t15,'TILE Q = ',f8.2,' MM'/                               &
     &        t15,'GROUNDWATER (SHAL AQ) Q = ',f8.2,' MM'/              &
     &        t15,'REVAP (SHAL AQ => SOIL/PLANTS) =',f8.2,' MM'/        &
     &        t15,'DEEP AQ RECHARGE = ',f8.2,' MM'/                     &
     &        t15,'TOTAL AQ RECHARGE =',f8.2,' MM'/                     &
     &        t15,'TOTAL WATER YLD = ',f8.2,' MM'/                      &
     &        t15,'PERCOLATION OUT OF SOIL =',f8.2,' MM'/               &
     &        t15,'ET = ',f8.1,' MM'/                                   &
     &        t15,'PET = ',f8.1,'MM'/                                   &
     &        t15,'TRANSMISSION LOSSES = ',f8.2,' MM'/                  &
     &        t15,'TOTAL SEDIMENT LOADING =  ',f8.3,' T/HA'/            &
     &        t15,'POND BUDGET'/                                        &
     &        t20,'EVAPORATION = ',f8.3,' MM'/                          &
     &        t20,'SEEPAGE = ',f8.3,' MM'/                              &
     &        t20,'RAINFALL ON POOL =',f8.3,' MM'/                      &
     &        t20,'INFLOW'/t25,'WATER = ',f8.3,' MM'/                   &
     &        t25,'SEDIMENT = ',f8.3,' T/HA'/                           &
     &        t20,'OUTFLOW'/t25,'WATER = ',f8.3,' MM'/                  &
     &        t25,'SEDIMENT = ',f8.3,' T/HA'/                           &
     &        t15,'RESERVOIR BUDGET'/                                   &
     &        t20,'EVAPORATION = ',f8.3,' MM'/                          &
     &        t20,'SEEPAGE = ',f8.3,' MM'/                              &
     &        t20,'RAINFALL ON RESERVOIR = ',f8.3,' MM'/                &
     &        t20,'INFLOW'/t25,'WATER = ',f8.3,' MM'/                   &
     &        t25,'SEDIMENT = ',f8.3,' T/HA'/                           &
     &        t20,'OUTFLOW'/t25,'WATER = ',f8.3,' MM'/                  &
     &        t25,'SEDIMENT = ',f8.3,' T/HA')
 2400 format (t15,'YIELD LOSS FROM PONDS'/t20,'WATER = ',f7.3,' MM'/t20,&
     &    'SEDIMENT = ',f7.3,' T/HA'/t15,'YIELD LOSS FROM RESERVOIRS'/  &
     &    t20,'WATER = ',f7.3,' MM'/t20,'SEDIMENT = ',f7.3,' T/HA')
 2500 format (t15,'OUTFLOW FROM IMPOUNDED WATER =  ',f8.3,' (MM)',/,t15,&
     &    'EVAPORATION FROM IMPOUNDED WATER =  ',f8.3,' (MM)',/,t15,    &
     &    'SEEPAGE INTO SOIL FROM IMPOUNDED WATER = ',f8.3,' (MM)',/,t15&
     &    ,'RAINFALL ON IMPOUNDED WATER = ',f8.3,' (MM)')
 2600 format (t15,'AVE ANNUAL BASIN VALUES')
 2700 format (//,t15,'NUTRIENTS',/,t20,'ORGANIC N =  ',f8.3,' (KG/HA)', &
     &    /,t20,'ORGANIC P =  ',f8.3,' (KG/HA)',/,t20,                  &
     &    'NO3 YIELD (SQ) =  ',f8.3,' (KG/HA)',/,t20,                   &
     &    'NO3 YIELD (SSQ) = ',f8.3,' (KG/HA)',/,t20,'SOL P YIELD = ',f8&
     &    .3,' (KG/HA)',/,t20,'NO3 LEACHED =  ',f8.3,' (KG/HA)',/,t20,  &
     &    'P LEACHED =  ',f8.3,' (KG/HA)',/,t20,                        &
     &    'N UPTAKE =  ',f8.3,' (KG/HA)',/,t20,'P UPTAKE = ',f8.3,      &
     &    ' (KG/HA)',/,t20,'NO3 YIELD (GWQ) =  ',f8.3,' (KG/HA)')
 2800 format (t20,'ACTIVE TO SOLUTION P FLOW =  ',f8.3,' (KG/HA)',/,t20,&
     &    'ACTIVE TO STABLE P FLOW =  ',f8.3,' (KG/HA)',/,t20,          &
     &    'N FERTILIZER APPLIED = ',f8.3,' (KG/HA)',/,t20,              &
     &    'P FERTILIZER APPLIED = ',f8.3,' (KG/HA)',/,t20,              &
     &    'N FIXATION = ',f8.3,' (KG/HA)',/,t20,'DENITRIFICATION = ',f8 &
     &    .3,' (KG/HA)')
 2900 format (t20,'HUMUS MIN ON ACTIVE ORG N =  ',f8.3,' (KG/HA)',/,t20,&
     &    'ACTIVE TO STABLE ORG N =  ',f8.3,' (KG/HA)',/,t20,           &
     &    'HUMUS MIN ON ACTIVE ORG P = ',f8.3,' (KG/HA)',/,t20,         &
     &    'MIN FROM FRESH ORG N = ',f8.3,' (KG/HA)',/,t20,              &
     &    'MIN FROM FRESH ORG P = ',f8.3,' (KG/HA)')
 3000 format (t20,'NO3 IN RAINFALL =  ',f8.3,' (KG/HA)',/,t20,          &
     &    'INITIAL NO3 IN SOIL =  ',f10.3,' (KG/HA)',/,t20,             &
     &    'FINAL NO3 IN SOIL =  ',f10.3,' (KG/HA)',/,t20,               &
     &    'INITIAL ORG N IN SOIL =  ',f10.3,' (KG/HA)',/,t20,           &
     &    'FINAL ORG N IN SOIL =  ',f10.3,' (KG/HA)',/,t20,             &
     &    'INITIAL MIN P IN SOIL =  ',f10.3,' (KG/HA)',/,t20,           &
     &    'FINAL MIN P IN SOIL =  ',f10.3,' (KG/HA)',/,t20,             &
     &    'INITIAL ORG P IN SOIL =  ',f10.3,' (KG/HA)',/,t20,           &
     &    'FINAL ORG P IN SOIL =  ',f10.3,' (KG/HA)',/,t20,             &
     &    'NO3 IN FERT =  ',f8.3,' (KG/HA)',/,t20,                      &
     &    'AMMONIA IN FERT =  ',f8.3,' (KG/HA)',/,t20,                  &
     &    'ORG N IN FERT =  ',f8.3,' (KG/HA)',/,t20,                    &
     &    'MINERAL P IN FERT =  ',f8.3,' (KG/HA)',/,t20,                &
     &    'ORG P IN FERT =  ',f8.3,' (KG/HA)',/,t20,                    &
     &    'N REMOVED IN YIELD =  ',f8.3,' (KG/HA)',/,t20,               &
     &    'P REMOVED IN YIELD =  ',f8.3,' (KG/HA)',/,t20,               &
     &    'AMMONIA VOLATILIZATION =  ',f8.3,' (KG/HA)',/,t20,           &
     &    'AMMONIA NITRIFICATION = ',f8.3,' (KG/HA)',/,t20,             &
     &    'NO3 EVAP-LAYER 2 TO 1 = ',f8.3)
 3100 format (/,t20,'DIE-GRO P Q =  ',f8.1,' (No/HA)',/,t20,            &
     &    'DIE-GRO LP Q =  ',f8.1,' (No/HA)',/,t20,                     &
     &    'DIE-GRO P SED = ',f8.1,' (No/HA)',/,t20,                     &
     &    'DIE-GRO LP SED = ',f8.1,' (No/HA)',/,t20,                    &
     &    'BACT P RUNOFF = ',f8.1,' (No/HA)',/,t20,                     &
     &    'BACT LP RUNOFF = ',f8.1,' (No/HA)',/,t20,                    &
     &    'BACT P SEDIMENT = ',f8.1,' (No/HA)',/,t20,                   &
     &    'BACT LP SEDIMENT = ',f8.1,' (No/HA)',/,t20,                  &
     &    'BACT P INCORP = ',f8.1,' (No/HA)',/,t20,                     &
     &    'BACT LP INCORP = ',f8.1,' (No/HA)',/)
      end

⌨️ 快捷键说明

复制代码 Ctrl + C
搜索代码 Ctrl + F
全屏模式 F11
切换主题 Ctrl + Shift + D
显示快捷键 ?
增大字号 Ctrl + =
减小字号 Ctrl + -