fid = ncdf_open('spectral_irradiance_Lean_1610-2009_ann_c170912.nc',/write) ncdf_varget,fid,'ssi',ssi ncdf_varget,fid,'ssi_ref',ssi_ref ncdf_varget,fid,'tsi_ref',tsi_ref ncdf_varget,fid,'tsi',tsi restore, 'scaled_ssi_tsi_extrmin_HD.sav' FOR iTime = 0, 399 DO ssi[*,iTime] = ssi_exmin[*] ssi_ref = ssi_exmin tsi[*] = tsi_exmin tsi_ref = tsi_exmin NCDF_VARPUT, fid,'ssi',ssi NCDF_VARPUT, fid,'ssi_ref',ssi_ref NCDF_VARPUT, fid,'tsi_ref',tsi_ref NCDF_VARPUT, fid,'tsi',tsi