Eduardo Eneas de Figueiredo

7.3.Produtos Técnicos

RAVINA

10 '-------------------------------------------------------

20 '------ PROGRAMA: ALARGAMENTO LATERAL DE RAVINAS -------

30 ' AUTOR: EDUARDO ENEAS DE FIGUEIREDO

40 '-------------------------------------------------------

50 '

60 ' DEFINE VARIAVEIS STRINGS

70 '

80 CLEAR:DEFSTR Z:KEY OFF

90 Z1=CHR$(12):Z2=CHR$(7):Z3="ALARGAMENTO LATERAL DE RAVINAS"

100 Z4=STRING$(70,"-"):Z5=STRING$(30,"-"):Z6="ALARGAMENTO TOTAL INFILTRACAO"

110 Z7="ALARGAMENTO TOTAL LENCOL E INFILTRACAO"

120 Z8="ALARGAMENTO PARCIAL FENDA DE TRACAO NA SUPERF. LENCOL E INFILT."

130 Z9=CHR$(13):Z0=CHR$(32)

140 Z11="ALARGAMENTO PARCIAL FENDA DE TRACAO NA PAREDE, LENCOL E INFILT."

150 Z12="ALARGAMENTO PARCIAL FENDA DE TRACAO NA SUPERFICIE E INFILTRACAO"

160 '

170 ' DEFINE CONSTANTES

180 '

190 C=0:FC=3.1416/180:HMA=2000:HMB=1000

200 '

210 ' APRESENTA OPCOES

220 '

230 ' MUDA PAGINA E APRESENTA CABECARIO 1

240 '

250 GOSUB 6220

260 PRINT TAB(10)" TIPO DO ALARGAMENTO CAUSA (S) DO ALARGAMENTO":PRINT TAB(10)Z4

270 PRINT TAB(10)"1 - ";Z6

280 PRINT TAB(10)"2 - ";Z7

290 PRINT TAB(10)"3 - ";Z8

300 PRINT TAB(10)"4 - ";Z11

310 PRINT TAB(10)"5 - ";Z12

320 PRINT TAB(10)Z4

330 PRINT TAB(10)"6 - FINALIZA"

340 PRINT TAB(10)Z4:PRINT:PRINT

350 PRINT TAB(10)"INDIQUE OPERACAO DESEJADA -> ";

360 Z=INPUT$(1):PRINT Z:IF Z<"1" OR Z>"6" THEN PRINT Z2;:GOTO 360

370 GOSUB 7170

380 IF ZA=Z0 THEN GOTO 200

390 ON VAL(Z) GOTO 400, 1880, 3260, 3290, 4590, 7380

400 '

410 ' ALARGAMENTO LATERAL PELA INFILTRACAO

420 '

430 ' MUDA PAGINA E APRESENTA ENTRADA DE DADOS

440 '

450 GOSUB 6220

460 PRINT TAB(10) "ALARGAMENTO TOTAL PELA INFILTRACAO ENTRADA DE DADOS DO SOLO"

470 PRINT TAB(10)Z4

480 INPUT" Capacidade de infiltracao maxima em mm/hora 'qo =";A1

490 INPUT" Capacidade de infiltracao minima em mm/hora 'qf =";A2

500 INPUT" Fator de amortecimento de Horton em (1/hora) ' k =";A3

510 INPUT" Angulo da parede da ravina em graus ' i =";A4

520 INPUT" Umidade inicial do solo em % do peso 'Uo =";A6

530 INPUT" Umidade do solo na saturacao em % do peso 'Us =";A7

540 INPUT" Coesao do solo na umidade inicial em ton/m2 'Co =";A8

550 INPUT" Coesao do solo na saturacao em ton/m2 'Cs =";A9

560 INPUT" Densidade dos graos ' G =";A10

570 INPUT" Indice de vazios do solo ' e =";A11

580 INPUT" Peso especifico do solo seco em ton/m3 'Ps =";A12

590 INPUT" Angulo de atrito na umidade inicial em graus 'Fo =";A13

600 INPUT" Angulo de atrito na saturacao em graus 'Fs =";A14

610 PRINT TAB(10)Z4:GOSUB 7170

620 IF ZA=Z0 THEN GOTO 400

630 '

640 ' IMPRESSAO NA TELA OU NO PAPEL

650 '

660 GOSUB 7240

670 IF ZA=Z9 THEN 680 ELSE 930

680 '

690 ' IMPRESSAO NA TELA

700 '

710 GOSUB 6220

720 PRINT TAB(10)"ALARGAMENTO TOTAL PELA INFILTRACAO DADOS FISICOS DO SOLO"

730 PRINT TAB(10)Z4

740 PRINT TAB(10)"Capacidade de infiltracao maxima qo = ";:PRINT USING "##.##"; A1;: PRINT " mm/h"

750 PRINT TAB(10)"Capacidade de infiltracao minima qf = ";:PRINT USING "##.##"; A2;: PRINT " mm/h"

760 PRINT TAB(10)"Fator de amortecimento de Horton k = ";:PRINT USING "##.##"; A3;: PRINT " 1/h"

770 PRINT TAB(10)"Angulo da parede da ravina i = ";:PRINT USING "##.##"; A4;: PRINT " graus"

780 PRINT TAB(10)"Umidade inicial do solo Uo = ";:PRINT USING "##.##"; A6;: PRINT " %"

790 PRINT TAB(10)"Umidade do solo na saturacao Us = ";:PRINT USING "##.##"; A7;: PRINT " %"

800 PRINT TAB(10)"Coesao do solo na umidade inicial Co = ";:PRINT USING "##.##"; A8;: PRINT " t/m2"

810 PRINT TAB(10)"Coesao do solo na saturacao Cs = ";:PRINT USING "##.##"; A9;: PRINT " t/m2"

820 PRINT TAB(10)"Densidade dos graos G = ";:PRINT USING "##.##"; A10

830 PRINT TAB(10)"Indice de vazios do solo e = ";:PRINT USING "##.##"; A11

840 PRINT TAB(10)"Peso especifico do solo seco Ps = ";:PRINT USING "##.##"; A12;: PRINT " t/m3"

850 PRINT TAB(10)"Angulo de atrito na umidade inicial Fo = ";:PRINT USING "##.##"; A13;: PRINT " graus"

860 PRINT TAB(10)"Angulo de atrito na saturacao Fs = ";:PRINT USING "##.##"; A14;: PRINT " graus"

870 PRINT TAB(10)Z4:GOSUB 7170

880 IF ZA=Z0 THEN GOTO 400

890 GOTO 1270

900 '

910 ' IMPRESSAO NO PAPEL

920 '

930 LPRINT:LPRINT:LPRINT:LPRINT:LPRINT:LPRINT:LPRINT:LPRINT:LPRINT:LPRINT:LPRINT:LPRINT:LPRINT

940 LPRINT TAB(30)Z5:LPRINT TAB(30)Z3:LPRINT TAB(30)Z5

950 LPRINT:LPRINT:LPRINT TAB(10)Z4

960 LPRINT TAB(10)"ALARGAMENTO TOTAL PELA INFILTRACAO DADOS FISICOS DO SOLO"

970 LPRINT TAB(10)Z4

980 LPRINT TAB(10)"Capacidade de infiltracao maxima q = ";:LPRINT USING "##.##"; A1;: LPRINT " mm/h"

990 LPRINT

1000 LPRINT TAB(10)"Capacidade de infiltracao minima q = ";:LPRINT USING "##.##"; A2;: LPRINT " mm/h"

1010 LPRINT

1020 LPRINT TAB(10)"Fator de amortecimento de Horton k = ";:LPRINT USING "##.##"; A3;: LPRINT " h "

1030 LPRINT

1040 LPRINT TAB(10)"Umidade inicial do solo U = ";:LPRINT USING "##.##"; A6;: LPRINT " %"

1050 LPRINT

1060 LPRINT TAB(10)"Umidade do solo na saturacao U = ";:LPRINT USING "##.##"; A7;: LPRINT " %"

1070 LPRINT

1080 LPRINT TAB(10)"Coesao do solo na umidade inicial C = ";:LPRINT USING "##.##"; A8;: LPRINT " t/m "

1090 LPRINT

1100 LPRINT TAB(10)"Coesao do solo na saturacao C = ";:LPRINT USING "##.##"; A9;: LPRINT " t/m "

1110 LPRINT

1120 LPRINT TAB(10)"Densidade dos graos G = ";:LPRINT USING "##.##"; A10

1130 LPRINT

1140 LPRINT TAB(10)"Indice de vazios do solo = ";:LPRINT USING "##.##"; A11

1150 LPRINT

1160 LPRINT TAB(10)"Peso especifico do solo seco = ";:LPRINT USING "##.##"; A12;: LPRINT " t/m "

1170 LPRINT

1180 LPRINT TAB(10)"Angulo de atrito na umidade inicial = ";:LPRINT USING "##.##"; A13;: LPRINT " graus"

1190 LPRINT

1200 LPRINT TAB(10)"Angulo de atrito na saturacao = ";:LPRINT USING "##.##"; A14;: LPRINT " graus"

1210 LPRINT TAB(10)Z4:GOSUB 7170

1220 IF ZA=Z0 THEN 400

1230 GOTO 1590

1240 '

1250 ' IMPRESSAO NA TELA

1260 '

1270 GOSUB 6290

1280 FOR I=A4 TO 90 STEP 15

1290 FOR FI=10 TO 30 STEP 5

1300 FA=(1-COS((I-FI)*FC))*(A12/1000)/(4*SIN(I*FC)*COS(FI*FC))

1310 FFI=(1-COS((I-A14)*FC))*(A12/1000)/(4*SIN(I*FC)*COS(A14*FC))

1320 HMI=(A9/10)/(FFI*(1+A7/100)):HM1=(A9/10)/(FFI*(1+A7/100))

1330 HMI=100*INT(HMI/100+1)

1340 FOR H=HMI TO HMB STEP 50

1350 C1=H*FA

1360 C2=10*(A8-A9)/(A7-A6)

1370 UC=((A8/10+C2*A6/100-C1)/(C1+C2))*100

1380 GOSUB 6960

1390 FB=SIN((I-FI)*FC/2)/(SIN(I*FC)*SIN((I+FI)*FC/2))

1400 BC=H*FB

1410 IF I=90 GOTO 1430

1420 FD=TAN(I*FC)/TAN((1/2)*((I+FI)*FC)):GOTO 1440

1430 Q=(UC/100)*A10*H/(2*(1+A11)):GOTO 1450

1440 Q=(UC/100)*A10*H*(1-1/FD)/(2*(1+A11))

1450 T=10*Q/A2 - (A1-A2)/(A3*A2)

1460 GOSUB 7080

1470 PRINT TAB(13);:PRINT USING "##.##";I;:PRINT TAB(23);:PRINT USING "##.##";FI;:PRINT TAB(32);:PRINT USING "####.#";H;:PRINT TAB(43);:PRINT USING "##.##";UC;:PRINT TAB(52);:PRINT USING "####.#";BC;: GOTO 1480

1480 PRINT TAB(63);:PRINT USING "###.#";Q;:PRINT TAB(73);:PRINT USING "####.#";B3:C=C+1

1490 IF I=90 AND H=HMB THEN C=0:PRINT TAB(10)Z4:PRINT TAB(10)"Altura Minima = ";:PRINT USING "####.#"; HM1;: PRINT " cm ":PRINT TAB(10)Z4:GOTO 1540

1500 IF H=HMB THEN C=0:PRINT TAB(10)Z4:PRINT TAB(10)"Altura Minima = ";:PRINT USING "####.#"; HM1;:PRINT " cm ":PRINT TAB(10)Z4:GOSUB 6290:GOTO 1530:C=C+1

1510 NEXT H

1520 NEXT FI

1530 NEXT I

1540 GOSUB 7310

1550 GOTO 10

1560 '

1570 ' IMPRESSAO NO PAPEL

1580 '

1590 GOSUB 6370

1600 FOR I=A4 TO 90 STEP 15

1610 FOR FI=10 TO 30 STEP 5

1620 FA=(1-COS((I-FI)*FC))*(A12/1000)/(4*SIN(I*FC)*COS(FI*FC))

1630 FFI=(1-COS((I-A14)*FC))*(A12/1000)/(4*SIN(I*FC)*COS(A14*FC))

1640 HMI=(A9/10)/(FFI*(1+A7/100)):HM1=(A9/10)/(FFI*(1+A7/100))

1650 HMI=100*INT(HMI/100+1)

1660 FOR H=HMI TO HMA STEP 50

1670 C1=H*FA

1680 C2=10*(A8-A9)/(A7-A6)

1690 UC=((A8/10+C2*A6/100-C1)/(C1+C2))*100

1700 GOSUB 6960

1710 FB=SIN((I-FI)*FC/2)/(SIN(I*FC)*SIN((I+FI)*FC/2))

1720 BC=H*FB

1730 IF I=90 GOTO 1750

1740 FD=TAN(I*FC)/TAN((1/2)*((I+FI)*FC)):GOTO 1760

1750 Q=(UC/100)*A10*H/(2*(1+A11)):GOTO 1770

1760 Q=(UC/100)*A10*H*(1-1/FD)/(2*(1+A11))

1770 T=10*Q/A2 - (A1-A2)/(A3*A2)

1780 GOSUB 7080

1790 LPRINT TAB(13);:LPRINT USING "##.##";I;:LPRINT TAB(23);:LPRINT USING "##.##";FI;:LPRINT TAB(32);:LPRINT USING "####.#";H;:LPRINT TAB(43);:LPRINT USING "##.##";UC;:LPRINT TAB(52);:LPRINT USING "####.#";BC;:GOTO 1800

1800 LPRINT TAB(63);:LPRINT USING "###.#";Q;:LPRINT TAB(73);:LPRINT USING "####.#";B3:C=C+1

1810 IF I=90 AND H=HMA THEN C=0:LPRINT TAB(10)Z4:LPRINT TAB(10)"Altura Minima = ";:LPRINT USING "####.#"; HM1;: LPRINT " cm ":LPRINT TAB(10)Z4:GOTO 1860

1820 IF H=HMA THEN C=0:LPRINT TAB(10)Z4:LPRINT TAB(10)"Altura Minima = ";:LPRINT USING "####.#"; HM1;:LPRINT " cm ":LPRINT TAB(10)Z4:GOSUB 6370:GOTO 1850:C=C+1

1830 NEXT H

1840 NEXT FI

1850 NEXT I

1860 GOSUB 7310

1870 GOTO 10

1880 '

1890 ' ALARGAMENTO PELO LENCOL: SISTEMA TOTAL (Z=2); SISTEMA PARCIAL (Z=3)

1900 '

1910 ' MUDA PAGINA E APRESENTA ENTRADA DE DADOS

1920 '

1930 GOSUB 6220

1940 IF VAL(Z)=3 THEN 1960 ELSE 1950

1950 PRINT TAB(10)"ALARGAMENTO TOTAL PELO LENCOL DADOS FISICOS DO SOLO":GOTO 1970

1960 PRINT TAB(10)"ALARGAMENTO PARCIAL:LENCOL+FENDA NA SUP. DADOS FISICOS DO SOLO"

1970 PRINT TAB(10)Z4

1980 INPUT" Angulo da parede da ravina em graus ' i =";A4

1990 INPUT" Umidade inicial do solo em % do peso 'Uo =";A6

2000 INPUT" Umidade do solo na saturacao em % do peso 'Us =";A7

2010 INPUT" Coesao do solo na umidade inicial em ton/m2 'Co =";A8

2020 INPUT" Coesao do solo na saturacao em ton/m2 'Cs =";A9

2030 INPUT" Peso especifico do solo seco em ton/m3 'Ps =";A12

2040 INPUT" Angulo de atrito na umidade inicial em graus 'Fo =";A13

2050 INPUT" Angulo de atrito na saturacao em graus 'Fs =";A14

2060 PRINT TAB(10)Z4:GOSUB 7170

2070 IF ZA=Z0 THEN GOTO 1880

2080 '

2090 ' IMPRESSAO NA TELA OU NO PAPEL

2100 '

2110 GOSUB 7240

2120 IF ZA=Z9 THEN 2130 ELSE 2350

2130 '

2140 ' IMPRESSAO NA TELA

2150 '

2160 GOSUB 6220

2170 IF VAL(Z)=3 THEN 2190 ELSE 2180

2180 PRINT TAB(10)"ALARGAMENTO TOTAL PELO LENCOL DADOS FISICOS DO SOLO":GOTO 2200

2190 PRINT TAB(10)"ALARGAMENTO PARCIAL:LENCOL+FENDA NA SUP. DADOS FISICOS DO SOLO":GOTO 2200

2200 PRINT TAB(10)Z4

2210 PRINT TAB(10)"Angulo da parede da ravina i = ";:PRINT USING "##.##"; A4;: PRINT " graus"

2220 PRINT TAB(10)"Umidade inicial do solo Uo = ";:PRINT USING "##.##"; A6;: PRINT " %"

2230 PRINT TAB(10)"Umidade do solo na saturacao Us = ";:PRINT USING "##.##"; A7;: PRINT " %"

2240 PRINT TAB(10)"Coesao do solo na umidade inicial Co = ";:PRINT USING "##.##"; A8;: PRINT " t/m2"

2250 PRINT TAB(10)"Coesao do solo na saturacao Cs = ";:PRINT USING "##.##"; A9;: PRINT " t/m2"

2260 PRINT TAB(10)"Peso especifico do solo seco Ps = ";:PRINT USING "##.##"; A12;: PRINT " t/m3"

2270 PRINT TAB(10)"Angulo de atrito na umidade inicial Fo = ";:PRINT USING "##.##"; A13;: PRINT " graus"

2280 PRINT TAB(10)"Angulo de atrito na saturacao Fs = ";:PRINT USING "##.##"; A14;: PRINT " graus"

2290 PRINT TAB(10)Z4:GOSUB 7170

2300 IF ZA=Z0 THEN GOTO 1880

2310 GOTO 2590

2320 '

2330 ' IMPRESSAO NO PAPEL

2340 '

2350 LPRINT:LPRINT:LPRINT:LPRINT:LPRINT:LPRINT:LPRINT:LPRINT:LPRINT:LPRINT:LPRINT:LPRINT:LPRINT:LPRINT:LPRINT

2360 LPRINT TAB(30)Z5:LPRINT TAB(30)Z3:LPRINT TAB(30)Z5

2370 LPRINT:LPRINT:LPRINT TAB(10)Z4

2380 IF VAL(Z)=2 THEN 2400 ELSE 2390

2390 LPRINT TAB(10)"ALARGAMENTO PARCIAL:LENCOL+FENDA NA SUP. DADOS FISICOS DO SOLO":GOTO 2410

2400 LPRINT TAB(10)"ALARGAMENTO TOTAL PELO LENCOL DADOS FISICOS DO SOLO"

2410 LPRINT TAB(10)Z4

2420 LPRINT TAB(10)"Umidade inicial do solo U = ";:LPRINT USING "##.##"; A6;: LPRINT " %"

2430 LPRINT

2440 LPRINT TAB(10)"Umidade do solo na saturacao U = ";:LPRINT USING "##.##"; A7;: LPRINT " %"

2450 LPRINT

2460 LPRINT TAB(10)"Coesao do solo na umidade inicial C = ";:LPRINT USING "##.##"; A8;: LPRINT " t/m "

2470 LPRINT

2480 LPRINT TAB(10)"Coesao do solo na saturacao C = ";:LPRINT USING "##.##"; A9;: LPRINT " t/m "

2490 LPRINT

2500 LPRINT TAB(10)"Peso especifico do solo seco = ";:LPRINT USING "##.##"; A12;: LPRINT " t/m "

2510 LPRINT

2520 LPRINT TAB(10)"Angulo de atrito na umidade inicial = ";:LPRINT USING "##.##"; A13;: LPRINT " graus"

2530 LPRINT

2540 LPRINT TAB(10)"Angulo de atrito na saturacao = ";:LPRINT USING "##.##"; A14;: LPRINT " graus"

2550 LPRINT TAB(10)Z4:GOSUB 7170

2560 IF ZA=Z0 THEN 1880

2570 GOTO 2950

2580 '

2590 ' IMPRESSAO NA TELA

2600 '

2610 GOSUB 6460

2620 FOR I=A4 TO 90 STEP 15

2630 FOR U=A6 TO A7 STEP 5.6

2640 FI=A13-((A13-A14)*100/(A7-A6))*(U/100-A6/100)

2650 FA=(1-COS((I-FI)*FC))*(A12/1000)/(4*SIN(I*FC)*COS(FI*FC))

2660 HMI=(A9/10)/(FA*(1+A7/100))

2670 HMI=100*INT(HMI/100+1)

2680 FOR H=HMI TO HMB STEP 50

2690 IF I=90 THEN H1=(1/2)*H:GOTO 2710

2700 H1=H*(1-SQR(TAN((1/2)*((I+FI)*FC))/TAN(I*FC)))

2710 A=(A7/100-U/100)*FA

2720 B=(A8/10-A9/10)-10*(U/100-A6/100)*(A8-A9)/(A7-A6)

2730 IF VAL(Z) = 2 GOTO 2750

2740 C = -H1^2*.5*A12*(1+U/100)*COS((1/2)*((I+FI)*FC))*SIN((1/2)*((I-FI)*FC))/(1000*COS(FI*FC))-(A8/10-10*(U/100-A6/100)*(A8-A9)/(A7-A6))*(H-H1)+FA*(1+U/100)*H^2:GOTO 2760

2750 C = FA*H^2*(1+U/100)-(A8/10-10*(U/100-A6/100)*(A8-A9)/(A7-A6))*H

2760 IF U=A7 GOTO 2810

2770 D=B^2-4*A*C

2780 E=((D)^.5-B)/(2*A)

2790 IF E<=0 THEN E=0

2800 PRINT TAB(14);:PRINT USING "##.##";I;:PRINT TAB(28);:PRINT USING "##.##";FI;:PRINT TAB(41);:PRINT USING "####.#";H;:PRINT TAB(56);:PRINT USING "##.##";U;:PRINT TAB(69);:PRINT USING "#####.#";E;:GOTO 2840:C=C+1

2810 PRINT " SE A UMIDADE INICIAL E A UMIDADE SATURADA,A RUPTURA OCORERA PARA QUAL- QUER ALTURA MAIOR OU IGUAL A ALTURA MINIMA INDEPENDENTEMENTE DA ZONA SATURADA":PRINT TAB(10)Z4

2820 IF I=90 GOTO 2900

2830 GOSUB 6460:GOTO 2890

2840 IF E<=0 THEN C=0:PRINT TAB(10)Z4:GOSUB 6460:GOTO 2880:C=C+1

2850 IF I=90 AND U=A7 AND H=HMB THEN C=0:PRINT TAB(10)Z4:GOTO 2900

2860 IF H=HMB THEN C=0:PRINT TAB(10)Z4:GOSUB 6460:C=C+1

2870 NEXT H

2880 NEXT U

2890 NEXT I

2900 GOSUB 7310

2910 GOTO 10

2920 '

2930 ' IMPRESSAO NO PAPEL

2940 '

2950 GOSUB 6540

2960 FOR I=A4 TO 90 STEP 15

2970 FOR U=A6 TO A7 STEP 5.6

2980 FI=A13-((A13-A14)*100/(A7-A6))*(U/100-A6/100)

2990 FA=(1-COS((I-FI)*FC))*(A12/1000)/(4*SIN(I*FC)*COS(FI*FC))

3000 HMI=(A9/10)/(FA*(1+A7/100))

3010 HMI=100*INT(HMI/100+1)

3020 FOR H=HMI TO HMA STEP 50

3030 IF I=90 THEN H1=(1/2)*H:GOTO 3050

3040 H1=H*(1-SQR(TAN((1/2)*((I+FI)*FC))/TAN(I*FC)))

3050 A=(A7/100-U/100)*FA

3060 B=(A8/10-A9/10)-10*(U/100-A6/100)*(A8-A9)/(A7-A6)

3070 IF VAL(Z) = 2 GOTO 3090

3080 C = -H1^2*.5*A12*(1+U/100)*COS((1/2)*((I+FI)*FC))*SIN((1/2)*((I-FI)*FC))/(1000*COS(FI*FC))-(A8/10-10*(U/100-A6/100)*(A8-A9)/(A7-A6))*(H-H1)+FA*(1+U/100)*H^2:GOTO 3100

3090 C = FA*H^2*(1+U/100)-(A8/10-10*(U/100-A6/100)*(A8-A9)/(A7-A6))*H

3100 IF U=A7 GOTO 3150

3110 D=B^2-4*A*C

3120 E=((D)^.5-B)/(2*A)

3130 IF E<=0 THEN E=0

3140 LPRINT TAB(14);:LPRINT USING "##.##";I;:LPRINT TAB(28);:LPRINT USING "##.##";FI;:LPRINT TAB(41);:LPRINT USING "####.#";H;:LPRINT TAB(56);:LPRINT USING "##.##";U;:LPRINT TAB(69);:LPRINT USING "#####.#";E:GOTO 3180:C=C+1

3150 LPRINT " SE A UMIDADE INICIAL E A UMIDADE SATURADA,A RUPTURA OCORERA PARA QUAL- QUER ALTURA MAIOR OU IGUAL A ALTURA MINIMA INDEPENDENTEMENTE DA ZONA SATURADA":LPRINT TAB(10)Z4

3160 IF I=90 GOTO 3240

3170 GOSUB 6540:GOTO 3230

3180 IF E<=0 THEN C=0:LPRINT TAB(10)Z4:GOSUB 6540:GOTO 3220:C=C+1

3190 IF I=90 AND U=A7 AND H=HMA THEN C=0:LPRINT TAB(10)Z4:GOTO 3240

3200 IF H=HMA THEN C=0:LPRINT TAB(10)Z4:GOSUB 6540:C=C+1

3210 NEXT H

3220 NEXT U

3230 NEXT I

3240 GOSUB 7310

3250 GOTO 10

3260 '

3270 ' ALARGAMENTO PARCIAL PELO LENCOL - FENDA NA SUPERFICIE (Z=3)

MOFIPE

PROGRAM MOFIP

c---Matriz coluna para arquivos de saída do fluxo

Dimension

* qt(15), qs(15),

* qss(15), qtCo(15),

* qtMu1(15), qtMu2(15),

* qtMu3(15), qtMu4(15),

* flowt(15), flowt9(15),

* flowt1(15), flowt2(15),

* flowt3(15), flowt4(15),

* flowt5(15), flowt6(15),

* flowt7(15), flowt8(15)

c---VARIÁVEIS DO SOLO: Saxton et al(1986), R.& Brakensiek(1989), B. & Corey (1964)

real*8

* Klamb1, Klamb2, Klamb3, Ksat1, Ksat2, Ksat3, Kwp1,

* Kwp2, Kwp3, Ki1, Ki2, Ki3, kii1,

* kii2, kii3, k10a, k10b, k10c,

* Lx1,Lx2,Lx3,Lx4,Lx5,Lx6,Lx7,Lx8,Lx9,Lx10,Lx11,Lx12,

* Ly1,Ly2,Ly3,Ly4,Ly5,Ly6,Ly7,Ly8,Ly9,Ly10,Ly11,Ly12

real*8

* dtet1, tet1, tetex, hss, hs, qss, qs, qt,

* rain, raint, ftet1, ftet2, ftet3, psii1, psii2,

* psii3, psie1, psie2, psie3, EVP, eaep1, eaep2, eaep3,

* ea1, ea2, ea3, psi10, dtet2, tet2, por1, por2, por3,

* tetfc1, tetfc2, tetfc3, tets1, tets2, tets3,

* tetre1, tetre2, tetre3, teti1, teti2, teti3,

* tet10a, tet10b, tet10c, tetwp1, tetwp2, tetwp3,

* A1, B1, A2, B2, A3, B3, C, D, E, T, U, V, R,

* fsolo, hc, hp, dtet3, tet3, bc, rd, ae, pm, AT, RL,

* HEB, HSB, BL, X,

* R1, R2, R3, R4, R5, R6, R7, R8, R9, R10, R11, R12,

* He1, He2, He3, He4, He5, He6, He7, He8, He9, He10,

* He11, He12,

* far1, far2, far3, h1, h2, h3, ftetA, ftetB, ftetC,

* dtet11, dtet12, dtet13,

* dtet21, dtet22, dtet23, dtet31, dtet32, dtet33,

* S1, S2, S3, YQ, Cg, Cr, Draina, cLd, sk, CW, Sw,

* CKr, CKf, Ds

real*8

* Ds1, Ds2, Ds3, Ds4, Ds5, Ds6, Ds7, Ds8, Ds9,

* Ps1, Ps2, Ps3, Ps4, Ps5, Ps6, Ps7, Ps8, Ps9,

* Em1, Em2, Em3, Em4, Em5, Em6, Em7, Em8, Em9,

* Em10, Em11, Em12,

* Tm1, Tm2, Tm3, Tm4, Tm5, Tm6, Tm7, Tm8, Tm9,

* Tm10, Tm11, Tm12,

* Rm1, Rm2, Rm3, Rm4, Rm5, Rm6, Rm7, Rm8, Rm9,

* Rm10, Rm11, Rm12,

* Vm1, Vm2, Vm3, Vm4, Vm5, Vm6, Vm7, Vm8, Vm9,

* Vm10, Vm11, Vm12,

* Um1, Um2, Um3, Um4, Um5, Um6, Um7, Um8, Um9,

* Um10, Um11, Um12,

* timeol, time

c----Leitura do Arquivo de Dados (Input.txt)

open (8,file='Input.txt',status='unknown')

2 format(40A)

c---(1)Ano inicial do banco de dados de chuva e Ano bisexto antecedente

read (8,2) ANOSIB

read (8,*) iyear0, iyearb, N

c---(2)Área da bacia(km2), comprimento rio (km), altitudes de entrada e saída (m)

read (8,2) BACICO

read (8,*) AT, RL, HEB, HSB, BL

c---(2)Áreas para modelo distribuído (km2)

read (8,2) AREASM

read (8,*)

* Ar1, Ar2, Ar3, Ar4, Ar5, Ar6, Ar7, Ar8, Ar9, Ar10, Ar11, Ar12

c---(2)Largura da célula na direção x (Lx1,2,...12)

read (8,2) DIREXM

read (8,*)

* Lx1, Lx2, Lx3, Lx4, Lx5, Lx6, Lx7, Lx8, Lx9, Lx10, Lx11, Lx12

c---(2)Largura da célula na direção x (Cy1,2,...12)

read (8,2) DIREYM

read (8,*)

* Ly1, Ly2, Ly3, Ly4, Ly5, Ly6, Ly7, Ly8, Ly9, Ly10, Ly11, Ly12

c---(2)Comprimento de trechos de rio para modelo distribuído (km)

read (8,2) COMPMD

read (8,*)

* R1, R2, R3, R4, R5, R6, R7, R8, R9, R10, R11, R12

c---(2)Altitudes de entrada das áreas para modelo distribuído (m)

read (8,2) ALTIMD

read (8,*)

* He1, He2, He3, He4, He5, He6, He7, He8, He9, He10, He11, He12

c---(3)Tensões do solo: a 10kpa, na cap.campo (33kPa) e no ponto de murcha(1500 kpa)

read (8,2) TENSAO

read (8,*) psi10, psifc, psiwp

c---(3)Fatores de ar no solo para cálculo da porosidade efetiva (-)

read (8,2) APRISA

read (8,*) far1, far2, far3

c---(3)Fatores que definem a umidade inicial dos solos p/cálculo da porosidade efetiva

read (8,2) UMIINI

read (8,*) ftetA, ftetB, ftetC

c---(3)Fatores C e D da função da capacidade de infiltração [C.f(tet)**D] (-)

c--- f(tet) é a função da cap.infiltração em termos da umidade do solo

read (8,2) CDTETA

read (8,*) C, D

c---(3)Teores de areia dos solos nos horizontes A, B, C ou 1, 2 e 3 (%)

read (8,2) TEORAR

read (8,*) S1, S2, S3

c---(3)Teores de argila dos solos nos horizontes A, B e C ou 1, 2 e 3 (%)

read (8,2) TEORAG

read (8,*) C1, C2, C3

c---(3)Profundidades dos horizontes dos solos (cm)

read (8,2) PROFSO

read (8,*) h1, h2, h3

c---(3)Fator da função razão evap.atual/evap.potencial [eaep=f(psi)] (-)

read (8,2) FATOEA

read (8,*) fsolo

c---(4)Altura de queda da chuva, oriunda da vegetação (m)

read (8,2) ALTQVG

read (8,*) YQ

c---(4)Proporções de áreas vegetadas Cg e de rochas Cr (-) [0 <= Cg,Cr <= 1]

read (8,2) FTCVGR

read (8,*) Cg, Cr

c---(4)Drenagem da intercepção Draina (excesso sobre a cap.intercepcao) (mm)

c--- Proporcao cLd da intercepcao que cai das folhas (mm) e (-)

read (8,2) DRNCLD

read (8,*) Draina, cLd

c---(5)Evaporações médias diárias dos meses,oriundas do tanque classe A(mm/dia)

read (8,2) EVAPTQ

read (8,*)

* Em1, Em2, Em3, Em4, Em5, Em6, Em7, Em8, Em9, Em10, Em11, Em12

c---(5)Temperaturas médias diarias dos meses do ano (oC)

read (8,2) TEMPME

read (8,*)

* Tm1, Tm2, Tm3, Tm4, Tm5, Tm6, Tm7, Tm8, Tm9, Tm10, Tm11, Tm12

c---(5)Radiações médias diárias globais (incidentes) dos meses do ano (cal/cm2)

read (8,2) RADIME

read (8,*)

* Rm1, Rm2, Rm3, Rm4, Rm5, Rm6, Rm7, Rm8, Rm9, Rm10, Rm11, Rm12

c---(5)Velocidades médias diárias do vento, dos meses do ano (m/s)

read (8,2) VELOME

read (8,*)

* Vm1, Vm2, Vm3, Vm4, Vm5, Vm6, Vm7, Vm8, Vm9, Vm10, Vm11, Vm12

c---(5)Umidades relativas do ar médias diárias, dos meses do ano (%)

read (8,2) UMIDME

read (8,*)

* Um1, Um2, Um3, Um4, Um5, Um6, Um7, Um8, Um9, Um10, Um11, Um12

c---(5)Opção p/a EVP: 1(EVP=Em); 2(EVP=bc); 3(EVP=rd); 4(EVP=ae); 5(EVP=pm)

read (8,2) OPTEVA

read (8,*) MVAP

c---(6)Intervalo de tempo p/propagação inttMD (min) e fator forma Muskingum

read (8,2) MUSKIN

read (8,*) inttMD, X

c---(6)Dados do canal: n de Manning(m**(1/3)/s),Largura (m),declividade (m/m)

read (8,2) CANALM

read (8,*) sk, CW, Sw

c---(7)Fatores de erosividade da chuva e do escoamento (-)

read (8,2) FATERF

read (8,*) CKr, CKf

c---(7)Diâmetro representativo do sedimento D50 (mm)

read (8,2) DIMSED

read (8,*) Ds

Ds=Ds/1000.

c---(7)Diâmetros dos sedimentos p/modelo distribuido (mm)

read (8,2) DIMSMD

read (8,*)

* Ds1, Ds2, Ds3, Ds4, Ds5, Ds6, Ds7, Ds8, Ds9

Ds1=Ds1/1000.

Ds2=Ds2/1000.

Ds3=Ds3/1000.

Ds4=Ds4/1000.

Ds5=Ds5/1000.

Ds6=Ds6/1000.

Ds7=Ds7/1000.

Ds8=Ds8/1000.

Ds9=Ds9/1000.

c---(7)Percentual de cada Diâmetro de sedimento p/modelo distribuído (%)

read (8,2) PORCDS

read (8,*)

* Ps1, Ps2, Ps3, Ps4, Ps5, Ps6, Ps7, Ps8, Ps9

MODFIPE

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc

c c

c Programa para Calculo do escoamento e erosao do solo (MOFIPE DISTRIBUIDO)

c Codigo escrito a posteriori: iniciado 26/09/06; encerrado: em construcao c

c c

c por E.E. de Figueiredo c

c c

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc

c----------------------------------------------------------------------------------------------------------c

c Implementacoes: c

c 1-Modelos para Evaporacao Potencial (EP) c

c Blaney e Criddle c

c Aerodinamico c

c Penman c

c 2-Modelo de Laursen com nove fracoes de sedimentos (distribuido)

c 3-Modelo distribuído na fase canal c

c Métodos: SCS(convexo) e Muskingum c

c c

c OBSERVACOES c

c 1-Lê dados de EP do tanque (médias diárias dos meses do ano c

c 2-Versao interativa com 3 horizontes de solo e arquivo de dados c

c c

c----------------------------------------------------------------------------c

C---VARIÁVEIS DO SOLO: SAXTON ET AL(1986), R.&BRAKENSIEK(1989), B&COREY (1964)

c---Variáveis e dados de entrada

real Klamb1, Klamb2, Klamb3, Ksat1, Ksat2, Ksat3, Kwp1, Kwp2, Kwp3, Kfc1,

* Kfc2, Kfc3, Ki1, Ki2, Ki3, kii1, kii2, kii3, k10a, k10b, k10c, dteta1,

* dteta2, dteta3

c----Leitura de Dados

open (1,file='Data_In1.txt',status='unknown')

2 format(40A4)

c---Tensoes do solo: 10Kpa, cap.campo (33kPa) e murcha(1500 kpa)

read (1,2) TENSOES

read (1,*) psi10, psifc, psiwp

c---Fatores de aprisionamento do ar que reduzem a porosidade total para a porosidade efetiva (-)

read (1,2) APRIS_AR

read (1,*) far1, far2, far3

c---Fatores que definem a umidade inicial dos horizontes A, B e C em termos da umidade de saturacao (porosidade efetiva)

read (1,2) UMI_INI

read (1,*) ftetaA, ftetaB, ftetaC

c---Fatores C e D da funcao da capacidade de infiltracao em termos da umidade [C.f(teta)^D] (-)

read (1,2) CD_TETA

read (1,*) C, D

c---Teor de areia dos horizontes A, B e C ou horizontes 1, 2 e 3 (%)

read (1,2) TEOR_ARE

read (1,*) S1, S2, S3

c---Teor de argila dos horizontes A, B e C ou horizontes 1, 2 e 3 (%)

read (1,2) TEOR_ARG

read (1,*) C1, C2, C3

c---Evaporacoes medias diarias dos meses do ano oriundas do tanque classe A (mm/dia)

read (1,2) EVAP_TQ

read (1,*) Em1, Em2, Em3, Em4, Em5, Em6, Em7, Em8, Em9, Em10, Em11, Em12

c---Temperaturas medias diarias dos meses do ano (oC)

read (1,2) TEMP_MES

read (1,*) Tm1, Tm2, Tm3, Tm4, Tm5, Tm6, Tm7, Tm8, Tm9, Tm10, Tm11, Tm12

c---Radiacoes medias diarias globais dos meses do ano (cal/cm2)

read (1,2) RADI_MES

read (1,*) Rm1, Rm2, Rm3, Rm4, Rm5, Rm6, Rm7, Rm8, Rm9, Rm10, Rm11, Rm12

c---Velocidades medias diarias do vento, dos meses do ano (m/s)

read (1,2) VELO_MES

read (1,*) Vm1, Vm2, Vm3, Vm4, Vm5, Vm6, Vm7, Vm8, Vm9, Vm10, Vm11, Vm12

c---Umidades relativas do ar medias diaria, dos meses do ano (%)

read (1,2) UMID_MES

read (1,*) Um1, Um2, Um3, Um4, Um5, Um6, Um7, Um8, Um9, Um10, Um11, Um12

c---Opcao para o modelo de EVP: 1(EVP=Em); 2(EVP=bc); 3(EVP=rd); 4(EVP=ae); 5(EVP=pm)

read (1,2) OPCAO_EV

read (1,*) EVAP

c---Ano Inicial do banco de dados de chuva e Ano bisexto antecedente ao ano inicial

read (1,2) ANOS_IBI

read (1,*) iyear0, iyearbi

c---Profundidades dos horizontes dos solos (cm)

read (1,2) PROF_SOL

read (1,*) h1, h2, h3

c---Fator da funcao eaep=f(psi)

read (1,2) FATO_EAP

read (1,*) fsolo

c---Dados da bacia: area(km2), comprimento do rio (km), altitudes de entrada e saida (m)

read (1,2) BACI_CON

read (1,*) area, CL, AltE, AltS

c---Areas para modelo distribuido (km2)

read (1,2) AREAS_MD

read (1,*) area1, area2, area3, area4, area5, area6, area7, area8, area9, area10,

* area11, area12

c---Comprimento de trechos de rio para modelo distribuido (km)

read (1,2) COMP_MD

read (1,*) CL1, CL2, CL3, CL4, CL5, CL6, CL7, CL8, CL9, CL10, CL11, CL12

c---Altitudes de entrada das areas para modelo distribuido (m)

read (1,2) ALTI_MD

read (1,*) Alt1, Alt2, Alt3, Alt4, Alt5, Alt6, Alt7, Alt8, Alt9, Alt10, Alt11, Alt12

c---Intervalo de tempo para calculo da propagacao no modelo distribuido inttMD (min) e fator de forma de Muskingum (X)

read (1,2) INTT_MU

read (1,*) inttMD, X

c---Dados do canal: n de Manning (-), Largura (m) e declividade (-)

read (1,2) CANAL_MD

read (1,*) sk, CW, Sw

c---Altura de queda da chuva, oriunda da vegetacao (m)

read (1,2) ALTQ_VG

read (1,*) YQ

c---Fatores de cobertura (-)

read (1,2) FATC_VGR

read (1,*) Cg, Cr

c---Intercepcao: Drenagem da intercepcao (excesso em relacao a cap de intercepcao) e proporcao da intercepcao que cai das folhas (mm) e (-)

read (1,2) INTER_VG

read (1,*) Draina, cLd

c---Fatores de erosividade da chuva e do escoamento (-)

read (1,2) FATE_RF

read (1,*) CKr, CKf

c---Diametro do sedimento D50 (mm)

read (1,2) DIAM_SED

read (1,*) Ds

c---Diametros dos sedimentos (distribuidos) para modelo distribuido (mm)

read (1,2) DIAMS_MD

read (1,*) Ds1, Ds2, Ds3, Ds4, Ds5, Ds6, Ds7, Ds8, Ds9

c---Porcentagem de cada Diametro de sedimento (distribuidos) para modelo distribuido (%)

read (1,2) PORCT_DS

read (1,*) Ps1, Ps2, Ps3, Ps4, Ps5, Ps6, Ps7, Ps8, Ps9

CCC data psi10, psifc, psiwp

CCC * /10., 33.3, 1500./

CCC data far1, far2, far3, ftetaA, ftetaB, ftetaC, C, D

CCC * /0.91, 0.8, 0.8, 0.24729, 0.381, 0.3, 2.5, 0.9/

CCC data S1, C1, S2, C2, S3, C3

CCC * /50.2, 15.77, 50.2, 32.5, 50.2, 32.5/

MODFIPE ENCERRADO EM DEZ 2006

c----------------------------------------------------------------------------c

c c

c MODFIPE c

c c

c Programa para Cálculo do Escoamento e Erosão do Solo c

c iniciado: 26/09/06; encerrado: 22/12/2006 c

c c

c por E.E. de Figueiredo c

c c

c----------------------------------------------------------------------------c

c c

c Implementações feitas no modelo MOFIPE desenvolvido em EXCEL por DAVI (2004): c

c 1-Modelos para Evaporação Potencial (EP) c

c Uso de dados de evaporações médias diárias dos meses c

c Blaney e Criddle c

c Radiação c

c Aerodinâmico c

c Penman c

c 2-Modelo de Laursen com nove frações de sedimentos (distribuído) c

c 3-Modelo distribuído na fase canal c

c SCS(convexo) c

c Muskingum c

c 4-Caracteristicas da Versao Atual: c

c Com 3 horizontes de solo c

c Com 2 arquivos de leitura de dados c

c Arquivo de dados da bacia e parametros do solo c

c Arquivo de dados de chuva c

c----------------------------------------------------------------------------c

c---VARIÁVEIS DO SOLO: Saxton et al(1986), R.& Brakensiek(1989), B. & Corey (1964)

c---Declaração de Variáveis reais

real*8

* Klamb1, Klamb2, Klamb3, Ksat1, Ksat2, Ksat3, Kwp1, Kwp2, Kwp3, Kfc1,

* Kfc2, Kfc3, Ki1, Ki2, Ki3, kii1, kii2, kii3, k10a, k10b, k10c,

* Lx, Ly, Lx1, Lx2, Lx3, Lx4, Lx5, Lx6, Lx7, Lx8, Lx9,

* Lx10, Lx11, Lx12, Ly1, Ly2, Ly3, Ly4, Ly5, Ly6, Ly7, Ly8, Ly9, Ly10,

* Ly11, Ly12

real*8

* dteta1, teta1, tetex, hss, hs, qss, qs, qt, rain, raint,

* fteta1, fteta2, fteta3, psii1, psii2, psii3, psie1, psie2, psie3

* EVP, eaep1, eaep2, eaep3, ea1, ea2, ea3, psi10, teta3

* dteta2, teta2, por1, por2, por3, tetafc1, tetafc2, tetafc3,

* tetas1, tetas2, tetas3, tetawp1, tetawp2, tetawp3, tetawp, tetafc,

* tetare1, tetare2, tetare3, tetai1, tetai2, tetai3, teta10a, teta10b,

* teta10c, A1, B1, A2, B2, A3, B3, C, D, fsolo, hc, hp, dteta3, teta3,

* E, bc, rd, ae, pm, AT, RL, H_E, H_S, BL, Cx1, Cx2, Cx3, Cx4, Cx5,

* Cx6, Cx7, Cx8, Cx9, Cx10, Cx11, Cx12, Cy1, Cy2, Cy3, Cy4, Cy5, Cy6,

* Cy7, Cy8, Cy9, Cy10, Cy11, Cy12, R1, R2, R3, R4, R5, R6, R7, R8, R9,

* R10, R11, R12, He1, He2, He3, He4, He5, He6, He7, He8, He9, He10,

* He11, He12, far1, far2, far3, h1, h2, h3, ftetaA, ftetaB, ftetaC,

* ftetai1, ftetai2, ftetai3, dteta11, dteta12, dteta13, dteta21, dteta22,

* dteta23, dteta31, dteta32, dteta33,

* Ar1, Ar2, Ar3, Ar4, Ar5, Ar6, Ar7, Ar8, Ar9, Ar10, Ar11, Ar12,

* S1, S2, S3, YQ, Cg, Cr, Draina, cLd, sk, CW, Sw, CKr, CKf, Ds,

* Ds1, Ds2, Ds3, Ds4, Ds5, Ds6, Ds7, Ds8, Ds9,

* Ps1, Ps2, Ps3, Ps4, Ps5, Ps6, Ps7, Ps8, Ps9,

* Em1, Em2, Em3, Em4, Em5, Em6, Em7, Em8, Em9, Em10, Em11, Em12,

* Tm1, Tm2, Tm3, Tm4, Tm5, Tm6, Tm7, Tm8, Tm9, Tm10, Tm11, Tm12,

* Rm1, Rm2, Rm3, Rm4, Rm5, Rm6, Rm7, Rm8, Rm9, Rm10, Rm11, Rm12,

* Vm1, Vm2, Vm3, Vm4, Vm5, Vm6, Vm7, Vm8, Vm9, Vm10, Vm11, Vm12,

* Um1, Um2, Um3, Um4, Um5, Um6, Um7, Um8, Um9, Um10, Um11, Um12

c---Matriz coluna para arquivos de saída do fluxo

Dimension qt(20000), qs(20000), qss(20000), qtCo(20000),

* qtMu1(20000), qtMu2(20000), qtMu3(20000), qtMu4(20000),

* flowtot(20000), flowtot9(20000),

* flowtot1(20000), flowtot2(20000), flowtot3(20000), flowtot4(20000),

* flowtot5(20000), flowtot6(20000), flowtot7(20000), flowtot8(20000)

c----Leitura do Arquivo de Dados (Data_input.txt)

open (1,file='Data_Input.txt',status='unknown')

2 format(40A4)

c---(1)Ano inicial do banco de dados de chuva e Ano bisexto antecedente

read (1,2) ANOS_IBI

read (1,*) iyear0, iyearbi, N

c---(2)Área da bacia(km2), comprimento rio (km), altitudes de entrada e saída (m)

read (1,2) BACI_CON

read (1,*) AT, RL, H_E, H_S, BL

c---(2)Áreas para modelo distribuído (km2)

read (1,2) AREAS_MD

read (1,*) Ar1, Ar2, Ar3, Ar4, Ar5, Ar6, Ar7, Ar8, Ar9, Ar10, Ar11, Ar12

c---(2)Largura da célula na direção x (Cx1,2,...12)

read (1,2) DIREX_MD

read (1,*) Cx1, Cx2, Cx3, Cx4, Cx5, Cx6, Cx7, Cx8, Cx9, Cx10, Cx11, Cx12

c---(2)Largura da célula na direção x (Cy1,2,...12)

read (1,2) DIREY_MD

read (1,*) Cy1, Cy2, Cy3, Cy4, Cy5, Cy6, Cy7, Cy8, Cy9, Cy10, Cy11, Cy12

c---(2)Comprimento de trechos de rio para modelo distribuído (km)

read (1,2) COMP_MD

read (1,*) R1, R2, R3, R4, R5, R6, R7, R8, R9, R10, R11, R12

c---(2)Altitudes de entrada das áreas para modelo distribuído (m)

read (1,2) ALTI_MD

read (1,*) He1, He2, He3, He4, He5, He6, He7, He8, He9, He10, He11, He12

c---(3)Tensões do solo: a 10kpa, na cap.campo (33kPa) e no ponto de murcha(1500 kpa)

read (1,2) TENSOES

read (1,*) psi10, psifc, psiwp

c---(3)Fatores de ar no solo para cálculo da porosidade efetiva (-)

read (1,2) APRIS_AR

read (1,*) far1, far2, far3

c---(3)Fatores que definem a umidade inicial dos solos p/cálculo da porosidade efetiva

read (1,2) UMI_INI

read (1,*) ftetaA, ftetaB, ftetaC

c---(3)Fatores C e D da função da capacidade de infiltração [C.f(teta)**D] (-)

c--- f(teta) é a função da cap.infiltração em termos da umidade do solo

read (1,2) CD_TETA

read (1,*) C, D

c---(3)Teores de areia dos solos nos horizontes A, B, C ou 1, 2 e 3 (%)

read (1,2) TEOR_ARE

read (1,*) S1, S2, S3

c---(3)Teores de argila dos solos nos horizontes A, B e C ou 1, 2 e 3 (%)

read (1,2) TEOR_ARG

read (1,*) C1, C2, C3

c---(3)Profundidades dos horizontes dos solos (cm)

read (1,2) PROF_SOL

read (1,*) h1, h2, h3

c---(3)Fator da função razão evap.atual/evap.potencial [eaep=f(psi)] (-)

read (1,2) FATO_EAP

read (1,*) fsolo

c---(4)Altura de queda da chuva, oriunda da vegetação (m)

read (1,2) ALTQ_VG

read (1,*) YQ

c---(4)Proporções de áreas vegetadas Cg e de rochas Cr (-) [0 <= Cg,Cr <= 1]

read (1,2) FATC_VGR

read (1,*) Cg, Cr

c---(4)Drenagem da intercepção Draina (excesso sobre a cap.intercepcao) (mm)

c--- Proporcao cLd da intercepcao que cai das folhas (mm) e (-)

read (1,2) INTER_VG

read (1,*) Draina, cLd

c---(5)Evaporações médias diárias dos meses,oriundas do tanque classe A(mm/dia)

read (1,2) EVAP_TQ

read (1,*) Em1, Em2, Em3, Em4, Em5, Em6, Em7, Em8, Em9, Em10, Em11, Em12

c---(5)Temperaturas médias diarias dos meses do ano (oC)

read (1,2) TEMP_MES

read (1,*) Tm1, Tm2, Tm3, Tm4, Tm5, Tm6, Tm7, Tm8, Tm9, Tm10, Tm11, Tm12

c---(5)Radiações médias diárias globais (incidentes) dos meses do ano (cal/cm2)

read (1,2) RADI_MES

read (1,*) Rm1, Rm2, Rm3, Rm4, Rm5, Rm6, Rm7, Rm8, Rm9, Rm10, Rm11, Rm12

c---(5)Velocidades médias diárias do vento, dos meses do ano (m/s)

read (1,2) VELO_MES

read (1,*) Vm1, Vm2, Vm3, Vm4, Vm5, Vm6, Vm7, Vm8, Vm9, Vm10, Vm11, Vm12

c---(5)Umidades relativas do ar médias diárias, dos meses do ano (%)

read (1,2) UMID_MES

read (1,*) Um1, Um2, Um3, Um4, Um5, Um6, Um7, Um8, Um9, Um10, Um11, Um12

c---(5)Opção para escolha do modelo de Evapotranspiração EVP:

c--- 1(EVP=Em); 2(EVP=bc); 3(EVP=rd); 4(EVP=ae); 5(EVP=pm)

read (1,2) OPCAO_EV

read (1,*) MVAP

c---(6)Intervalo de tempo p/ cálculo da propagação do escoamento inttMD (min)

c--- Fator de forma de Muskingum (X)

read (1,2) INTT_MU

read (1,*) inttMD, X

c---(6)Dados do canal: n de Manning(m**(1/3)/s),Largura (m),declividade (m/m)

read (1,2) CANAL_MD

read (1,*) sk, CW, Sw

c---(7)Fatores de erosividade da chuva e do escoamento (-)

read (1,2) FATE_RF

read (1,*) CKr, CKf

c---(7)Diâmetro representativo do sedimento D50 (mm)

read (1,2) DIAM_SED

read (1,*) Ds

Ds=Ds/1000.

c---(7)Diâmetros dos sedimentos p/modelo distribuido (mm)

read (1,2) DIAMS_MD

read (1,*) Ds1, Ds2, Ds3, Ds4, Ds5, Ds6, Ds7, Ds8, Ds9

Ds1=Ds1/1000.

Ds2=Ds2/1000.

Ds3=Ds3/1000.

Ds4=Ds4/1000.

Ds5=Ds5/1000.

Ds6=Ds6/1000.

Ds7=Ds7/1000.

Ds8=Ds8/1000.

Ds9=Ds9/1000.

c---(7)Percentual de cada Diâmetro de sedimento p/modelo distribuído (%)

read (1,2) PORCT_DS

read (1,*) Ps1, Ps2, Ps3, Ps4, Ps5, Ps6, Ps7, Ps8, Ps9

c---Parâmetros de Saxton et al

A1 = 100.*EXP(-4.396-0.0715*C1-(0.000488)*S1**2.-

* (0.00004285)*(S1**2.)*C1)

B1 = -3.14-(0.00222)*C1**2.-(0.00003484)*(S1**2.)*C1

A2 = 100.*EXP(-4.396-0.0715*C2-(0.000488)*S2**2.-

*(0.00004285)*(S2**2.)*C2)

B2 = -3.14-(0.00222)*C2**2.-(0.00003484)*(S2**2.)*C2

A3 = 100.*EXP(-4.396-0.0715*C3-(0.000488)*S3**2.-

*(0.00004285)*(S3**2.)*C3)

B3 = -3.14-(0.00222)*C3**2.-(0.00003484)*(S3**2.)*C3

c---Parâmetro de Brooks & Corey

Klamb1 = 1/B1

Klamb2 = 1/B2

Klamb3 = 1/B3

c---Variáveis de Saxton et al

c---Porosidades totais dos solos (m3/m3)

por1 = 0.332-0.0007251*S1+0.1276*ALOG10(C1)

por2 = 0.332-0.0007251*S2+0.1276*ALOG10(C2)

por3 = 0.332-0.0007251*S3+0.1276*ALOG10(C3)

c---Tensões dos solos (kPa)

psie1 = 100.*(-0.108+0.341*por1)

psie2 = 100.*(-0.108+0.341*por2)

psie3 = 100.*(-0.108+0.341*por3)

c---Umidades (m3/m3) correspondentes a saturação dos solos (porosidade efetiva)

c---far1,2,3 são fatores (-) para o aprisionamento de ar (0.8 < far1,2,3 < 1)

tetas1 = far1*por1

tetas2 = far2*por2

tetas3 = far3*por3

c---Umidades (m3/m3) dos solos correspondentes a capacidade de campo (33.3 kPa)

tetafc1 = (psifc/A1)**(1/B1)

tetafc2 = (psifc/A2)**(1/B2)

tetafc3 = (psifc/A3)**(1/B3)

c---Umidades (m3/m3) dos solos correspondentes ao ponto de murcha (1500 kPa)

tetawp1 = (psiwp/A1)**(1/B1)

tetawp2 = (psiwp/A2)**(1/B2)

tetawp3 = (psiwp/A3)**(1/B3)

c---Umidade (m3/m3) correspondentes a 10kPa (solo considerado saturado)

teta10a = (psi10/A1)**(1/B1)

teta10b = (psi10/A2)**(1/B2)

teta10c = (psi10/A3)**(1/B3)

c---Umidades residuais (Rawls and Brakensiek) dos solos (m3/m3)

tetare1 = -0.0182482+0.00087269*S1+0.00513488*C1+

* 0.02939286*por1-0.00015395*C1**2.-0.0010827*S1*por1-

* 0.00018233*C1**2.*por1**2.+0.00030703*C1**2.*por1-

* 0.0023584*por1**2.*C1

tetare2 = -0.0182482+0.00087269*S2+0.00513488*C2+

* 0.02939286*por2-0.00015395*C2**2.-0.0010827*S2*por2-

* 0.00018233*C2**2.*por2**2.+0.00030703*C2**2.*por2-

* 0.0023584*por2**2.*C2

tetare3 = -0.0182482+0.00087269*S3+0.00513488*C3+

* 0.02939286*por3-0.00015395*C3**2.-0.0010827*S3*por3-

* 0.00018233*C3**2.*por3**2.+0.00030703*C3**2.*por3-

* 0.0023584*por3**2.*C3

c---Umidades iniciais dos solos (= uma parcela da porosidade total)

c---ftetai1,2,3 são os fatores definidores (0 < ftetai1,2,3 < 1)

tetai1 = ftetaA*por1

tetai2 = ftetaB*por2

tetai3 = ftetaC*por3

MOSESS 2D

c----------------------------------------------------------------------------c

c c

c MOSESS_2D c

c c

c 2 Dimensional model for runoff and soil erosion prediction c

c at the catchment scale: Oct 2008-Sept 2009 c

c c

c Author: E.E. de Figueiredo c

c UAEC/UFCG/CTRN c

c Brazil c

c c

c Program coded during Post-doctoral research in c

c the University of Sheffield/UK with the c

c support of CAPES/BRASIL c

c c

c----------------------------------------------------------------------------c

c c

c Program code in Fortran g77 which includes two components: c

c c

c A) FLOW COMPONENT c

c c

c -Rainfall intensities based on the rainfall event duration c

c -Potential evapotranspiration between events of rainfall based on c

c the following alternative evaporation methods: c

c .Pan evaporation data c

c .Blaney e Criddle equation c

c .Radiation c

c .Aerodynamic c

c .Penman equation for potential evaporation c

c -Surface infiltration can be calculated by: c

c .Potential function c

c .Horton c

c -Percolation is based on the hydraulic conductivities of the soils c

c layers (limited to 3 soils layers) c

c -Direct runoff is a result of excess rainfall over infiltration c

c capacity, saturation of the surface layer, or both c

c -Sub-surface flow takes place after the surface soil layer is c

c saturated, and is based on the soil hydraulic conductivity c

c -Base flow (not modelled) c

c -Surface runoff routing can be calculated by: c

c .Muskingum equation c

c .Convex method of the SCS c

c Obs: Both methods takes the runoff travel time as the time of c

c concentration, which can be calculated alternatively by: c

c .Kirpich c

c .Kerby c

c c

c B) SOIL EROSION COMPONENT c

c c

c -Detachment by raindrop is based on the rainfall squared moment c

c -Detachment by runoff is based on the critical shear stress by c

c Shields c

c -Transport/Deposition is based on the flow transport capacity by: c

c .Engelund-Hansen (total load) with D50 c

c .Yalin (bed load) with D50 c

c .Laursen (bed and/or total load) with D50 c

c .Laursen (bed and/or total load) with Sed Size Distributionc

c c

c NOTES: c

c -Catchment is divided into sub-catchments related to the main c

c river (limited to 20) c

c -channels reaches areas are not taken into account c

c -Reaches length and height are set for the routing procedures c

c -Reaches flow velocity is based on Manning's equation that requires c

c roughness coefficient, channel width, and slope c

c -Each reach define a sub-catchment for which data are assigned c

c ex.: rainfall, areas, relief, soils parameters, etc c

c -Application of the model must be for relatively uniform basins c

c----------------------------------------------------------------------------c

PROGRAM MOSESS2D

c---Matrix variables

c---1,2,3,,,,20(sub-basin number or reach number)

c---qt,qs,qsb,qss (total,surface,sub-surface,and excess saturation runoff mm)........calculated

c---qtMu1,,,20,flowtot1,,,,20 (propagated runoff at reache1,,,20)....................calculated

Dimension

* qt(20000), qs(20000),

* qsb(20000), qss(20000),

* qtMu1(20000), qtMu2(20000),

* qtMu3(20000), qtMu4(20000),

* qtMu5(20000), qtMu6(20000),

* qtMu7(20000), qtMu8(20000),

* qtMu9(20000), qtMu10(20000),

* qtMu11(20000), qtMu12(20000),

* qtMu13(20000), qtMu14(20000),

* qtMu15(20000), qtMu16(20000),

* qtMu17(20000), qtMu18(20000),

* qtMu19(20000), qtMu20(20000),

* flowtot1(20000), flowtot2(20000),

* flowtot3(20000), flowtot4(20000),

* flowtot5(20000), flowtot6(20000),

* flowtot7(20000), flowtot8(20000),

* flowtot9(20000), flowtot10(20000),

* flowtot11(20000), flowtot12(20000),

* flowtot13(20000), flowtot14(20000),

* flowtot15(20000), flowtot16(20000),

* flowtot17(20000), flowtot18(20000),

* flowtot19(20000), flowtot20(20000),

* hpt(20000), set(20000), sload(20000),

* qm(20000), ftet(20000), KtetA(20000),

* KtetB(20000), KtetC(20000), psiA(20000),

* tday(20000), tmon(20000), tyea(20000),

* tche(20000), tetatA(20000), tetatB(20000),

* tetatC(20000), qgw(20000)

c---Real variables

c---? means horizons(A, B, C or soils layers A, B or C)

c---1,2,3,,,,20(sub-basin number or reach number)

c---KlamA?(exponent of B&Corey function) ............................................calculated

c---Ksat?(saturated hydraulic conductivity - mm/h)...................................calculated

c---ki?=kwp?(hydraulic conductivity at wilting point - mm/h).........................calculated

c---kii?(hydraulic conductivity as a function of soil moisture - mm/h)...............calculated

c---k10?(hydraulic conductivity at a tension of 10 kPa - mm/h).......................calculated

c---ktet,kt1,,,20(infiltration parameters[f=Ct(tet)**Dt f=FC+(FO-FC)exp(-ktet.tetA)].input data

c---Real variables

real*8

* KlamA, KlamB, KlamC,

* KlamA1, KlamA2, KlamA3, KlamA4, KlamA5, KlamA6, KlamA7,

* KlamA8, KlamA9, KlamA10, KlamA11, KlamA12, KlamA13,KlamA14,

* KlamA15, KlamA16,KlamA17, KlamA18, KlamA19, KlamA20,

* KlamB1, KlamB2, KlamB3, KlamB4, KlamB5, KlamB6, KlamB7,

* KlamB8, KlamB9, KlamB10, KlamB11, KlamB12, KlamB13,KlamB14,

* KlamB15, KlamB16,KlamB17, KlamB18, KlamB19, KlamB20,

* KlamC1, KlamC2, KlamC3, KlamC4, KlamC5, KlamC6, KlamC7,

* KlamC8, KlamC9, KlamC10, KlamC11, KlamC12, KlamC13,KlamC14,

* KlamC15, KlamC16,KlamC17, KlamC18, KlamC19, KlamC20,

* KwpA, KwpB, KwpC,

* KwpA1, KwpA2, KwpA3, KwpA4, KwpA5, KwpA6, KwpA7,

* KwpA8, KwpA9, KwpA10, KwpA11, KwpA12, KwpA13,KwpA14,

* KwpA15, KwpA16,KwpA17, KwpA18, KwpA19, KwpA20,

* KwpB1, KwpB2, KwpB3, KwpB4, KwpB5, KwpB6, KwpB7,

* KwpB8, KwpB9, KwpB10, KwpB11, KwpB12, KwpB13,KwpB14,

* KwpB15, KwpB16,KwpB17, KwpB18, KwpB19, KwpB20,

* KwpC1, KwpC2, KwpC3, KwpC4, KwpC5, KwpC6, KwpC7,

* KwpC8, KwpC9, KwpC10, KwpC11, KwpC12, KwpC13,KwpC14,

* KwpC15, KwpC16,KwpC17, KwpC18, KwpC19, KwpC20,

* KiiA, KiiB, KiiC,

* KiiA1, KiiA2, KiiA3, KiiA4, KiiA5, KiiA6, KiiA7,

* KiiA8, KiiA9, KiiA10, KiiA11, KiiA12, KiiA13,KiiA14,

* KiiA15, KiiA16,KiiA17, KiiA18, KiiA19, KiiA20,

* KiiB1, KiiB2, KiiB3, KiiB4, KiiB5, KiiB6, KiiB7,

* KiiB8, KiiB9, KiiB10, KiiB11, KiiB12, KiiB13,KiiB14,

* KiiB15, KiiB16,KiiB17, KiiB18, KiiB19, KiiB20,

* KiiC1, KiiC2, KiiC3, KiiC4, KiiC5, KiiC6, KiiC7,

* KiiC8, KiiC9, KiiC10, KiiC11, KiiC12, KiiC13,KiiC14,

* KiiC15, KiiC16,KiiC17, KiiC18, KiiC19, KiiC20,

* KtetAm, KtetBm, KtetCm,Kmed,

* kt1, kt2, kt3, kt4, kt5, kt6, kt7, kt8, kt9, kt10,

* kt11, kt12, kt13, kt14, kt15, kt16, kt17, kt18, kt19, kt20

c---Other real variables

c---? means horizons A, B or C (or soils layers A, B or C)

c---1,2,3,,,,20(sub-basin number or reach number)

c---AT(area of catchment - km2)......................................................input data

c---RL(main river lenght - km).......................................................input data

c---HEB,HSB(upstream and downstream catchment elevations - m)........................input data

c---AR1,,,20(sub-areas 1,,,20 - km2).................................................input data

c---RL1,,,20(reaches lenght - km)....................................................input data

c---RH1,,,20(sub-basin elevations - m)...............................................input data

c---rain1,,,20(input data of rainfall for sub-basins 1,,,20 - mm)....................input data

c---fr1,,,20(rainfall,slope,routing distibution factor)..............................calculated

c---fch1,,,20(rainfall and infiltration distribution factor).........................calculated

c---fs1,,,20(slope distribution factor)..............................................calculated

c---fk1,,,20(routing distribution factor)............................................calculated

c---tc,tc1,,,tc20(time of concentration - min).......................................calculated

c---X(shape parameter of Muskingum)..................................................input data

c---BD?,BD1,,,20(soil bulk density - g/cm3)..........................................input data

c---por?,por?1,,,20(soil porosity=1 - BD/SPD; SPD=2.65 g/cm3)........................calculated

c---far?far?1,,,20(air entrapment factor)............................................input data

c---tets?tets?1,,,20(moisture content at saturation m3/m3)...........................calculated

c---fumin?,fumin?1,,,20(initial moisture content calculation factor).................input data

c---teti?,teti?1,,,20(initial moisture content m3/m3)................................calculated

c---Ct,Ct1,,,20;Dt,Dt1,,,Dt20 (infiltration capacity parameters f = Ct.tetA.**Dt)....input data

c---FO1,,,FO20 ; FC1,,,FC20 (intiltration capacity parameters - mm/h)................input data

c---ftet?(infiltration capacity mm/h)................................................calculated

c---fte0,ftetF,ktet(infiltration parameters of f=fC+(fo-fc)exp(-k.tetA) - mm/h)......input data

c---ftet?1,,,,ftet20(infiltration capacity distribution - mm/h)......................calculated

c---S?(sand percentage - %)..........................................................input data

c---S?1,,,,,S?20 (sand percentage distributed - %)...................................input data

c---C?(clay percentage - %)..........................................................input data

c---C?1,,,,,C?20 (clay percentage distributed - %)...................................input data

c---BB?(parameters of the Rawls et al equation for Ksat?)............................calculated

c---Ksat?(parameters of the Rawls et al equation for Ksat?)..........................calculated

c---Ksat?1,,,20(parameters of the Rawls et al equation for Ksat? distributed)........calculated

c---h?(soils layers depth - cm)......................................................input data

c---h?1,,,h?20 (soils layers depth distributed - cm).................................input data

c---A?,B?(Saxton's parameters).......................................................calculated

c---Em1,,,Em12(mean daily evaporation for months 1,,,12 - mm)........................input data

c---Rm1,,,Rm12(mean daily radiation for months 1,,,12 - cal/cm2).....................input data

c---Tm1,,,Tm12(mean daily temperature for months 1,,,12 - oC)........................input data

c---Vm1,,,Vm12(mean daily wind velocity for months 1,,,12 - m/s).....................input data

c---Um1,,,Um12(mean daily moisture content of the air for months 1,,,12 - %).........input data

c---MEVAPT(method to calculate potential evaporation - mm/day: see as follows).......calculated

c---fsolo(factor for the relation eaep?=f(soil tension)).............................input data

c---bcd(Blanney Criddle method for evaporation - mm).................................calculated

c---rd(radiation method for evaporation - cal/cm2)...................................calculated

c---ae(aerodynamic method for evaporation - mm)......................................calculated

c---pm(Penman method for evaporation - mm)...........................................calculated

c---E(mean daily evaporation - mm)...................................................calculated

c---T(mean daily temperature - oC)...................................................calculated

c---R(mean daily radiation - cal/cm2)................................................calculated

c---V(mean daily wind velocity - m/s)................................................calculated

c---U(mean daily moisture content of the air - %)....................................calculated

c---sk(Manning's coefficient)........................................................input data

c---CW(channel width - m)............................................................input data

c---Sw(channel slope - fraction).....................................................input data

c---sk?1,,,20(Manning's coefficient distributed......................................input data

c---CW?1,,,20(channel width distributed - m).........................................input data

c---Sw?1,,,20(channel slope distributed - fraction)..................................input data

c---Cg,Cr(canopy and rock cover proportions).........................................input data

c---Cg?1,,,20,Cr?1,,,20(canopy and rock cover proportions distributed)...............input data

c---Draina(drainage proportion from the canopy)......................................input data

c---cLd(proportion of rain that fall as leaf drip)...................................input data

c---CKr,CKf(rainfall and runoff erosivity coefficients)..............................input data

c---CKr?1,,,20,CKf?1,,,200(rainfall and runoff erosivity coefficients distributed)...input data

c---Ds,Ds1,,,Ds9(sediment size diameters - mm).......................................input data

c---Ps,Ps1,,,Ps9(sediment percentages % for sediment sizes Ds,Ds1,,,Ds20 mm).........input data

c---YQ,YV(vegetation heights - m)....................................................fixed data

c---qtMu1,,,20,flowtot1,,,,20 (propagated runoff at reache1,,,20)....................calculated

c---AA,BA,AB,BB,AC,BC (Saxton et al or B-Corey parameters) = f(SA,SB,SC,CA,CB,CC)....calculated

c---hcc(channel water depth - m).....................................................calculated

c---hp(overland flow depth - m)......................................................calculated

c---BB?(parameters of the Rawls et al equation for Ksat?)............................calculated

c---dtet1?;dtet2?;dtet3?(moisture content variation m3/m3)...........................calculated

c---dtet?,tet?(moisture content variation and actual moisture content m3/m3).........calculated

c---tetex? (excess of soil saturation m3/m3).........................................calculated

c---hss(lamina from the excess of soil saturation - mm)..............................calculated

c---hs(lamina from the excess of rainfall intensity over infiltration capacity mm)...calculated

c---qs,qss,qt,raint (runoff laminas, rainfall mm)....................................calculated

c---psii?(soil tension at the initial soil moisture - kPa)...........................calculated

c---psie?(soil tension at the air entry - kPa).......................................calculated

c---eaep?(ration of actual to potential evapotranspiration)..........................calculated

c---ea?(actual evapotranspirarion - mm)..............................................calculated

c---EVP (evaporation mm).............................................................calculated

c---psi10(tension of 10 kPa).........................................................fixed data

c---tetfc?(moisture content at field capacity m3/m3).................................calculated

c---tetre?(residual moisture content m3/m3)..........................................calculated

c---tet10?(moisture content at soil tension of 10 kPa - m3/m3).......................calculated

c---tetwp?(moisture content at wilting point m3/m3)..................................calculated

c Double Precision

real*8

* AT, RL, HEB, HSB,

* rain1, rain2, rain3, rain4, rain5, rain6, rain7, rain8,

* rain9, rain10, rain11, rain12, rain13, rain14, rain15, rain16,

* rain17,rain18, rain19, rain20, rain, raint, rainM,

* fr1, fr2, fr3, fr4, fr5, fr6, fr7, fr8, fr9, fr10,

* fr11, fr12, fr13, fr14, fr15, fr16, fr17, fr18, fr19, fr20,

* fch1, fch2, fch3, fch4, fch5, fch6, fch7, fch8, fch9,

* fch10,fch11,fch12,fch13,fch14,fch15,fch16,fch17,fch18,

* fch19,fch20,

* fs1, fs2, fs3, fs4, fs5, fs6, fs7, fs8, fs9, fs10,

* fs11, fs12, fs13, fs14, fs15, fs16, fs17, fs18, fs19, fs20,

* fk1, fk2, fk3, fk4, fk5, fk6, fk7, fk8, fk9, fk10,

* fk11, fk12, fk13, fk14, fk15, fk16, fk17, fk18, fk19, fk20,

* AR1, AR2, AR3, AR4, AR5, AR6, AR7, AR8, AR9, AR10,

* AR11, AR12, AR13, AR14, AR15, AR16, AR17, AR18, AR19, AR20,

* RL1, RL2, RL3, RL4, RL5, RL6, RL7, RL8, RL9, RL10,

* RL11, RL12, RL13, RL14, RL15, RL16, RL17, RL18, RL19, RL20,

* RH1, RH2, RH3, RH4, RH5, RH6, RH7, RH8, RH9, RH10,

* RH11, RH12, RH13, RH14, RH15, RH16, RH17, RH18, RH19, RH20,

* tc, X,

* tc1, tc2, tc3, tc4, tc5, tc6, tc7, tc8, tc9, tc10,

* tc11, tc12, tc13, tc14, tc15, tc16, tc17, tc18, tc19, tc20,

* BDA, BDB, BDC,

* BDA1, BDA2, BDA3, BDA4, BDA5, BDA6, BDA7, BDA8, BDA9,

* BDA10,BDA11,BDA12,BDA13,BDA14,BDA15,BDA16,BDA17,BDA18,

* BDA19,BDA20,

* BDB1, BDB2, BDB3, BDB4, BDB5, BDB6, BDB7, BDB8, BDB9,

* BDB10,BDB11,BDB12,BDB13,BDB14,BDB15,BDB16,BDB17,BDB18,

* BDB19,BDB20,

* BDC1, BDC2, BDC3, BDC4, BDC5, BDC6, BDC7, BDC8, BDC9,

* BDC10,BDC11,BDC12,BDC13,BDC14,BDC15,BDC16,BDC17,BDC18,

* BDC19,BDC20,

* porA, porB, porC,

* porA1, porA2, porA3, porA4, porA5, porA6, porA7, porA8, porA9,

* porA10,porA11,porA12,porA13,porA14,porA15,porA16,porA17,porA18,

* porA19,porA20,

* porB1, porB2, porB3, porB4, porB5, porB6, porB7, porB8, porB9,

* porB10,porB11,porB12,porB13,porB14,porB15,porB16,porB17,porB18,

* porB19,porB20,

* porC1, porC2, porC3, porC4, porC5, porC6, porC7, porC8, porC9,

* porC10,porC11,porC12,porC13,porC14,porC15,porC16,porC17,porC18,

* porC19,porC20,

* farA, farB, farC,

* farA1, farA2, farA3, farA4, farA5, farA6, farA7, farA8, farA9,

* farA10,farA11,farA12,farA13,farA14,farA15,farA16,farA17,farA18,

* farA19,farA20,

* farB1, farB2, farB3, farB4, farB5, farB6, farB7, farB8, farB9,

* farB10,farB11,farB12,farB13,farB14,farB15,farB16,farB17,farB18,

* farB19,farB20,

* farC1, farC2, farC3, farC4, farC5, farC6, farC7, farC8, farC9,

* farC10,farC11,farC12,farC13,farC14,farC15,farC16,farC17,farC18,

* farC19,farC20,

* tetsA, tetsB, tetsC,

* tetsA1, tetsA2, tetsA3, tetsA4, tetsA5, tetsA6, tetsA7,

* tetsA8, tetsA9, tetsA10, tetsA11, tetsA12, tetsA13,tetsA14,

* tetsA15, tetsA16,tetsA17, tetsA18, tetsA19, tetsA20,

* tetsB1, tetsB2, tetsB3, tetsB4, tetsB5, tetsB6, tetsB7,

* tetsB8, tetsB9, tetsB10, tetsB11, tetsB12, tetsB13,tetsB14,

* tetsB15, tetsB16,tetsB17, tetsB18, tetsB19, tetsB20,

* tetsC1, tetsC2, tetsC3, tetsC4, tetsC5, tetsC6, tetsC7,

* tetsC8, tetsC9, tetsC10, tetsC11, tetsC12, tetsC13,tetsC14,

* tetsC15, tetsC16,tetsC17, tetsC18, tetsC19, tetsC20,

* fuminA, fuminB, fuminC,

* fuminA1, fuminA2, fuminA3, fuminA4, fuminA5, fuminA6,

* fuminA7, fuminA8, fuminA9, fuminA10, fuminA11, fuminA12,

* fuminA13, fuminA14, fuminA15, fuminA16, fuminA17, fuminA18,

* fuminA19, fuminA20,

* fuminB1, fuminB2, fuminB3, fuminB4, fuminB5, fuminB6,

* fuminB7, fuminB8, fuminB9, fuminB10, fuminB11, fuminB12,

* fuminB13, fuminB14, fuminB15, fuminB16, fuminB17, fuminB18,

* fuminB19, fuminB20,

* fuminC1, fuminC2, fuminC3, fuminC4, fuminC5, fuminC6,

* fuminC7, fuminC8, fuminC9, fuminC10, fuminC11, fuminC12,

* fuminC13, fuminC14, fuminC15, fuminC16, fuminC17, fuminC18,

* fuminC19, fuminC20,

* tetiA, tetiB, tetiC,

* tetiA1, tetiA2, tetiA3, tetiA4, tetiA5, tetiA6, tetiA7,

* tetiA8, tetiA9, tetiA10, tetiA11, tetiA12, tetiA13, tetiA14,

* tetiA15, tetiA16, tetiA17, tetiA18, tetiA19, tetiA20,

* tetiB1, tetiB2, tetiB3, tetiB4, tetiB5, tetiB6, tetiB7,

* tetiB8, tetiB9, tetiB10, tetiB11, tetiB12, tetiB13, tetiB14,

* tetiB15, tetiB16, tetiB17, tetiB18, tetiB19, tetiB20,

* tetiC1, tetiC2, tetiC3, tetiC4, tetiC5, tetiC6, tetiC7,

* tetiC8, tetiC9, tetiC10, tetiC11, tetiC12, tetiC13, tetiC14,

* tetiC15, tetiC16, tetiC17, tetiC18, tetiC19, tetiC20,

* Ct, Dt,

* Ct1, Ct2, Ct3, Ct4, Ct5, Ct6, Ct7, Ct8, Ct9, Ct10,

* Ct11, Ct12, Ct13, Ct14, Ct15, Ct16, Ct17, Ct18, Ct19, Ct20,

* Dt1, Dt2, Dt3, Dt4, Dt5, Dt6, Dt7, Dt8, Dt9, Dt10,

* Dt11, Dt12, Dt13, Dt14, Dt15, Dt16, Dt17, Dt18, Dt19, Dt20,

* FO1, FO2, FO3, FO4, FO5, FO6, FO7, FO8, FO9, FO10,

* FO11, FO12, FO13, FO14, FO15, FO16, FO17, FO18, FO19, FO20,

* FC1, FC2, FC3, FC4, FC5, FC6, FC7, FC8, FC9, FC10,

* FC11, FC12, FC13, FC14, FC15, FC16, FC17, FC18, FC19, FC20,

* ftetA, PHI, ftetm,

* ftet0, ftetF, ktet,

* ftetA1, ftetA2, ftetA3, ftetA4, ftetA5, ftetA6, ftetA7,

* ftetA8, ftetA9, ftetA10, ftetA11, ftetA12, ftetA13, ftetA14,

* ftetA15, ftetA16, ftetA17, ftetA18, ftetA19, ftetA20,

* SA, SB, SC,

* SA1, SA2, SA3, SA4, SA5, SA6, SA7, SA8, SA9, SA10,

* SA11, SA12, SA13, SA14, SA15, SA16, SA17, SA18, SA19, SA20,

* SB1, SB2, SB3, SB4, SB5, SB6, SB7, SB8, SB9, SB10,

* SB11, SB12, SB13, SB14, SB15, SB16, SB17, SB18, SB19, SB20,

* SC1, SC2, SC3, SC4, SC5, SC6, SC7, SC8, SC9, SC10,

* SC11, SC12, SC13, SC14, SC15, SC16, SC17, SC18, SC19, SC20,

* CA, CB, CC,

* CA1, CA2, CA3, CA4, CA5, CA6, CA7, CA8, CA9, CA10,

* CA11, CA12, CA13, CA14, CA15, CA16, CA17, CA18, CA19, CA20,

* CB1, CB2, CB3, CB4, CB5, CB6, CB7, CB8, CB9, CB10,

* CB11, CB12, CB13, CB14, CB15, CB16, CB17, CB18, CB19, CB20,

* CC1, CC2, CC3, CC4, CC5, CC6, CC7, CC8, CC9, CC10,

* CC11, CC12, CC13, CC14, CC15, CC16, CC17, CC18, CC19, CC20,

* KsatA, KsatB, KsatC,

* KsatA1, KsatA2, KsatA3, KsatA4, KsatA5, KsatA6, KsatA7,

* KsatA8, KsatA9, KsatA10, KsatA11, KsatA12, KsatA13,KsatA14,

* KsatA15, KsatA16,KsatA17, KsatA18, KsatA19, KsatA20,

* KsatB1, KsatB2, KsatB3, KsatB4, KsatB5, KsatB6, KsatB7,

* KsatB8, KsatB9, KsatB10, KsatB11, KsatB12, KsatB13,KsatB14,

* KsatB15, KsatB16,KsatB17, KsatB18, KsatB19, KsatB20,

* KsatC1, KsatC2, KsatC3, KsatC4, KsatC5, KsatC6, KsatC7,

* KsatC8, KsatC9, KsatC10, KsatC11, KsatC12, KsatC13,KsatC14,

* KsatC15, KsatC16,KsatC17, KsatC18, KsatC19, KsatC20,

* AvA, AvB, AvC,

* hA, hB, hC,

* hA1, hA2, hA3, hA4, hA5, hA6, hA7, hA8, hA9, hA10,

* hA11, hA12, hA13, hA14, hA15, hA16, hA17, hA18, hA19, hA20,

* hB1, hB2, hB3, hB4, hB5, hB6, hB7, hB8, hB9, hB10,

* hB11, hB12, hB13, hB14, hB15, hB16, hB17, hB18, hB19, hB20,

* hC1, hC2, hC3, hC4, hC5, hC6, hC7, hC8, hC9, hC10,

* hC11, hC12, hC13, hC14, hC15, hC16, hC17, hC18, hC19, hC20,

* Em1, Em2, Em3, Em4, Em5, Em6, Em7, Em8, Em9,

* Em10, Em11, Em12,

* Tm1, Tm2, Tm3, Tm4, Tm5, Tm6, Tm7, Tm8, Tm9,

* Tm10, Tm11, Tm12,

* Rm1, Rm2, Rm3, Rm4, Rm5, Rm6, Rm7, Rm8, Rm9,

* Rm10, Rm11, Rm12,

* Vm1, Vm2, Vm3, Vm4, Vm5, Vm6, Vm7, Vm8, Vm9,

* Vm10, Vm11, Vm12,

* Um1, Um2, Um3, Um4, Um5, Um6, Um7, Um8, Um9,

* Um10, Um11, Um12,

* Ra1, Ra2, Ra3, Ra4, Ra5, Ra6, Ra7, Ra8, Ra9,

* Ra10, Ra11, Ra12,

* Vg1, Vg2, Vg3, Vg4, Vg5, Vg6, Vg7, Vg8, Vg9,

* Vg10, Vg11, Vg12,

* fsolo,

* bcd, rd, ae, pm, z0, ROar, ROag, RaR, VgR, RaR1, VgR1,

* vK, cp, Pbar, pmt, pmte,

* E, T, U, V, R,

* sk, CW, Sw,

* sk1, sk2, sk3, sk4, sk5, sk6, sk7, sk8, sk9, sk10,

* sk11, sk12, sk13, sk14, sk15, sk16, sk17, sk18, sk19, sk20,

* CW1, CW2, CW3, CW4, CW5, CW6, CW7, CW8, CW9, CW10,

* CW11, CW12, CW13, CW14, CW15, CW16, CW17, CW18, CW19, CW20,

* Sw1, Sw2, Sw3, Sw4, Sw5, Sw6, Sw7, Sw8, Sw9, Sw10,

* Sw11, Sw12, Sw13, Sw14, Sw15, Sw16, Sw17, Sw18, Sw19, Sw20,

* Dit, Cg, Cr, cLD,

* Dit1, Dit2, Dit3, Dit4, Dit5,

* Dit6, Dit7, Dit8, Dit9, Dit10,

* Dit11, Dit12, Dit13, Dit14, Dit15,

* Dit16, Dit17, Dit18, Dit19, Dit20, DitM, Draina,

* Cg1, Cg2, Cg3, Cg4, Cg5, Cg6, Cg7, Cg8, Cg9, Cg10,

* Cg11, Cg12, Cg13, Cg14, Cg15, Cg16, Cg17, Cg18, Cg19, Cg20,

* Cr1, Cr2, Cr3, Cr4, Cr5, Cr6, Cr7, Cr8, Cr9, Cr10,

* Cr11, Cr12, Cr13, Cr14, Cr15, Cr16, Cr17, Cr18, Cr19, Cr20,

* cLd1, cLd2, cLd3, cLd4, cLd5, cLd6, cLd7, cLd8, cLd9,

* cLd10, cLd11, cLd12, cLd13, cLd14, cLd15, cLd16, cLd17, cLd18,

* cLd19, cLd20,

* CKr, CKf,

* CKr1, CKr2, CKr3, CKr4, CKr5, CKr6, CKr7, CKr8, CKr9,

* CKr10, CKr11, CKr12, CKr13, CKr14, CKr15, CKr16, CKr17, CKr18,

* CKr19, CKr20,

* CKf1, CKf2, CKf3, CKf4, CKf5, CKf6, CKf7, CKf8, CKf9,

* CKf10, CKf11, CKf12, CKf13, CKf14, CKf15, CKf16, CKf17, CKf18,

* CKf19, CKf20,

* Ds, Ds1, Ds2, Ds3, Ds4, Ds5, Ds6, Ds7, Ds8, Ds9,

* Ps, Ps1, Ps2, Ps3, Ps4, Ps5, Ps6, Ps7, Ps8, Ps9,

* YV, YQ,

* qtMu1, qtMu2, qtMu3, qtMu4,

* qtMu5, qtMu6, qtMu7, qtMu8,

* qtMu9, qtMu10, qtMu11, qtMu12,

* qtMu13, qtMu14, qtMu15, qtMu16,

* qtMu17, qtMu18, qtMu19, qtMu20,

* flowtot1, flowtot2, flowtot3, flowtot4,

* flowtot5, flowtot6, flowtot7, flowtot8,

* flowtot9, flowtot10, flowtot11, flowtot12,

* flowtot13, flowtot14, flowtot15, flowtot16,

* flowtot17, flowtot18, flowtot19, flowtot20,

* AA, BA, AB, BB, AC, BC,

* AA1, AA2, AA3, AA4, AA5, AA6, AA7, AA8, AA9, AA10,

* AA11, AA12, AA13, AA14, AA15, AA16, AA17, AA18, AA19, AA20,

* BA1, BA2, BA3, BA4, BA5, BA6, BA7, BA8, BA9, BA10,

* BA11, BA12, BA13, BA14, BA15, BA16, BA17, BA18, BA19, BA20,

* AB1, AB2, AB3, AB4, AB5, AB6, AB7, AB8, AB9, AB10,

* AB11, AB12, AB13, AB14, AB15, AB16, AB17, AB18, AB19, AB20,

* BB1, BB2, BB3, BB4, BB5, BB6, BB7, BB8, BB9, BB10,

* BB11, BB12, BB13, BB14, BB15, BB16, BB17, BB18, BB19, BB20,

* AC1, AC2, AC3, AC4, AC5, AC6, AC7, AC8, AC9, AC10,

* AC11, AC12, AC13, AC14, AC15, AC16, AC17, AC18, AC19, AC20,

* BC1, BC2, BC3, BC4, BC5, BC6, BC7, BC8, BC9, BC10,

* BC11, BC12, BC13, BC14, BC15, BC16, BC17, BC18, BC19, BC20,

* hcc, hp,

* BBA, BBB, BBC,

* dtetA, tetA, dtetB, tetB, dtetC, tetC,

* dtet1A,

* dtet1B, dtet2B, dtet3B, dtet1C, dtet2C, dtet3C,

* tetexA, tetexB, tetexC, hss, hs,

* qss, qs, qt,

* psiiA, psiiB, psiiC, psiAm,

* eaepA, eaepB, eaepC, eaA, eaB, eaC, EVP, psi10,

* tetfcA, tetfcB, tetfcC,

* tetfcA1, tetfcA2, tetfcA3, tetfcA4, tetfcA5, tetfcA6,

* tetfcA7, tetfcA8, tetfcA9, tetfcA10, tetfcA11, tetfcA12,

* tetfcA13, tetfcA14, tetfcA15, tetfcA16, tetfcA17, tetfcA18,

* tetfcA19, tetfcA20,

* tetfcB1, tetfcB2, tetfcB3, tetfcB4, tetfcB5, tetfcB6,

* tetfcB7, tetfcB8, tetfcB9, tetfcB10, tetfcB11, tetfcB12,

* tetfcB13, tetfcB14, tetfcB15, tetfcB16, tetfcB17, tetfcB18,

* tetfcB19, tetfcB20,

* tetreA, tetreB, tetreC,

* tetreA1, tetreA2, tetreA3, tetreA4, tetreA5, tetreA6,

* tetreA7, tetreA8, tetreA9, tetreA10, tetreA11, tetreA12,

* tetreA13, tetreA14, tetreA15, tetreA16, tetreA17, tetreA18,

* tetreA19, tetreA20,

* tetreB1, tetreB2, tetreB3, tetreB4, tetreB5, tetreB6,

* tetreB7, tetreB8, tetreB9, tetreB10, tetreB11, tetreB12,

* tetreB13, tetreB14, tetreB15, tetreB16, tetreB17, tetreB18,

* tetreB19, tetreB20,

* tet10A, tet10B, tet10C,

* tet10A1, tet10A2, tet10A3, tet10A4, tet10A5, tet10A6,

* tet10A7, tet10A8, tet10A9, tet10A10, tet10A11, tet10A12,

* tet10A13, tet10A14, tet10A15, tet10A16, tet10A17, tet10A18,

* tet10A19, tet10A20,

* tet10B1, tet10B2, tet10B3, tet10B4, tet10B5, tet10B6,

* tet10B7, tet10B8, tet10B9, tet10B10, tet10B11, tet10B12,

* tet10B13, tet10B14, tet10B15, tet10B16, tet10B17, tet10B18,

* tet10B19, tet10B20,

* tetwpA, tetwpB, tetwpC,

* tetwpA1, tetwpA2, tetwpA3, tetwpA4, tetwpA5, tetwpA6,

* tetwpA7, tetwpA8, tetwpA9, tetwpA10, tetwpA11, tetwpA12,

* tetwpA13, tetwpA14, tetwpA15, tetwpA16, tetwpA17, tetwpA18,

* tetwpA19, tetwpA20,

* tetwpB1, tetwpB2, tetwpB3, tetwpB4, tetwpB5, tetwpB6,

* tetwpB7, tetwpB8, tetwpB9, tetwpB10, tetwpB11, tetwpB12,

* tetwpB13, tetwpB14, tetwpB15, tetwpB16, tetwpB17, tetwpB18,

* tetwpB19, tetwpB20, timeold,

* FA1, FA2, FA3, FA4, FA5, FA6, FA7, FA8, FA9, FA10,

* FA11, FA12, FA13, FA14, FA15, FA16, FA17, FA18, FA19, FA20,

* setm, hpche, hpday, hpmon, hpyea, qche, qday, qmon, qyea,

* VPI, VPIR, FCOS,

* zriv, zwt, ztet, hgw, htetsC, qgwi

c----Open Input Data File (Input.txt)

open (8,file='Input_Data.txt',status='unknown')

2 format(40A4)

c---(1)Read Initial year, antecedent leap year, number of rainfall data

read (8,2) ANOSIB

read (8,*) iyear0, iyearbi, N

c---(2)Read Basin area (km2), river lenght (km) and elevation (m)

read (8,2) BASIN

read (8,*) AT, RL, HEB, HSB

c---(3)Read Method of propagation MPROP(0=global; 1=Convex; 2=Muskingum) and No of Areas MDIST (1 to 20)

read (8,2) PROPAG

read (8,*) MPROP, MDIST

c---(3)Read Method of Time of Concentration MTCON(1=Kirpich; 2=Kerby, 3=obs) (min); Shape Muskingum factor (X)

read (8,2) TCONT

read (8,*) MTCON, tc, tinttMD, X

c---(4)Read Method to calculate porosity MPOROS (1=Saxton; 2=(2.65-BDA)/2.65; 3=obs)

read (8,2) POROSITY

read (8,*) MPOROS, BDA, BDB, BDC, porA, porB, porC

c---(4)Read Method to calculate Teta Saturated MTETAS(1=farA*por; 2=obs)

read (8,2) TETASAT

read (8,*) MTETAS, farA, farB, farC, tetsA, tetsB, tetsC

c---(4)Read Method to calculate initial moisture MTETAI(1=fuminA*tetasA; 2=obs) (0 < fumin < 1)

read (8,2) UMIINI

read (8,*) MTETAI, fuminA, fuminB, fuminC, tetiA, tetiB, tetiC

c---(4)Read Method to calculate Infiltration MCITET(1=Ct.(tet)**Dt;] (-); 2=ftetF+(ftet0-ftetF)exp(-ktet(tet))

read (8,2) CITETA

read (8,*) MCITET, Ct, Dt, ftet0, ftetF, ktet, PHI

c---(4)Read Soils texture (%Sand; %Clay)

read (8,2) TEORAR

read (8,*) SA, SB, SC, CA, CB, CC

c---(4)Read Method to calculate Ksat: MKSAT(1=Saxton; 2=Rawls; 3=Brak; 4=Cosby; 5=obs) (mm/h)

c---(4)Method to calculate K(tet):MKTET(1=B-Cor;2=Camp;3=Sax;4=VanG) (mm/h)

read (8,2) SATHYC

read (8,*) MKSAT, MKTET, KsatA, KsatB, KsatC, AvA, AvB, AvC

c---(4)Read Soils depth (cm)

read (8,2) PROFSO

read (8,*) hA, hB, hC

c---(5)Read data of interception (mm)

read (8,2) DINTCEP

read (8,*) Dit, Cg, Cr, cLd

c---(6)Mean daily pan evaporation rate(mm/dia)

read (8,2) EVAPTQ

read (8,*) Em1, Em2, Em3, Em4, Em5, Em6, Em7, Em8, Em9, Em10,

* Em11, Em12

c---(6)Read Mean daily temperatures of each month (oC)

read (8,2) TEMPME

read (8,*) Tm1, Tm2, Tm3, Tm4, Tm5, Tm6, Tm7, Tm8, Tm9, Tm10,

* Tm11, Tm12

c---(6)Read Mean daily radiations of each month (cal/cm2)

read (8,2) RADIME

read (8,*) Rm1, Rm2, Rm3, Rm4, Rm5, Rm6, Rm7, Rm8, Rm9, Rm10,

* Rm11, Rm12

c---(6)Read Mean daily wind velocities of each month (m/s)

read (8,2) VELOME

read (8,*) Vm1, Vm2, Vm3, Vm4, Vm5, Vm6, Vm7, Vm8, Vm9, Vm10,

* Vm11, Vm12

c---(6)Read Mean daily relative moisture content of the air (%)

read (8,2) UMIDME

read (8,*) Um1, Um2, Um3, Um4, Um5, Um6, Um7, Um8, Um9, Um10,

* Um11, Um12

c---(6)Read Mean aerodynamic resistance to vapour transpot (s/m)

read (8,2) AERODR

read (8,*) Ra1, Ra2, Ra3, Ra4, Ra5, Ra6, Ra7, Ra8, Ra9, Ra10,

* Ra11, Ra12

c---(6)Read Mean vegetation resistance to water transport (s/m)

read (8,2) VEGETR

read (8,*) Vg1, Vg2, Vg3, Vg4, Vg5, Vg6, Vg7, Vg8, Vg9, Vg10,

* Vg11, Vg12

c---(6)Read Method to calculate potential evaporation MEVAPT(1=obs; 2=Bl-Criddle; 3=radiation; 4=aerodynamic; 5=Penman)

c--- fsolo=factor for the actual evapotranspiration function (eaep)

read (8,2) OPTEVA

read (8,*) MEVAPT, fsolo, z0

c---(7)Read Channel data: Manning's n(=sk) (m**(1/3)/s), channel width CW (m), channel slope sw (m/m)

read (8,2) CANALM

read (8,*) sk, CW, Sw

c---(7)Read Canopy cover Cg and rock cover Cr (-) [0 <= Cg,Cr <= 1] - Soil erosion component

c read (8,2) FTCVGR

c read (8,*) Cg, Cr

c---(7)Read Drainage (DRAINA) from the vegetation and proportion (cLd) that falls as leaf drip (mm)

c read (8,2) INTRVG

c read (8,*) Draina, cLd

c---(7)Read Rainfall and runoff erosivity factors (-)

read (8,2) FATERF

read (8,*) CKr, CKf

c---(7)Read Method to calculate flow transport capacity MSYTC(1=E-hansen; 2=Yalin; 3=Laursen_Ds; 4=Laursen_Distributed)

read (8,2) CAPTRA

read (8,*) MSYTC

c---(7)Read Representative sediment diameter D50 (mm); Sediment size distributions (mm)

read (8,2) DIMSED

read (8,*) Ds, Ds1, Ds2, Ds3, Ds4, Ds5, Ds6, Ds7, Ds8, Ds9

c---Sediment size (m)

Ds=Ds/1000.

Ds1=Ds1/1000.

Ds2=Ds2/1000.

Ds3=Ds3/1000.

Ds4=Ds4/1000.

Ds5=Ds5/1000.

Ds6=Ds6/1000.

Ds7=Ds7/1000.

Ds8=Ds8/1000.

Ds9=Ds9/1000.

c---(7)Read Percentage of each sediment size (%)

read (8,2) PORCDS

read (8,*) Ps, Ps1, Ps2, Ps3, Ps4, Ps5, Ps6, Ps7, Ps8, Ps9

c---(8)Read Water table position in relation to ground surface (cm) and initial ground water (l/s)

read (8,2) GWATER

read (8,*) zwt, qgwi

c---DISTRIBUTED DATA

c---(4)Open Distributed data File

open (9,file='Dist_Data.txt',status='unknown')

3 format(40A4)

c---(2)Read Subbasins areas (km2)

read (9,3) AREAS

read (9,*) AR1, AR2, AR3, AR4, AR5, AR6, AR7, AR8, AR9, AR10,

* AR11, AR12, AR13, AR14, AR15, AR16, AR17, AR18, AR19, AR20

c---(2)Read Reaches lenght (m)

read (9,3) RLENGHT1

read (9,*) RL1, RL2, RL3, RL4, RL5, RL6, RL7, RL8, RL9, RL10,

* RL11, RL12, RL13, RL14, RL15, RL16, RL17, RL18, RL19, RL20

c---(2)Read Reaches elevation (m)

read (9,3) RHEIGHT1

read (9,*) RH1, RH2, RH3, RH4, RH5, RH6, RH7, RH8, RH9, RH10,

* RH11, RH12, RH13, RH14, RH15, RH16, RH17, RH18, RH19, RH20

c---(4)Read Bulk Density (BDA)

read (9,3) BULKDA

read (9,*) BDA1, BDA2, BDA3, BDA4, BDA5,

* BDA6, BDA7, BDA8, BDA9, BDA10,

* BDA11, BDA12, BDA13, BDA14, BDA15,

* BDA16, BDA17, BDA18, BDA19, BDA20

c---(4)Read Bulk Density (BDB)

read (9,3) BULKDB

read (9,*) BDB1, BDB2, BDB3, BDB4, BDB5,

* BDB6, BDB7, BDB8, BDB9, BDB10,

* BDB11, BDB12, BDB13, BDB14, BDB15,

* BDB16, BDB17, BDB18, BDB19, BDB20

c---(4)Read Bulk Density (BDC)

read (9,3) BULKDC

read (9,*) BDC1, BDC2, BDC3, BDC4, BDC5,

* BDC6, BDC7, BDC8, BDC9, BDC10,

* BDC11, BDC12, BDC13, BDC14, BDC15,

* BDC16, BDC17, BDC18, BDC19, BDC20

c---(4)Read Soil Porosity (porA)

read (9,3) POROSA

read (9,*) porA1, porA2, porA3, porA4, porA5,

* porA6, porA7, porA8, porA9, porA10,

* porA11, porA12, porA13, porA14, porA15,

* porA16, porA17, porA18, porA19, porA20

c---(4)Read Soil Porosity (porB)

read (9,3) POROSB

read (9,*) porB1, porB2, porB3, porB4, porB5,

* porB6, porB7, porB8, porB9, porB10,

* porB11, porB12, porB13, porB14, porB15,

* porB16, porB17, porB18, porB19, porB20

c---(4)Read Soil Porosity (porC)

read (9,3) POROSC

read (9,*) porC1, porC2, porC3, porC4, porC5,

* porC6, porC7, porC8, porC9, porC10,

* porC11, porC12, porC13, porC14, porC15,

* porC16, porC17, porC18, porC19, porC20

c---(4)Read Air entrapment factor (farA): soil A

read (9,3) FAIRENTA

read (9,*) farA1, farA2, farA3, farA4, farA5,

* farA6, farA7, farA8, farA9, farA10,

* farA11, farA12, farA13, farA14, farA15,

* farA16, farA17, farA18, farA19, farA20

c---(4)Read Air entrapment factor (farB): soil B

read (9,3) FAIRENTB

read (9,*) farB1, farB2, farB3, farB4, farB5,

* farB6, farB7, farB8, farB9, farB10,

* farB11, farB12, farB13, farB14, farB15,

* farB16, farB17, farB18, farB19, farB20

c---(4)Read Air entrapment factor (farC): soil C

read (9,3) FAIRENTC

read (9,*) farC1, farC2, farC3, farC4, farC5,

* farC6, farC7, farC8, farC9, farC10,

* farC11, farC12, farC13, farC14, farC15,

* farC16, farC17, farC18, farC19, farC20

c---(4)Read saturated moisture content (DTETASA): soil A

read (9,3) DTETASA

read (9,*) tetsA1, tetsA2, tetsA3, tetsA4, tetsA5,

* tetsA6, tetsA7, tetsA8, tetsA9, tetsA10,

* tetsA11, tetsA12, tetsA13, tetsA14, tetsA15,

* tetsA16, tetsA17, tetsA18, tetsA19, tetsA20

c---(4)Read saturated moisture content (DTETASB): soil B

read (9,3) DTETASB

read (9,*) tetsB1, tetsB2, tetsB3, tetsB4, tetsB5,

* tetsB6, tetsB7, tetsB8, tetsB9, tetsB10,

* tetsB11, tetsB12, tetsB13, tetsB14, tetsB15,

* tetsB16, tetsB17, tetsB18, tetsB19, tetsB20

c---(4)Read saturated moisture content (DTETASC): soil C

read (9,3) DTETASC

read (9,*) tetsC1, tetsC2, tetsC3, tetsC4, tetsC5,

* tetsC6, tetsC7, tetsC8, tetsC9, tetsC10,

* tetsC11, tetsC12, tetsC13, tetsC14, tetsC15,

* tetsC16, tetsC17, tetsC18, tetsC19, tetsC20

c---(4)Read Factor for initial moisture content (fuminA): soil A

read (9,3) FUMINA

read (9,*) fuminA1, fuminA2, fuminA3, fuminA4, fuminA5,

* fuminA6, fuminA7, fuminA8, fuminA9, fuminA10,

* fuminA11, fuminA12, fuminA13, fuminA14, fuminA15,

* fuminA16, fuminA17, fuminA18, fuminA19, fuminA20

c---(4)Read Factor for initial moisture content (fuminB): soil B

read (9,3) FUMINB

read (9,*) fuminB1, fuminB2, fuminB3, fuminB4, fuminB5,

* fuminB6, fuminB7, fuminB8, fuminB9, fuminB10,

* fuminB11, fuminB12, fuminB13, fuminB14, fuminB15,

* fuminB16, fuminB17, fuminB18, fuminB19, fuminB20

c---(4)Read Factor for initial moisture content (fuminC): soil C

read (9,3) FUMINC

read (9,*) fuminC1, fuminC2, fuminC3, fuminC4, fuminC5,

* fuminC6, fuminC7, fuminC8, fuminC9, fuminC10,

* fuminC11, fuminC12, fuminC13, fuminC14, fuminC15,

* fuminC16, fuminC17, fuminC18, fuminC19, fuminC20

c---(4)Read Initial soil moisture content (tetaiA): soil A

read (9,3) TETAI_A

read (9,*) tetiA1, tetiA2, tetiA3, tetiA4, tetiA5,

* tetiA6, tetiA7, tetiA8, tetiA9, tetiA10,

* tetiA11, tetiA12, tetiA13, tetiA14, tetiA15,

* tetiA16, tetiA17, tetiA18, tetiA19, tetiA20

c---(4)Read Initial soil moisture content (tetaiB): soil B

read (9,3) TETAI_B

read (9,*) tetiB1, tetiB2, tetiB3, tetiB4, tetiB5,

* tetiB6, tetiB7, tetiB8, tetiB9, tetiB10,

* tetiB11, tetiB12, tetiB13, tetiB14, tetiB15,

* tetiB16, tetiB17, tetiB18, tetiB19, tetiB20

c---(4)Read Initial soil moisture content (tetaiC): soil C

read (9,3) TETAI_C

read (9,*) tetiC1, tetiC2, tetiC3, tetiC4, tetiC5,

* tetiC6, tetiC7, tetiC8, tetiC9, tetiC10,

* tetiC11, tetiC12, tetiC13, tetiC14, tetiC15,

* tetiC16, tetiC17, tetiC18, tetiC19, tetiC20

c---(4)Read Parameters of the infiltration function f=C(Teta)**(-D)

read (9,3) CTTETA

read (9,*) Ct1, Ct2, Ct3, Ct4, Ct5, Ct6, Ct7, Ct8, Ct9, Ct10,

* Ct11, Ct12, Ct13, Ct14, Ct15, Ct16, Ct17, Ct18, Ct19, Ct20

c---(4)Read Parameters of the infiltration function f=C(Teta)**(-D)

read (9,3) DTTETA

read (9,*) Dt1, Dt2, Dt3, Dt4, Dt5, Dt6, Dt7, Dt8, Dt9, Dt10,

* Dt11, Dt12, Dt13, Dt14, Dt15, Dt16, Dt17, Dt18, Dt19, Dt20

c---(4)Read Parameters of the infiltration function f=f0+(f0-fC)e**(-kteta)

read (9,3) FO

read (9,*) FO1, FO2, FO3, FO4, FO5, FO6, FO7, FO8, FO9, FO10,

* FO11, FO12, FO13, FO14, FO15, FO16, FO17, FO18, FO19, FO20

c---(4)Read Parameters of the infiltration function f=f0+(f0-fC)e**(-kteta)

read (9,3) FC

read (9,*) FC1, FC2, FC3, FC4, FC5, FC6, FC7, FC8, FC9, FC10,

* FC11, FC12, FC13, FC14, FC15, FC16, FC17, FC18, FC19, FC20

c---(4)Read Parameters of the infiltration function f=f0+(f0-fC)e**(-kteta)

read (9,3) KT

read (9,*) kt1, kt2, kt3, kt4, kt5, kt6, kt7, kt8, kt9, kt10,

* kt11, kt12, kt13, kt14, kt15, kt16, kt17, kt18, kt19, kt20

c---(4)Read Soils Texture (SANDA): soil A

read (9,3) SANDA

read (9,*) SA1, SA2, SA3, SA4, SA5,

* SA6, SA7, SA8, SA9, SA10,

* SA11, SA12, SA13, SA14, SA15,

* SA16, SA17, SA18, SA19, SA20

c---(4)Read Soils Texture (SANDB): soil B

read (9,3) SANDB

read (9,*) SB1, SB2, SB3, SB4, SB5,

* SB6, SB7, SB8, SB9, SB10,

* SB11, SB12, SB13, SB14, SB15,

* SB16, SB17, SB18, SB19, SB20

c---(4)Read Soils Texture (SANDC): soil C

read (9,3) SANDC

read (9,*) SC1, SC2, SC3, SC4, SC5,

* SC6, SC7, SC8, SC9, SC10,

* SC11, SC12, SC13, SC14, SC15,

* SC16, SC17, SC18, SC19, SC20

c---(4)Read Soils Texture (CLAYA): soil A

read (9,3) CLAYA

read (9,*) CA1, CA2, CA3, CA4, CA5,

* CA6, CA7, CA8, CA9, CA10,

* CA11, CA12, CA13, CA14, CA15,

* CA16, CA17, CA18, CA19, CA20

c---(4)Read Soils Texture (CLAYB): soil B

read (9,3) CLAYB

read (9,*) CB1, CB2, CB3, CB4, CB5,

* CB6, CB7, CB8, CB9, CB10,

* CB11, CB12, CB13, CB14, CB15,

* CB16, CB17, CB18, CB19, CB20

c---(4)Read Soils Texture (CLAYC): soil C

read (9,3) CLAYC

read (9,*) CC1, CC2, CC3, CC4, CC5,

* CC6, CC7, CC8, CC9, CC10,

* CC11, CC12, CC13, CC14, CC15,

* CC16, CC17, CC18, CC19, CC20

c---(4)Read Saturated conductivity (SATKA): soil A

read (9,3) KSAT_A

read (9,*) KsatA1, KsatA2, KsatA3, KsatA4, KsatA5,

* KsatA6, KsatA7, KsatA8, KsatA9, KsatA10,

* KsatA11, KsatA12, KsatA13, KsatA14, KsatA15,

* KsatA16, KsatA17, KsatA18, KsatA19, KsatA20

c---(4)Read Saturated conductivity (SATKB): soil B

read (9,3) KSAT_B

read (9,*) KsatB1, KsatB2, KsatB3, KsatB4, KsatB5,

* KsatB6, KsatB7, KsatB8, KsatB9, KsatB10,

* KsatB11, KsatB12, KsatB13, KsatB14, KsatB15,

* KsatB16, KsatB17, KsatB18, KsatB19, KsatB20

c---(4)Read Saturated conductivity (SATKC): soil C

read (9,3) KSAT_C

read (9,*) KsatC1, KsatC2, KsatC3, KsatC4, KsatC5,

* KsatC6, KsatC7, KsatC8, KsatC9, KsatC10,

* KsatC11, KsatC12, KsatC13, KsatC14, KsatC15,

* KsatC16, KsatC17, KsatC18, KsatC19, KsatC20

c---(4)Read Soils Depth (hA): soil A

read (9,3) DEPTHA

read (9,*) hA1, hA2, hA3, hA4, hA5,

* hA6, hA7, hA8, hA9, hA10,

* hA11, hA12, hA13, hA14, hA15,

* hA16, hA17, hA18, hA19, hA20

c---(4)Read Soils Depth (hB): soil B

read (9,3) DEPTHB

read (9,*) hB1, hB2, hB3, hB4, hB5,

* hB6, hB7, hB8, hB9, hB10,

* hB11, hB12, hB13, hB14, hB15,

* hB16, hB17, hB18, hB19, hB20

c---(4)Read Soils Depth (hC): soil C

read (9,3) DEPTHC

read (9,*) hC1, hC2, hC3, hC4, hC5,

* hC6, hC7, hC8, hC9, hC10,

* hC11, hC12, hC13, hC14, hC15,

* hC16, hC17, hC18, hC19, hC20

c---(5)Read data of interception (Dit)

read (9,3) DITCP

read (9,*) Dit1, Dit2, Dit3, Dit4, Dit5,

* Dit6, Dit7, Dit8, Dit9, Dit10,

* Dit11, Dit12, Dit13, Dit14, Dit15,

* Dit16, Dit17, Dit18, Dit19, Dit20

c---(5)Read Groud Cover (Cg)

read (9,3) GCOVER1

read (9,*) Cg1, Cg2, Cg3, Cg4, Cg5,

* Cg6, Cg7, Cg8, Cg9, Cg10,

* Cg11, Cg12, Cg13, Cg14, Cg15,

* Cg16, Cg17, Cg18, Cg19, Cg20

c---(5)Read Ground Cover (Cr)

read (9,3) GCOVER2

read (9,*) Cr1, Cr2, Cr3, Cr4, Cr5,

* Cr6, Cr7, Cr8, Cr9, Cr10,

* Cr11, Cr12, Cr13, Cr14, Cr15,

* Cr16, Cr17, Cr18, Cr19, Cr20

c---(5)Read Canopy Cover (CLD)

read (9,3) CLDD1

read (9,*) cLd1, cLd2, cLd3, cLd4, cLd5,

* cLd6, cLd7, cLd8, cLd9, cLd10,

* cLd11, cLd12, cLd13, cLd14, cLd15,

* cLd16, cLd17, cLd18, cLd19, cLd20

c---(7)Read Manning-Strickler coefficient (sk)

read (9,3) ROUGHNSS

read (9,*) sk1, sk2, sk3, sk4, sk5, sk6, sk7, sk8, sk9, sk10,

* sk11, sk12, sk13, sk14, sk15, sk16, sk17, sk18, sk19, sk20

c---(7)Read Channel Width (CW)

read (9,3) CHANNWID

read (9,*) CW1, CW2, CW3, CW4, CW5, CW6, CW7, CW8, CW9, CW10,

* CW11, CW12, CW13, CW14, CW15, CW16, CW17, CW18, CW19, CW20

c---(7)Read Channel Slope (sw)

read (9,3) CHANNSLP

read (9,*) sw1, sw2, sw3, sw4, sw5, sw6, sw7, sw8, sw9, sw10,

* sw11, sw12, sw13, sw14, sw15, sw16, sw17, sw18, sw19, sw20

c---(7)Read Erositivy Coefficients (CKr)

read (9,3) EROSIV1

read (9,*) CKr1, CKr2, CKr3, CKr4, CKr5,

* CKr6, CKr7, CKr8, CKr9, CKr10,

* CKr11, CKr12, CKr13, CKr14, CKr15,

* CKr16, CKr17, CKr18, CKr19, CKr20

c---(7)Read Erositivy Coefficients (CKf)

read (9,3) EROSIV2

read (9,*) CKf1, CKf2, CKf3, CKf4, CKf5,

* CKf6, CKf7, CKf8, CKf9, CKf10,

* CKf11, CKf12, CKf13, CKf14, CKf15,

* CKf16, CKf17, CKf18, CKf19, CKf20