Changeset 97 for trunk/src/paper01/fig1
- Timestamp:
- 08/30/11 14:38:29 (13 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/src/paper01/fig1/fig1_no_obs_map.pro
r94 r97 13 13 ; has been produced by :ref:`oaflux_mask_30N30S.pro`. 14 14 ; 15 ; Figure of 15 ; Figure of 16 16 ; ++ 17 17 ; is saved in … … 107 107 reinitplt, /z,/invert 108 108 key_portrait = 1 109 109 ; 110 110 openps, FILENAME = project_od_env + 'fig1_no_obs_map.ps' 111 111 close,/all 112 112 ; 113 113 yr=365. 114 lh_mi=365/yr 115 lh_ma=3285/yr 114 lh_mi=365/yr 115 lh_ma=3285/yr 116 116 lh_int=0.25 117 sw_mi=365/yr 118 sw_ma=3285/yr 117 sw_mi=365/yr 118 sw_ma=3285/yr 119 119 sw_int=0.25 120 lw_mi=365/yr 121 lw_ma=1460/yr 120 lw_mi=365/yr 121 lw_ma=1460/yr 122 122 lw_int=0.125 123 123 lct=65 124 124 ; 125 125 file=project_id_env + 'mask_oaflux_30N30S.nc' 126 126 initncdf, file 127 127 domdef, 30,390,-30,30 128 128 msk=ncdf_lec(file,var='msk') 129 129 ; 130 130 fi=project_id_env + 'lhf_no_obs.txt' 131 131 res=read_ascii(fi,data_start=1) … … 135 135 no=reform(ff(2,*)) 136 136 no=no/yr 137 ;ind=where(no ge lh_ma) 137 ;ind=where(no ge lh_ma) 138 138 ; no(ind)=no-0.5 139 139 ; 140 140 marge1=[0,0,-5,0] 141 141 plt, .4+msk*1.e-5,realcont=2,/nocont,/nofill,xminor=1,yminor=1,box=box,lct=lct, lh_mi, lh_ma, int=lh_int, $ 142 142 title='a) Turbulent Flux validation sites', subtitle='', small=[1,4,1],/rempl, marge=marge1, format='(i5)' 143 143 NN=n_elements(lat) 144 144 ; 145 145 usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5],/fill 146 147 for n=0,NN-1 do begin 148 x=lon[n] 149 y=lat[n] 150 c=no[n] 151 cmi=lh_mi 152 cma=lh_ma 146 ; 147 for n=0,NN-1 do begin 148 x=lon[n] 149 y=lat[n] 150 c=no[n] 151 cmi=lh_mi 152 cma=lh_ma 153 153 dc=cma-cmi 154 154 col=((10+244*(c-cmi)/dc) > 10) < 254 155 155 plots, x,y,psym=8,symsize=1.5,color=col 156 156 endfor 157 157 ; 158 158 usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5] 159 160 for n=0,NN-1 do begin 161 x=lon[n] 162 y=lat[n] 163 c=no[n] 164 159 ; 160 for n=0,NN-1 do begin 161 x=lon[n] 162 y=lat[n] 163 c=no[n] 164 ; 165 165 plots, x,y,psym=8,symsize=1.5,color=0 166 166 endfor 167 167 ; 168 168 fi=project_id_env + 'sw_no_obs.txt' 169 169 res=read_ascii(fi,data_start=1) … … 173 173 no=reform(ff(2,*)) 174 174 no=no/yr 175 ;ind=where(no ge sw_ma) 175 ;ind=where(no ge sw_ma) 176 176 ;no(ind)=sw_ma-0.5 177 177 marge1=[0,0,-5,0] … … 179 179 title='b) SWR validation sites', subtitle='', small=[1,4,2],/noer,/rempl, marge=marge1, format='(i5)' 180 180 NN=n_elements(lat) 181 181 ; 182 182 usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5],/fill 183 184 for n=0,NN-1 do begin 185 x=lon[n] 186 y=lat[n] 187 c=no[n] 188 cmi=sw_mi 189 cma=sw_ma 183 ; 184 for n=0,NN-1 do begin 185 x=lon[n] 186 y=lat[n] 187 c=no[n] 188 cmi=sw_mi 189 cma=sw_ma 190 190 dc=cma-cmi 191 191 col=((10+244*(c-cmi)/dc) > 10) < 254 192 192 plots, x,y,psym=8,symsize=1.5,color=col 193 193 endfor 194 194 ; 195 195 usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5] 196 197 for n=0,NN-1 do begin 198 x=lon[n] 199 y=lat[n] 200 c=no[n] 201 196 ; 197 for n=0,NN-1 do begin 198 x=lon[n] 199 y=lat[n] 200 c=no[n] 201 ; 202 202 plots, x,y,psym=8,symsize=1.5,color=0 203 203 endfor 204 204 ; 205 205 fi=project_id_env + 'lw_no_obs.txt' 206 206 res=read_ascii(fi,data_start=1) … … 210 210 no=reform(ff(2,*)) 211 211 no=no/yr 212 ;ind=where(no ge lw_ma) 212 ;ind=where(no ge lw_ma) 213 213 ;no(ind)=lw_ma-0.25 214 214 ; 215 215 marge1=[0,0,-5,0] 216 216 plt, .4+msk*1.e-5,realcont=2,/nocont,/nofill,/noer,xminor=1,yminor=1,box=box,lct=lct,lw_mi, lw_ma, int=lw_int, $ 217 217 title='c) LWR validation sites', subtitle='', small=[1,4,3],/rempl, marge=marge1, format='(i5)' 218 218 NN=n_elements(lat) 219 219 ; 220 220 usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5],/fill 221 222 for n=0,NN-1 do begin 223 x=lon[n] 224 y=lat[n] 225 c=no[n] 226 cmi=lw_mi 227 cma=lw_ma 221 ; 222 for n=0,NN-1 do begin 223 x=lon[n] 224 y=lat[n] 225 c=no[n] 226 cmi=lw_mi 227 cma=lw_ma 228 228 dc=cma-cmi 229 229 col=((10+244*(c-cmi)/dc) > 10) < 254 230 230 plots, x,y,psym=8,symsize=1.5,color=col 231 231 endfor 232 232 ; 233 233 usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5] 234 235 for n=0,NN-1 do begin 236 x=lon[n] 237 y=lat[n] 238 c=no[n] 239 234 ; 235 for n=0,NN-1 do begin 236 x=lon[n] 237 y=lat[n] 238 c=no[n] 239 ; 240 240 plots, x,y,psym=8,symsize=1.5,color=0 241 241 endfor 242 243 ; ;independent validation sites244 245 plots, 61.5, 15.5,psym=2 ; ;as location246 plots, 331, 25.5,psym=2 ; ;subduction (central)247 plots, 326, 18,psym=2 ; ;subduction (sw)248 plots, 338, 18,psym=2 ; ;subduction (se)249 plots, 156, -1.7 , psym=2 ; ;coare250 plots, 274.2, -21, psym=2 ; ;stratus251 plots, 202.46, 22.45, psym=2 ; ;WHOTS252 242 ; 243 ; independent validation sites 244 ; 245 plots, 61.5, 15.5,psym=2 ; as location 246 plots, 331, 25.5,psym=2 ; subduction (central) 247 plots, 326, 18,psym=2 ; subduction (sw) 248 plots, 338, 18,psym=2 ; subduction (se) 249 plots, 156, -1.7 , psym=2 ; coare 250 plots, 274.2, -21, psym=2 ; stratus 251 plots, 202.46, 22.45, psym=2 ; WHOTS 252 ; 253 253 closeps 254 254 ; 255 255 end
Note: See TracChangeset
for help on using the changeset viewer.