TRMM 3B42 자료 읽기

fn_in = 'hdfname'
print, fn_in

sd_id=HDF_SD_START(fn_in,/read)

sds_id=HDF_SD_SELECT(sd_id,0) ; rain
HDF_SD_GETINFO,sds_id,name=n,ndims=r,type=t,natts=nats,dims=dims
rain = fltarr(dims(0), dims(1))
HDF_SD_GETDATA, sds_id, rain

sds_id=HDF_SD_SELECT(sd_id,1) ; err
HDF_SD_GETINFO,sds_id,name=n,ndims=r,type=t,natts=nats,dims=dims
err = fltarr(dims(0), dims(1))
HDF_SD_GETDATA, sds_id, err
rain1 = transpose(rain)
err1 = transpose(err)
HDF_SD_END, sd_id