c********************************************************************** c c********************************************************************** subroutine saeinp c c +++ PURPOSE +++ c print out input file for stand alone erosion c c +++ ARGUMENT DECLARATIONS +++ c +++ ARGUMENT DEFINITIONS +++ c c +++ PARAMETER +++ c c + + + GLOBAL COMMON BLOCKS + + + c compiler instr. - no warn of unreferenced symbols in include files *$noreference include 'p1werm.inc' include 'm1subr.inc' include 'p1const.inc' include 'b1glob.inc' c include 'b1geom.inc' include 'm1geo.inc' include 'w1wind.inc' include 'w1pavg.inc' include 's1dbh.inc' include 's1layr.inc' include 's1phys.inc' include 's1agg.inc' include 's1surf.inc' include 's1sgeo.inc' include 'h1db1.inc' include 'm1flag.inc' include 'm1sim.inc' include 'wpath.inc' *$reference c c + + + LOCAL COMMON BLOCKS + + + *$noreference include 'erosion/e2grid.inc' include 'erosion/e3grid.inc' include 'erosion/m2geo.inc' include 'erosion/s2sgeo.inc' include 'erosion/s2agg.inc' include 'erosion/s2surf.inc' include 'erosion/e2erod.inc' c include 'erosion/w2wind.inc' *$reference c c +++ LOCAL VARIABLES +++ integer i,j,k,l, wustfl, icsr, outfl integer n integer b, day, mon, yr real wuref, a real wzzo, wzzov real tempegt c real wr c c + + + LOCAL VARIABLE DEFINITIONS + + + c i = c wustfl= c icsr = c outfl = c nhill = c wuref = c a = c wzzo = c wus = c wust = c wusp = c wr = c c c +++ SAVES +++ c c +++ SUBROUTINES CALLED +++ c sbzo c sbwus c sbwust c sbinit c sbigrd c sbwind c sberod c c sbout c sbsfdi c c +++ FUNCTION DECLARATIONS +++ integer lentrm c c +++ END SPECIFICATIONS +++ c c code for output call caldat (am0jd,day,mon,yr) if ((day.eq.31).and.(mon.eq.8).and.(yr.eq.29)) then call fopenk (42, rootp(1:lentrm(rootp)) // 'saeros.inp', * 'unknown') write(42,2101) day, mon, yr 2101 format('# day mon yr',3(i4)) do 23 j = jmax, 0, -1 write (40,10) (egt(i,j), i = 0, imax) 23 continue do 24 j = jmax, 0, -1 write (40,10) (egtss(i,j), i = 0, imax) 24 continue 10 format (1x, 31f14.8) tempegt = 0. do 550 j = 1, jmax-1 do 545 i = 1, imax-1 tempegt = tempegt+egt(i,j) 545 continue 550 continue write(42,2102) tempegt 2102 format('#tempegt', f15.8) write(42,*) '0' write(42,*) '.TRUE.' write(42,*) '1' write(42,2100) 2100 format ('# amxsim(1,1),amxsim(2,1),amxsim(1,2),amxsim(2,2)') write(42,*) amxsim(1,1), amxsim(2,1), amxsim(1,2),amxsim(2,2) write(42,2105) 2105 format('# amasim') write(42,*) amasim write(42,2110) 2110 format('# nacctr') write(42,*) nacctr write(42,2115) 2115 format('# amxar(1,1), amxar(2,1), amxar(1,2),amxar(2,2)') write(42,*) amxar(1,1,nacctr), amxar(2,1,nacctr), & amxar(1,2,nacctr),amxar(2,2,nacctr) c barriers write(42,2120) 2120 format ('# nbr') write(42,*) nbr do 33 b = 1,nbr write(42,2125) 2125 format('# amxar(1,1,b), amxar(2,1), amxar(1,2),amxar(2,2)') write(42,*) amxbr(1,1,b),amxbr(2,1,b),amxbr(1,2,b),amxbr(2,2,b) write (42,2130) 2130 format('# amzbr(b), ampbr(b), amxbrw(b)') write (42,*) amzbr(b), ampbr(b), amxbrw(b) 33 continue c subregions write(42,2135) 2135 format('# nsubr') write(42,*) '1' write(42,2140) 2140 format('# amxsr(1,1), amxsr(2,1), amxsr(1,2),amxsr(2,2)') write(42,*) amxsr(1,1,1),amxsr(2,1,1),amxsr(1,2,1),amxsr(2,2,1) c crop write(42,2145) 2145 format('# abzht(1)') write(42,*) abzht(1) write(42,2150) 2150 format('# abrsai(1), abrlai(1)') write(42,*) abrsai(1), abrlai(1) write(42,2155) 2155 format ('# abffcv(1)') write(42,*) abffcv(1) c soil write(42,2160) 2160 format('# nslay(1)') write(42,*) nslay(1) write(42,2165) 2165 format( '# aszlyt(l,1)') write(42,*) (aszlyt(l,1), l=1,nslay(1)) write(42,2170) 2170 format('# asdblk(l,1)') write(42,*) (asdblk(l,1), l=1,nslay(1)) write(42,2175) 2175 format('# asfsan(l,1)') write(42,*) (asfsan(l,1), l=1,nslay(1)) write(42,2180) 2180 format('# asfsil(l,1)') write(42,*) (asfsil(l,1), l=1,nslay(1)) write(42,2185) 2185 format('# asfcla(l,1)') write(42,*) (asfcla(l,1), l=1,nslay(1)) write(42,2190) 2190 format('# asvroc(l,1)') write(42,*) (asvroc(l,1), l=1,nslay(1)) write(42,2195) 2195 format('# asdagd(l,1)') write(42,*) (asdagd(l,1), l=1,nslay(1)) write(42,2200) 2200 format('# aseags(l,1)') write(42,*) (aseags(l,1), l=1,nslay(1)) write(42,2205) 2205 format('# aslagm(l,1)') write(42,*) (aslagm(l,1), l=1,nslay(1)) write(42,2210) 2210 format('#aslagn(l,1)') write(42,*) (aslagn(l,1), l=1,nslay(1)) write(42,2215) 2215 format('# aslagx(l,1)') write(42,*) (aslagx(l,1), l=1,nslay(1)) write(42,2220) 2220 format('# as0ags(l,1)') write(42,*) (as0ags(l,1), l=1,nslay(1)) write(42,2225) 2225 format('# asfcr(1) aszcr(1) asflos(1) asfmos(1) asdcr(1)', & ' asecr(1)') write(42,*) asfcr(1), aszcr(1), asflos(1), asmlos(1), asdcr(1), & asecr(1) write(42,2230) 2230 format('#aslrr(1)') write(42,*) aslrr(1) write(42,2235) 2235 format('# aszrgh(1) asxrgs(1) asxrgw(1) asargo(1)') write(42,*) aszrgh(1), asxrgs(1), asxrgw(1), asargo(1) write(42,2240) 2240 format('# asxdks(1)') write(42,*) asxdks(1) write(42,2245) c hydrology 2245 format('# ahzsnd(1)') write(42,*) ahzsnd(1) write(42,2250) 2250 format('# ahrwcw(l,1), l=1,nslay(1)') write(42,*) (ahrwcw(l,1), l=1,nslay(1)) write(42,2255) 2255 format('# ahrwca(l,1), l=1,nslay(1)') write(42,*) (ahrwca(l,1), l=1,nslay(1)) write(42,2260) 2260 format('# ahrwc0(l,1), l=1,12') write(42,22) (ahrwc0(l,1), l=1,12) write(42,2265) 2265 format('# ahrwc0(l,1), l=13,24') 22 format(12(1x,f10.8)) write(42,22) (ahrwc0(l,1), l=13,24) c weather write(42,2270) 2270 format('# awdair') write(42,*) awdair write(42,2275) 2275 format('# awadir') write(42,*) awadir write(42,2280) 2280 format('# ntstep') write(42,*) ntstep write(42,2285) 2285 format('# anemht awzzo wzoflg') write(42,*) anemht, awzzo, wzoflg write(42,2290) 2290 format('#wflg (read wind speeds only, no Weibull params)') write(42,*) '1' write(42,2295) 2295 format('# awu(k), k=1,24') write(42,*) (awu(k), k=1,6) write(42,*) (awu(k), k=7,12) write(42,*) (awu(k), k=13,18) write(42,*) (awu(k), k=19,24) end if c return end c++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++