source: trunk/src/paper01/fig14/as_validation_net_1994_95_v10.pro @ 73

Last change on this file since 73 was 50, checked in by pinsard, 13 years ago

replace TROPFLUX by PROJECT

File size: 9.9 KB
Line 
1;+
2; .. _as_validation_net_1994_95_v10.pro:
3;
4; =================================
5; as_validation_net_1994_95_v10.pro
6; =================================
7;
8; DESCRIPTION
9; ===========
10;
11; .. graphviz::
12;
13;    digraph as_validation_net_1994_95_v10 {
14;       graph [
15;       rankdir="LR",
16;       ]
17;       met_data [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/met_data_as_expt_19*.txt"];
18;       swr_erai [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/erai_swr_19910101_20091231_oafluxgrid.nc"];
19;       lwr_erai [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/erai_lwr_19910101_20091231_oafluxgrid.nc"];
20;       lhf_erai [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/erai_lhf_19910101_20091231_oafluxgrid.nc"];
21;       tropflux_general [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/TropFlux_19890101_20091231.nc"];
22;       swr_tropflux_nrt [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/TropFlux_swr_19890101_20091231_NRT_v50.nc"];
23;       swr_oafluxgrid [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/swr_oafluxgrid_1985_2007.nc"];
24;       lwr_oafluxgrid [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/lwr_oafluxgrid_1985_2007.nc"];
25;       lhf_oafluxgrid [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/lhf_oafluxgrid_1985_2007.nc"];
26;       shf_oafluxgrid [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/shf_oafluxgrid_1985_2007.nc"];
27;       swr_ncep2 [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/swr_ncep2_oaflxgrid_19890101_20091231.nc"];
28;       lwr_ncep2 [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/NCEP2_flux_19890101_20090729.nc"];
29;       lhf_ncep2 [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/lhf_ncep2_oafluxgrid_19890101_20091231.nc"];
30;       shf_ncep2 [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/shf_ncep2_oafluxgrid_19890101_20091231.nc"];
31;
32;       ncep1 [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/ncep1_flux_19890101_20091231.nc"];
33;
34;       figure [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/as_validation_net_1994_95_v10.ps"];
35;
36;       as_validation_net_1994_95_v10 [shape=box,
37;       fontname=Courier,
38;       color=blue,
39;       URL="http://forge.ipsl.jussieu.fr/tropflux/broswrer/trunk/src/paper01/fig14/as_validation_net_1994_95_v10.pro",
40;       label="${TROPFLUX}/src/paper01/fig14/as_validation_net_1994_95_v10.pro"];
41;       {met_data swr_erai lwr_erai lhf_erai tropflux_general swr_tropflux_nrt swr_oafluxgrid lwr_oafluxgrid lhf_oafluxgrid shf_oafluxgrid swr_ncep2 lwr_ncep2 lhf_ncep2 shf_ncep2 ncep1} -> {as_validation_net_1994_95_v10} -> {figure}
42;    }
43;
44; SEE ALSO
45; ========
46;
47; :ref:`project_profile.sh`
48; :ref:`project_init.pro`
49; :ref:`cm_project.pro`
50;
51; :ref:`statistics.pro`
52;
53; EXAMPLES
54; ========
55;
56; ::
57;
58;  IDL> as_validation_net_1994_95_v10
59;
60;
61; TODO
62; ====
63;
64; make it work on cratos : missing data erai_swr_19910101_20091231_oafluxgrid.nc, etc
65;
66; coding rules
67;
68; get rid of spwan
69;
70; complete description
71;
72; handle IO error
73;
74; EVOLUTIONS
75; ==========
76;
77; $Id$
78;
79; - fplod 20110420T115926Z aedon.locean-ipsl.upmc.fr (Darwin)
80;
81;   * remove hard coding path
82;   * add graphviz
83;
84; - fplod 20110411T142955Z aedon.locean-ipsl.upmc.fr (Darwin)
85;
86;   * minimal header
87;
88;-
89pro as_validation_net_1994_95_v10
90@cm_general
91@cm_project
92reinitplt, /z,/invert
93key_portrait = 1
94coefpalit=.9
95
96openps, FILENAME = 'idl.ps'
97; partie a changer
98asbox=[61,62,15,16]
99date1=19941016 & date2=19951019
100
101fi=project_id_env+"met_data_as_expt_1994.txt"
102res=read_ascii(fi,data_start=1)
103ff=res.field1
104swr94=reform(ff(7,*))*0.93
105lwr94=reform(ff(8,*))
106lhf94=reform(ff(5,*))
107shf94=reform(ff(6,*))
108net94=swr94+lwr94+lhf94+shf94
109
110fi=project_id_env+"met_data_as_expt_1995.txt"
111res=read_ascii(fi,data_start=1)
112ff=res.field1
113swr95=reform(ff(7,*))*0.93
114lwr95=reform(ff(8,*))
115lhf95=reform(ff(5,*))
116shf95=reform(ff(6,*))
117net95=swr95+lwr95+lhf95+shf95
118
119net=[net94, net95]
120ind=where(finite(net))
121net=net(ind)
122
123;; start data - 16/10/94 and end date 31/12/94
124;; start data - 01/01/95 and end date 19/10/95
125
126help, net
127;; reading other data sets
128file=project_id_env+"erai_swr_19910101_20091231_oafluxgrid.nc"
129initncdf, file
130swr=read_ncdf('swr',date1,date2,file=file,/nostr, box=asbox) & swr=reform(swr)
131file=project_id_env+'erai_lwr_19890101_20091231_oafluxgrid.nc'
132initncdf, file
133lwr=read_ncdf('lwr',date1-1,date2,file=file,/nostr, box=asbox) & lwr=reform(lwr)
134file=project_id_env+'erai_shf_19890101_20091231_oafluxgrid.nc'
135initncdf, file
136shf=-1*read_ncdf('shf',date1,date2,file=file,/nostr, box=asbox) & shf=reform(shf)
137file=project_id_env+'erai_lhf_19890101_20091231_oafluxgrid.nc'
138initncdf, file
139lhf=-1*read_ncdf('lhf',date1,date2,file=file,/nostr, box=asbox) & lhf=reform(lhf)
140
141era=swr+lwr+lhf+shf
142help, era
143era=era(ind)
144help, era
145
146swr=0. & lwr=0. & shf=0. & lhf=0.
147
148file=project_id_env+"TropFlux_19890101_20091231.nc"
149initncdf, file
150swr=read_ncdf('swr',date1,date2,file=file,/nostr, box=asbox) & swr=reform(swr)
151lwr=read_ncdf('lwr',date1,date2,file=file,/nostr, box=asbox) & lwr=reform(lwr)
152shf=read_ncdf('shf',date1,date2,file=file,/nostr, box=asbox) & shf=reform(shf)
153lhf=read_ncdf('lhf',date1,date2,file=file,/nostr, box=asbox) & lhf=reform(lhf)
154trop=swr+lwr+shf+lhf
155help, trop
156
157file=project_id_env+"TropFlux_swr_19890101_20091231_NRT_v50.nc"
158initncdf, file
159olr1=read_ncdf("sw", date1-1, date2, file=file,/nostr, box=asbox) ;; & olr1=grossemoyenne(olr1, "xy",/nan, box=asbox) &
160olr1=reform(olr1)
161
162olr=olr1+lwr+shf+lhf
163
164swr=0. & lwr=0. & shf=0. & lhf=0. & olr1=0.
165
166trop=trop(ind) & olr=olr(ind)
167help, trop, olr
168
169file=project_id_env+'swr_oafluxgrid_1985_2007.nc'
170initncdf, file
171swr=read_ncdf("swr", date1, date2, file=file,/nostr, box=asbox) & swr=reform(swr)
172
173file=project_id_env+'lwr_oafluxgrid_1985_2007.nc'
174initncdf, file
175lwr=read_ncdf("lwr", date1, date2, file=file,/nostr, box=asbox) & lwr=-1*reform(lwr)
176
177file=project_id_env+'lhf_oafluxgrid_1985_2009.nc'
178initncdf, file
179lhf=read_ncdf("lhf", date1, date2, file=file,/nostr, box=asbox) & lhf=-1*reform(lhf)
180
181file=project_id_env+'shf_oafluxgrid_1985_2009.nc'
182initncdf, file
183shf=read_ncdf("shf", date1, date2, file=file,/nostr, box=asbox) & shf=-1*reform(shf)
184
185oaf=swr+lwr+shf+lhf
186oaf=oaf(ind)
187help, oaf
188
189swr=0. & lwr=0. & shf=0. & lhf=0.
190
191fi=project_id_env+'swr_ncep2_oaflxgrid_19890101_20091231.nc'
192initncdf, fi
193swr=read_ncdf("swr", date1, date2, file=fi,/nostr, box=asbox) & swr=reform(swr)
194
195fi=project_id_env+"NCEP2_flux_19890101_20090729.nc"
196initncdf, fi
197lwr=read_ncdf("lwr", date1, date2, file=fi,/nostr, box=asbox) & lwr=grossemoyenne(lwr, "xy",/nan, box=asbox) & lwr=reform(lwr)
198
199lwr=-1*lwr
200
201fi=project_id_env+'lhf_ncep2_oafluxgrid_19890101_20091231.nc'
202initncdf, fi
203lhf=read_ncdf("lhf", date1, date2, file=fi,/nostr, box=asbox) & lhf=-1*reform(lhf)
204
205fi=project_id_env+'shf_ncep2_oafluxgrid_19890101_20091231.nc'
206initncdf, fi
207shf=read_ncdf("shf", date1, date2, file=fi,/nostr, box=asbox) & shf=-1*reform(shf)
208
209nce2=swr+lwr+shf+lhf
210help, nce2
211
212nce2=nce2(ind)
213help, nce2
214swr=0. & lwr=0. & shf=0. & lhf=0.
215
216file=project_id_env+"ncep1_flux_19890101_20091231.nc"
217initncdf, file
218swr=read_ncdf("swr", date1, date2, file=file,/nostr, box=asbox) & swr=reform(swr)
219lwr=read_ncdf("lwr", date1, date2, file=file,/nostr, box=asbox) & lwr=reform(lwr)
220lhf=read_ncdf("lhf", date1, date2, file=file,/nostr, box=asbox) & lhf=reform(lhf)
221shf=read_ncdf("shf", date1, date2, file=file,/nostr, box=asbox) & shf=reform(shf)
222nce1=swr+lwr+shf+lhf
223
224help, nce1
225
226nce1=nce1(ind)
227help, nce1
228;ind=where(finite(net))
229;net=net(ind) & era=era(ind) & trop=trop(ind) & oaf=oaf(ind) & nce1=nce1(ind) & nce2=nce2(ind)
230;olr=olr(ind)
231
232statistics, net, era, $
233                 cor, bias, std, rmsd
234print, cor, bias, std,  rmsd
235cstat_era=string(cor, bias, std, rmsd, format='(f4.2,2x,f6.2,2x,f4.2,1x,f6.2)')
236
237statistics, net, trop, $
238                 cor, bias, std, rmsd
239print, cor, bias, std,  rmsd
240cstat_trop=string(cor, bias, std, rmsd, format='(f4.2,2x,f6.2,2x,f4.2,1x,f6.2)')
241
242statistics, net, oaf, $
243                 cor, bias, std, rmsd
244print, cor, bias, std,  rmsd
245cstat_oaf=string(cor, bias, std, rmsd, format='(f4.2,2x,f6.2,2x,f4.2,1x,f6.2)')
246
247statistics, net, nce1, $
248                 cor, bias, std, rmsd
249print, cor, bias, std,  rmsd
250cstat_nce1=string(cor, bias, std, rmsd, format='(f4.2,2x,f6.2,2x,f4.2,1x,f6.2)')
251
252statistics, net, nce2, $
253                 cor, bias, std, rmsd
254print, cor, bias, std,  rmsd
255cstat_nce2=string(cor, bias, std, rmsd, format='(f4.2,2x,f6.2,2x,f4.2,1x,f6.2)')
256
257statistics, net, olr, $
258                 cor, bias, std, rmsd
259print, cor, bias, std,  rmsd
260cstat_olr=string(cor, bias, std, rmsd, format='(f4.2,2x,f6.2,2x,f4.2,1x,f6.2)')
261
262splot, net, era, charsize=1.1, title='Net Flux - AS Vs ERAI', $
263     xrange=[-200,200], yrange=[-200,200], small=[2,3,1], psym=2, xmin=1, ymin=1
264oplot, [-200,200], [-200,200]
265xyouts, -150,-175, cstat_era, charsize=1.
266
267splot, net, trop, charsize=1.1, title='Net Flux - AS Vs TropFlux',/noer, $
268     xrange=[-200,200], yrange=[-200,200], small=[2,3,2], psym=2, xmin=1, ymin=1
269oplot, [-200,200], [-200,200]
270xyouts, -150,-175, cstat_trop, charsize=1.
271
272splot, net, oaf, charsize=1.1, title='Net Flux - AS Vs OAFlux',/noer, $
273     xrange=[-200,200], yrange=[-200,200], small=[2,3,3], psym=2, xmin=1, ymin=1
274oplot, [-200,200], [-200,200]
275xyouts, -150,-175, cstat_oaf, charsize=1.
276
277splot, net, nce1, charsize=1.1, title='Net Flux - AS Vs NCEP1',/noer, $
278     xrange=[-200,200], yrange=[-200,200], small=[2,3,4], psym=2, xmin=1, ymin=1
279oplot, [-200,200], [-200,200]
280xyouts, -150,-175, cstat_nce1, charsize=1.
281
282splot, net, nce2, charsize=1.1, title='Net Flux - AS Vs NCEP2',/noer, $
283     xrange=[-200,200], yrange=[-200,200], small=[2,3,5], psym=2, xmin=1, ymin=1
284oplot, [-200,200], [-200,200]
285xyouts, -150,-175, cstat_nce2, charsize=1.
286
287splot, net, olr, charsize=1.1, title='Net Flux - AS Vs TropFlux_NRT',/noer, $
288     xrange=[-200,200], yrange=[-200,200], small=[2,3,6], psym=2, xmin=1, ymin=1
289oplot, [-200,200], [-200,200]
290xyouts, -150,-175, cstat_olr, charsize=1.
291
292closeps
293fig=project_od_env+'as_validation_net_1994_95_v10.ps'
294spawn, 'mv '+psdir+'idl.ps '+fig
295spawn, 'gv '+fig
296end
Note: See TracBrowser for help on using the repository browser.