Changeset 94 for trunk/src/paper01/fig11
- Timestamp:
- 08/25/11 13:02:20 (13 years ago)
- Location:
- trunk/src/paper01/fig11
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/src/paper01/fig11/swr_isccp_tropflux_correlation_fig10.pro
r85 r94 105 105 106 106 openps, FILENAME = project_od_env+'correlation_rmsd_olra_swra_smooth.ps' 107 st=19890101 & en=20071231 107 st=19890101 108 en=20071231 108 109 109 110 restore, project_od_env + "correlation_isccp_olr.idl" -
trunk/src/paper01/fig11/swr_isccp_tropflux_new_v1.pro
r85 r94 104 104 105 105 openps, FILENAME = project_od_env+'correlation_rmsd_olra_swra_smooth.ps' 106 st=19890101 & en=20071231 106 st=19890101 107 en=20071231 107 108 108 109 file=project_od_env+'swr_oafluxgrid_1985_2007.nc' … … 132 133 nsmooth=45 133 134 134 swi=[[[swi]],[[swi]],[[swi]]] & olri=[[[olri]],[[olri]],[[olri]]] 135 si=size(sw) & n1=si(1) & n2=si(2) 135 swi=[[[swi]],[[swi]],[[swi]]] 136 olri=[[[olri]],[[olri]],[[olri]]] 137 si=size(sw) 138 n1=si(1) 139 n2=si(2) 136 140 137 141 for j=0,n2-1 do begin … … 144 148 endfor 145 149 146 swi=swi[*,*,365:365+364] & olri=olri[*,*,365:365+364] 150 swi=swi[*,*,365:365+364] 151 olri=olri[*,*,365:365+364] 147 152 help, swi, olri 148 153 149 154 caldat, time,mon,day,yea 150 swm=olr*0. & olrm=olr*0. 155 swm=olr*0. 156 olrm=olr*0. 151 157 for jt=0,jpt-1 do begin 152 158 jtt=(time(jt)-julday(1,1,yea(jt))) < 364 … … 160 166 m=reform(sw_trop(*,*,0)) 161 167 tmask=finite(m) 162 cor_olr=fltarr(jpi,jpj)+!values.f_nan & rms_olr=fltarr(jpi,jpj)+!values.f_nan 168 cor_olr=fltarr(jpi,jpj)+!values.f_nan 169 rms_olr=fltarr(jpi,jpj)+!values.f_nan 163 170 164 171 for jj=0,jpj-1 do begin … … 177 184 endfor 178 185 179 cor=fltarr(jpi,jpj)+!values.f_nan & rms=fltarr(jpi,jpj)+!values.f_nan 186 cor=fltarr(jpi,jpj)+!values.f_nan 187 rms=fltarr(jpi,jpj)+!values.f_nan 180 188 for jj=1,jpj-2 do begin 181 189 for ji=1,jpi-2 do begin … … 188 196 save, cor, cor_olr, filename=project_od_env + 'correlation_isccp_olr.idl' 189 197 190 lct=64 & cs=0.9 198 lct=64 199 cs=0.9 191 200 marge1=[0,0,-5,0] 192 201 help, cor
Note: See TracChangeset
for help on using the changeset viewer.