file="skt.mon.mean.nc"
        file_id=ncdf_open(file)
        file_info=ncdf_inquire(file_id)
        nvars=file_info.nvars
        print, "No. of variables", nvars
        for var_id=0,nvars-1 do begin
          varinfo=ncdf_varinq(file_id, var_id)
          print, var_id, "  ", varinfo.name
        endfor

        var_id=3
        ncdf_attget, file_id, var_id, "missing_value", missing_value
        ncdf_attget, file_id, var_id, "scale_factor", scale_factor
        ncdf_attget, file_id, var_id, "add_offset", add_offset
        print,"missing_value", missing_value
        print, "scale_factor", scale_factor
        print, "add_offset", add_offset

        ncdf_varget,file_id, "lon", lon
        ncdf_varget,file_id, "lat", lat
        ncdf_varget,file_id, "time", time
        ncdf_varget,file_id, "skt", skt

        ncdf_close,file_id

	d=reform(skt(120,47,*))
	d=d-mean(d)

	nt=n_elements(time)
	fs_base,nt,nf,f,cosu,sinu
	fs,nt,cosu,sinu,d,a,b,psd

	dt=1.0/12.0
	t=1948.0+findgen(nt)*dt
        wave = WAVELET(d,dt,PERIOD=period,COI=coi,/PAD,SIGNIF=signif)
        nscale = N_ELEMENTS(period)
        CONTOUR,ABS(wave)^2,t,period, charsize=siz1, $
          XSTYLE=1,XTITLE="Time",YTITLE="Period (year)",$
          TITLE="Wavelet Spectrum:  Total", $
          YRANGE=[MAX(period),MIN(period)], $
          /YTYPE, $
          NLEVELS=25,/FILL

        signif = REBIN(TRANSPOSE(signif),nt,nscale)
        CONTOUR,ABS(wave)^2/signif,t,period, $
          /OVERPLOT,LEVEL=1.0,C_ANNOT='95%'
        PLOTS,t,coi,NOCLIP=0



	end