,

 

1.     .

 

2.     (, , , )

 

3.    

 

4.    

 

5.     ,

 

6.    

 

7.    

 

8.    

 

9.     ʻ

 

 

 

 

1.     .

 

 

 

 

2.     (, , , )

 

 

 

 

 

3.    

 

 

 

4.    

 

' Programma rascheta zagrjaznenija

' Avtor prof.dtn. A.I.GOLOVANOV

' M G M I

' 2 zagrjaznitelja:Cd i Pb

' w=wm+(m-wm)*exp(-(psi/mju*hk)); 4-tochechanja shema po vode i soljam

' Dk=w*(Dm+L*v)

' d(wC/dt)=d(Dk*dC/dx-vC)/dx-dN/dt-q

' dN/dt=-bet(alf*S-w*C)

' Mnogoletie,Smax - neogranich.

' Konechn rezultati

' Zagr. vesnoj; Zagr s osadkami; Variazija pogodi

' Otbor Cd i Pb

' Oroshenie presnoj vodoj

' Razmernosti: Konz[gr/m^3 vodi]; Soderganje[gr/v^3 pochvi]

' F A R!

 

Dim p(40), mju(40), wm(40), kf(40), hk(40), h(40), w0(40), x(40), w1(40), w2(40)

Dim w3(40), h1(40), h2(40), r(40), pe(40), q(40), er(40), s(20), s1(20)

Dim oc(100), toc(100), ef(40), ek(40), pol(20), tpol(20), tsu(20), teu(20)

Dim tets(40), hgs(40), hgm(40), gs(40), kw(40), al(40), plt(40), vs(40)

Dim ca0(40), ca1(40), cn0(40), cn1(40), na0(40), nn0(40), na1(40), nn1(40)

Dim sora0(40), sora1(40), pn(40), qn(40), qam(40), qni(40), za(20), zn(20)

Dim za1(20), zn1(20), kk1(40), kk2(40), smax(40), zmn(20, 50), szmn(20)

Dim tetopt(40), ga(40), gum1(40), gum2(40), kk3(40), kisp(50), kosad(50)

Dim sorn0(40), sorn1(40), smax1(40), kfar(50)

DefInt I

'****************************************************************************

k1$ = "df1_05w.bas" 'imja faila rezultatov po vode

k1s$ = "df1_05s.bas" 'imja faila rezultatov po azotu

k2s$ = "df1_05e.bas"

k3s$ = "df1_05j.bas"

 

k2$ = " UFA, 50%(1969), Mn.travi" 'zagolovok

 

' Uslovija

l1 = 0 '-1: hg=peremen, 0: hg=const, +1: hg=gluboko

l2 = -1 '-1: 1-var progr, 0: 2-var progr, +1: 3-var progr

l3 = -1 '-1: net vivoda cherez dtv, 0: vivod w2(I), +1: vivod cherez dtv

l4 = -1 '-1: w0(I)= ravnovesnaja, s vesennimi zapasami

'+1: w0(I)= ljubaja

l6 = 1 '-1: bez prognoza Oc; +1: prognoz Oc

 

n = 40 'k-vo sloev

h(1) = 0: h(2) = 0.02: h(3) = 0.03: h(4) = 0.05: h(5) = 0.05: h(6) = 0.05: h(7) = 0.05: h(8) = 0.05

h(9) = 0.05: h(10) = 0.05: h(11) = 0.05: h(12) = 0.05: h(13) = 0.05: h(14) = 0.05: h(15) = 0.05

h(16) = 0.05: h(17) = 0.05: h(18) = 0.05: h(19) = 0.1: h(20) = 0.1: h(21) = 0.1: h(22) = 0.1

h(23) = 0.2: h(24) = 0.2: h(25) = 0.2: h(26) = 0.2: h(27) = 0.2: h(28) = 0.25: h(29) = 0.25

h(30) = 0.25: h(31) = 0.25: h(32) = 0.25: h(33) = 0.25: h(34) = 0.25: h(35) = 0.25: h(36) = 0.25

h(37) = 0.25: h(38) = 0.25: h(39) = 0.25: h(40) = 0

horA = 0.5: horB1 = 0.3: horB2 = 0.2: horC = 5 'tolschina genet.gorizontov,m

 

pA = 0.52: wmA = 0.11: mjuA = 1: hkA = 2#: kfA = 0.4

pB1 = 0.48: wmB1 = 0.09: mjuB1 = 1: hkB1 = 1.8: kfB1 = 0.2

pB2 = 0.47: wmB2 = 0.08: mjuB2 = 1: hkB2 = 1.7: kfB2 = 0.15

pC = 0.5: wmC = 0.12: mjuC = 1: hkC = 1.7: kfC = 0.1

 

naA = 0.5: nnA = 20: pltA = 1.2

naB1 = 0.5: nnB1 = 20: pltB1 = 1.4

naB2 = 0.3: nnB2 = 10: pltB2 = 1.3

naC = 0.05: nnC = 0.5: pltC = 1.3

 

'naA=Cd,mg/kg.pochvi, nachaln.

'nnA=Pb,mg/kg.pochvi, nachaln.

'pltA=pflotnost pochvi,gr/cm^3

ppv = 0.6 'v dolljah ot poristosti

wz = 1.4 * wmA 'vlagn.zavjadanija

reg = 0.85: mpolmin = 40

hw0 = 0.4: hwm = 0.4: hwk = 0.4: he0 = 0.5: hem = 0.5

hg0 = 2.5: kf0 = 0: t0 = 3: hp = 0.2: fk = 0.0001: d0 = 0.002: hdr = 3

Uvlves = 50

tgod = 365: tveg = 170: dte = 0.75: dtu = 0.01: dtef = 10: dtv = 10

ngod = 10: nbezN = 0

L = 0.2: Dm = 0.0001

Alfa = 0.2: Beta = 0.2

Alfa1 = 0.1: Beta1 = 0.2

catza = 1 'zagr Cd s osadkami,mg/litr

catzn = 30 'zagr Pb s osadkami,mg/litr

urmax = 12000 'max ur.trav,kg/ha

capv = 0: cnpv = 0 'Cd i Pb v polivn.vode,mg/litr

horsol = 0.5: horg = 1#

al(1) = 0.04: al(2) = 0.05: al(3) = 0.06: al(4) = 0.08: al(5) = 0.09: al(6) = 0.1

al(7) = 0.09: al(8) = 0.03: al(9) = 0.06: al(10) = 0.08: al(11) = 0.08: al(12) = 0.06

al(13) = 0.05: al(14) = 0.04: al(15) = 0.03: al(16) = 0.03: al(17) = 0.03

ga(1) = 6.6: ga(2) = 6.6: ga(3) = 6.6: ga(4) = 6.6: ga(5) = 6.6: ga(6) = 6.6

ga(7) = 6.6: ga(8) = 6.6: ga(9) = 6.6: ga(10) = 6.6: ga(11) = 6.6: ga(12) = 6.6

ga(13) = 6.6: ga(14) = 6.6: ga(15) = 6.6: ga(16) = 6.6: ga(17) = 6.6

tetopt(1) = 0.26: tetopt(2) = 0.26: tetopt(3) = 0.26: tetopt(4) = 0.26

tetopt(5) = 0.26: tetopt(6) = 0.26: tetopt(7) = 0.26: tetopt(8) = 0.26

tetopt(9) = 0.26: tetopt(10) = 0.26: tetopt(11) = 0.26: tetopt(12) = 0.26

tetopt(13) = 0.26: tetopt(14) = 0.26: tetopt(15) = 0.26: tetopt(16) = 0.26

tetopt(17) = 0.26

noc = 43

oc(1) = 6: oc(2) = 3: oc(3) = 2: oc(4) = 2: oc(5) = 5: oc(6) = 23: oc(7) = 6: oc(8) = 4

oc(9) = 2: oc(10) = 3: oc(11) = 3: oc(12) = 19: oc(13) = 2: oc(14) = 33: oc(15) = 5

oc(16) = 7: oc(17) = 26: oc(18) = 17: oc(19) = 2: oc(20) = 3: oc(21) = 2: oc(22) = 7: oc(23) = 2

oc(24) = 4: oc(25) = 16: oc(26) = 7: oc(27) = 27: oc(28) = 4: oc(29) = 4: oc(30) = 4

oc(31) = 4: oc(32) = 5: oc(33) = 2: oc(34) = 2: oc(35) = 10: oc(36) = 4: oc(37) = 4

oc(38) = 10: oc(39) = 4: oc(40) = 6: oc(41) = 7: oc(42) = 10: oc(43) = 4

 

toc(1) = 6: toc(2) = 7: toc(3) = 8: toc(4) = 11: toc(5) = 16: toc(6) = 26: toc(7) = 27

toc(8) = 29: toc(9) = 32: toc(10) = 44: toc(11) = 53: toc(12) = 54: toc(13) = 59:

toc(14) = 64: toc(15) = 65: toc(16) = 66: toc(17) = 67: toc(18) = 75: toc(19) = 81

toc(20) = 85: toc(21) = 92: toc(22) = 95: toc(23) = 96: toc(24) = 97: toc(25) = 98

toc(26) = 102: toc(27) = 110: toc(28) = 112: toc(29) = 113: toc(30) = 118: toc(31) = 119

toc(32) = 120: toc(33) = 121: toc(34) = 125: toc(35) = 126: toc(36) = 135: toc(37) = 136

toc(38) = 139: toc(39) = 144: toc(40) = 149: toc(41) = 156: toc(42) = 159: toc(43) = 164

 

'ek=Esum, mm/day !

ek(1) = 2.3: ek(2) = 2.7: ek(3) = 3.4: ek(4) = 4.9: ek(5) = 4.6: ek(6) = 3.9: ek(7) = 4.6

ek(8) = 4.6: ek(9) = 3.1: ek(10) = 4.2: ek(11) = 3.5: ek(12) = 4.1: ek(13) = 3.1

ek(14) = 2.8: ek(15) = 2.6: ek(16) = 1.3: ek(17) = 1.2

 

'Efis, doli ot Esum !

ef(1) = 0.2: ef(2) = 0.2: ef(3) = 0.2: ef(4) = 0.15: ef(5) = 0.1: ef(6) = 0.1: ef(7) = 0.1

ef(8) = 0.1: ef(9) = 0.4: ef(10) = 0.3: ef(11) = 0.2: ef(12) = 0.1: ef(13) = 0.1: ef(14) = 0.1

ef(15) = 0.1: ef(16) = 0.1: ef(17) = 0.1: ef(18) = 0.1: ef(19) = 0.1: ef(20) = 0.1: ef(21) = 0.1

ef(22) = 0.1: ef(23) = 0.1

'Moduln.koef dlja osadkov kosad, m=1.027;sigma=.194

kosad(1) = 0.92: kosad(2) = 0.84: kosad(3) = 1.1: kosad(4) = 1.35: kosad(5) = 1.3

kosad(6) = 0.75: kosad(7) = 1.04: kosad(8) = 1.16: kosad(9) = 1.2: kosad(10) = 0.92

kosad(11) = 1.06: kosad(12) = 0.65: kosad(13) = 1.06: kosad(14) = 0.95: kosad(15) = 0.96

kosad(16) = 0.93: kosad(17) = 1.14: kosad(18) = 0.76: kosad(19) = 1.09: kosad(20) = 1.12

kosad(21) = 1.15: kosad(22) = 1.2: kosad(23) = 0.9: kosad(24) = 0.54: kosad(25) = 1.34

kosad(26) = 1.32: kosad(27) = 0.91: kosad(28) = 1.02: kosad(29) = 1.08: kosad(30) = 1.05

'Moduln.koef dlja isparenija kisp, m=1.031;sigma=.129

kisp(1) = 0.89: kisp(2) = 0.92: kisp(3) = 1.16: kisp(4) = 1.28: kisp(5) = 0.96

kisp(6) = 0.92: kisp(7) = 1.06: kisp(8) = 1.1: kisp(9) = 1.1: kisp(10) = 0.84

kisp(11) = 1.04: kisp(12) = 0.96: kisp(13) = 1.01: kisp(14) = 0.79: kisp(15) = 1.27

kisp(16) = 1.06: kisp(17) = 0.96: kisp(18) = 1.05: kisp(19) = 1.01: kisp(20) = 1.26

kisp(21) = 1.02: kisp(22) = 1.05: kisp(23) = 1.14: kisp(24) = 0.86: kisp(25) = 0.91

kisp(26) = 0.91: kisp(27) = 1.29: kisp(28) = 0.99: kisp(29) = 1.06: kisp(30) = 1.09

'Koef FAR

kfar(1) = 0.87: kfar(2) = 1.07: kfar(3) = 1.06: kfar(4) = 1.15: kfar(5) = 0.88

kfar(6) = 1.02: kfar(7) = 1.03: kfar(8) = 0.95: kfar(9) = 1.02: kfar(10) = 1.12

kfar(11) = 0.94: kfar(12) = 1.07: kfar(13) = 1.15: kfar(14) = 1.03: kfar(15) = 0.89

kfar(16) = 1.02: kfar(17) = 0.99: kfar(18) = 0.93: kfar(19) = 1.26: kfar(20) = 1.12

kfar(21) = 1.09: kfar(22) = 0.97: kfar(23) = 0.95: kfar(24) = 0.9: kfar(25) = 1.1

kfar(26) = 1.04: kfar(27) = 1.06: kfar(28) = 0.97: kfar(29) = 1.05: kfar(30) = 0.78

 

npol = 11

pol(1) = 40: pol(2) = 40: pol(3) = 40: pol(4) = 40: pol(5) = 40: pol(6) = 40

pol(7) = 40: pol(8) = 40: pol(9) = 40: pol(10) = 40: pol(11) = 40

tpol(1) = 5555: tpol(2) = 25: tpol(3) = 45: tpol(4) = 65: tpol(5) = 85: tpol(6) = 105

tpol(7) = 125: tpol(8) = 145: tpol(9) = 165: tpol(10) = 185: tpol(11) = 205

'***************************************************************************

 

'for i=1 to ngod:kisp(i)=1:kosad(i)=1:kfar(i)=1:next

 

360 k3$ = ", Bez solej": k4$ = ", Cd i Pb ": k5$ = ", Soli s gipsom"

open "o",#1,k1$

open "o",#2,k1s$

open "o",#3,k2s$

open "o",#4,k3s$

JJ = 1

nev = Int(tveg / dtef): nd = Int(tgod / 10)

tetmax = (ppv * pA - wz) / (pA - wz)

l5 = -1 '-1: Cw=analit, +1: Cw=dif

For i = 1 To noc: soc = soc + oc(i): Next

For i = 1 To nd: sec = sec + ek(i) * dtef: Next

z$ = k4$

Print #1,: Print #1,: Print #1, k2$ + z$: Print #1,

Print #2,: Print #2,: Print #2, k2$ + z$: Print #1,

Print #1, " Ishodnaja informazija"

Print #1, " K-vo sloev: N=";: Print #1, using; "##"; n

If l1 = -1 Then Print #1, " hg peremen.,drenag"

If l1 = 0 Then Print #1, " hg const"

If l1 = 1 Then Print #1, " hg gluboko"

If l2 = -1 Then Print #1, " Schet regima polivov"

If l2 = 0 Then Print #1, " Ozenka regimov poliva"

If l2 = 1 Then Print #1, " Kapitalnaja promivka"

If l4 = -1 Then Print #1, " W0(I)-ravnovesnaja" Else Print #1, " W0(I) - ljubaja"

Print #1, " Vlagn.zav.=";: Print #1, using; "##.###"; wz

Print #1, " Tetmin/tetopt=";: Print #1, using; "##.##"; reg;

'print#1," Kfar=";:print#1,using "##.##";kfar

Print #1, " Hw0=";: Print #1, using; "##.##"; hw0;

Print #1, " Hwmax=";: Print #1, using; "##.##"; hwm;

Print #1, " Hwk=";: Print #1, using; "##.##"; hwk

Print #1, " He0=";: Print #1, using; "##.##"; he0;

Print #1, " Hemax=";: Print #1, using; "##.##"; hem

Print #1, " Hg0=";: Print #1, using; "##.##"; hg0;

Print #1, " Kf0=";: Print #1, using; "##.###"; kf0;

Print #1, " Hp=";: Print #1, using; "##.##"; hp;

Print #1, " Fk=";: Print #1, using; "##.####"; fk;

Print #1, " D0=";: Print #1, using; "##.####"; d0;

Print #1, " Hdr=";: Print #1, using; "##.##"; hdr

Print #1, " Tveg=";: Print #1, using; "###"; tveg;

Print #1, " Dte=";: Print #1, using; "##.##"; dte;

Print #1, " Dtu=";: Print #1, using; "##.##"; dtu;


5.     ,

 

 

 

 

 

6.

 

 

 

7.

 

 

 

8.

 

 

 

9. ʻ

 

' Programma rascheta poliva i urodgaja

' na zasolennyh zemljah

' Avtor prof.dtn. A.I.GOLOVANOV

' M G M I

' w=wm+(p-wm)*exp(-psi/mju/hk)

' S grafikoj, Vs=w*D*dc/dx, PPK

' C - koncentracija r-ra,kgekv/m^3 vody

' S - k-vo iona v PPK,kgekv/m^3 pochvy

' Isotermy: Sna/Sca^.5=Knc*w*Cna/(w*Cca)^.5

' Sca/Smg=Kmc*w*Cca/(w*Cmg)

 

' M N O G O L E T I J E

' Kfzimn=Kf*exp(-4*(t-tveg)/(tgod-tveg)); Kf zimn.=.001Kf v sloje 1.2m

' Rezultat: Gzim/G=.67, Dzim/D=.72

' Gipsovanie s promivkoj

DefInt I-J

Dim p(40), mju(40), wm(40), kf(40), hk(40), h(40), w0(40), x(40), w1(40)

Dim w2(40), w3(40), h1(40), h2(40), r(40), pe(40), q(40), er(40), s(25)

Dim s1(25), oc(20, 75), toc(20, 75), noc(75), ef(40), ek(20, 37), pol(20), tpol(20), tsu(20)

Dim teu(20), tets(40), hgs(40), hgm(40), gs(40), kw(40), al(40), plt(40)

Dim n0(40), sna0(40), sca0(40), smg0(40), n1(40), sna1(40), sca1(40)

Dim smg1(40), cna0(40), cca0(40), cmg0(40), cna1(40), cca1(40), cmg1(40)

Dim vs1(40), vs2(40), sppk(40), pec(40), pem(40), qc(40), qm(40), ga(40)

Dim tetopt(40), alfs(40), bets(40), gams(40), sol(9, 400)

Dim z(10), cna2(40), cca2(40), cmg2(40), ys0(15), ys1(15), zmn(19, 50), zmn1(9, 50)

Dim sumzmn(20), sum(20)

'****************************************************************************

k1$ = "rez_07.bas"

k2$ = " ELISTA, Luzerna, 10 let "

open "o",#1,k1$

open "i",#4,"tvvolg.bas"

open "i",#5,"ocvolg.bas"

 

' Uslovija

l1 = -1 '-1: hg=peremen, 0: hg=const, +1: hg=gluboko

l2 = -1 '-1: 1-var progr, 0: 2-var progr, +1: 3-var progr

l3 = -1 '-1: net vivoda cherez dtv, 0: vivod w2(I), +1: vivod cherez dtv

l4 = -1 '-1: w0(I)= ravnovesnaja, +1: w0(I)= ljubaja

l6 = -1 '-1: bez prognoza Oc; +1: prognoz Oc

n = 35

h(1) = 0: h(2) = 0.05: h(3) = 0.05: h(4) = 0.05: h(5) = 0.05: h(6) = 0.05: h(7) = 0.05: h(8) = 0.05

h(9) = 0.05: h(10) = 0.05: h(11) = 0.05: h(12) = 0.05: h(13) = 0.05: h(14) = 0.05: h(15) = 0.05

h(16) = 0.05: h(17) = 0.05: h(18) = 0.05: h(19) = 0.05: h(20) = 0.05: h(21) = 0.05: h(22) = 0.05

h(23) = 0.1: h(24) = 0.1: h(25) = 0.15: h(26) = 0.2: h(27) = 0.2: h(28) = 0.2: h(29) = 0.2

h(30) = 0.2: h(31) = 0.2: h(32) = 0.2: h(33) = 0.2: h(34) = 0.2: h(35) = 0

'h(36)=.2:h(37)=.2:h(38)=.2:h(39)=0

horA = 0.5: horB1 = 0.3: horB2 = 0.2: horC = 5

pA = 0.52: wmA = 0.11: mjuA = 1: hkA = 2#: kfA = 0.4: pltA = 1.2

pB1 = 0.48: wmB1 = 0.1: mjuB1 = 1: hkB1 = 1.8: kfB1 = 0.25: pltB1 = 1.4

pB2 = 0.47: wmB2 = 0.09: mjuB2 = 1: hkB2 = 1.5: kfB2 = 0.1: pltB2 = 1.35

pC = 0.5: wmC = 0.12: mjuC = 1: hkC = 1.4: kfC = 0.25: pltC = 1.3

nA = 3.5: snaA = 1#: scaA = 27: smgA = 7#

nB1 = 1.5: snaB1 = 0.5: scaB1 = 21: smgB1 = 6.5

nB2 = 1#: snaB2 = 0.2: scaB2 = 13.8: smgB2 = 3#

nC = 0.5: snaC = 0.1: scaC = 6#: smgC = 0.9

'nA i t.d. - (Na+Ca+Mg),mgekv/100gr pochvi v vodnoj vytjagke

'snaA i t.d.-Na,Ca,Mg,mgekv/100gr.pochvi v PPK

'sna+sca+smg=emkoct pogloschenija

wz = 1.5 * wmA: ppv = 0.62 * pA

mpolmin = 50

reg = 0.6: Kbiol = 0.8

hw0 = 0.5: hwm = 0.5: hwk = 0.5: he0 = 0.9: hem = 0.9

hg0 = 3#: kf0 = 0#: t0 = 10: hp = 1.5: d0 = 0.002: hdr = 3#

fk = 0 '0.0005

tgod = 365: tveg = 200: ttepl = 250

dte = 0.25: dtu = 0.02: dtef = 10: dtv = 5

Uvlves = 50 'Vesen.uvl,mm

ngod = 20: npresn = 0

L = 0.5: Dm = 0.0004: Knca = 0.6: Kmca = 1.2: Beta = 0.05

cpvnz = 0.6: cpvcz = 0.4: cpvmz = 0.2 'mineraliz. polivn.vodi,gr/litp

cpvnpres = 0.1: cpvcpres = 0.5: cpvmpres = 0.1 'mineraliz. presn.polivn.vodi,gr/litp

horsol = 0.3 'dlja grafika solej

horg = 1 'dlja vodoobmena

ekna = 23: ekca = 20: ekmg = 12

sosdop = 5 'dop. stepen osolonzev,%

Dozagipsa = 20 't/ha

hgi = horsol '.4 'sloj gipsivanija,m

tgips = 171

 

 

 

al(1) = 0.03: al(2) = 0.04: al(3) = 0.05: al(4) = 0.08: al(5) = 0.09: al(6) = 0.1

al(7) = 0.09: al(8) = 0.03: al(9) = 0.06: al(10) = 0.08: al(11) = 0.08: al(12) = 0.06

al(13) = 0.05: al(14) = 0.04: al(15) = 0.03: al(16) = 0.03: al(17) = 0.03

al(18) = 0.03: al(19) = 0: al(20) = 0

al(21) = 0: al(22) = 0: al(23) = 0

 

ga(1) = 6.6: ga(2) = 6.6: ga(3) = 6.6: ga(4) = 6.6: ga(5) = 6.6: ga(6) = 6.6

ga(7) = 6.6: ga(8) = 6.6: ga(9) = 6.6: ga(10) = 6.6: ga(11) = 6.6: ga(12) = 6.6

ga(13) = 6.6: ga(14) = 6.6: ga(15) = 6.6: ga(16) = 6.6: ga(17) = 6.6

ga(18) = 6.6: ga(19) = 6.6: ga(20) = 6.6: ga(21) = 6.6

'ga(22)=6.6:ga(23)=6.6

tetopt(1) = 0.26: tetopt(2) = 0.26: tetopt(3) = 0.26: tetopt(4) = 0.26

tetopt(5) = 0.26: tetopt(6) = 0.26: tetopt(7) = 0.26: tetopt(8) = 0.26

tetopt(9) = 0.26: tetopt(10) = 0.26: tetopt(11) = 0.26: tetopt(12) = 0.26

tetopt(13) = 0.26: tetopt(14) = 0.26: tetopt(15) = 0.26: tetopt(16) = 0.26

tetopt(17) = 0.26: tetopt(18) = 0.26: tetopt(19) = 0.26: tetopt(20) = 0.26

tetopt(21) = 0.26: tetopt(22) = 0.26: tetopt(23) = 0.26

 

'Efis, doli ot Esum !

ef(1) = 0.2: ef(2) = 0.2: ef(3) = 0.2: ef(4) = 0.15: ef(5) = 0.1: ef(6) = 0.1: ef(7) = 0.1

ef(8) = 0.1: ef(9) = 0.4: ef(10) = 0.3: ef(11) = 0.2: ef(12) = 0.1: ef(13) = 0.1: ef(14) = 0.1

ef(15) = 0.15: ef(16) = 0.2: ef(17) = 0.25: ef(18) = 0.25: ef(19) = 0.3: ef(20) = 0.3: ef(21) = 0.3

ef(22) = 0.3: ef(23) = 0.4

kfar = 1

 

npol = 1

pol(1) = 300: pol(2) = 40: pol(3) = 40: pol(4) = 40: pol(5) = 40

'tpol(1)=240

tpol(1) = 5555

'tpol(2)=39:tpol(3)=51:tpol(4)=86:tpol(5)=107

 

ntakt = 1 'dlja kapitalnyh promyvok

tsu(1) = 200: tsu(2) = 500: tsu(3) = 80

teu(1) = 201.2: teu(2) = 600: teu(3) = 90

'Nachaln.epjura wlagnosti

w0(1) = 0.186: w0(2) = 0.187: w0(3) = 0.189: w0(4) = 0.191: w0(5) = 0.193: w0(6) = 0.194

w0(7) = 0.196: w0(8) = 0.198: w0(9) = 0.201: w0(10) = 0.203: w0(11) = 0.205: w0(12) = 0.176

w0(13) = 0.177: w0(14) = 0.178: w0(15) = 0.179: w0(16) = 0.18: w0(17) = 0.181

w0(18) = 0.157: w0(19) = 0.157: w0(20) = 0.156: w0(21) = 0.155

w0(22) = 0.192: w0(23) = 0.193: w0(24) = 0.193: w0(25) = 0.194: w0(26) = 0.195

w0(27) = 0.196: w0(28) = 0.197: w0(29) = 0.197: w0(30) = 0.197: w0(31) = 0.197

'w0(32)=.14:w0(33)=.14

'***************************************************************************

360 k3$ = ", Bez oroschenija": k4$ = ", Oroschenie": k5$ = ", Prognoz os."

k5s$ = ", Soli s PPK"

tcal0 = Timer

l5 = -1 '-1: Cw=analit, +1: Cw=dif

tetmax = (ppv - wz) / (pA - wz)

m = -1

nev = Int(tveg / dtef): nd = Int(tgod / 10) + 1

If l2 = 0 Then z$ = k3$

If l2 = -1 Then z$ = k4$

If l6 = 1 And l2 = -1 Then z1$ = k5$

zs$ = k5s$

Print #1,: Print #1,: Print #1, k2$ + z$ + z1$ + zs$: Print #1,

Print #1, " Ishodnaja informazija"

'Vvod Oc(j,i), Toc(j,i), Ek(j,i)

nlet = 20: nisp = 20

Print #1, " Ngod God Osadki Isparen Defizit"

Print " Ngod God Osadki Isparen Defizit"

For j = 1 To nlet

For i = 1 To nisp: Input #4, tvoz, vlvoz

ek(j, i) = 0.0061 * (25 + tvoz) ^ 2 * (1 - 0.01 * vlvoz) * Kbiol

sum(j) = sum(j) + 10 * ek(j, i)

Next

sumisp = sumisp + sum(j)

Next

For j = 1 To nlet

sum1 = 0

Input #5, noc(j)

For i = 1 To noc(j)

Input #5, oc(j, i), toc(j, i)

oc(j, i) = 1 * oc(j, i) '1.274

sum1 = sum1 + oc(j, i)

Next

god = 1967 + j

Print using; "######"; j; god;: Print using; "#########"; sum1; sum(j); sum(j) - sum1

Print #1, using; "######"; j; god;

Print #1, using; "#########"; sum1; sum(j); sum(j) - sum1

sumosadk = sumosadk + sum1

Next

Print #1, " Summa ";

Print #1, using; "#########"; sumosadk; sumisp; sumisp - sumosadk

Print #1, " Srednee ";

Print #1, using; "#########"; sumosadk / nlet; sumisp / nlet; (sumisp - sumosadk) / nlet