Changeset 44
- Timestamp:
- 04/12/11 18:14:21 (13 years ago)
- Location:
- trunk
- Files:
-
- 13 edited
- 3 copied
Legend:
- Unmodified
- Added
- Removed
-
trunk/adm/guides/data_content.rst
r38 r44 95 95 - :ref:`interp_olr_30n30s_1989_2009.pro` 96 96 97 - :ref:`fig1_no_obs_map.pro` 98 97 99 .. _data_in_erai: 98 100 -
trunk/src/paper01/fig1/fig1_no_obs_map.pro
r43 r44 1 1 ;+ 2 ; .. _fig1_no_obs_map.pro 2 ; .. _fig1_no_obs_map.pro: 3 3 ; 4 4 ; =================== … … 9 9 ; =========== 10 10 ; 11 ; .. graphviz:: 12 ; 13 ; digraph fig1_no_obs_map { 14 ; graph [ 15 ; rankdir="LR", 16 ; ] 17 ; mask [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/mask_oaflux_30N30S.nc"]; 18 ; 19 ; lhf_no_obs [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/lhf_no_obs.txt"]; 20 ; 21 ; sw_no_obs [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/sw_no_obs.txt"]; 22 ; 23 ; lw_no_obs [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/lw_no_obs.txt"]; 24 ; 25 ; figure [shape=ellipse,fontname=Courier,label="${TROPFLUX_OD}/fig1_no_obs_map.ps"]; 26 ; 27 ; fig1_no_obs_map [shape=box, 28 ; fontname=Courier, 29 ; color=blue, 30 ; URL="http://forge.ipsl.jussieu.fr/tropflux/browser/trunk/src/paper01/fig1/fig1_no_obs_map.pro", 31 ; label="${TROPFLUX}/src/paper01/fig1/fig1_no_obs_map.pro"]; 32 ; 33 ; {mask lhf_no_obs sw_no_obs lw_no_obs} -> {fig1_no_obs_map} -> {figure} 34 ; } 35 ; 11 36 ; SEE ALSO 12 37 ; ======== 13 38 ; 39 ; :ref:`guide data OAFLUX <data_in_oaflux>` 40 ; 14 41 ; :ref:`tropflux_profile.sh` 42 ; :ref:`tropflux_init.pro` 43 ; 44 ; :ref:`oaflux_mask_30N30S.pro` 15 45 ; 16 46 ; EXAMPLES … … 22 52 ; IDL> fig1_no_obs_map 23 53 ; 54 ; TODO 55 ; ==== 56 ; 57 ; coding rules 58 ; 59 ; get rid of spwan 60 ; 61 ; complete description 62 ; 63 ; externalize function 64 ; 65 ; handle IO error 66 ; 24 67 ; EVOLUTIONS 25 68 ; ========== … … 27 70 ; $Id$ 28 71 ; 72 ; - fplod 20110412T073818Z aedon.locean-ipsl.upmc.fr (Darwin) 73 ; 74 ; * make it work on cratos 75 ; * remove hard coding path 76 ; * add graphviz 77 ; 29 78 ; - fplod 20110411T142955Z aedon.locean-ipsl.upmc.fr (Darwin) 30 79 ; … … 33 82 ;- 34 83 pro fig1_no_obs_map 35 @common 36 84 ; 85 @cm_tropflux 86 ; 37 87 reinitplt, /z,/invert 38 88 key_portrait = 1 … … 47 97 lct=65 48 98 49 file= '/Volumes/Iomega_HDD/work/flux_reconstruction/gridded_data/mask_oaflux_30N30S.nc'99 file=tropflux_id_env + 'mask_oaflux_30N30S.nc' 50 100 initncdf, file 51 101 domdef, 30,390,-30,30 52 102 msk=ncdf_lec(file,var='msk') 53 103 54 fi= '/Users/pkb/work/MY_SAXO/flux_automat/lhf_no_obs.txt'104 fi=tropflux_id_env + 'lhf_no_obs.txt' 55 105 res=read_ascii(fi,data_start=1) 56 106 ff=res.field1 … … 69 119 70 120 for n=0,NN-1 do begin 71 x=lon (n)72 y=lat (n)73 c=no (n)121 x=lon[n] 122 y=lat[n] 123 c=no[n] 74 124 cmi=lh_mi & cma=lh_ma & dc=cma-cmi 75 125 col=((10+244*(c-cmi)/dc) > 10) < 254 … … 80 130 81 131 for n=0,NN-1 do begin 82 x=lon (n)83 y=lat (n)84 c=no (n)132 x=lon[n] 133 y=lat[n] 134 c=no[n] 85 135 86 136 plots, x,y,psym=8,symsize=1.5,color=0 87 137 endfor 88 138 89 fi= '/Users/pkb/work/MY_SAXO/flux_automat/sw_no_obs.txt'139 fi=tropflux_id_env + 'sw_no_obs.txt' 90 140 res=read_ascii(fi,data_start=1) 91 141 ff=res.field1 … … 103 153 104 154 for n=0,NN-1 do begin 105 x=lon (n)106 y=lat (n)107 c=no (n)155 x=lon[n] 156 y=lat[n] 157 c=no[n] 108 158 cmi=sw_mi & cma=sw_ma & dc=cma-cmi 109 159 col=((10+244*(c-cmi)/dc) > 10) < 254 … … 114 164 115 165 for n=0,NN-1 do begin 116 x=lon (n)117 y=lat (n)118 c=no (n)166 x=lon[n] 167 y=lat[n] 168 c=no[n] 119 169 120 170 plots, x,y,psym=8,symsize=1.5,color=0 121 171 endfor 122 172 123 fi= '/Users/pkb/work/MY_SAXO/flux_automat/lw_no_obs.txt'173 fi=tropflux_id_env + 'lw_no_obs.txt' 124 174 res=read_ascii(fi,data_start=1) 125 175 ff=res.field1 … … 138 188 139 189 for n=0,NN-1 do begin 140 x=lon (n)141 y=lat (n)142 c=no (n)190 x=lon[n] 191 y=lat[n] 192 c=no[n] 143 193 cmi=lw_mi & cma=lw_ma & dc=cma-cmi 144 194 col=((10+244*(c-cmi)/dc) > 10) < 254 … … 149 199 150 200 for n=0,NN-1 do begin 151 x=lon (n)152 y=lat (n)153 c=no (n)201 x=lon[n] 202 y=lat[n] 203 c=no[n] 154 204 155 205 plots, x,y,psym=8,symsize=1.5,color=0 … … 168 218 closeps 169 219 170 fig='fig1_no_obs_map.ps' 171 spawn, 'mv '+psdir+'idl.ps '+cpsdir+fig 172 spawn, 'gv '+cpsdir+fig 173 return 220 fig=tropflux_od_env + 'fig1_no_obs_map.ps' 221 spawn, 'mv '+psdir+'idl.ps '+fig 222 spawn, 'gv '+fig 174 223 end 175 224 -
trunk/src/paper01/fig2/fig2_timeline_diagram.pro
r43 r44 9 9 ; =========== 10 10 ; 11 ; .. graphviz:: 12 ; 13 ; digraph fig2_timeline_diagram { 14 ; graph [ 15 ; rankdir="LR", 16 ; ] 17 ; 18 ; qlat [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/qlat*_dy.cdf"]; 19 ; lwnet [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/lwnet*_dy.cdf"]; 20 ; swnet [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/swnet*_dy.cdf"]; 21 ; fig2_timeline_diagram [shape=box, 22 ; fontname=Courier, 23 ; color=blue, 24 ; URL="http://forge.ipsl.jussieu.fr/tropflux/browser/trunk/src/paper01/fig2/fig2_timeline_diagram.pro", 25 ; label="${TROPFLUX}/src/paper01/fig2/fig2_timeline_diagram.pro"]; 26 ; 27 ; figure [shape=ellipse,fontname=Courier,label="${TROPFLUX_OD}/fig2_timeline_diagram.ps"]; 28 ; 29 ; {qlat lwnet swnet} -> {fig2_timeline_diagram} -> {figure} 30 ; } 31 ; 11 32 ; SEE ALSO 12 33 ; ======== 13 34 ; 14 35 ; :ref:`tropflux_profile.sh` 36 ; :ref:`tropflux_init.pro` 37 ; 38 ; :ref:`read_lh.pro` 39 ; :ref:`read_lw.pro` 40 ; :ref:`read_sw.pro` 41 ; 42 ; :func:`x_site_location` 43 ; :func:`y_site_location` 15 44 ; 16 45 ; EXAMPLES … … 22 51 ; IDL> fig2_timeline_diagram 23 52 ; 53 ; TODO 54 ; ==== 55 ; 56 ; make it work on cratos :: 57 ; % Error: 58 ; % the vector size (7670) is incompatible 59 ; % with the domain dimensions 60 ; % [jpi/nx, jpj/ny, jpk/nz, jpt] = [1/1, 1/1, 1/1, 1] 61 ; 62 ; coding rules 63 ; 64 ; get rid of spwan 65 ; 66 ; complete description 67 ; 68 ; handle IO error 69 ; 24 70 ; EVOLUTIONS 25 71 ; ========== … … 27 73 ; $Id$ 28 74 ; 75 ; - fplod 20110412T084618Z aedon.locean-ipsl.upmc.fr (Darwin) 76 ; 77 ; * remove hard coding path 78 ; * add graphviz 79 ; * externalize function 80 ; 29 81 ; - fplod 20110411T142955Z aedon.locean-ipsl.upmc.fr (Darwin) 30 82 ; … … 33 85 ;- 34 86 pro fig2_timeline_diagram 35 @common 36 ;----------------------------------------------- 87 88 @cm_tropflux 89 37 90 reinitplt, /z,/invert 38 91 key_portrait = 1 … … 40 93 41 94 openps, FILENAME = 'idl.ps' 42 ;------------------------------------------------------------ 43 ; partie a changer 44 ;------------------------------------------------------------ 45 ;; site locations 46 95 96 ; site locations 47 97 sitelist=['8s67e','12s55e', '8s55e', '8s80.5e', '1.5s80.5e', '0n80.5e', '1.5n80.5e', '1.5s90e', $ 48 98 '0n90e', '1.5n90e', '4n90e','8n90e','12n90e', '15n90e', '5s95e', $ … … 56 106 '21n23w', '4n23w', '4n38w', '6s10w', '8n38w', '8s30w'] 57 107 58 59 108 nn=n_elements(sitelist) 60 109 date1=19890101 … … 77 126 78 127 nsmooth=1. 79 read_lh, csite,date1,date2,nsmooth, $ 80 lh 81 82 ind=where(finite(lh)) & lh(ind)=1. 83 ind=where(finite(lh,/nan)) & lh(ind)=0. 128 print, csite 129 print, date1 130 print, date2 131 read_lh, csite, date1, date2, nsmooth, lh 132 133 ind=where(finite(lh)) & lh[ind]=1. 134 ind=where(finite(lh,/nan)) & lh[ind]=0. 84 135 no_obs_lh=no_obs_lh+lh 85 136 … … 113 164 114 165 nsmooth=1. 115 read_sw, csite, date1, date2, box, $ 116 sw 166 read_sw, csite, date1, date2, box, sw 117 167 swr=sw & ind=where(finite(swr)) & valid=n_elements(ind) 118 168 if (valid ge 10) then begin 119 ind=where(finite(sw)) & sw (ind)=1.120 ind=where(finite(sw,/nan)) & sw (ind)=0.169 ind=where(finite(sw)) & sw[ind]=1. 170 ind=where(finite(sw,/nan)) & sw[ind]=0. 121 171 no_obs_sw=no_obs_sw+sw 122 172 endif … … 128 178 129 179 nn=n_elements(sitelist) 130 date1=date1131 date2=date2132 133 180 134 181 for n=0, nn-1 do begin … … 144 191 145 192 nsmooth=1. 146 read_lw, csite, date1, date2, box, $ 147 lw 193 read_lw, csite, date1, date2, box, lw 148 194 lwr=lw & ind=where(finite(lwr)) & valid=n_elements(ind) 149 195 if (valid ge 10) then begin 150 151 ind=where(finite(lw)) & lw(ind)=1. 152 ind=where(finite(lw,/nan)) & lw(ind)=0. 196 ind=where(finite(lw)) & lw[ind]=1. 197 ind=where(finite(lw,/nan)) & lw[ind]=0. 153 198 no_obs_lw=no_obs_lw+lw 154 199 endif … … 163 208 pltt, no_obs_lw, "t",/ov1d, color=50 164 209 pltt, no_obs_lh, "t",/ov1d 165 ;------------------------------------------------------------166 210 closeps 167 fig='fig2_timeline_diagram.ps' 168 spawn, 'mv '+psdir+'idl.ps '+cpsdir+fig 169 spawn, 'gv '+cpsdir+fig 170 ;------------------------------------------------------------ 171 return 211 fig=tropflux_od_env + 'fig2_timeline_diagram.ps' 212 spawn, 'mv '+psdir+'idl.ps '+fig 213 spawn, 'gv '+fig 172 214 end 173 174 175 function x_site_location, site176 n1=strpos(site, 's')177 if (n1 gt -1) then begin178 ns=-1.179 x=strmid(site, 0, n1)180 x=float(x)*ns181 endif else begin182 n1=strpos(site, 'n')183 x=strmid(site, 0, n1)184 ny=1.185 x=float(x)*ny186 endelse187 return, float(x)188 end189 ;--------------------------------------------------------------------------190 function y_site_location, site191 n1=strpos(site, 'e')192 if (n1 gt -1) then begin193 n=strpos(site, 's')194 if (n gt -1) then begin195 y=strmid(site, n+1, n1-n-1)196 endif else begin197 n=strpos(site, 'n')198 y=strmid(site, n+1, n1-n-1)199 endelse200 201 endif else begin202 n1=strpos(site, 'w')203 n=strpos(site, 's')204 if (n gt -1) then begin205 y=strmid(site, n+1, n1-n-1)206 y=180+(180-float(y))207 endif else begin208 n=strpos(site, 'n')209 y=strmid(site, n+1, n1-n-1)210 y=180+(180-float(y))211 endelse212 endelse213 return,float(y)214 end -
trunk/src/paper01/fig2/read_lh.pro
r43 r44 1 1 ;+ 2 ; .. _read_lh.pro 2 ; .. _read_lh.pro: 3 3 ; 4 4 ; =========== … … 9 9 ; =========== 10 10 ; 11 ; .. graphviz:: 12 ; 13 ; digraph read_lh { 14 ; graph [ 15 ; rankdir="LR", 16 ; ] 17 ; qlat [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/qlat*_dy.cdf"]; 18 ; 19 ; lh [shape=tab,fontname=Courier,label="lh"] 20 ; 21 ; read_lh [shape=box, 22 ; fontname=Courier, 23 ; color=blue, 24 ; URL="http://forge.ipsl.jussieu.fr/tropflux/browser/trunk/src/paper01/fig2/read_lh.pro", 25 ; label="${TROPFLUX}/src/paper01/fig2/read_lh.pro"]; 26 ; 27 ; {qlat} -> {read_lh} -> {lh} 28 ; } 29 ; 11 30 ; SEE ALSO 12 31 ; ======== 13 32 ; 14 33 ; :ref:`tropflux_profile.sh` 34 ; :ref:`tropflux_init.pro` 35 ; :ref:`cm_tropflux.pro` 36 ; 37 ; :func:`time_lec` 38 ; 39 ; :ref:`fig2_timeline_diagram.pro` 15 40 ; 16 41 ; EXAMPLES … … 20 45 ; 21 46 ; IDL> @tropflux_init 22 ; IDL> csite=++ 23 ; IDL> ++ 47 ; IDL> csite='8s67e' 48 ; IDL> date1=19890101L 49 ; IDL> date2=20091231L 50 ; IDL> nsmooth=1. 24 51 ; IDL> read_lh, csite, date1, date2, nsmooth, lh 52 ; 53 ; TODO 54 ; ==== 55 ; 56 ; coding rules 25 57 ; 26 58 ; EVOLUTIONS … … 28 60 ; 29 61 ; $Id$ 62 ; 63 ; - fplod 20110412T093201Z aedon.locean-ipsl.upmc.fr (Darwin) 64 ; 65 ; * externalize function 66 ; * remove hard coding path 67 ; * add graphviz 30 68 ; 31 69 ; - fplod 20110411T142955Z aedon.locean-ipsl.upmc.fr (Darwin) … … 36 74 pro read_lh, csite, date1, date2, nsmooth, lh 37 75 ; 38 @c ommon76 @cm_tropflux 39 77 40 dir='/Users/pkb/data/TPR/'41 ;42 78 ;; DEFINE THE OUTPUT TIME AXIS 43 79 ; … … 55 91 56 92 ;;LHF 57 fi= dir+'qlat'+csite+'_dy.cdf'93 fi=tropflux_id_env+'qlat'+csite+'_dy.cdf' 58 94 f=file_test(fi) 59 95 lh=fltarr(jpt)+!values.f_nan … … 66 102 ind=where((tt0 ge -1e-5) and (tt0 le jpt-1+1e-5)) 67 103 if (ind(0) ne -1) then begin 68 lh(tt0(ind))=lh0(ind)69 lh_q(tt0(ind))=lh_q0(ind)104 lh(tt0(ind))=lh0(ind) 105 lh_q(tt0(ind))=lh_q0(ind) 70 106 endif 71 endif 107 endif else begin 108 print, fi, ' not found' 109 endelse 72 110 73 111 ind=where(lh_q ne 1 and lh_q ne 2) … … 93 131 94 132 end 95 96 97 function time_lec, fi98 tt=ncdf_lec(fi,var='time')99 fid=ncdf_open(fi) & vid=ncdf_varid(fid,'time')100 ncdf_attget, fid,vid,'units',orig101 orig=string(orig)102 yy=long(strmid(orig,11,4))103 mm=long(strmid(orig,16,2))104 dd=long(strmid(orig,19,2))105 tt=julday(mm,dd,yy,12,00)+tt106 return, tt107 end -
trunk/src/paper01/fig2/read_lw.pro
r43 r44 1 1 ;+ 2 ; .. _read_lw.pro 2 ; .. _read_lw.pro: 3 3 ; 4 4 ; =========== … … 9 9 ; =========== 10 10 ; 11 ; .. graphviz:: 12 ; 13 ; digraph read_lw { 14 ; graph [ 15 ; rankdir="LR", 16 ; ] 17 ; lwnet [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/lwnet*_dy.cdf"]; 18 ; lw [shape=tab,fontname=Courier,label="lw"] 19 ; 20 ; read_lw [shape=box, 21 ; fontname=Courier, 22 ; color=blue, 23 ; URL="http://forge.ipsl.jussieu.fr/tropflux/browser/trunk/src/paper01/fig2/read_lw.pro", 24 ; label="${TROPFLUX}/src/paper01/fig2/read_lw.pro"]; 25 ; 26 ; {lwnet} -> {read_lw} -> {lw} 27 ; } 28 ; 29 ; 11 30 ; SEE ALSO 12 31 ; ======== 13 32 ; 14 33 ; :ref:`tropflux_profile.sh` 34 ; 35 ; :func:`time_lec` 36 ; 37 ; :ref:`fig2_timeline_diagram.pro` 15 38 ; 16 39 ; EXAMPLES … … 20 43 ; 21 44 ; IDL> @tropflux_init 22 ; IDL> csite=++ 23 ; IDL> ++ 45 ; IDL> csite='8s67e' 46 ; IDL> date1=19890101L 47 ; IDL> date2=20091231L 48 ; IDL> box=[ 249.500,250.500,-0.500000,0.500000] 24 49 ; IDL> read_lw, csite, date1, date2, box, lw 50 ; 51 ; TODO 52 ; ==== 53 ; 54 ; make it work with a file : now nsmooth pb 55 ; 56 ; confirm example parameters values 25 57 ; 26 58 ; EVOLUTIONS … … 28 60 ; 29 61 ; $Id$ 62 ; 63 ; - fplod 20110412T093250Z aedon.locean-ipsl.upmc.fr (Darwin) 64 ; 65 ; * externalize function 66 ; * remove hard coding path 67 ; * add graphviz 68 ; * remove "smooth" lines (no nsmooth parameter here) 30 69 ; 31 70 ; - fplod 20110411T142955Z aedon.locean-ipsl.upmc.fr (Darwin) … … 36 75 pro read_lw, csite, date1, date2, box, lw 37 76 ; 38 @c ommon77 @cm_tropflux 39 78 40 79 dir='/Users/pkb/data/TPR/' … … 54 93 55 94 ;; LWR 56 fi= dir+'lwnet'+csite+'_dy.cdf'95 fi=tropflux_id_env+'lwnet'+csite+'_dy.cdf' 57 96 f=file_test(fi) 58 97 lw=fltarr(jpt)+!values.f_nan … … 65 104 lw(tt0(ind))=lw0(ind) 66 105 endif 67 endif 106 endif else begin 107 print, fi, ' not found' 108 endelse 68 109 69 110 … … 83 124 endfor 84 125 85 nsmooth=nsmooth86 lw=smooth(lw,nsmooth,/nan)126 ;++nsmooth=nsmooth 127 ;++lw=smooth(lw,nsmooth,/nan) 87 128 88 129 end 89 90 91 function time_lec, fi92 tt=ncdf_lec(fi,var='time')93 fid=ncdf_open(fi) & vid=ncdf_varid(fid,'time')94 ncdf_attget, fid,vid,'units',orig95 orig=string(orig)96 yy=long(strmid(orig,11,4))97 mm=long(strmid(orig,16,2))98 dd=long(strmid(orig,19,2))99 tt=julday(mm,dd,yy,12,00)+tt100 return, tt101 end -
trunk/src/paper01/fig2/read_sw.pro
r43 r44 1 1 ;+ 2 ; .. _read_sw.pro 2 ; 3 ; .. _read_sw.pro: 3 4 ; 4 5 ; =========== … … 9 10 ; =========== 10 11 ; 12 ; .. graphviz:: 13 ; 14 ; digraph read_sw { 15 ; graph [ 16 ; rankdir="LR", 17 ; ] 18 ; swnet [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/swnet*_dy.cdf"]; 19 ; 20 ; sw [shape=tab,fontname=Courier,label="sw"] 21 ; 22 ; read_sw [shape=box, 23 ; fontname=Courier, 24 ; color=blue, 25 ; URL="http://forge.ipsl.jussieu.fr/tropflux/browser/trunk/src/paper01/fig2/read_sw.pro", 26 ; label="${TROPFLUX}/src/paper01/fig2/read_sw.pro"]; 27 ; 28 ; {swnet} -> {read_sw} -> {sw} 29 ; } 30 11 31 ; SEE ALSO 12 32 ; ======== 13 33 ; 14 34 ; :ref:`tropflux_profile.sh` 35 ; 36 : :func:`time_lec` 37 ; 38 ; :ref:`fig2_timeline_diagram.pro` 15 39 ; 16 40 ; EXAMPLES … … 20 44 ; 21 45 ; IDL> @tropflux_init 22 ; IDL> csite=++ 23 ; IDL> ++ 46 ; IDL> csite='8s67e' 47 ; IDL> date1=19890101L 48 ; IDL> date2=20091231L 49 ; IDL> box=++ 24 50 ; IDL> read_sw, csite, date1, date2, box, sw 51 ; 52 ; TODO 53 ; ==== 54 ; 55 ; make it work with a file 56 ; 25 57 ; 26 58 ; EVOLUTIONS … … 28 60 ; 29 61 ; $Id$ 62 ; 63 ; - fplod 20110412T093325Z aedon.locean-ipsl.upmc.fr (Darwin) 64 ; 65 ; * externalize function 66 ; * remove hard coding path 67 ; * add graphviz 68 ; * remove smooth call (unreacheable statement after end !!) 30 69 ; 31 70 ; - fplod 20110411T142955Z aedon.locean-ipsl.upmc.fr (Darwin) … … 36 75 pro read_sw, csite, date1, date2, box, sw 37 76 ; 38 @common 39 40 dir='/Users/pkb/data/TPR/' 41 77 @cm_tropflux 42 78 ; 43 79 ;; DEFINE THE OUTPUT TIME AXIS … … 55 91 56 92 ;; SWR 57 fi= dir+'swnet'+csite+'_dy.cdf'93 fi=tropflux_id_env+'swnet'+csite+'_dy.cdf' 58 94 f=file_test(fi) 59 95 sw=fltarr(jpt)+!values.f_nan ;; swr … … 70 106 sw_q(tt0(ind))=sw_q0(ind) 71 107 endif 72 endif 108 endif else begin 109 print, fi, ' not found' 110 endelse 73 111 74 112 ind1=where((sw_q ne 1) and (sw_q ne 2)) … … 90 128 r=execute(com) 91 129 endfor 130 92 131 end 93 132 94 nsmooth=nsmooth95 sw=smooth(sw,nsmooth,/nan)96 97 98 99 function time_lec, fi100 tt=ncdf_lec(fi,var='time')101 fid=ncdf_open(fi) & vid=ncdf_varid(fid,'time')102 ncdf_attget, fid,vid,'units',orig103 orig=string(orig)104 yy=long(strmid(orig,11,4))105 mm=long(strmid(orig,16,2))106 dd=long(strmid(orig,19,2))107 tt=julday(mm,dd,yy,12,00)+tt108 return, tt109 end -
trunk/src/paper01/fig2/time_lec.pro
r43 r44 1 1 ;+ 2 ; .. _read_lh.pro 2 ; 3 ; ============ 4 ; time_lec.pro 5 ; ============ 6 ; 7 ; .. function:: time_lec(fi) 3 8 ; 4 ; =========== 5 ; read_lh.pro 6 ; =========== 7 ; 8 ; DESCRIPTION 9 ; =========== 10 ; 11 ; SEE ALSO 12 ; ======== 13 ; 14 ; :ref:`tropflux_profile.sh` 9 ; :param fi: file 15 10 ; 16 11 ; EXAMPLES … … 19 14 ; :: 20 15 ; 21 ; IDL> @tropflux_init 22 ; IDL> csite=++ 23 ; IDL> ++ 24 ; IDL> read_lh, csite, date1, date2, nsmooth, lh 16 ; IDL> fi=++ 17 ; IDL> result=time_lec(fi) 18 ; 19 ; SEE ALSO 20 ; ======== 21 ; 22 ; :ref:`read_lh.pro` 23 ; :ref:`read_lw.pro` 24 ; :ref:`read_sw.pro` 25 ; 26 ; TODO 27 ; ==== 28 ; 29 ; complete header 30 ; 31 ; coding rules 25 32 ; 26 33 ; EVOLUTIONS … … 29 36 ; $Id$ 30 37 ; 31 ; - fplod 2011041 1T142955Z aedon.locean-ipsl.upmc.fr (Darwin)38 ; - fplod 20110412T160108Z aedon.locean-ipsl.upmc.fr (Darwin) 32 39 ; 40 ; * externalized from read_lh.pro 33 41 ; * minimal header 34 42 ; 35 43 ;- 36 pro read_lh, csite, date1, date2, nsmooth, lh37 ;38 @common39 40 dir='/Users/pkb/data/TPR/'41 ;42 ;; DEFINE THE OUTPUT TIME AXIS43 ;44 jda1=date2jul(date1)45 jda2=date2jul(date2)46 jpt=(jda2-jda1+1l)47 time=jda1+dindgen(jpt)48 49 ;50 ;; FIRST READ ALL DATA FROM SITE51 ;; - if file exists, extract correct time axis52 ;; - if not, fill variable with missing values53 ;54 55 56 ;;LHF57 fi=dir+'qlat'+csite+'_dy.cdf'58 f=file_test(fi)59 lh=fltarr(jpt)+!values.f_nan60 lh_q=fltarr(jpt)+!values.f_nan61 if (f) then begin62 tt0=time_lec(fi)63 lh0=ncdf_lec(fi,var='QL_137')64 lh_q0=ncdf_lec(fi,var='QQL_5137')65 tt0=tt0-time(0)66 ind=where((tt0 ge -1e-5) and (tt0 le jpt-1+1e-5))67 if (ind(0) ne -1) then begin68 lh(tt0(ind))=lh0(ind)69 lh_q(tt0(ind))=lh_q0(ind)70 endif71 endif72 73 ind=where(lh_q ne 1 and lh_q ne 2)74 75 76 if (ind(0) ne -1) then begin77 lh(ind)=!Values.f_nan78 endif79 ;80 ;; Replace missing values by "NaN"81 ;82 tsvars=['lh']83 vars=[tsvars]84 nn=n_elements(vars)85 for n=0,nn-1 do begin86 var=vars(n)87 com='ind=where('+var+' ge 1.e20) & if (ind(0) ne -1) then '+var+'(ind)=!values.f_nan'88 r=execute(com)89 endfor90 91 nsmooth=nsmooth92 lh=smooth(lh,nsmooth,/nan)93 94 end95 96 97 44 function time_lec, fi 98 45 tt=ncdf_lec(fi,var='time') -
trunk/src/paper01/fig2/x_site_location.pro
r43 r44 1 1 ;+ 2 ; .. _fig2_timeline_diagram.pro:3 2 ; 4 ; =================== ======5 ; fig2_timeline_diagram.pro6 ; =================== ======3 ; =================== 4 ; x_site_location.pro 5 ; =================== 7 6 ; 8 ; DESCRIPTION 7 ; .. function:: x_site_location(site) 8 ; 9 ; DESCRITPION 9 10 ; =========== 11 ; 12 ; :param site: ++ 13 ; EXAMPLES 14 ; ======== 15 ; 16 ; :: 17 ; 18 ; IDL> site=++ 19 ; IDL> result=x_site_location(site) 10 20 ; 11 21 ; SEE ALSO 12 22 ; ======== 13 23 ; 14 ; :ref:` tropflux_profile.sh`24 ; :ref:`fig2_timeline_diagram.pro` 15 25 ; 16 ; EXAMPLES17 ; ==== ====26 ; TODO 27 ; ==== 18 28 ; 19 ; :: 20 ; 21 ; IDL> @tropflux_init 22 ; IDL> fig2_timeline_diagram 23 ; 29 ; coding rules 30 ; 24 31 ; EVOLUTIONS 25 32 ; ========== … … 27 34 ; $Id$ 28 35 ; 29 ; - fplod 2011041 1T142955Z aedon.locean-ipsl.upmc.fr (Darwin)30 ; 31 ; * minimal header36 ; - fplod 20110412T083358Z aedon.locean-ipsl.upmc.fr (Darwin) 37 ; 38 ; * creation from fig2_timeline_diagram.pro 32 39 ; 33 40 ;- 34 pro fig2_timeline_diagram 35 @common 36 ;----------------------------------------------- 37 reinitplt, /z,/invert 38 key_portrait = 1 39 coefpalit=.9 41 function x_site_location, site 40 42 41 openps, FILENAME = 'idl.ps' 42 ;------------------------------------------------------------ 43 ; partie a changer 44 ;------------------------------------------------------------ 45 ;; site locations 46 47 sitelist=['8s67e','12s55e', '8s55e', '8s80.5e', '1.5s80.5e', '0n80.5e', '1.5n80.5e', '1.5s90e', $ 48 '0n90e', '1.5n90e', '4n90e','8n90e','12n90e', '15n90e', '5s95e', $ 49 '8s165e', '8s180w', '8s155w', '8s125w', '8s110w', '8s95w', '5s156e', '5s165e', '5s180w', '5s170w', $ 50 '5s155w', '5s140w', '5s125w', '5s110w', '5s95w', '2s156e', '2s165e', '2s180w', '2s170w', '2s155w', '2s140w', $ 51 '2s125w', '2s110w', '2s95w', '0n147e', '0n156e', '0n165e', '0n180w', '0n170w', '0n155w', '0n140w', '0n125w', $ 52 '0n110w', '0n95w', '2n147e', '2n156e', '2n165e', '2n180w', '2n170w', '2n155w', '2n140w', '2n125w', '2n110w', $ 53 '2n95w', '5n147e', '5n156e', '5n165e', '5n170w', '5n155w', '5n140w', '5n125w', '5n110w', '5n95w', $ 54 '8n156e', '8n165e', '8n180w', '8n170w', '9n140w', '8n125w', '8n110w', '8n95w', $ 55 '0n0e', '0n10w', '0n23w', '0n35w', '10s10w', '12n23w', '12n38w', '14s32w', '15n38w', '19s34w', '20n38w', $ 56 '21n23w', '4n23w', '4n38w', '6s10w', '8n38w', '8s30w'] 57 58 59 nn=n_elements(sitelist) 60 date1=19890101 61 date2=20091231 62 63 no_obs_lh=findgen(7670.)*0. 64 no_obs_lw=findgen(7670.)*0. 65 no_obs_sw=findgen(7670.)*0. 66 67 for n=0, nn-1 do begin 68 69 ;; reading data from mooring 70 71 site=sitelist(n) & csite=site 72 print, csite 73 x=x_site_location(site) 74 y=y_site_location(site) 75 if (y ge 0. and y le 30.) then y=y+360. 76 dx=0.5 & dy=0.5 & box=[y-dy, y+dy, x-dx, x+dx] 77 78 nsmooth=1. 79 read_lh, csite,date1,date2,nsmooth, $ 80 lh 81 82 ind=where(finite(lh)) & lh(ind)=1. 83 ind=where(finite(lh,/nan)) & lh(ind)=0. 84 no_obs_lh=no_obs_lh+lh 85 86 endfor 87 88 sitelist=['5n165e','8s67e','12s55e', '8s55e', '8s80.5e', '0n80.5e', '1.5n80.5e', '1.5s90e', $ 89 '0n90e', '1.5n90e', '4n90e','8n90e','12n90e', '15n90e', '5s95e', $ 90 '8s165e', '8s125w', '8s110w', '8s95w', '5s156e', '5s165e', '5s180w', '5s170w', $ 91 '5s155w', '5s140w', '5s125w', '5s110w', '5s95w', '2s156e', '2s165e', '2s180w', '2s170w', '2s155w', '2s140w', $ 92 '2s125w', '2s110w', '2s95w', '0n147e', '0n156e', '0n165e', '0n180w', '0n170w', '0n155w', '0n140w', '0n125w', $ 93 '0n110w', '0n95w', '2n147e', '2n156e', '2n165e', '2n180w', '2n170w', '2n155w', '2n140w', '2n125w', '2n110w', $ 94 '2n95w', '5n147e', '5n156e', '5n170w', '5n155w', '5n140w', '5n125w', '5n110w', '5n95w', $ 95 '8n156e', '8n165e', '8n180w', '8n170w', '9n140w', '8n125w', '8n110w', '8n95w', $ 96 '0n0e', '0n10w', '0n23w', '0n35w', '10s10w', '12n23w', '12n38w', '14s32w', '15n38w', '19s34w', '20n38w', $ 97 '21n23w', '4n23w', '4n38w', '6s10w', '8n38w', '8s30w'] 98 99 nn=n_elements(sitelist) 100 date1=date1 101 date2=date2 102 103 for n=0, nn-1 do begin 104 105 ;; reading data from mooring 106 107 site=sitelist(n) & csite=site 108 print, csite 109 x=x_site_location(site) 110 y=y_site_location(site) 111 if (y ge 0. and y le 30.) then y=y+360. 112 dx=0.5 & dy=0.5 & box=[y-dy, y+dy, x-dx, x+dx] 113 114 nsmooth=1. 115 read_sw, csite, date1, date2, box, $ 116 sw 117 swr=sw & ind=where(finite(swr)) & valid=n_elements(ind) 118 if (valid ge 10) then begin 119 ind=where(finite(sw)) & sw(ind)=1. 120 ind=where(finite(sw,/nan)) & sw(ind)=0. 121 no_obs_sw=no_obs_sw+sw 122 endif 123 endfor 124 help, no_obs_sw 125 126 sitelist=['0n110w', '0n140w', '0n165e', '0n170w', '0n23w', '0n80.5e', '0n95w', '10s10w', '12n23w', '15n38w', $ 127 '15n90e', '2n95w', '2s95w', '5n95w', '5s95w', '8n95w', '8s67e', '8s95w'] 128 129 nn=n_elements(sitelist) 130 date1=date1 131 date2=date2 132 133 134 for n=0, nn-1 do begin 135 136 ;; reading data from mooring 137 138 site=sitelist(n) & csite=site 139 print, csite 140 x=x_site_location(site) 141 y=y_site_location(site) 142 if (y ge 0. and y le 30.) then y=y+360. 143 dx=0.5 & dy=0.5 & box=[y-dy, y+dy, x-dx, x+dx] 144 145 nsmooth=1. 146 read_lw, csite, date1, date2, box, $ 147 lw 148 lwr=lw & ind=where(finite(lwr)) & valid=n_elements(ind) 149 if (valid ge 10) then begin 150 151 ind=where(finite(lw)) & lw(ind)=1. 152 ind=where(finite(lw,/nan)) & lw(ind)=0. 153 no_obs_lw=no_obs_lw+lw 154 endif 155 156 endfor 157 158 help,no_obs_lw 159 pltt, no_obs_lh, "t", small=[1,3,1],/rempl,0,80, $ 160 title='No of valid observations for turbulent flux validation', charsize=1., $ 161 ytitle="No. of observations",subtitle='black(turbulent), red(SWR), blue(LWR)' 162 pltt, no_obs_sw, "t",/ov1d, color=250 163 pltt, no_obs_lw, "t",/ov1d, color=50 164 pltt, no_obs_lh, "t",/ov1d 165 ;------------------------------------------------------------ 166 closeps 167 fig='fig2_timeline_diagram.ps' 168 spawn, 'mv '+psdir+'idl.ps '+cpsdir+fig 169 spawn, 'gv '+cpsdir+fig 170 ;------------------------------------------------------------ 171 return 172 end 173 174 175 function x_site_location, site 176 n1=strpos(site, 's') 43 n1=strpos(site, 's') 177 44 if (n1 gt -1) then begin 178 45 ns=-1. … … 187 54 return, float(x) 188 55 end 189 ;--------------------------------------------------------------------------190 function y_site_location, site191 n1=strpos(site, 'e')192 if (n1 gt -1) then begin193 n=strpos(site, 's')194 if (n gt -1) then begin195 y=strmid(site, n+1, n1-n-1)196 endif else begin197 n=strpos(site, 'n')198 y=strmid(site, n+1, n1-n-1)199 endelse200 201 endif else begin202 n1=strpos(site, 'w')203 n=strpos(site, 's')204 if (n gt -1) then begin205 y=strmid(site, n+1, n1-n-1)206 y=180+(180-float(y))207 endif else begin208 n=strpos(site, 'n')209 y=strmid(site, n+1, n1-n-1)210 y=180+(180-float(y))211 endelse212 endelse213 return,float(y)214 end -
trunk/src/paper01/fig2/y_site_location.pro
r43 r44 1 1 ;+ 2 ; .. _fig2_timeline_diagram.pro:3 2 ; 4 ; =================== ======5 ; fig2_timeline_diagram.pro6 ; =================== ======3 ; =================== 4 ; y_site_location.pro 5 ; =================== 7 6 ; 8 ; DESCRIPTION 7 ; .. function:: y_site_location(site) 8 ; 9 ; DESCRITPION 9 10 ; =========== 10 11 ; 11 ; SEE ALSO 12 ; ======== 13 ; 14 ; :ref:`tropflux_profile.sh` 12 ; :param site: ++ 15 13 ; 16 14 ; EXAMPLES … … 18 16 ; 19 17 ; :: 18 ; 19 ; IDL> site=++ 20 ; IDL> result=y_site_location(site) 20 21 ; 21 ; IDL> @tropflux_init22 ; IDL> fig2_timeline_diagram22 ; SEE ALSO 23 ; ======== 23 24 ; 25 ; :ref:`fig2_timeline_diagram.pro` 26 ; 27 ; TODO 28 ; ==== 29 ; 30 ; coding rules 31 ; 24 32 ; EVOLUTIONS 25 33 ; ========== … … 27 35 ; $Id$ 28 36 ; 29 ; - fplod 2011041 1T142955Z aedon.locean-ipsl.upmc.fr (Darwin)30 ; 31 ; * minimal header37 ; - fplod 20110412T083358Z aedon.locean-ipsl.upmc.fr (Darwin) 38 ; 39 ; * creation from fig2_timeline_diagram.pro 32 40 ; 33 41 ;- 34 pro fig2_timeline_diagram35 @common36 ;-----------------------------------------------37 reinitplt, /z,/invert38 key_portrait = 139 coefpalit=.940 41 openps, FILENAME = 'idl.ps'42 ;------------------------------------------------------------43 ; partie a changer44 ;------------------------------------------------------------45 ;; site locations46 47 sitelist=['8s67e','12s55e', '8s55e', '8s80.5e', '1.5s80.5e', '0n80.5e', '1.5n80.5e', '1.5s90e', $48 '0n90e', '1.5n90e', '4n90e','8n90e','12n90e', '15n90e', '5s95e', $49 '8s165e', '8s180w', '8s155w', '8s125w', '8s110w', '8s95w', '5s156e', '5s165e', '5s180w', '5s170w', $50 '5s155w', '5s140w', '5s125w', '5s110w', '5s95w', '2s156e', '2s165e', '2s180w', '2s170w', '2s155w', '2s140w', $51 '2s125w', '2s110w', '2s95w', '0n147e', '0n156e', '0n165e', '0n180w', '0n170w', '0n155w', '0n140w', '0n125w', $52 '0n110w', '0n95w', '2n147e', '2n156e', '2n165e', '2n180w', '2n170w', '2n155w', '2n140w', '2n125w', '2n110w', $53 '2n95w', '5n147e', '5n156e', '5n165e', '5n170w', '5n155w', '5n140w', '5n125w', '5n110w', '5n95w', $54 '8n156e', '8n165e', '8n180w', '8n170w', '9n140w', '8n125w', '8n110w', '8n95w', $55 '0n0e', '0n10w', '0n23w', '0n35w', '10s10w', '12n23w', '12n38w', '14s32w', '15n38w', '19s34w', '20n38w', $56 '21n23w', '4n23w', '4n38w', '6s10w', '8n38w', '8s30w']57 58 59 nn=n_elements(sitelist)60 date1=1989010161 date2=2009123162 63 no_obs_lh=findgen(7670.)*0.64 no_obs_lw=findgen(7670.)*0.65 no_obs_sw=findgen(7670.)*0.66 67 for n=0, nn-1 do begin68 69 ;; reading data from mooring70 71 site=sitelist(n) & csite=site72 print, csite73 x=x_site_location(site)74 y=y_site_location(site)75 if (y ge 0. and y le 30.) then y=y+360.76 dx=0.5 & dy=0.5 & box=[y-dy, y+dy, x-dx, x+dx]77 78 nsmooth=1.79 read_lh, csite,date1,date2,nsmooth, $80 lh81 82 ind=where(finite(lh)) & lh(ind)=1.83 ind=where(finite(lh,/nan)) & lh(ind)=0.84 no_obs_lh=no_obs_lh+lh85 86 endfor87 88 sitelist=['5n165e','8s67e','12s55e', '8s55e', '8s80.5e', '0n80.5e', '1.5n80.5e', '1.5s90e', $89 '0n90e', '1.5n90e', '4n90e','8n90e','12n90e', '15n90e', '5s95e', $90 '8s165e', '8s125w', '8s110w', '8s95w', '5s156e', '5s165e', '5s180w', '5s170w', $91 '5s155w', '5s140w', '5s125w', '5s110w', '5s95w', '2s156e', '2s165e', '2s180w', '2s170w', '2s155w', '2s140w', $92 '2s125w', '2s110w', '2s95w', '0n147e', '0n156e', '0n165e', '0n180w', '0n170w', '0n155w', '0n140w', '0n125w', $93 '0n110w', '0n95w', '2n147e', '2n156e', '2n165e', '2n180w', '2n170w', '2n155w', '2n140w', '2n125w', '2n110w', $94 '2n95w', '5n147e', '5n156e', '5n170w', '5n155w', '5n140w', '5n125w', '5n110w', '5n95w', $95 '8n156e', '8n165e', '8n180w', '8n170w', '9n140w', '8n125w', '8n110w', '8n95w', $96 '0n0e', '0n10w', '0n23w', '0n35w', '10s10w', '12n23w', '12n38w', '14s32w', '15n38w', '19s34w', '20n38w', $97 '21n23w', '4n23w', '4n38w', '6s10w', '8n38w', '8s30w']98 99 nn=n_elements(sitelist)100 date1=date1101 date2=date2102 103 for n=0, nn-1 do begin104 105 ;; reading data from mooring106 107 site=sitelist(n) & csite=site108 print, csite109 x=x_site_location(site)110 y=y_site_location(site)111 if (y ge 0. and y le 30.) then y=y+360.112 dx=0.5 & dy=0.5 & box=[y-dy, y+dy, x-dx, x+dx]113 114 nsmooth=1.115 read_sw, csite, date1, date2, box, $116 sw117 swr=sw & ind=where(finite(swr)) & valid=n_elements(ind)118 if (valid ge 10) then begin119 ind=where(finite(sw)) & sw(ind)=1.120 ind=where(finite(sw,/nan)) & sw(ind)=0.121 no_obs_sw=no_obs_sw+sw122 endif123 endfor124 help, no_obs_sw125 126 sitelist=['0n110w', '0n140w', '0n165e', '0n170w', '0n23w', '0n80.5e', '0n95w', '10s10w', '12n23w', '15n38w', $127 '15n90e', '2n95w', '2s95w', '5n95w', '5s95w', '8n95w', '8s67e', '8s95w']128 129 nn=n_elements(sitelist)130 date1=date1131 date2=date2132 133 134 for n=0, nn-1 do begin135 136 ;; reading data from mooring137 138 site=sitelist(n) & csite=site139 print, csite140 x=x_site_location(site)141 y=y_site_location(site)142 if (y ge 0. and y le 30.) then y=y+360.143 dx=0.5 & dy=0.5 & box=[y-dy, y+dy, x-dx, x+dx]144 145 nsmooth=1.146 read_lw, csite, date1, date2, box, $147 lw148 lwr=lw & ind=where(finite(lwr)) & valid=n_elements(ind)149 if (valid ge 10) then begin150 151 ind=where(finite(lw)) & lw(ind)=1.152 ind=where(finite(lw,/nan)) & lw(ind)=0.153 no_obs_lw=no_obs_lw+lw154 endif155 156 endfor157 158 help,no_obs_lw159 pltt, no_obs_lh, "t", small=[1,3,1],/rempl,0,80, $160 title='No of valid observations for turbulent flux validation', charsize=1., $161 ytitle="No. of observations",subtitle='black(turbulent), red(SWR), blue(LWR)'162 pltt, no_obs_sw, "t",/ov1d, color=250163 pltt, no_obs_lw, "t",/ov1d, color=50164 pltt, no_obs_lh, "t",/ov1d165 ;------------------------------------------------------------166 closeps167 fig='fig2_timeline_diagram.ps'168 spawn, 'mv '+psdir+'idl.ps '+cpsdir+fig169 spawn, 'gv '+cpsdir+fig170 ;------------------------------------------------------------171 return172 end173 174 175 function x_site_location, site176 n1=strpos(site, 's')177 if (n1 gt -1) then begin178 ns=-1.179 x=strmid(site, 0, n1)180 x=float(x)*ns181 endif else begin182 n1=strpos(site, 'n')183 x=strmid(site, 0, n1)184 ny=1.185 x=float(x)*ny186 endelse187 return, float(x)188 end189 ;--------------------------------------------------------------------------190 42 function y_site_location, site 191 43 n1=strpos(site, 'e') 192 44 if (n1 gt -1) then begin 193 45 n=strpos(site, 's') … … 198 50 y=strmid(site, n+1, n1-n-1) 199 51 endelse 200 201 52 endif else begin 202 53 n1=strpos(site, 'w') -
trunk/src/paper01/fig3/air_validation_scatter_2000_2009_v50.pro
r43 r44 1 1 ;+ 2 ; _air_validation_scatter_2000_2009_v50.pro: 2 ; 3 ; .. _air_validation_scatter_2000_2009_v50.pro: 3 4 ; 4 5 ; ======================================== … … 9 10 ; =========== 10 11 ; 12 ; .. graphviz:: 13 ; 14 ; digraph air_validation_scatter_2000_2009_v50 { 15 ; graph [ 16 ; rankdir="LR", 17 ; ] 18 ; 19 ; air_erai [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/air_2000_2009_erai_v50.txt"]; 20 ; 21 ; air_trop [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/air_2000_2009_trop_v50.txt"]; 22 ; air_oaflux [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/air_2000_2009_oaflx_v50.txt"]; 23 ; air_ncep [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/air_2000_2009_ncep_v50.txt"]; 24 ; air_ncep1 [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/air_2000_2009_ncep1_v50.txt"]; 25 ; air_ncep2 [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/air_2000_2009_ncep2_v50.txt"]; 26 ; air_ncep2_oafluxgrid [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/air_2m_ncep2_oafluxgrid_19890101_20091231.nc"]; 27 ; 28 ; t2m_ncep1 [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/t2m_ncep1_19890101_20091231.nc"]; 29 ; figure [shape=ellipse,fontname=Courier,label="${TROPFLUX_OD}/air_validation_scatter_2000_2009_v50.ps"]; 30 ; 31 ; erai_t2m [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/erai_t2m_19890101_20091231_oafluxgrid.nc"]; 32 ; 33 ; tropflux_t2m [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/TropFlux_t2m_19890101_20091231_v50.nc"]; 34 ; 35 ; oaflux_basic [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/OAFlux_basic_variables_1985_2009.nc"]; 36 ; 37 ; air_validation_scatter_2000_2009_v50 [shape=box, 38 ; fontname=Courier, 39 ; color=blue, 40 ; URL="http://forge.ipsl.jussieu.fr/tropflux/browser/trunk/src/paper01/fig3/air_validation_scatter_2000_2009_v50.pro", 41 ; label="${TROPFLUX}/src/paper01/fig3/air_validation_scatter_2000_2009_v50.pro"]; 42 ; 43 ; {air_erai air_trop air_oaflux air_ncep air_ncep1 air_ncep2 erai_t2m tropflux_t2m oaflux_basic air_ncep2_oafluxgrid t2m_ncep1} -> {air_validation_scatter_2000_2009_v50} -> {figure} 44 ; 45 ; } 46 ; 11 47 ; SEE ALSO 12 48 ; ======== … … 14 50 ; :ref:`tropflux_profile.sh` 15 51 ; 52 ; :ref:`read_variables_v2.pro` 53 ; 16 54 ; EXAMPLES 17 55 ; ======== … … 20 58 ; 21 59 ; IDL> @tropflux_init 22 ; IDL> date1= ++23 ; IDL> date2= ++60 ; IDL> date1=20001001L 61 ; IDL> date2=20091231L 24 62 ; IDL> air_validation_scatter_2000_2009_v50, date1, date2 25 63 ; 64 ; TODO 65 ; ==== 66 ; 67 ; make it work (missing data) 68 ; 69 ; Praveen confirm date selection 70 ; 71 ; coding rules 72 ; 73 ; get rid of spwan 74 ; 75 ; complete description 76 ; 77 ; externalize function 78 ; 79 ; handle IO error 26 80 ; 27 81 ; EVOLUTIONS … … 30 84 ; $Id$ 31 85 ; 86 ; - fplod 20110412T101323Z aedon.locean-ipsl.upmc.fr (Darwin) 87 ; 88 ; * remove hard coding path 89 ; * add graphviz 90 ; 32 91 ; - fplod 20110411T142955Z aedon.locean-ipsl.upmc.fr (Darwin) 33 92 ; … … 35 94 ; 36 95 ;- 37 38 96 pro air_validation_scatter_2000_2009_v50, date1, date2 39 @common 97 98 @cm_tropflux 40 99 41 100 reinitplt, /z,/invert … … 68 127 close,/all 69 128 70 fi_air_erai= '/Users/pkb/work/MY_SAXO/flux_automat/air_2000_2009_erai_v50.txt'129 fi_air_erai=tropflux_id_env+'air_2000_2009_erai_v50.txt' 71 130 openw,1,fi_air_erai 72 fi_air_trop= '/Users/pkb/work/MY_SAXO/flux_automat/air_2000_2009_trop_v50.txt'131 fi_air_trop=tropflux_id_env+'air_2000_2009_trop_v50.txt' 73 132 openw,2,fi_air_trop 74 fi_air_oaflx= '/Users/pkb/work/MY_SAXO/flux_automat/air_2000_2009_oaflx_v50.txt'133 fi_air_oaflx=tropflux_id_env+'air_2000_2009_oaflx_v50.txt' 75 134 openw,3,fi_air_oaflx 76 fi_air_ncep2= '/Users/pkb/work/MY_SAXO/flux_automat/air_2000_2009_ncep2_v50.txt'135 fi_air_ncep2=tropflux_id_env+'air_2000_2009_ncep2_v50.txt' 77 136 openw,4,fi_air_ncep2 78 fi_air_ncep1= '/Users/pkb/work/MY_SAXO/flux_automat/air_2000_2009_ncep1_v50.txt'137 fi_air_ncep1=tropflux_id_env+'air_2000_2009_ncep1_v50.txt' 79 138 openw,5,fi_air_ncep1 80 139 … … 88 147 ;; first reading the whole ERAI uncorrected and corrected data 89 148 90 file= '/Volumes/Iomega_HDD/TropFlux/input_uncor/erai_t2m_19890101_20091231_oafluxgrid.nc'149 file=tropflux_id_env+'erai_t2m_19890101_20091231_oafluxgrid.nc' 91 150 initncdf, file 92 151 unc=read_ncdf('t2m',date1,date2,file=file,/nostr) & unc=unc-273.15 … … 94 153 95 154 ;file='/Volumes/Iomega_HDD/TropFlux/input_cor/full_cor/TropFlux_t2m_19890101_20091231_v3.nc' 96 file= '/Users/pkb/data/TropFlux/TropFlux_t2m_19890101_20091231_v50.nc'155 file=tropflux_id_env+'TropFlux_t2m_19890101_20091231_v50.nc' 97 156 initncdf, file 98 157 cor=read_ncdf('t2m',date1,date2,file=file,/nostr) & cor=cor-273.15 99 158 help, cor 100 159 101 file= '/Volumes/Iomega_HDD/work/flux_reconstruction/OAFLX_GRID/OAFlux_basic_variables_1985_2009.nc'160 file=tropflux_id_env+'OAFlux_basic_variables_1985_2009.nc' 102 161 initncdf, file 103 162 oaf=read_ncdf("air", date1, date2, file=file,/nostr) 104 163 help, oaf 105 164 106 fi= '/Volumes/Iomega_HDD/flux_reconstruction/ncep2/air_2m_ncep2_oafluxgrid_19890101_20091231.nc'165 fi=tropflux_id_env+'air_2m_ncep2_oafluxgrid_19890101_20091231.nc' 107 166 initncdf, fi 108 167 nce=read_ncdf("air", date1-1, date2, file=fi,/nostr) … … 110 169 help, nce 111 170 112 file= '/Volumes/Iomega_HDD/flux_reconstruction/ncep/t2m_ncep1_19890101_20091231.nc'171 file=tropflux_id_env+'t2m_ncep1_19890101_20091231.nc' 113 172 initncdf, file 114 173 nce1=read_ncdf("t2m", date1, date2, file=file,/nostr) … … 187 246 close,/all 188 247 189 fi_air_erai= '/Users/pkb/work/MY_SAXO/flux_automat/air_2000_2009_erai_v50.txt'248 fi_air_erai=tropflux_id_env+'air_2000_2009_erai_v50.txt' 190 249 res=read_ascii(fi_air_erai,data_start=1) 191 250 ff=res.field1 … … 213 272 oplot, mean_tao, yfit, color=250, thick=2 214 273 215 fi_air_trop= '/Users/pkb/work/MY_SAXO/flux_automat/air_2000_2009_trop_v50.txt'274 fi_air_trop=tropflux_id_env+'air_2000_2009_trop_v50.txt' 216 275 res=read_ascii(fi_air_trop,data_start=1) 217 276 ff=res.field1 … … 240 299 oplot, mean_tao, yfit, color=250, thick=2 241 300 242 fi_air_oaflx= '/Users/pkb/work/MY_SAXO/flux_automat/air_2000_2009_oaflx_v50.txt'301 fi_air_oaflx=tropflux_id_env+'air_2000_2009_oaflx_v50.txt' 243 302 res=read_ascii(fi_air_oaflx,data_start=1) 244 303 ff=res.field1 … … 269 328 270 329 271 fi_air_ncep= '/Users/pkb/work/MY_SAXO/flux_automat/air_2000_2009_ncep_v50.txt'330 fi_air_ncep=tropflux_id_env+'air_2000_2009_ncep_v50.txt' 272 331 res=read_ascii(fi_air_ncep,data_start=1) 273 332 ff=res.field1 … … 296 355 oplot, mean_tao, yfit, color=250, thick=2 297 356 298 fi_air_ncep1= '/Users/pkb/work/MY_SAXO/flux_automat/air_2000_2009_ncep1_v50.txt'357 fi_air_ncep1=tropflux_id_env+'air_2000_2009_ncep1_v50.txt' 299 358 res=read_ascii(fi_air_ncep1,data_start=1) 300 359 ff=res.field1 … … 325 384 closeps 326 385 327 fig= 'air_validation_scatter_2000_2009_v50.ps'328 spawn, 'mv '+psdir+'idl.ps '+ cpsdir+fig329 spawn, 'gv '+ cpsdir+fig386 fig=tropflux_od_env+'air_validation_scatter_2000_2009_v50.ps' 387 spawn, 'mv '+psdir+'idl.ps '+fig 388 spawn, 'gv '+fig 330 389 return 331 390 end -
trunk/src/paper01/fig3/q2m_validation_scatter_2000_2009_v50.pro
r43 r44 1 1 ;+ 2 ; .. _q2m_validation_scatter_2000_2009_v50.pro 2 ; 3 ; .. _q2m_validation_scatter_2000_2009_v50.pro: 3 4 ; 4 5 ; ======================================== … … 9 10 ; =========== 10 11 ; 12 ; .. graphviz:: 13 ; 14 ; digraph q2m_validation_scatter_2000_2009_v50 { 15 ; graph [ 16 ; rankdir="LR", 17 ; ] 18 ; 19 ; q2m_erai [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/q2m_2000_2009_erai_v50.txt"]; 20 ; 21 ; q2m_trop [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/q2m_2000_2009_trop_v50.txt"]; 22 ; q2m_oaflux [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/q2m_2000_2009_oaflx_v50.txt"]; 23 ; q2m_ncep [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/q2m_2000_2009_ncep_v50.txt"]; 24 ; q2m_ncep1 [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/q2m_2000_2009_ncep1_v50.txt"]; 25 ; 26 ; erai_q2m_oafluxgrid [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/erai_t2m_19890101_20091231_oafluxgrid.nc"]; 27 ; 28 ; tropflux_q2m [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/TropFlux_q2m_19890101_20091231_v50.nc"]; 29 ; 30 ; oaflux_basic [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/OAFlux_basic_variables_1985_2009.nc"]; 31 ; 32 ; sphum_ncep2_oafluxgrid [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/sphum_ncep2_oafluxgrid_19890101_20091231.nc"]; 33 ; 34 ; q2m_validation_scatter_2000_2009_v50 [shape=box, 35 ; fontname=Courier, 36 ; color=blue, 37 ; URL="http://forge.ipsl.jussieu.fr/tropflux/browser/trunk/src/paper01/fig3/q2m_validation_scatter_2000_2009_v50.pro", 38 ; label="${TROPFLUX}/src/paper01/fig3/q2m_validation_scatter_2000_2009_v50.pro"]; 39 ; 40 ; figure [shape=ellipse,fontname=Courier,label="${TROPFLUX_OD}/q2m_validation_scatter_2000_2009_v50.ps"]; 41 ; 42 ; {q2m_erai q2m_trop q2m_oaflux q2m_ncep q2m_ncep1 erai_q2m_oafluxgrid tropflux_q2m oaflux_basic sphum_ncep2_oafluxgrid} -> {q2m_validation_scatter_2000_2009_v50} -> {figure} 43 ; 44 ; } 45 ; 11 46 ; SEE ALSO 12 47 ; ======== … … 20 55 ; 21 56 ; IDL> @tropflux_init 22 ; IDL> date1= ++23 ; IDL> date2= ++57 ; IDL> date1=20001001L 58 ; IDL> date2=20091231L 24 59 ; IDL> q2m_validation_scatter_2000_2009_v50, date1, date2 60 ; 61 ; TODO 62 ; ==== 63 ; 64 ; make it work (missing data) 65 ; 66 ; Praveen confirm date selection 67 ; 68 ; coding rules 69 ; 70 ; get rid of spwan 71 ; 72 ; complete description 73 ; 74 ; externalize function 75 ; 76 ; handle IO error 25 77 ; 26 78 ; EVOLUTIONS … … 29 81 ; $Id$ 30 82 ; 83 ; - fplod 20110412T104125Z aedon.locean-ipsl.upmc.fr (Darwin) 84 ; 85 ; * remove hard coding path 86 ; * add graphviz 87 ; 31 88 ; - fplod 20110411T142955Z aedon.locean-ipsl.upmc.fr (Darwin) 32 89 ; … … 35 92 ;- 36 93 pro q2m_validation_scatter_2000_2009_v50, date1, date2 37 @common 94 95 @cm_tropflux 38 96 39 97 reinitplt, /z,/invert … … 66 124 close,/all 67 125 68 fi_q2m_erai=' /Users/pkb/work/MY_SAXO/flux_automat/q2m_2000_2009_erai_v50.txt'126 fi_q2m_erai='tropflux_id_env+'q2m_2000_2009_erai_v50.txt' 69 127 openw,1,fi_q2m_erai 70 fi_q2m_trop=' /Users/pkb/work/MY_SAXO/flux_automat/q2m_2000_2009_trop_v50.txt'128 fi_q2m_trop='tropflux_id_env+'q2m_2000_2009_trop_v50.txt' 71 129 openw,2,fi_q2m_trop 72 fi_q2m_oaflx=' /Users/pkb/work/MY_SAXO/flux_automat/q2m_2000_2009_oaflx_v50.txt'130 fi_q2m_oaflx='tropflux_id_env+'q2m_2000_2009_oaflx_v50.txt' 73 131 openw,3,fi_q2m_oaflx 74 fi_q2m_ncep=' /Users/pkb/work/MY_SAXO/flux_automat/q2m_2000_2009_ncep_v50.txt'132 fi_q2m_ncep='tropflux_id_env+'q2m_2000_2009_ncep_v50.txt' 75 133 openw,4,fi_q2m_ncep 76 fi_q2m_ncep1=' /Users/pkb/work/MY_SAXO/flux_automat/q2m_2000_2009_ncep1_v50.txt'134 fi_q2m_ncep1='tropflux_id_env+'q2m_2000_2009_ncep1_v50.txt' 77 135 openw,5,fi_q2m_ncep1 78 136 … … 86 144 ;; first reading the whole ERAI uncorrected and corrected data 87 145 88 file= '/Volumes/Iomega_HDD/TropFlux/input_uncor/erai_q2m_19890101_20091231_oafluxgrid.nc'146 file=tropflux_id_env+'erai_q2m_19890101_20091231_oafluxgrid.nc' 89 147 initncdf, file 90 148 unc=read_ncdf('q2m',date1,date2,file=file,/nostr) 91 149 help, unc 92 150 93 file= '/Volumes/Iomega_HDD/TropFlux/input_cor/full_cor/TropFlux_q2m_19890101_20091231_v20.nc'151 file=tropflux_id_env+'TropFlux_q2m_19890101_20091231_v20.nc' 94 152 initncdf, file 95 153 cor=read_ncdf('q2m',date1,date2,file=file,/nostr) 96 154 help, cor 97 155 98 file= '/Volumes/Iomega_HDD/work/flux_reconstruction/OAFLX_GRID/OAFlux_basic_variables_1985_2009.nc'156 file=tropflux_id_env+'OAFlux_basic_variables_1985_2009.nc' 99 157 initncdf, file 100 158 oaf=read_ncdf("sphum", date1, date2, file=file,/nostr) 101 159 help, oaf 102 160 103 fi= '/Volumes/Iomega_HDD/flux_reconstruction/ncep2/sphum_ncep2_oafluxgrid_19890101_20091231.nc'161 fi=tropflux_id_env+'sphum_ncep2_oafluxgrid_19890101_20091231.nc' 104 162 initncdf, fi 105 163 nce=read_ncdf("shum", date1-1, date2, file=fi,/nostr) … … 320 378 closeps 321 379 322 fig='q2m_validation_scatter_2000_2009_v50.ps' 323 spawn, 'mv '+psdir+'idl.ps '+cpsdir+fig 324 spawn, 'gv '+cpsdir+fig 325 return 380 fig=tropflux_od_env+'q2m_validation_scatter_2000_2009_v50.ps' 381 spawn, 'mv '+psdir+'idl.ps '+fig 382 spawn, 'gv '+fig 326 383 end 384 327 385 function x_site_location, site 328 386 n1=strpos(site, 's') -
trunk/src/paper01/fig3/read_variables_v2.pro
r43 r44 1 1 ;+ 2 ; .. _read_variables_v2.pro 2 ; 3 ; .. _read_variables_v2.pro: 3 4 ; 4 5 ; ===================== … … 9 10 ; =========== 10 11 ; 12 ; .. graphviz:: 13 ; 14 ; digraph read_variables_v2 { 15 ; graph [ 16 ; rankdir="LR", 17 ; ] 18 ; rad [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/rad*_dy.cdf"]; 19 ; qlat[shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/qlat*_dy.cdf"]; 20 ; pos[shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/pos*_dy.cdf"]; 21 ; met[shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/met*_dy.cdf"]; 22 ; 23 ; at [shape=tab,fontname=Courier,label="at"] 24 ; sw [shape=tab,fontname=Courier,label="sw"] 25 ; rh [shape=tab,fontname=Courier,label="rh"] 26 ; sst [shape=tab,fontname=Courier,label="sst"] 27 ; wu [shape=tab,fontname=Courier,label="wu"] 28 ; wv [shape=tab,fontname=Courier,label="wv"] 29 ; ws [shape=tab,fontname=Courier,label="ws"] 30 ; lh [shape=tab,fontname=Courier,label="lh"] 31 ; 32 ; read_variables_v2 [shape=box, 33 ; fontname=Courier, 34 ; color=blue, 35 ; URL="http://forge.ipsl.jussieu.fr/tropflux/browser/trunk/src/paper01/fig2/read_variables_v2.pro", 36 ; label="${TROPFLUX}/src/paper01/fig2/read_variables_v2.pro"]; 37 ; 38 ; {qlat rad pos met} -> {read_variables_v2} -> {at sw rh sst wu wv ws lh} 39 ; } 40 ; 11 41 ; SEE ALSO 12 42 ; ======== 13 43 ; 14 44 ; :ref:`tropflux_profile.sh` 45 ; 46 ; :ref:`air_validation_scatter_2000_2009_v50.pro` 47 ; :ref:`q2m_validation_scatter_2000_2009_v50.pro` 48 ; :ref:`sst_validation_scatter_2000_2009_v50.pro` 49 ; :ref:`ws_validation_scatter_2000_2009_v50.pro` 15 50 ; 16 51 ; EXAMPLES … … 22 57 ; IDL> csite=++ 23 58 ; IDL> ++ 24 ; IDL> read_variables_v2, csite,date1,date2,nsmooth, at, sw,rh,sst,wu,wv,ws, lh 59 ; IDL> read_variables_v2, csite, date1, date2, nsmooth, at, sw, rh, sst, wu, wv, ws, lh 60 ; 61 ; TODO 62 ; ==== 63 ; 64 ; externalize function 65 ; 66 ; make it work (missing data) 25 67 ; 26 68 ; EVOLUTIONS … … 29 71 ; $Id$ 30 72 ; 73 ; - fplod 20110412T110534Z aedon.locean-ipsl.upmc.fr (Darwin) 74 ; 75 ; * remove hard coding path 76 ; * add graphviz 77 ; 31 78 ; - fplod 20110411T142955Z aedon.locean-ipsl.upmc.fr (Darwin) 32 79 ; … … 34 81 ; 35 82 ;- 36 pro read_variables_v2, csite,date1,date2,nsmooth, $ 37 at, sw,rh,sst,wu,wv,ws, lh 38 39 ; 40 @common 41 42 dir='/Volumes/Iomega_HDD/work/flux_reconstruction/mooriing_all/' 83 pro read_variables_v2, csite, date1, date2, nsmooth, $ 84 at, sw, rh, sst, wu, wv, ws, lh 85 @cm_tropflux 43 86 ; 44 87 ;; DEFINE THE OUTPUT TIME AXIS … … 54 97 ;; - if not, fill variable with missing values 55 98 ; 56 57 58 99 ;;SHORTWAVE 59 fi= dir+'rad'+csite+'_dy.cdf'100 fi=tropflux_id_env+'rad'+csite+'_dy.cdf' 60 101 f=file_test(fi) 61 102 sw=fltarr(jpt)+!values.f_nan … … 69 110 ind=where((tt0 ge -1e-5) and (tt0 le jpt-1+1e-5)) 70 111 if (ind(0) ne -1) then begin 71 sw(tt0(ind))=sw0(ind) 72 sw_q(tt0(ind))=sw_q0(ind) 73 endif 74 endif 112 sw(tt0(ind))=sw0(ind) 113 sw_q(tt0(ind))=sw_q0(ind) 114 endif 115 endif else begin 116 print, fi, ' not found' 117 endelse 75 118 76 119 ind_sw=where(sw_q ne 1. and sw_q ne 2.) 77 120 78 121 ;;LHF 79 fi= dir+'qlat'+csite+'_dy.cdf'122 fi=tropflux_id_env+'qlat'+csite+'_dy.cdf' 80 123 f=file_test(fi) 81 124 lh=fltarr(jpt)+!values.f_nan … … 88 131 ind=where((tt0 ge -1e-5) and (tt0 le jpt-1+1e-5)) 89 132 if (ind(0) ne -1) then begin 90 lh(tt0(ind))=lh0(ind) 91 lh_q(tt0(ind))=lh_q0(ind) 92 endif 93 endif 133 lh(tt0(ind))=lh0(ind) 134 lh_q(tt0(ind))=lh_q0(ind) 135 endif 136 endif else begin 137 print, fi, ' not found' 138 endelse 94 139 95 140 ;;POSITION 96 fi= dir+'pos'+csite+'_dy.cdf'141 fi=tropflux_id_env+'pos'+csite+'_dy.cdf' 97 142 f=file_test(fi) 98 143 lat=fltarr(jpt)+!values.f_nan … … 108 153 lon(tt0(ind))=lon0(ind) 109 154 endif 110 endif 155 endif else begin 156 print, fi, ' not found' 157 endelse 111 158 112 159 113 160 ;;MET PARAMETERS 114 fi= dir+'met'+csite+'_dy.cdf'161 fi=tropflux_id_env+'met'+csite+'_dy.cdf' 115 162 f=file_test(fi) 116 163 wu=fltarr(jpt)+!values.f_nan … … 151 198 sst_q(tt0(ind))=sst_q0(ind) 152 199 endif 153 endif 200 endif else begin 201 print, fi, ' not found' 202 endelse 154 203 155 204 ind_at=where(at_q ne 1 and at_q ne 2) & ind_ws=where(ws_q ne 1 and ws_q ne 2) -
trunk/src/paper01/fig3/sst_validation_scatter_2000_2009_v50.pro
r43 r44 1 1 ;+ 2 ; .. _sst_validation_scatter_2000_2009_v50.pro 2 ; .. _sst_validation_scatter_2000_2009_v50.pro: 3 3 ; 4 4 ; ======================================== … … 362 362 closeps 363 363 364 fig= 'sst_validation_scatter_2000_2009_v50.ps'365 spawn, 'mv '+psdir+'idl.ps '+ cpsdir+fig366 spawn, 'gv '+ cpsdir+fig364 fig=tropflux_od_env+'sst_validation_scatter_2000_2009_v50.ps' 365 spawn, 'mv '+psdir+'idl.ps '+fig 366 spawn, 'gv '+fig 367 367 return 368 368 end -
trunk/src/paper01/fig3/statistics_3var_v1.pro
r43 r44 1 1 ;+ 2 ; .. _statistics_3var_v1.pro 2 ; .. _statistics_3var_v1.pro: 3 3 ; 4 4 ; ====================== -
trunk/src/paper01/fig3/ws_validation_scatter_2000_2009_v50.pro
r43 r44 1 1 ;+ 2 ; .. _ws_validation_scatter_2000_2009_v50.pro 2 ; .. _ws_validation_scatter_2000_2009_v50.pro: 3 3 ; 4 4 ; ======================================= … … 117 117 help, nce 118 118 119 dir='/Volumes/Iomega_HDD/work/flux_reconstruction/gridded_data/' 120 fi=dir+'zonal_wind_speed_oafluxgrid_30N30S.nc' 119 fi=tropflux_id_env+'zonal_wind_speed_oafluxgrid_30N30S.nc' 121 120 initncdf, fi 122 121 u=read_ncdf("u", date1, date2, file=fi,/nostr) 123 fi= dir+'meridional_wind_speed_oafluxgrid_30N30S.nc'122 fi=tropflux_id_env+'meridional_wind_speed_oafluxgrid_30N30S.nc' 124 123 initncdf, fi 125 124 v=read_ncdf("v", date1, date2, file=fi,/nostr) … … 368 367 closeps 369 368 370 fig= 'ws_validation_scatter_2000_2009_v50.ps'371 spawn, 'mv '+psdir+'idl.ps '+ cpsdir+fig372 spawn, 'gv '+ cpsdir+fig369 fig=tropflux_od_env+'ws_validation_scatter_2000_2009_v50.ps' 370 spawn, 'mv '+psdir+'idl.ps '+fig 371 spawn, 'gv '+fig 373 372 return 374 373 end -
trunk/src/tropflux_init.pro
r40 r44 77 77 @cm_4mesh 78 78 @cm_4data 79 @cm_4ps 80 @cm_general 79 81 ; 80 82 ; define default directories
Note: See TracChangeset
for help on using the changeset viewer.