Discussion :

[arthy]

Bonjour

Je suis en train de travailler sur des données acquises selon des intervalles de dates différents (par huitaine, par décade, ...). Plus tard je vais lire ces dates pour pouvoir traiter mes données

Je souhaiterai écrire un programme qui édite les différentes dates ligne par ligne (cf fichier joint) Ici les 4 premiers caractères correspondent à l'année et les 3 derniers au jour julien correspondant à la date d'acquisition et l'on a 1 acquisition tous les 8 jours. Ne voulant pas avoir à saisir manuellement toutes ces dates surtout que c'est très fastidieux vu le nombre d'années sur lesquels je travaille et en outre la variabilité du cycle d'acquisition, comment procéder pour parvenir à "fusionner" l'année et le jour?

Voici un exp de code qui "bugg" en considérant que je travaille sur 7années (de 2001 à 2007), les données étant acquises toutes les huit jours

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
 
implicit none
 
integer,parameter::n_year = 7, n_day = 8 , n_dec = 1+ 365/n_day
! n_dec = nb d'acquisition par an, n_day = durée d'uncycle
 
integer:: year(n_year), jour(n_dec), i, j
 
year = (/(2001 +i,i=0,n_year-1)/)
 
jour =(/(i*n_day + 1,i=0,n_dec -1)/)
 
open(10, file='exp', form = 'formatted')
 
do i = 1, n_year
 
  do j=1,n_dec
 
   write(10,fmt ='(a4,a3)')year(i),jour(j)
 
  enddo
 
enddo
 
close(10)
 
end

Merci