Escoamento do Solo

Boa tarde,

Gostaria de saber onde eu posso encontrar os dados de escoamento intermediário e do solo depois de feita a simulação. Seria no arquivo "NOSOLO" da pasta "Output"? E o que significa cada uma das colunas deste arquivo?

Obrigado.

Olá

Você tem interesse em saber as variáveis 'fluxo sub-superficial (QINT)' e 'fluxo subterrâneo (QBAS)' estimadas pelo modelo, para cada minibacia em cada intervalo de tempo?

[As variáveis QINT e QBAS estão descritas no artigo "The MGB-IPH model for large-scale rainfall—runoff modelling" (Collischonn et al., 2007, Hydrological Sciences Journal).]

No momento as saídas do modelo relacionadas a isto são os arquivos binários 'QBASTUDO_bin.MGB' (QBAS em todas minibacias e todos intervalos de tempo) e 'QITUDO_bin.MGB' (vazão gerada em cada minibacia em cada intervalo de tempo, isto é, somatório de QINT, QBAS e QSUP (fluxo superficial)) na pasta 'output'.

Se tiver interesse em ler os arquivos 'QBASTUDO_bin.MGB' e 'QITUDO_bin.MGB', posso te mostrar como fazer.

O arquivo 'NOSOLO.txt' não tem informações relevantes para isto, ele é apenas usado para checar algumas informações pelo modelo.

Att
Ayan

Ayan,

Tenho interesse, sim, em ler estes arquivos.

Obrigado.

Ok. Tu tem acesso a Matlab? Pois tenho uma rotina bem fácil de rodar no Matlab pra ler estes arquivos.
Outra possibilidade seria te passar um código em Fortran.

Tenho Matlab

Aí pode usar a seguinte rotina no matlab (exemplo para abrir o arquivo 'QBASTUDO_bin.MGB', para os outros arquivos binários é igual):


MiniB=XXX; %número de minibacias
Tempo=XXX; %número de intervalos de tempo da simulação (isto é, o número de dias de simulação; pode ser visto no arquivo 'infomgb.sim' na pasta 'input')

fileID = fopen('QBASTUDO_bin.MGB'); %abre o arquivo
BIN = fread(fileID,'single'); %lê o arquivo
p=0;qbas=zeros(Tempo,MiniB);

%obtém para cada dia de simulação o valor da variavel qbas para todas as minibacias
for k=1:Tempo
for l=1:MiniB;p=p+1;qbas(k,l)=BIN(p);
end
end
fclose(fileID);

%para plotar a serie temporal da variavel qbas para a minibacia 10:
figure;plot(qbas(:,10))


Abraço
Ayan

Obrigado