Changeset 94 for trunk/src/paper01/fig6/icoads_ws_stats_paper.pro
- Timestamp:
- 08/25/11 13:02:20 (13 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/src/paper01/fig6/icoads_ws_stats_paper.pro
r85 r94 86 86 ; partie a changer 87 87 marge=[-2,-2, -4,2] 88 st=19890101 & en=20091231 88 st=19890101 89 en=20091231 89 90 domdef, box 90 91 … … 124 125 help, noc 125 126 126 si=size(trop) & nt=si(3) 127 si=size(trop) 128 nt=si(3) 127 129 128 130 lat_rama=[67, 55, 55, 80.5, 80.5, 80.5, 80.5, 90, 90, 90, 90, 90, 90, 90, 95] … … 142 144 lon=[lon_rama, lon_pirata, lon_tao] 143 145 144 lat_tao=lat & lon_tao=lon 146 lat_tao=lat 147 lon_tao=lon 145 148 nn=n_elements(lat) 146 149 … … 149 152 x=lat_tao(n) 150 153 y=lon_tao(n) 151 dx=abs(reform(glamt-x)) & dy=abs(reform(gphit)-y) 154 dx=abs(reform(glamt-x)) 155 dy=abs(reform(gphit)-y) 152 156 ind=where((dx le 0.5) and (dy le 0.5)) 153 157 if (total(ind) ge 0.) then begin … … 166 170 167 171 param=icoads*mask 168 trop=trop*mask & oaf=oaf*mask & erai=erai*mask 169 ncep2=ncep2*mask & ncep1=ncep1*mask 172 trop=trop*mask 173 oaf=oaf*mask 174 erai=erai*mask 175 ncep2=ncep2*mask 176 ncep1=ncep1*mask 170 177 noc=noc*mask 171 178 172 179 ind=where(finite(param) and finite(noc) and finite(erai) and finite(trop) and finite(oaf) and finite(ncep2) and finite(ncep1)) 173 x=param(ind) & y=trop(ind) 174 param=param(ind) & trop=trop(ind) & erai=erai(ind) & oaf=oaf(ind) & ncep2=ncep2(ind) & ncep1=ncep1(ind) 180 x=param(ind) 181 y=trop(ind) 182 param=param(ind) 183 trop=trop(ind) 184 erai=erai(ind) 185 oaf=oaf(ind) 186 ncep2=ncep2(ind) 187 ncep1=ncep1(ind) 175 188 noc=noc(ind) 176 189 … … 219 232 abs_noc=string(abs_diff, format='(f5.2)') 220 233 ; 221 mio=0 & mao=15 & inx=1 & nx=(mao-mio)/inx+1l & xx=mio+indgen(nx)*inx 222 mis=0 & mas=15 & iny=1 & ny=(mas-mis)/iny+1l & yy=mis+indgen(ny)*iny 223 xp=4 & yp=13 224 x1=0 & x2=15 & y1=x1 & y2=x2 234 mio=0 235 mao=15 236 inx=1 237 nx=(mao-mio)/inx+1l 238 xx=mio+indgen(nx)*inx 239 mis=0 240 mas=15 241 iny=1 242 ny=(mas-mis)/iny+1l 243 yy=mis+indgen(ny)*iny 244 xp=4 245 yp=13 246 x1=0 247 x2=15 248 y1=x1 249 y2=x2 225 250 ; 226 251 pdf, param,trop,mio,mao,inx,mis,mas,iny, $
Note: See TracChangeset
for help on using the changeset viewer.