README 3.45 KB
Newer Older
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
28
Ce repertoire contient les scripts necessaires à l'exploitation des données safrans.

1) La premiere etape consiste à extraire des fichiers de donnees SAFRAN les champs nous interessant.
   C'est le role du programme SRC/importSafran.py (beaucoup trop lent)
   input: dans le repertoire (PATH_INPUT_TO_DATA_SAFRAN="/home/olivierb/solvers/trunk/SandBox/FF/AedesAlbopictus/intraAnnuel/donnees_SAFRAN/SRC/SAFRAN/") les fichiers SAFRAN_2008080107_2009080106.nc
   output: sont écrites dans le repertoire (PATH_OUTPUT_TO_DATA_SAFRAN="/home/olivierb/solvers/trunk/SandBox/FF/AedesAlbopictus/intraAnnuel/donnees_SAFRAN/AUTOROUTES/DATA/numA et numA+1"). Il y a un fichier par pas de temps ("TAIR_%i",step) où step est le pas de temps, ce fichier contient les 9892 valeurs de la grille (incompléte) SAFRAN
   

2) La seconde étape consiste à contruire un raster à partir des données safran extraites, en vue de construire les champs sur des maillage.
   Pour cela, c'est le code rasterToSafranIndex.c.
   input : #define PATH_TO_DATA_SAFRAN "/home/olivierb/solvers/trunk/SandBox/FF/AedesAlbopictus/intraAnnuel/donnees_SAFRAN/AUTOROUTES/DATA/[2010-2017]" 
   output : dans le repertoire #define PATH_TO_DATA_RASTER "/home/olivierb/solvers/trunk/SandBox/FF/AedesAlbopictus/intraAnnuel/donnees_SAFRAN/AUTOROUTES/DATA_RASTER/[2010-2017]", il y a un fichier (TAIR_%i,step) contenant le raster
   
3) La troisiéme etape consiste à construire les champs sur les maillages.
   3.1) Cas des maillages 2D/1D (ie sous repertoire AUTOROUTE)
   	
	3.1.1) extraire les coordonnees en lambertIIE (epsg:27572) des noeuds du maillage.
   	       Il s'agit d'extraire les coordonnées de tous les noeuds des 16 zones, c'est le programme buildXY.edp:
	       	  input: les 16 fichiers th0.msh,....,th15.msh
		  output: les deux fichiers X Y.
	3.1.2) construire les champs correspondants aux 16 zones
	       C'est le programme preparTemp.c:
	       	  input: les fichiers X,Y et les rasters DATA_RASTER/TAIR_200x_200(x+1)/TAIR*
		  output: les fichiers des repertoire DATA_200x/TAIRxsy
		Le fichier view.edp permet de visualiser les champs obtenues.
   3.2) Cas d'un seul maillage representant la france (ie sous repertoire FRANCE)


29
30
//fichier computeIjxy.c
permet de tracer le nuage de point (Ijxy,Tj)
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63




// commentaire obsolete, je les garde pour me souvenir du lien entre le compteur et le temps
  //chaque serie_numA contient 1 an de donnee commencant le 1 aout à 6h:
  //    annee numA:
  //    -----------
  //     nbr heure du 1 aout à 6h aux 31 decembre à minuit=152 j + 6 h=152*24 + (24-6)=3666
  //     correspondnat aux indices {8759-3666+1,...,8759}
  //     annee numA+1:
  //    -------------
  //     nbr heure du 1 janvier à 0h aux 1 aout à 6h= 8760-3666
  //     correspondnat aux indices {0,...,8759-3666}


contenu du repertoire PATH_INPUT_TO_DATA_SAFRAN:
------------------------------------------------

 SAFRAN_2010080107_2011080106.nc - ... -SAFRAN_(numA+2010)080107_(numA+2011)080106.nc- ... - SAFRAN_2016080107_2017080106.nc

 chaque fichier SAFRAN_(numA)080107_(numA+1)080106.nc contient :
 - 24*365=8760 pas de tempsles donnees horaires à partir du 1 aout à 6h
 - 9892 relevé par pas de temps (pour chaque champs)

 les données de l'année numA correspondent aux pas de temps [(numA-1)*8760 +3666, numA*8760+3666]


calcul de la photoperiode
-------------------------

1) precacul.r permet le calcul de la photoperiode par jour en fct de la latitude
2) pour tout les points du maillages, on  fait la regle de 3.