PDA

View Full Version : Input dati [fortran]


gtr84
19-04-2012, 11:19
ciao a tutti

ho dei dati strutturati in questo modo

858,100858,19510101, 73, 0
858,100858,19510102, 84, 0
858,100858,19510103, 122, 0
858,100858,19510104, 116, 0
858,100858,19510105, 120, 0
858,100858,19510106, 115, 0
858,100858,19510107, 108, 0
858,100858,19510108, 105, 0
858,100858,19510109, 93, 0
858,100858,19510110, 128, 0
858,100858,19510111, 128, 0
858,100858,19510112, 110, 0
858,100858,19510113, 136, 0
eccetera.....


vorrei inserirli in una matrice N x 5

con il comando read però leggo solo la prima
colonna

per l'acquisizione ho scritto questo codice



program mediatemp

implicit none

!parameter(N=25000)
integer, dimension(23000,2) :: br
integer ::i,k,j




open(1, FILE='br19512009.txt', STATUS='old')

do i=1,1000
do j=1,6
read(1,*) br(i,j)
enddo
enddo

do i=1,10
write(*,*) br(i,3)
enddo


close(1)

end program mediatemp


dove sbaglio?

grazie

gtr84
19-04-2012, 15:25
ho risolto cambiando il codice
e senza usare il comando format


do i=1,22340
read(1,*) br(i,1),br(i,2),br(i,3),br(i,4),br(i,5)

end do


br(i,j) è la matrice in cui inserisco i dati



c'è un altro problema. la terza colonna rappresenta
la data (aaaammgg) e la quarta la temperatura media
di tale giorno.

voglio creare un vettore con le medie mensili ma
la parte difficile è che non tutti i mesi sono di 30
giorni. vorrei scrivere un codice che riesce a calcolare
le medie mensili usando precisamente i valori del
mese (tenendo anche conto che ci sono gli anni
bisestili, perché ho dati anche di fine '800)

qualcuno ha un suggerimento?