Changeset 94 for trunk/src/paper01/fig16/fig16_timeseries_nino3.pro
- Timestamp:
- 08/25/11 13:02:20 (13 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/src/paper01/fig16/fig16_timeseries_nino3.pro
r85 r94 102 102 103 103 openps, FILENAME = project_od_env+'fig16_timeseries_nino3.ps' 104 date1=19940101 & date2=20071231 104 date1=19940101 105 date2=20071231 105 106 box=[240,330,-5,5] 106 107 … … 123 124 sst_nclim=total(reform(sst(0:nt-1),365,ny),2)/ny 124 125 125 trop_nclim=[trop_nclim ,trop_nclim, trop_nclim] & trop_nclim=smooth(trop_nclim,nsmooth,/nan) & trop_nclim=trop_nclim[365:365+364] 126 oafl_nclim=[oafl_nclim ,oafl_nclim, oafl_nclim] & oafl_nclim=smooth(oafl_nclim,nsmooth,/nan) & oafl_nclim=oafl_nclim[365:365+364] 127 nce2_nclim=[nce2_nclim ,nce2_nclim ,nce2_nclim] & nce2_nclim=smooth(nce2_nclim,nsmooth,/nan) & nce2_nclim=nce2_nclim[365:365+364] 128 era_nclim=[era_nclim ,era_nclim, era_nclim] & era_nclim=smooth(era_nclim,nsmooth,/nan) & era_nclim=era_nclim[365:365+364] 129 nce_nclim=[nce_nclim ,nce_nclim ,nce_nclim] & nce_nclim=smooth(nce_nclim,nsmooth,/nan) & nce_nclim=nce_nclim[365:365+364] 130 sst_nclim=[sst_nclim ,sst_nclim ,sst_nclim] & sst_nclim=smooth(sst_nclim,nsmooth,/nan) & sst_nclim=sst_nclim[365:365+364] 131 132 time=time & jpt=n_elements(time) 126 trop_nclim=[trop_nclim ,trop_nclim, trop_nclim] 127 trop_nclim=smooth(trop_nclim,nsmooth,/nan) 128 trop_nclim=trop_nclim[365:365+364] 129 oafl_nclim=[oafl_nclim ,oafl_nclim, oafl_nclim] 130 oafl_nclim=smooth(oafl_nclim,nsmooth,/nan) 131 oafl_nclim=oafl_nclim[365:365+364] 132 nce2_nclim=[nce2_nclim ,nce2_nclim ,nce2_nclim] 133 nce2_nclim=smooth(nce2_nclim,nsmooth,/nan) 134 nce2_nclim=nce2_nclim[365:365+364] 135 era_nclim=[era_nclim ,era_nclim, era_nclim] 136 era_nclim=smooth(era_nclim,nsmooth,/nan) 137 era_nclim=era_nclim[365:365+364] 138 nce_nclim=[nce_nclim ,nce_nclim ,nce_nclim] 139 nce_nclim=smooth(nce_nclim,nsmooth,/nan) 140 nce_nclim=nce_nclim[365:365+364] 141 sst_nclim=[sst_nclim ,sst_nclim ,sst_nclim] 142 sst_nclim=smooth(sst_nclim,nsmooth,/nan) 143 sst_nclim=sst_nclim[365:365+364] 144 145 time=time 146 jpt=n_elements(time) 133 147 134 148 caldat, time,mon,day,yea 135 trop_clim=trop*0. & era_clim=trop*0. 136 oafl_clim=trop*0. & nce_clim=trop*0. 137 nce2_clim=trop*0. & sst_clim=trop*0. 149 trop_clim=trop*0. 150 era_clim=trop*0. 151 oafl_clim=trop*0. 152 nce_clim=trop*0. 153 nce2_clim=trop*0. 154 sst_clim=trop*0. 138 155 139 156 for jt=0,jpt-1 do begin … … 147 164 endfor 148 165 149 trop_ano=trop-trop_clim & oafl_ano=oafl-oafl_clim 150 nce2_ano=nce2-nce2_clim & era_ano=era-era_clim 151 nce_ano=nce-nce_clim & sst_ano=sst-sst_clim 152 153 trop_hf=trop*0. & nce2_hf=trop*0. 154 oafl_hf=trop*0. & era_hf=trop*0. 166 trop_ano=trop-trop_clim 167 oafl_ano=oafl-oafl_clim 168 nce2_ano=nce2-nce2_clim 169 era_ano=era-era_clim 170 nce_ano=nce-nce_clim 171 sst_ano=sst-sst_clim 172 173 trop_hf=trop*0. 174 nce2_hf=trop*0. 175 oafl_hf=trop*0. 176 era_hf=trop*0. 155 177 nce_hf=trop*0. 156 178 157 Thf=1 & Tlf=90 179 Thf=1 180 Tlf=90 158 181 tsamp=1. ; monthly 159 182 freq=findgen(jpt/2+1)/(jpt*tsamp) … … 189 212 help, trop_hf, nce2_hf, era_hf, nce_hf 190 213 191 trop_ano=trop_ano-trop_hf & oafl_ano=oafl_ano-oafl_hf 192 nce2_ano=nce2_ano-nce2_hf & era_ano=era_ano-era_hf 214 trop_ano=trop_ano-trop_hf 215 oafl_ano=oafl_ano-oafl_hf 216 nce2_ano=nce2_ano-nce2_hf 217 era_ano=era_ano-era_hf 193 218 nce_ano=nce_ano-nce_hf 194 219
Note: See TracChangeset
for help on using the changeset viewer.