PROGRAMMI DI
VALUTAZIONE DATI
I.Cionni, M.Mircea
Scopo: Creare tools di diagnostica e visualizzazione dati
Strumento: NCL, NCAR Comand Language
Lista dei desideri
Time
Series
Mappe
2D Mod
Mappe 2D
Mod Obs
Vertical Scatter
profiles plots
Var vs.
time
contours Contour e
“palline”
XZ
Var1 vs.
var2
Add
Standa
rd dev.
Contours
on
contours
YZ
regressi
on
Wind
labels
Time
vs. Z
Bar
Plots
Tort
e
Histograms Graphic
Scores
Taylor
diagram
Time series Modello Osservazioni
1) Leggere nel database del modello e delle osservazioni
2) Estrarre dal modello valori orari in corrispondenza dei
punti geografici in cui c’è la stazione (interpolazione
bilineare
3) Calcolo delle medie giornaliere
4) Plot dei valori orari e giornalieri
;**************namelist*************************
date_min=2005020800
date_max=2005021200
year_s="2005"
SPEC="PM10"
data="BRACE"
cspec="pm10"
sim="IT0“
………………….
linint2_points(xo,yo,O3_dummy, False, X,Y, 0)
Mappa 2D dati da modello
•
•
•
legge nell’archivio i dati e media solo i valori includi nell’intervallo temporale
scelto dall’utente
Crea un file netcdf con la media calcolata
Visualizza graficamente i dati su una mappa amministrativa
print("
NAMELIST
")
spec="c_NO2"
typ1=(/"CI0","NI0","SA0","SC0","SI0"/)
lev=20
data_typ="conc"
date_start1=2005010100
date_end1=2005123123
mod1="2005_farm.2.13.5_saprc90noTUV" ;name of dir
model1="2005REF"
;name used in
titles and labels
dir1=("/gpor_proj/minni/briganti/minnifarm/farm/out/"
) ; directory raw data
mean_file1=new(1,"string")
;dir+file already
processed
work_dir=
"/afs/enea.it/bol/user/cionni/PROGRAMMI/Tools/Reposit
ory/"
Mappa 2D
dati da
modello
Creare un’unica mappa
dai 4 domini4x4
Considera i punti di
sovrapposizione 2
realizzazioni in uno stesso
ensemble statistico.
Mappa 2D differenze
spec="c_O3"
typ1=(/"SC0","SA0","SI0","CI0","NI0"/)
typ2=(/"SC2","SA2","SI2","CI2","NI2"/)
lev=20
data_typ="conc”
date_start1=2005010100
date_end1=2005123123
date_start2=2005010100
date_end2=2005123123
mod1="2005_farm.2.13.5_saprc90noTUV"
;name
used in titles and labelsr
mod2="2005_farm.2.13.5_saprc90noTUV_GAINSnoCP
model1="2005REF"
;name used in
titles and labels
model2="2005GAINSnoCP"
dir1=("/gpor_proj/minni/briganti/minnifarm/farm/ou
t/") ; directory raw data
dir2=("/gpor_proj/minni/briganti/minnifarm/farm/ou
t/")
mean_file1=new(1,"string")
;dir+file already
processed
mean_file2=new(1,"string")
work_dir=
"/afs/enea.it/bol/user/cionni/PROGRAMMI/Tools/Repo
sitory/"
Mappa 2D Modello Osservazioni
•Legge nel database del modello
•Media sul periodo temporale selezionato
•Crea la mappa
•Leggere nel file degli scores
• Estrarre le variabili delle stazioni
• Graficare sulla mappa le stazioni con pallini
nella scala di colori della mappa
….
;********************Namelist
Osservazioni********************************
****
dom="IT0"
year="2005"
spec_o="PM10"
sel_score_mod="SIM_AVG"
sel_score_obs="OBS_AVG"
sel_station_type="Tutte"
sel_area_type="Tutte"
sel_elevation_min="Tutte"
sel_elevation_max="Tutte“
Scatter Plots
•Leggere nel file degli scores e estrae le
variabile modello e la variabile osservazioni
•Selezionare le stazioni
•Calcolare la regression line (minimi quadrati)e
la confidenza attraverso il ttest
•Graficare i valori dei dati (punti) ela linea di
regressione
; NAMELIST
dom="IT0"
year="2005"
spec="O3"
sel_score_mod="SIM_AVG"
sel_score_obs="OBS_AVG"
sel_station_type="Tutte"
sel_area_type="urbana"
sel_elevation_min="Tutte"
sel_elevation_max="Tutte"
sel_region="EMILIA_-_ROMAGNA"
Bar plot
•Leggere nel file degli scores ed estrarre le variabile modello e la variabile
osservazioni
•Selezionare le caratteristiche delle stazioni
•Calcolare la media e la deviazione standard su stazioni che hanno le stesse
caratteristiche
•Graficare i valori medi e le barre di errore.
sel_station_type=
;
;
;
sel_area_type=
;
;
;
sel_elevation_min=
;
sel_elevation_max=
;
sel_region=
;
(/\
"Tutte"\
"fondo",\
"industriale",\
"traffico"
\
/)
(/\
"Tutte"\
"rurale",\
"urbana",\
"suburbana"\
/)
(/\
"Tutte"\
"150",\
"250",\
"1000"\
/)
(/\
"Tutte"\
"249",\
"999",\
"2000"\
/)
(/\
"Tutte"\
End