Changeset 97 for trunk/src/paper01/fig14/subdctn_ses_validation_net_1991_93_all_v10.pro

Timestamp:

08/30/11 14:38:29 (13 years ago)

Author:

pinsard

Message:

suppress blank lines trailing blank

File:

• trunk/src/paper01/fig14/subdctn_ses_validation_net_1991_93_all_v10.pro (modified) (3 diffs)

trunk/src/paper01/fig14/subdctn_ses_validation_net_1991_93_all_v10.pro

@@ -92 +92 @@
 key_portrait = 1
 coefpalit=.9
-
+;
 openps, FILENAME = project_od_env+"subdctn_ses_validation_net_1991_93_all_v10.ps"
 ; partie a changer
-
-;; subduction experiment in the atlantic in 1991-93
+;
+; subduction experiment in the atlantic in 1991-93
 box=[337.5, 338.5, 17.5, 18.5]
-
-;; first leg
-
-;; second leg
-
+;
+; first leg
+;
+; second leg
+;
 fi=project_id_env+"met_data_subduct_ses_92.txt"
 res=read_ascii(fi,data_start=1)
@@ -112 +112 @@
 net_92=[swr_92+lwr_92+shf_92+lhf_92]
 help, net_92
-
-;; second leg
-
+;
+; second leg
+;
 fi=project_id_env+"met_data_subduct_ses_93.txt"
 res=read_ascii(fi,data_start=1)
@@ -124 +124 @@
 net_93=[swr_93+lwr_93+shf_93+lhf_93]
 help, net_93
-
+;
 net=[net_92, net_93]
 help, net
-
-;; reading other data sets
-
-st1=19920209 
+;
+; reading other data sets
+;
+st1=19920209
 en1=19920911
-st2=19920921 
+st2=19920921
 en2=19930618
-
+;
 file=project_id_env+"erai_swr_19910101_20091231_oafluxgrid.nc"
 initncdf, file
-swr1=read_ncdf('swr',st1,en1,file=file,/nostr, box=box) 
-swr1=grossemoyenne(swr1, "xy", box=box,/nan) 
+swr1=read_ncdf('swr',st1,en1,file=file,/nostr, box=box)
+swr1=grossemoyenne(swr1, "xy", box=box,/nan)
 swr1=reform(swr1)
-swr2=read_ncdf('swr',st2,en2,file=file,/nostr, box=box) 
-swr2=grossemoyenne(swr2, "xy", box=box,/nan) 
+swr2=read_ncdf('swr',st2,en2,file=file,/nostr, box=box)
+swr2=grossemoyenne(swr2, "xy", box=box,/nan)
 swr2=reform(swr2)
-
+;
 file=project_id_env+'erai_lhf_19890101_20091231_oafluxgrid.nc'
 initncdf, file
-lhf1=-1*read_ncdf('lhf',st1,en1,file=file,/nostr, box=box) 
-lhf1=grossemoyenne(lhf1, "xy", box=box,/nan) 
+lhf1=-1*read_ncdf('lhf',st1,en1,file=file,/nostr, box=box)
+lhf1=grossemoyenne(lhf1, "xy", box=box,/nan)
 lhf1=reform(lhf1)
-lhf2=-1*read_ncdf('lhf',st2,en2,file=file,/nostr, box=box) 
-lhf2=grossemoyenne(lhf2, "xy", box=box,/nan) 
+lhf2=-1*read_ncdf('lhf',st2,en2,file=file,/nostr, box=box)
+lhf2=grossemoyenne(lhf2, "xy", box=box,/nan)
 lhf2=reform(lhf2)
-
+;
 file=project_id_env+'erai_shf_19890101_20091231_oafluxgrid.nc'
 initncdf, file
-shf1=-1*read_ncdf('shf',st1,en1,file=file,/nostr, box=box) 
-shf1=grossemoyenne(shf1, "xy", box=box,/nan) 
+shf1=-1*read_ncdf('shf',st1,en1,file=file,/nostr, box=box)
+shf1=grossemoyenne(shf1, "xy", box=box,/nan)
 shf1=reform(shf1)
-shf2=-1*read_ncdf('shf',st2,en2,file=file,/nostr, box=box) 
-shf2=grossemoyenne(shf2, "xy", box=box,/nan) 
+shf2=-1*read_ncdf('shf',st2,en2,file=file,/nostr, box=box)
+shf2=grossemoyenne(shf2, "xy", box=box,/nan)
 shf2=reform(shf2)
-
+;
 file=project_id_env+'erai_lwr_19890101_20091231_oafluxgrid.nc'
 initncdf, file
-lwr1=read_ncdf('lwr',st1-1,en1,file=file,/nostr, box=box) 
-lwr1=grossemoyenne(lwr1, "xy", box=box,/nan) 
+lwr1=read_ncdf('lwr',st1-1,en1,file=file,/nostr, box=box)
+lwr1=grossemoyenne(lwr1, "xy", box=box,/nan)
 lwr1=reform(lwr1)
-lwr2=read_ncdf('lwr',st2-1,en2,file=file,/nostr, box=box) 
-lwr2=grossemoyenne(lwr2, "xy", box=box,/nan) 
+lwr2=read_ncdf('lwr',st2-1,en2,file=file,/nostr, box=box)
+lwr2=grossemoyenne(lwr2, "xy", box=box,/nan)
 lwr2=reform(lwr2)
-
-net1=swr1+(lwr1+shf1+lhf1) 
+;
+net1=swr1+(lwr1+shf1+lhf1)
 net2=swr2+(lwr2+shf2+lhf2)
 era=[net1,net2]
 help, era
-
-lhf1=0. 
-lhf2=0. 
-shf1=0. 
+;
+lhf1=0.
+lhf2=0.
+shf1=0.
 shf2=0.
-lwr1=0. 
-lwr2=0. 
-swr1=0. 
+lwr1=0.
+lwr2=0.
+swr1=0.
 swr2=0.
-net1=0. 
+net1=0.
 net2=0.
-
+;
 file=project_id_env+"TropFlux_19890101_20091231.nc"
 initncdf, file
-lhf1=read_ncdf('lhf',st1,en1,file=file,/nostr, box=box) 
-lhf1=grossemoyenne(lhf1, "xy", box=box,/nan) 
+lhf1=read_ncdf('lhf',st1,en1,file=file,/nostr, box=box)
+lhf1=grossemoyenne(lhf1, "xy", box=box,/nan)
 lhf1=reform(lhf1)
-lhf2=read_ncdf('lhf',st2,en2,file=file,/nostr, box=box) 
-lhf2=grossemoyenne(lhf2, "xy", box=box,/nan) 
+lhf2=read_ncdf('lhf',st2,en2,file=file,/nostr, box=box)
+lhf2=grossemoyenne(lhf2, "xy", box=box,/nan)
 lhf2=reform(lhf2)
-shf1=read_ncdf('shf',st1,en1,file=file,/nostr, box=box) 
-shf1=grossemoyenne(shf1, "xy", box=box,/nan) 
+shf1=read_ncdf('shf',st1,en1,file=file,/nostr, box=box)
+shf1=grossemoyenne(shf1, "xy", box=box,/nan)
 shf1=reform(shf1)
-shf2=read_ncdf('shf',st2,en2,file=file,/nostr, box=box) 
-shf2=grossemoyenne(shf2, "xy", box=box,/nan) 
+shf2=read_ncdf('shf',st2,en2,file=file,/nostr, box=box)
+shf2=grossemoyenne(shf2, "xy", box=box,/nan)
 shf2=reform(shf2)
-lwr1=read_ncdf('lwr',st1,en1,file=file,/nostr, box=box) 
-lwr1=grossemoyenne(lwr1, "xy", box=box,/nan) 
+lwr1=read_ncdf('lwr',st1,en1,file=file,/nostr, box=box)
+lwr1=grossemoyenne(lwr1, "xy", box=box,/nan)
 lwr1=reform(lwr1)
-lwr2=read_ncdf('lwr',st2,en2,file=file,/nostr, box=box) 
-lwr2=grossemoyenne(lwr2, "xy", box=box,/nan) 
+lwr2=read_ncdf('lwr',st2,en2,file=file,/nostr, box=box)
+lwr2=grossemoyenne(lwr2, "xy", box=box,/nan)
 lwr2=reform(lwr2)
-swr1=read_ncdf('swr',st1,en1,file=file,/nostr, box=box) 
-swr1=grossemoyenne(swr1, "xy", box=box,/nan) 
+swr1=read_ncdf('swr',st1,en1,file=file,/nostr, box=box)
+swr1=grossemoyenne(swr1, "xy", box=box,/nan)
 swr1=reform(swr1)
-swr2=read_ncdf('swr',st2,en2,file=file,/nostr, box=box) 
-swr2=grossemoyenne(swr2, "xy", box=box,/nan) 
+swr2=read_ncdf('swr',st2,en2,file=file,/nostr, box=box)
+swr2=grossemoyenne(swr2, "xy", box=box,/nan)
 swr2=reform(swr2)
-
+;
 file=project_id_env+"TropFlux_swr_19890101_20091231_NRT_v50.nc"
 initncdf, file
-olr1=read_ncdf("sw", st1-1, en1, file=file,/nostr, box=box) 
-olr1=grossemoyenne(olr1, "xy",/nan, box=box) 
+olr1=read_ncdf("sw", st1-1, en1, file=file,/nostr, box=box)
+olr1=grossemoyenne(olr1, "xy",/nan, box=box)
 olr1=reform(olr1)
-olr2=read_ncdf("sw", st2-1, en2, file=file,/nostr, box=box) 
-olr2=grossemoyenne(olr2, "xy",/nan, box=box) 
+olr2=read_ncdf("sw", st2-1, en2, file=file,/nostr, box=box)
+olr2=grossemoyenne(olr2, "xy",/nan, box=box)
 olr2=reform(olr2)
-
-net1=swr1+lwr1+(shf1+lhf1) 
+;
+net1=swr1+lwr1+(shf1+lhf1)
 net2=swr2+lwr2+(shf2+lhf2)
-neto_1=olr1+lwr1+lhf1+shf1 
+neto_1=olr1+lwr1+lhf1+shf1
 neto_2=olr2+lwr2+lhf2+shf2
-trop=[net1,net2] 
+trop=[net1,net2]
 olr=[neto_1,neto_2]
 help, trop, olr
-
+;
 file=project_id_env+'swr_oafluxgrid_1985_2007.nc'
 initncdf, file
-swr1=read_ncdf("swr", st1, en1, file=file,/nostr, box=box) 
-swr1=grossemoyenne(swr1, "xy",/nan, box=box) 
+swr1=read_ncdf("swr", st1, en1, file=file,/nostr, box=box)
+swr1=grossemoyenne(swr1, "xy",/nan, box=box)
 swr1=reform(swr1)
-swr2=read_ncdf("swr", st2, en2, file=file,/nostr, box=box) 
-swr2=grossemoyenne(swr2, "xy",/nan, box=box) 
+swr2=read_ncdf("swr", st2, en2, file=file,/nostr, box=box)
+swr2=grossemoyenne(swr2, "xy",/nan, box=box)
 swr2=reform(swr2)
-
+;
 file=project_id_env+'lwr_oafluxgrid_1985_2007.nc'
 initncdf, file
-lwr1=-1*read_ncdf("lwr", st1, en1, file=file,/nostr, box=box) 
-lwr1=grossemoyenne(lwr1, "xy",/nan, box=box) 
+lwr1=-1*read_ncdf("lwr", st1, en1, file=file,/nostr, box=box)
+lwr1=grossemoyenne(lwr1, "xy",/nan, box=box)
 lwr1=reform(lwr1)
-lwr2=-1*read_ncdf("lwr", st2, en2, file=file,/nostr, box=box) 
-lwr2=grossemoyenne(lwr2, "xy",/nan, box=box) 
+lwr2=-1*read_ncdf("lwr", st2, en2, file=file,/nostr, box=box)
+lwr2=grossemoyenne(lwr2, "xy",/nan, box=box)
 lwr2=reform(lwr2)
-
+;
 file=project_id_env+'lhf_oafluxgrid_1985_2009.nc'
 initncdf, file
-lhf1=-1*read_ncdf("lhf", st1, en1, file=file,/nostr, box=box) 
-lhf1=grossemoyenne(lhf1, "xy",/nan, box=box) 
+lhf1=-1*read_ncdf("lhf", st1, en1, file=file,/nostr, box=box)
+lhf1=grossemoyenne(lhf1, "xy",/nan, box=box)
 lhf1=reform(lhf1)
-lhf2=-1*read_ncdf("lhf", st2, en2, file=file,/nostr, box=box) 
-lhf2=grossemoyenne(lhf2, "xy",/nan, box=box) 
+lhf2=-1*read_ncdf("lhf", st2, en2, file=file,/nostr, box=box)
+lhf2=grossemoyenne(lhf2, "xy",/nan, box=box)
 lhf2=reform(lhf2)
-
+;
 file=project_id_env+'shf_oafluxgrid_1985_2009.nc'
 initncdf, file
-shf1=-1*read_ncdf("shf", st1, en1, file=file,/nostr, box=box) 
-shf1=grossemoyenne(shf1, "xy",/nan, box=box) 
+shf1=-1*read_ncdf("shf", st1, en1, file=file,/nostr, box=box)
+shf1=grossemoyenne(shf1, "xy",/nan, box=box)
 shf1=reform(shf1)
-shf2=-1*read_ncdf("shf", st2, en2, file=file,/nostr, box=box) 
-shf2=grossemoyenne(shf2, "xy",/nan, box=box) 
+shf2=-1*read_ncdf("shf", st2, en2, file=file,/nostr, box=box)
+shf2=grossemoyenne(shf2, "xy",/nan, box=box)
 shf2=reform(shf2)
-
-net1=swr1+lwr1+lhf1+shf1 
+;
+net1=swr1+lwr1+lhf1+shf1
 net2=swr2+lwr2+lhf2+shf2
 oaf=[net1,net2]
 help, oaf
-
+;
 fi=project_id_env+'swr_ncep2_oaflxgrid_19890101_20091231.nc'
 initncdf, fi
-swr1=read_ncdf("swr", st1, en1, file=fi,/nostr, box=box) 
-swr1=grossemoyenne(swr1, "xy",/nan, box=box) 
+swr1=read_ncdf("swr", st1, en1, file=fi,/nostr, box=box)
+swr1=grossemoyenne(swr1, "xy",/nan, box=box)
 swr1=reform(swr1)
-swr2=read_ncdf("swr", st2, en2, file=fi,/nostr, box=box) 
-swr2=grossemoyenne(swr2, "xy",/nan, box=box) 
+swr2=read_ncdf("swr", st2, en2, file=fi,/nostr, box=box)
+swr2=grossemoyenne(swr2, "xy",/nan, box=box)
 swr2=reform(swr2)
-
+;
 fi=project_id_env+"NCEP2_flux_19890101_20090729.nc"
 initncdf, fi
-lwr1=-1*read_ncdf("lwr", st1, en1, file=fi,/nostr, box=box) 
-lwr1=grossemoyenne(lwr1, "xy",/nan, box=box) 
+lwr1=-1*read_ncdf("lwr", st1, en1, file=fi,/nostr, box=box)
+lwr1=grossemoyenne(lwr1, "xy",/nan, box=box)
 lwr1=reform(lwr1)
-lwr2=-1*read_ncdf("lwr", st2, en2, file=fi,/nostr, box=box) 
-lwr2=grossemoyenne(lwr2, "xy",/nan, box=box) 
+lwr2=-1*read_ncdf("lwr", st2, en2, file=fi,/nostr, box=box)
+lwr2=grossemoyenne(lwr2, "xy",/nan, box=box)
 lwr2=reform(lwr2)
-
+;
 fi=project_id_env+'lhf_ncep2_oafluxgrid_19890101_20091231.nc'
 initncdf, fi
-lhf1=-1*read_ncdf("lhf", st1, en1, file=fi,/nostr, box=box) 
-lhf1=grossemoyenne(lhf1, "xy",/nan, box=box) 
+lhf1=-1*read_ncdf("lhf", st1, en1, file=fi,/nostr, box=box)
+lhf1=grossemoyenne(lhf1, "xy",/nan, box=box)
 lhf1=reform(lhf1)
-lhf2=-1*read_ncdf("lhf", st2, en2, file=fi,/nostr, box=box) 
-lhf2=grossemoyenne(lhf2, "xy",/nan, box=box) 
+lhf2=-1*read_ncdf("lhf", st2, en2, file=fi,/nostr, box=box)
+lhf2=grossemoyenne(lhf2, "xy",/nan, box=box)
 lhf2=reform(lhf2)
-
+;
 fi=project_id_env+'shf_ncep2_oafluxgrid_19890101_20091231.nc'
 initncdf, fi
-shf1=-1*read_ncdf("shf", st1, en1, file=fi,/nostr, box=box) 
-shf1=grossemoyenne(shf1, "xy",/nan, box=box) 
+shf1=-1*read_ncdf("shf", st1, en1, file=fi,/nostr, box=box)
+shf1=grossemoyenne(shf1, "xy",/nan, box=box)
 shf1=reform(shf1)
-shf2=-1*read_ncdf("shf", st2, en2, file=fi,/nostr, box=box) 
-shf2=grossemoyenne(shf2, "xy",/nan, box=box) 
+shf2=-1*read_ncdf("shf", st2, en2, file=fi,/nostr, box=box)
+shf2=grossemoyenne(shf2, "xy",/nan, box=box)
 shf2=reform(shf2)
-
-net1=swr1+lwr1+lhf1+shf1 
+;
+net1=swr1+lwr1+lhf1+shf1
 net2=swr2+lwr2+lhf2+shf2
 nce2=[net1,net2]
 help, nce2
-
+;
 file=project_id_env+"ncep1_flux_19890101_20091231.nc"
 initncdf, file
-lhf1=read_ncdf('lhf',st1,en1,file=file,/nostr, box=box) 
-lhf1=grossemoyenne(lhf1, "xy", box=box,/nan) 
+lhf1=read_ncdf('lhf',st1,en1,file=file,/nostr, box=box)
+lhf1=grossemoyenne(lhf1, "xy", box=box,/nan)
 lhf1=reform(lhf1)
-lhf2=read_ncdf('lhf',st2,en2,file=file,/nostr, box=box) 
-lhf2=grossemoyenne(lhf2, "xy", box=box,/nan) 
+lhf2=read_ncdf('lhf',st2,en2,file=file,/nostr, box=box)
+lhf2=grossemoyenne(lhf2, "xy", box=box,/nan)
 lhf2=reform(lhf2)
-shf1=read_ncdf('shf',st1,en1,file=file,/nostr, box=box) 
-shf1=grossemoyenne(shf1, "xy", box=box,/nan) 
+shf1=read_ncdf('shf',st1,en1,file=file,/nostr, box=box)
+shf1=grossemoyenne(shf1, "xy", box=box,/nan)
 shf1=reform(shf1)
-shf2=read_ncdf('shf',st2,en2,file=file,/nostr, box=box) 
-shf2=grossemoyenne(shf2, "xy", box=box,/nan) 
+shf2=read_ncdf('shf',st2,en2,file=file,/nostr, box=box)
+shf2=grossemoyenne(shf2, "xy", box=box,/nan)
 shf2=reform(shf2)
-lwr1=read_ncdf('lwr',st1,en1,file=file,/nostr, box=box) 
-lwr1=grossemoyenne(lwr1, "xy", box=box,/nan) 
+lwr1=read_ncdf('lwr',st1,en1,file=file,/nostr, box=box)
+lwr1=grossemoyenne(lwr1, "xy", box=box,/nan)
 lwr1=reform(lwr1)
-lwr2=read_ncdf('lwr',st2,en2,file=file,/nostr, box=box) 
-lwr2=grossemoyenne(lwr2, "xy", box=box,/nan) 
+lwr2=read_ncdf('lwr',st2,en2,file=file,/nostr, box=box)
+lwr2=grossemoyenne(lwr2, "xy", box=box,/nan)
 lwr2=reform(lwr2)
-swr1=read_ncdf('swr',st1,en1,file=file,/nostr, box=box) 
-swr1=grossemoyenne(swr1, "xy", box=box,/nan) 
+swr1=read_ncdf('swr',st1,en1,file=file,/nostr, box=box)
+swr1=grossemoyenne(swr1, "xy", box=box,/nan)
 swr1=reform(swr1)
-swr2=read_ncdf('swr',st2,en2,file=file,/nostr, box=box) 
-swr2=grossemoyenne(swr2, "xy", box=box,/nan) 
+swr2=read_ncdf('swr',st2,en2,file=file,/nostr, box=box)
+swr2=grossemoyenne(swr2, "xy", box=box,/nan)
 swr2=reform(swr2)
-net1=swr1+(lwr1+shf1+lhf1) 
+net1=swr1+(lwr1+shf1+lhf1)
 net2=swr2+(lwr2+shf2+lhf2)
 nce1=[net1,net2]
-
+;
 help, nce1
-
+;
 ind=where(finite(net))
-net=net(ind) 
-era=era(ind) 
-trop=trop(ind) 
-oaf=oaf(ind) 
-nce2=nce2(ind) 
+net=net(ind)
+era=era(ind)
+trop=trop(ind)
+oaf=oaf(ind)
+nce2=nce2(ind)
 nce1=nce1(ind)
 olr=olr(ind)
-
+;
 statistics, net, era, $
                  cor, bias, std, rmsd
 print, cor, bias, std,  rmsd
 cstat_era=string(cor, bias, std, rmsd, format='(f4.2,2x,f6.2,2x,f4.2,1x,f6.2)')
-
+;
 statistics, net, trop, $
                  cor, bias, std, rmsd
 print, cor, bias, std,  rmsd
 cstat_trop=string(cor, bias, std, rmsd, format='(f4.2,2x,f6.2,2x,f4.2,1x,f6.2)')
-
+;
 statistics, net, oaf, $
                  cor, bias, std, rmsd
 print, cor, bias, std,  rmsd
 cstat_oaf=string(cor, bias, std, rmsd, format='(f4.2,2x,f6.2,2x,f4.2,1x,f6.2)')
-
+;
 statistics, net, nce1, $
                  cor, bias, std, rmsd
 print, cor, bias, std,  rmsd
 cstat_nce1=string(cor, bias, std, rmsd, format='(f4.2,2x,f6.2,2x,f4.2,1x,f6.2)')
-
+;
 statistics, net, nce2, $
                  cor, bias, std, rmsd
 print, cor, bias, std,  rmsd
 cstat_nce2=string(cor, bias, std, rmsd, format='(f4.2,2x,f6.2,2x,f4.2,1x,f6.2)')
-
+;
 statistics, net, olr, $
                  cor, bias, std, rmsd
 print, cor, bias, std,  rmsd
 cstat_olr=string(cor, bias, std, rmsd, format='(f4.2,2x,f6.2,2x,f4.2,1x,f6.2)')
-
-
+;
+;
 splot, net, era, charsize=1.1, title='Net Flux - Subdn. ses Vs ERAI', $
      xrange=[-200,200], yrange=[-200,200], small=[2,3,1], psym=2, xmin=1, ymin=1
 oplot, [-200,200], [-200,200]
 xyouts, 0,-175, cstat_era, charsize=1.
-
+;
 splot, net, trop, charsize=1.1, title='Net Flux - Subdn. ses Vs TropFlux',/noer, $
      xrange=[-200,200], yrange=[-200,200], small=[2,3,2], psym=2, xmin=1, ymin=1
 oplot, [-200,200], [-200,200]
 xyouts, 0,-175, cstat_trop, charsize=1.
-
+;
 splot, net, oaf, charsize=1.1, title='Net Flux - Subdn. ses Vs OAFlux',/noer, $
      xrange=[-200,200], yrange=[-200,200], small=[2,3,3], psym=2, xmin=1, ymin=1
 oplot, [-200,200], [-200,200]
 xyouts, 0,-175, cstat_oaf, charsize=1.
-
+;
 splot, net, nce1, charsize=1.1, title='Net Flux - Subdn. ses Vs NCEP1',/noer, $
      xrange=[-200,200], yrange=[-200,200], small=[2,3,4], psym=2, xmin=1, ymin=1
 oplot, [-200,200], [-200,200]
 xyouts, 0,-175, cstat_nce1, charsize=1.
-
+;
 splot, net, nce2, charsize=1.1, title='Net Flux - Subdn. ses Vs NCEP2',/noer, $
      xrange=[-200,200], yrange=[-200,200], small=[2,3,5], psym=2, xmin=1, ymin=1
 oplot, [-200,200], [-200,200]
 xyouts, 0,-175, cstat_nce2, charsize=1.
-
+;
 splot, net, olr, charsize=1.1, title='Net Flux - Subdn. ses Vs TropFlux_NRT',/noer, $
      xrange=[-200,200], yrange=[-200,200], small=[2,3,6], psym=2, xmin=1, ymin=1
 oplot, [-200,200], [-200,200]
 xyouts, 0,-175, cstat_olr, charsize=1.
-
+;
 closeps
+;
 end