c file: 'tdbug.for' subroutine tdbug(isr,slay,output) c + + + PURPOSE + + + c This program prints out many of the global variables before c and after the call to various MANAGEMENT practices c + + + KEY WORDS + + + c wind, erosion, tillage, soil, crop, decomposition c management c + + + GLOBAL COMMON BLOCKS + + + *$noereference include 'p1werm.inc' include 's1layr.inc' include 's1surf.inc' include 's1phys.inc' include 's1agg.inc' include 's1dbh.inc' include 's1dbc.inc' include 's1sgeo.inc' include 'c1gen.inc' include 'c1glob.inc' include 'd1glob.inc' include 'h1hydro.inc' include 'h1scs.inc' include 'h1db1.inc' include 'decomp/decomp.inc' c + + + LOCAL COMMON BLOCKS + + + include 'main/main.inc' *$reference c + + + LOCAL VARIABLES + + + integer i,l,isr,slay,output c + + + LOCAL DEFINITIONS + + + c l - This variable is an index on soil layers. c + + + SUBROUTINES CALLED + + + c + + + FUNCTIONS CALLED + + + c + + + UNIT NUMBERS FOR INPUT/OUTPUT DEVICES + + + c * = screen and keyboard c 29 = debug MANAGement c + + + DATA INITIALIZATIONS + + + c + + + INPUT FORMATS + + + c + + + OUTPUT FORMATS + + + 2040 format(3x,'aslagn aslagx aslagm as0ags') 2041 format(3x,'asdblk asdsblk aszlyt') 2050 format (1x,4f7.2) 2051 format (1x,f7.2,2x,f7.2,2x,f7.2) 2060 format (1x,i4,1x,f7.2,1x,f7.2,f6.2,4f7.2,f6.2,3f7.2) 2061 format (4f7.3,f6.2,4f7.2,f6.2,3f7.2) 2062 format (f7.2) 2063 format (4x,i1,6(1x,f8.4)) 2064 format (3x,5f7.3) 2164 format (3x,3f7.3) 2065 format (3x,'layer asdblk aszlyt sfsan asfsil asfcla ', & 'as0ph ascmg ascna asfcce asfcec asfesp') 2066 format(3x,'asfom asfnoh asfpoh asfpsp asfsmb asdagd aseags ', & 'ahrwc aheaep ahrwcw ahrwcf ahrwca ahrwcs') 2067 format('aslrr') 2068 format(3x,'layer admrtz(1) admrtz(2) admrtz(3) ', & ' admbgz(1) admbgz(2) admbgz(3)') 2069 format(4x,'adfscv adffcv acmyld aczht aczrtd') 2169 format(4x,'acmyld aczht aczrtd') 2269 format(4x,'adfscv adffcv ') 2070 format(3x,'aszrgh asxrgw asxrgs asargo asxdks asxdkh') 2071 format (4x,6(2x,f7.3)) 2072 format(3x,'asfcr asflos') 2073 format (1x,2f7.3) 2074 format(3x,'admf(1) admf(2) admf(3) admst(1)', & ' admst(2) admst(3)') 2075 format (6(2x,f7.3)) c + + + END SPECIFICATIONS + + + c c crushing process if (output.eq.11) then write(29,2040) do 100 l = 1,slay write(29,2050) aslagn(l,isr),aslagx(l,isr),aslagm(l,isr), & as0ags(l,isr) 100 continue endif c loosening process if (output.eq.12) then write(29,2041) do 110 l = 1,slay write(29,2051) asdblk(l,isr),asdsblk(l,isr),aszlyt(l,isr) 110 continue endif c mixing process if (output.eq.13) then write(29,2065) do 200 l = 1,slay write(29,2060) l,asdblk(l,isr),aszlyt(l,isr),asfsan(l,isr), & asfsil(l,isr),asfcla(l,isr),as0ph(l,isr),ascmg(l,isr), & ascna(l,isr),asfcce(l,isr),asfcec(l,isr),asfesp(l,isr) 200 continue write(29,2066) do 201 l = 1,slay write(29,2061) asfom(l,isr),asfnoh(l,isr),asfpoh(l,isr), & asfpsp(l,isr),asfsmb(l,isr),asdagd(l,isr),aseags(l,isr), & ahrwc(l,isr),aheaep(l,isr),ahrwcw(l,isr),ahrwcf(l,isr), & ahrwca(l,isr),ahrwcs(l,isr) 201 continue write(29,2068) do 202 l = 1,slay write(29,2063) & l,admrtz(l,1,isr),admrtz(l,2,isr), admrtz(l,3,isr), & admbgz(l,1,isr),admbgz(l,2,isr),admbgz(l,3,isr) 202 continue endif c random roughness calculation if (output.eq.2) then write(29,2067) write(29,2062) aslrr(isr) endif c crust breakdown if (output.eq.5) then write(29,2072) write(29,2073) asfcr(isr), asflos(isr) endif c inversion if (output.eq.14) then do 120 l = 1,slay write(29,2060) l,asdblk(l,isr),aszlyt(l,isr),asfsan(l,isr), & asfsil(l,isr),asfcla(l,isr),as0ph(l,isr),ascmg(l,isr), & ascna(l,isr),asfcce(l,isr),asfcec(l,isr),asfesp(l,isr) 120 continue write(29,2066) do 125 l = 1,slay write(29,2061) asfom(l,isr),asfnoh(l,isr),asfpoh(l,isr), & asfpsp(l,isr),asfsmb(l,isr),asdagd(l,isr),aseags(l,isr), & ahrwc(l,isr),aheaep(l,isr),ahrwcw(l,isr),ahrwcf(l,isr), & ahrwca(l,isr),ahrwcs(l,isr) 125 continue endif c remove process if (output.eq.61) then write(29,2169) write(29,2164) acmyld(isr), aczht(isr),aczrtd(isr) write(29,2269) do 126 i=1,mnbpls write(29,2073) adfscv(i,isr),adffcv(i,isr) 126 continue endif c oriented roughness calculation if (output.eq.3 .or. output.eq.4 .or. output.eq.5) then write(29,2070) write(29,2071) aszrgh(isr),asxrgw(isr),asxrgs(isr),asargo(isr), & asxdks(isr),asxdkh(isr) endif c bury/lift process if (output.eq.34) then write(29,2068) do 230 l = 1,slay write(29,2063) & l,admrtz(l,1,isr),admrtz(l,2,isr),admrtz(l,3,isr), & admbgz(l,1,isr),admbgz(l,2,isr),admbgz(l,3,isr) 230 continue write(29,2074) write(29,2075) admf(1,isr),admf(2,isr),admf(3,isr), & admst(1,isr),admst(2,isr),admst(3,isr) endif return end