Blank 07 W2 Form

Document Sample
Blank 07 W2 Form Powered By Docstoc
					c SEDIMENT MODEL - Steady State Version
c
      include 'sedsfdim.for'
c-- Input variables
     integer monzz,dayzz,yearzz
     real jcinzz,sodzz,jnh4zz,jno3zz,jpo4zz,jsizz,
    .o2zz,nh4zz,no3zz,po4zz,sizz,tempzz
c- output variables
     real nh40z,nh41z,nh42z,nh4t2z,no30z,no31z,no32z,no3t2z
c-- function for zbrent root finder
     external sedf
c-- leave the file names blank - the file names come from the command line
     open(11,file=' ',form='formatted')                                      !Parameter input
     open(13,file=' ',form='formatted')                                      !Data input
     open(15,file=' ' ,form='formatted')                                     !Output file
     read(11,1001) m1,m2,Dp,w2,Dd,thtaDp,thtaDd                              !(A)
     read(11,1001) kappnh4,pienh4,thtanh4,kmnh4,thtakmnh4,kmnh4o2
     read(11,1001) kapp1no3,k2no3,thtano3
     read(11,1001) kappd1,kappp1,pie1s,pie2s,thtapd1,kmhso2
     read(11,1001) ksi,csisat,dpie1si,pie2si
     read(11,1001) h2,thtasi,kmpsi,o2critsi
     read(11,1001) dpie1po4,pie2po4,o2crit,kmo2Dp
     read(11,1001) frpon1,kpon1,thtapon1
     read(11,1001) frpon2,kpon2,thtapon2
     read(11,1001) frpon3,kpon3,thtapon3
     read(11,1001) frpoc1,kpoc1,thtapoc1
     read(11,1001) frpoc2,kpoc2,thtapoc2
     read(11,1001) frpoc3,kpoc3,thtapoc3
     read(11,1001) frpop1,kpop1,thtapop1
     read(11,1001) frpop2,kpop2,thtapop2
     read(11,1001) frpop3,kpop3,thtapop3
     read(11,1001) ratiocn, ratiocp, ratiocsi
1001 format(//7f12.2)
     read(13,*) ninp

c-- Main loop
     do 999 i=1,ninp
     read(13,*) monzz,dayzz,yearzz,sodzz,jnh4zz,jno3zz,jpo4zz,jsizz,
    . o2zz,nh4zz,no3zz,po4zz,sizz,tempzz                                     !(B)
c-- if missing data do not process
     if((jnh4zz.eq.0.0).or.(sodzz.eq.0.0).or.
    . (o2zz.eq.0.0).or.(tempzz.eq.0.0).or.
    . (jnh4zz.eq.-999.0).or.(sodzz.eq.-999.0).or.
    . (o2zz.eq.-999.0).or.(tempzz.eq.-999.0)) then
     goto 999
     endif
c-- Compute jn                                                               !(C)
     xknh4=kappnh4*thtanh4**(tempzz-20.)
c-- making the conversion factor for jn to sod a function of temperature
     ao2n=106./16.*12./14.*32./12.*1.068**(tempzz-20.)
     o2=o2zz
c-- convert from mg/m2-d to g/m2-d
     xjnh4=jnh4zz/1000.
     xnh4=nh4zz/1000.
c-- compute intermediate variables
     z1=xjnh4**2
     z2=xknh4**2
     z3=ao2n**2
     z4=1./z3
     z5=o2**2
     z5b=(-4.*ao2n*o2*xknh4**4*xnh4**3
    . -z3*z1*z2*xnh4**2+18.*ao2n*o2*z1*z2*xnh4+27.*z5*z1*z2+4.*z3
    . *xjnh4**4)
c-- check for negative root. If this occurs, assume that the olw nh4=0
     if(z5b.le.0.0) xnh4=0.0
     z5a=(z4*o2*xknh4*sqrt(-4.*ao2n*o2*xknh4**4*xnh4**3
    . -z3*z1*z2*xnh4**2+18.*ao2n*o2*z1*z2*xnh4+27.*z5*z1*z2+4.*z3
    . *xjnh4**4)/sqrt(3.)/6.0+z4*(9.*ao2n*o2*xjnh4*z2*xnh4
    . +27.*z5*xjnh4*z2
    . +2.*z3*xjnh4**3)/54.0)
     z6=z5a/abs(z5a)*(abs(z5a))**(1.0/3.0)
     xjn = z6+(3.*o2*z2*xnh4+ao2n*z1)/(ao2n*z6)/9.0+xjnh4/3.0
     acn=106./16.*12./14.
     jcinzz=acn*xjn*1000.                                                !(D)
c-- do not process missing data record
     if(jcinzz.le.0.0) goto 999
c- Set maximum no of iterations and error                                !(E)
     isteady=20
     esteady=0.001

c-- set min and max and accuracy for root finder
    sodmin = 1.0e-6
    sodmax = 1.0e6
    eps = 1.0e-5

c-- set overlying water concentrations
    nh40=nh4zz
    no30=no3zz
    o20=o2zz
    if(o20.eq.0.0) o20=1.0e-2
    po40=po4zz
    si0=sizz
    tempd=tempzz
    jcin=jcinzz/(frpoc1+frpoc2)                                          !(F)

c-- stoichiometric ratios for N,P and Si
    jnin=jcinzz/(frpon1+frpon2)/ratiocn
    jpin=jcinzz/(frpop1+frpop2)/ratiocp
    jsiin=jcin/ratiocsi

c-- evaluate the temperature dependent coefficients
    temp20=tempd-20.0
    temp202=temp20/2.0
    xappnh4=kappnh4*thtanh4**temp202
   xappd1=kappd1*thtapd1**temp202
   xappp1=kappp1*thtapd1**temp202
   xapp1no3=kapp1no3*thtano3**temp202
   xk2no3=k2no3*thtano3**temp20*h2
   xksi=ksi*thtasi**temp20*h2
   kl12=Dd/(h2/2.0)*thtaDd**temp20
   w12=Dp/(h2/2.0)*(thtaDp**temp20)*poc1/1.0e5
   xkpon1=kpon1*thtapon1**temp20*h2
   xkpon2=kpon2*thtapon2**temp20*h2
   xkpon3=kpon3*thtapon3**temp20*h2
   xkpoc1=kpoc1*thtapoc1**temp20*h2
   xkpoc2=kpoc2*thtapoc2**temp20*h2
   xkpoc3=kpoc3*thtapoc3**temp20*h2
   xkpop1=kpop1*thtapop1**temp20*h2
   xkpop2=kpop2*thtapop2**temp20*h2
   xkpop3=kpop3*thtapop3**temp20*h2
   xksi=ksi*thtasi**temp20*h2

c- set t minus 1 concentrations to anything       !(G)
    nh4t2tm1=0.0
    sit2tm1=0.0
    psitm1=0.0

c-- Loop to converge                              !(H)
    do 1113 istst=1,isteady
    iststsav=istst
c-- compute diagenesis                            !(I)
    pon1=(frpon1*jnin)/(xkpon1+w2)
    pon2=(frpon2*jnin)/(xkpon2+w2)
    pon3=(frpon3*jnin)/(xkpon3+w2)
    xjn=xkpon1*pon1+xkpon2*pon2+xkpon3*pon3
    jn=xkpon1*pon1+xkpon2*pon2+xkpon3*pon3
    poc1=(frpoc1*jcin)/(xkpoc1+w2)
    poc2=(frpoc2*jcin)/(xkpoc2+w2)
    poc3=(frpoc3*jcin)/(xkpoc3+w2)
    xjc=xkpoc1*poc1+xkpoc2*poc2+xkpoc3*poc3
    jc=xkpoc1*poc1+xkpoc2*poc2+xkpoc3*poc3
    pop1=(frpop1*jpin)/(xkpop1+w2)
    pop2=(frpop2*jpin)/(xkpop2+w2)
    pop3=(frpop3*jpin)/(xkpop3+w2)
    xjp=xkpop1*pop1+xkpop2*pop2+xkpop3*pop3
    jp=xkpop1*pop1+xkpop2*pop2+xkpop3*pop3
    fd2=1.0/(1.0+m2*pie2si)
    k3=xksi*(csisat-fd2*sit2tm1)/(psitm1+kmpsi)
    psi=jsiin/(k3+w2)

c-- solve the nh4,no3,sod, po4 and si equations
    sod=zbrent(sedf,sodmin,sodmax,eps)            !(J)

c-- replace the t minus 1 concentrations          !(K)
     nh41tm1=nh41
     psitm1=psi
     sit2tm1=sit2
c-- compare to the previous iteration                                       !(L)
     iagain=0
c- do at least one iteration
     if(istst.eq.1) iagain=1
     if(abs(nh4t2sav-nh4t2).gt.abs(esteady*nh4t2sav)) iagain=1
     if(abs(no3t2sav-no3t2).gt.abs(esteady*no3t2sav)) iagain=1
     if(abs(hst2sav-hst2).gt.abs(esteady*hst2sav)) iagain=1
     if(abs(sit2sav-sit2).gt.abs(esteady*sit2sav)) iagain=1
     if(abs(po4t2sav-po4t2).gt.abs(esteady*po4t2sav)) iagain=1
     if(abs(pon1sav-pon1).gt.abs(esteady*pon1sav)) iagain=1
     if(abs(pon2sav-pon2).gt.abs(esteady*pon2sav)) iagain=1
     if(abs(pon3sav-pon3).gt.abs(esteady*pon3sav)) iagain=1
     if(abs(poc1sav-poc1).gt.abs(esteady*poc1sav)) iagain=1
     if(abs(poc2sav-poc2).gt.abs(esteady*poc2sav)) iagain=1
     if(abs(poc3sav-poc3).gt.abs(esteady*poc3sav)) iagain=1
     if(abs(pop1sav-pop1).gt.abs(esteady*pop1sav)) iagain=1
     if(abs(pop2sav-pop2).gt.abs(esteady*pop2sav)) iagain=1
     if(abs(pop3sav-pop3).gt.abs(esteady*pop3sav)) iagain=1
     if(abs(psisav-psi).gt.abs(esteady*psisav)) iagain=1
c- iagain=0 => converged; otherwise save the ending conditions and repeat
     if(iagain.eq.0) go to 1114
     if(istst.ne.isteady) then
     nh4t2sav=nh4t2                                                         !(M)
     no3t2sav=no3t2
     hst2sav=hst2
     sit2sav=sit2
     po4t2sav=po4t2
     pon1sav=pon1
     pon2sav=pon2
     pon3sav=pon3
     poc1sav=poc1
     poc2sav=poc2
     poc3sav=poc3
     pop1sav=pop1
     pop2sav=pop2
     pop3sav=pop3
     psisav=psi
     end if
1113 continue
     print *,'Did not equilibrate in',isteady,' iterations'
1114 continue
c- output solid phase concentrations in mg/g (O2) mg/kg (po4) mg/g (Si)
     hst2r=hst2/1000./m2
     po4t2r=po4t2/1000./m2
     sit2r=(psi+sit2)/1.0e6/m2
     psir=(psi)/1.0e6/m2
c- output POC,PON,etc in mg/g
     pon1r=pon1/1.0e6/m2
     pon2r=pon2/1.0e6/m2
     pon3r=pon3/1.0e6/m2
     ponr=pon1r+pon2r+pon3r
    poc1r=poc1/1.0e6/m2
    poc2r=poc2/1.0e6/m2
    poc3r=poc3/1.0e6/m2
    pocr=poc1r+poc2r+poc3r
    pop1r=pop1/1.0e6/m2
    pop2r=pop2/1.0e6/m2
    pop3r=pop3/1.0e6/m2
    popr=pop1r+pop2r+pop3r
c- output dissolved concentrations in mg/L
    nh40z=nh40/1000.
    nh41z=nh41/1000.
    nh42z=nh42/1000.
    nh4t2z=nh4t2/1000.
    no30z=no30/1000.
    no31z=no31/1000.
    no32z=no32/1000.
    no3t2z=no3t2/1000.
    si0z=si0/1000.
    si1z=si1/1000.
    si2z=si2/1000.
    sit2z=sit2/1000.
    po40z=po40/1000.
    po41z=po41/1000.
    po42z=po42/1000.
    po4t2z=po4t2/1000.
    hs1z=hs1
    hs2z=hs2
    hst2z=hst2
     if(i.eq.1)
   . write(15,1003)'mon,','day,','year,','sodzz,','jnh4zz,','jno3zz,',
   . 'jpo4zz,','jsizz,', 'o2zz,','nh4zz,','no3zz,','po4zz,','sizz,',
   . 'tempzz,','jnin,','jn,','jnh4,','nh40z,','nh41z,','nh42z,',
   . 'pon1r,','pon2r,','pon3r,','jno3,','no30z,','no31z,','no32z,',
   . 'jpin,','jp,','jpo4,','po40z,','po41z,','po42z,','pop1r,',
   . 'pop2r,','pop3r,','jsiin,','jsi,','si0z,','si1z,','si2z,',
   . 'sit2r,','jcin,','jc,','jhs,','hs1z,','hs2z,','hst2r,','poc1r,',
   . 'poc2r,','poc3r,','o20,','sod,','s,','h1,','csod'
1003 format(60a11)
    h1=kl12*h2/s
    xmon=monzz
    xday=dayzz
    xyear=yearzz
    write(15,103) xmon,xday,xyear,sodzz,jnh4zz,jno3zz,jpo4zz,
   . jsizz, o2zz,nh4zz,no3zz,po4zz,sizz,tempzz,jnin,jn,jnh4,nh40z,
   . nh41z,nh42z,pon1r,pon2r,pon3r,jno3,no30z,no31z,no32z,jpin,jp,
   . jpo4,po40z,po41z,po42z,pop1r,pop2r,pop3r,jsiin,jsi,si0z,si1z,
   . si2z,sit2r,jcin,jc,jhs,hs1z,hs2z,hst2r,poc1r,poc2r,poc3r,
   . o20,sod,s,h1,csod
103 format(60(e10.3,','))

999 continue
   end
    function sedf(sod1)                                            !(N)
    include 'sedsfdim.for'
c- Ammonia
    s = sod1/o20                                                   !(O)
    k0h1p=0.
    k1h1p=0.
    k2h2d=0.
    k2h2p=0.
    if(kmnh4.ne.0.0) then
      k0h1d=xappnh4**2/s*kmnh4*(o20/(kmnh4o2+o20))
      k1h1d=s
    else
      k1h1d=xappnh4**2/s*(o20/(kmnh4o2+o20))+s
      k0h1d=0.
    endif
    j1=s*nh40
    k3=0.0
    j2=jn
    pie1=pienh4
    pie2=pienh4
    kmc1=kmnh4
    call sedsffnl(nh41,nh42,nh4t1,nh4t2,nh41tm1)
    jnh4=s*(nh41-nh40)
c- Oxygen consumed by nitrification
c a1 = 64/14 * 1/1000 - mole ratio and mg/m2-day to gm/m2-day
    a1 = 0.0045714
    if(kmnh4.ne.0.0) then
      jo2nh4=a1*k0h1d*nh41/(kmnh4+nh41tm1)
    else
      jo2nh4=a1*(k1h1d-s)*nh41
    endif

c- Nitrate
    k0h1d=0.
    k0h1p=0.
    kmc1=0.0
    k1h1d=xapp1no3**2/s+s
    k1h1p=0.
    k2h2d=xk2no3
    k2h2p=0.
    if(kmnh4.ne.0.0) then
      j1=s*no30+
   . xappnh4**2/s*kmnh4*(o20/(kmnh4o2+o20))*nh41/(kmnh4+nh41tm1)
    else
      j1=s*no30+xappnh4**2/s*(o20/(kmnh4o2+o20))*nh41
    endif
    k3=0.0
    j2=0.0
    pie1=0.
    pie2=0.
    call sedsffnl(no31,no32,no3t1,no3t2,0.0)
   jno3=s*(no31-no30)

c- Sulfide
c diagenesis consumed by denitrification
c a2 = 10/8 * 32/14 * 1/1000
    a2 = 0.00285714
    xjcno3=a2*xapp1no3**2/s*no31 + a2*xk2no3*no32
c- Convert carbon diagenesis flux to o2* units and decrement jc
c a0= 32./12./1000.
    a0=2.666666666e-3
    xjco2 = a0*jc
    xjc1 = amax1(xjco2-xjcno3,0.0)
203 continue
c-- Sulfide in O2 equivalents
    k0h1d=0.
    k0h1p=0.
    kmc1=0.0
    k1h1d=xappd1**2/s*(o20/kmhso2) + s
    k1h1p=xappp1**2/s*(o20/kmhso2)
    k2h2d=0.
    k2h2p=0.
    j1=0.
    k3=0.0
    j2=xjc1
    pie1=pie1s
    pie2=pie2s
    call sedsffnl(hs1,hs2,hst1,hst2,0.0)
    jhs=s*hs1
    csod=(xappd1**2/s*fd1 + xappp1**2/s*fp1)*(o20/kmhso2)*hst1

c- Silica
    k0h1d=0.
    k0h1p=0.
    kmc1=0.0
    k1h1d=s
    k1h1p=0.
    k2h2d=0.
    k2h2p=0.
    j1=s*si0
c-- Oxygen dependency of pie1
    if(o20.lt.o2crit) then
       pie1=pie2si*dpie1si**(o20/o2crit)
    else
       pie1=pie2si*dpie1si
    endif
    pie2=pie2si
c- Silica dissolution kinetics
    fd2=1./(1.+m2*pie2)
    k3=xksi*psi/(psi+kmpsi)*fd2
    j2=xksi*psi/(psi+kmpsi)*csisat
    call sedsffnl(si1,si2,sit1,sit2,0.0)
    jsi=s*(si1-si0)
c-- Phosphate
    k0h1d=0.
    k0h1p=0.
    kmc1=0.0
    k1h1d=s
    k1h1p=0.
    k2h2d=0.
    k2h2p=0.
    j1=s*po40
    k3=0.0
    j2=jp
c-- Oxygen dependency of pie1
    if(o20.lt.o2crit) then
       pie1=pie2po4*dpie1po4**(o20/o2crit)
    else
       pie1=pie2po4*dpie1po4
    endif
    pie2=pie2po4
    call sedsffnl(po41,po42,po4t1,po4t2,0.0)
    jpo4=s*(po41-po40)

c- SOD - evaluate the error
    sod=csod + jo2nh4
    sedf = sod - sod1                              !(P)
    return
    end

     subroutine sedsffnl(c1,c2,ct1,ct2,c1tm1)      !(Q)
     include 'sedsfdim.for'
c-   Notation
c     c1 = layer 1 dissolved conc.
c     c2 = layer 2 dissolved conc.
c     ct1 = layer 1 total conc.
c     ct2 = layer 2 total conc.
c     c1tm1 = c1 at time level t - 1

    fd1=1./(1.+m1*pie1)
    fp1=m1*pie1/(1.+m1*pie1)
    fd2=1./(1.+m2*pie2)
    fp2=m2*pie2/(1.+m2*pie2)
c- Transport and Decay terms
    f12 = w12*fp1 + kl12*fd1
    f21 = w12*fp2 + kl12*fd2
c evaluate the MM term at time level t-1           !(R)
    if(kmc1.ne.0.0) then
      xk0 = (k0h1d*fd1 + k0h1p*fp1)/(kmc1+c1tm1)
    else
      xk0=0.0
    endif
    xk1 = xk0 + k1h1d*fd1 + k1h1p*fp1
    xk2 = k2h2d*fd2 + k2h2p*fp2
c- Matrix and forcing function
    a11 = -f12 -xk1 - w2
    a21 = f12 + w2
    a12 = f21
    b1 = -j1
    a22=-f21-xk2-w2-k3
    b2=-j2
c-- Solve the 2x2 set of linear equations, stop if determinant=0
    det=a11*a22-a12*a21
    if(det.eq.0.0) stop                                                  !(S)
    ct1=(b1*a22-b2*a12)/det
    ct2=(b2*a11-b1*a21)/det
    c1=fd1*ct1
    c2=fd2*ct2
    return
    end

sedsfdim.for                                                             !(T)
------------
c- Input parameters
     real
    . m1,m2,kappnh4,kmnh4,kmnh4o2,kapp1no3,k2no3,kappd1,kappp1,
    . kmhso2,ksi,kmpsi,kmo2Dp,kmo2Dd,kpon1,kpon2,kpon3,kpoc1,
    . kpoc2,kpoc3,kpop1,kpop2,kpop3
     common /in/
    . m1,m2,Dp,w2,Dd,thtaDp,thtaDd,kappnh4,pienh4,thtanh4,kmnh4,
    . kmnh4o2,kapp1no3,k2no3,thtano3,kappd1,kappp1,pie1s,pie2s,
    . thtapd1,kmhso2,ksi,csisat,dpie1si,pie2si,h2,thtasi,kmpsi,
    . dpie1po4,pie2po4,o2crit,kmo2Dp,kmo2Dd,frpon1,kpon1,thtapon1,
    . frpon2,kpon2,thtapon2,frpon3,kpon3,thtapon3,frpoc1,kpoc1,
    . thtapoc1,frpoc2,kpoc2,thtapoc2,frpoc3,kpoc3,thtapoc3,frpop1,
    . kpop1,thtapop1,frpop2,kpop2,thtapop2,frpop3,kpop3,thtapop3,
    . ratiocn, ratiocp, ratiocsi
c- concentrations
     real
    . nh40,nh41,nh42,nh4t1,nh4t2,nh41tm1,jnh4,nh4t2sav,no30,no31,no32,
    . no3t1,no3t2,jno3,no3t2sav,jhs,jsi,jpo4
     common /conc/
    . nh40,nh41,nh42,nh4t1,nh4t2,nh41tm1,jnh4,nh4t2sav,no30,no31,no32,
    . no3t1,no3t2,jno3,no3t2sav,hs0,hs1,hs2,hst1,hst2,jhs,hst2sav,
    . si0,si1,si2,sit1,sit2,jsi,sit2sav,po40,po41,po42,po4t1,po4t2,
    . jpo4,po4t2sav
c- diagenesis
     real
    . jnin,xjn,jn,jcin,xjc,jc,jpin,xjp,jp,jsiin,jo2nh4
     common /diag/
    . jnin,xjn,jn,jcin,xjc,jc,jpin,xjp,jp,psi,psitm1,psisav,psiav,
    . jsiin,xjcno3,xjco2,xjc1,jo2nh4,xjco2av,xjc1av
c- sedsffnl storage
     real
    . k0h1d,k0h1p,k1h1d,k1h1p,k2h2d,k2h2p,k3,
    . pie1,pie2,j1,j2,kmc1,w12,kl12,tempd,o20
     common /pars/
    . k0h1d,k0h1p,k1h1d,k1h1p,k2h2d,k2h2p,k3,
    . pie1,pie2,j1,j2,kmc1,w12,kl12,tempd,o20
c- sedf storage
     real kl12nom
     common /store/
    . xappnh4,xappd1,xappp1,xapp1no3,xk2no3,xksi,temp20,temp202,
    . fd1,fp1,fd2,fp2,sod,csod,s,kl12nom
  Input file
   ----------
 Test Run
    m1            m2               Dp            w2             Dd
          0.5             0.5   6.00E-05       6.85E-06       5.00E-04
    kg/L          kg/L            m2/d           m/d           m2/d
 kappnh4         pienh4          thtanh4       kmnh4        thtakmnh4
  1.31E-01                  1        1.123          728           1.125
    m/d            L/kg               -          ug/L             -
kapp1no3         k2no3           thtano3
          0.1          0.25            1.08
    m/d             m/d               -
  kappd1        kappp1            pie1s          pie2s       thtapd1
          0.2             0.4           100          100           1.08
    m/d            m/d             L/kg          L/kg            -
     ksi          csisat          dpie1si        pie2si
          0.5       40000                 10         100
      /d         ug Si/L           L/kg          L/kg
     h2            thtasi        kmpsi           o2critsi
          0.1             1.1   5.00E+07                1
     m                 -          ug Si/L      mg O2/L
 dpie1po4       pie2po4            o2crit      kmo2Dp
         300             100               2            4
    L/kg           L/kg         mg O2/L         mg O2/L
   frpon1        kpon1          thtapon1
        0.65         0.035               1.1
  gN/gN               /d             -
   frpon2        kpon2          thtapon2
        0.25       0.0018              1.15
  gN/gN               /d             -
   frpon3        kpon3          thtapon3
          0.1               0          1.17
  gN/gN               /d             -
   frpoc1        kpoc1          thtapoc1
        0.65         0.035               1.1
  gC/gC               /d             -
   frpoc2        kpoc2          thtapoc2
          0.2      0.0018              1.15
  gC/gC               /d             -
   frpoc3        kpoc3          thtapoc3
        0.15                0          1.17
  gC/gC               /d             -
   frpop1        kpop1          thtapop1
        0.65         0.035               1.1
   gP/gP              /d             -
   frpop2        kpop2          thtapop2
          0.2      0.0018              1.15
   gP/gP              /d             -
   frpop3        kpop3          thtapop3
        0.15                0          1.17
   gP/gP              /d             -
   ratiocn        ratiocp         ratiocsi
    5.68        41       2.5
gC/gN      gC/gP     gC/gSi
              !(U)


 thtaDp     thtaDd
     1.15      1.08
     -            -
kmnh4o2
     0.74
 mg O2/L



kmhso2
      4
mg O2/L
Output file - One line, transposed                   !(V)
-----------
 sodzz         g/m2-d          1.17     jpo4   mg/m2-d       9.64
jnh4zz         mg/m2-d         62.8    po40z   mg/L          0.01
jno3zz         mg/m2-d        -14.1    po41z   mg/L         0.056
jpo4zz         mg/m2-d             0   po42z   mg/L          4.32
 jsizz         mg/m2-d          325    pop1r   mg/g         0.007
  o2zz         mg/L            5.57    pop2r   mg/g          0.05
 nh4zz         mg/m3            125    pop3r   mg/g         0.585
 no3zz         mg/m3            103    jsiin   mg/m2-d        219
 po4zz         mg/m3             9.6     jsi   mg/m2-d       94.2
  sizz         mg/m3            385     si0z   mg/L         0.385
tempzz         oC              15.8     si1z   mg/L         0.841
  jnin         mg/m2-d            91    si2z   mg/L          27.7
    jn         mg/m2-d         80.3    sit2r   mg/L          36.4
  jnh4         mg/m2-d         67.4     jcin   mg/m2-d        547
 nh40z         mg/L           0.125       jc   mg/m2-d        457
 nh41z         mg/L           0.451      jhs   mg/m2-d      0.006
 nh42z         mg/L            22.6     hs1z   mg/L         0.028
 pon1r         mg/g            0.05     hs2z   mg/L           303
 pon2r         mg/g           0.426    hst2r   mg/g            31
 pon3r         mg/g            2.66    poc1r   mg/g         0.303
  jno3         mg/m2-d         7.33    poc2r   mg/g          2.05
 no30z         mg/L           0.103    poc3r   mg/g            24
 no31z         mg/L           0.139      o20   mg/L          5.57
 no32z         mg/L           0.023      sod   g/m2-d        1.15
  jpin         mg/m2-d         13.4        s   m/d          0.206
    jp         mg/m2-d         11.2       h1   m            0.002
                                        csod   g/m2-d        1.09

				
DOCUMENT INFO
Shared By:
Categories:
Tags:
Stats:
views:26
posted:1/7/2011
language:English
pages:15
Description: Blank 07 W2 Form document sample