PRO GRS_RESP, EDGES_OUT=EOUT, EDGES_IN=EIN, AREA=AREA, DRM=GRSMATRIX, ERROR=ERROR

error=1

file = loc_file('grsmatrix.dat',path=[curdir(),'SSWDB_SMM'], count=count)
if count ge 1 then begin
    error=0 
    grsmatrix=0
    readme=0
    dmatrix=0
    ematrix=0
    wmatrix=0
    em=0
    restore,file
    edge_products, [ematrix(0)-dmatrix(0)/2, ematrix+dmatrix/2]*1e3, edges_2=ein
    edge_products, [em(0)-wmatrix(0)/2, em+wmatrix/2]*1e3, edges_2=eout
    area = 317.55 ;cm2
    grsmatrix= temporary( f_div(grsmatrix, rebin(wmatrix,n_elements(em),n_elements(ematrix)))/area)
    
    
endif

end