File name : GRS_spectra_V001.idl File description : A file containing integrated gamma ray spectra for the entire Moon. The individual gamma ray spectra have been integrated over the time and binned into 5 x 5 degree cells. Each integrated spectra has 512 channels. Values are in units of counts per 32-seconds. To restore these data, under IDL, run the command: restore,'GRS_spectra_V001.idl' You will get 1 structure : GRSPDS_SPECTRA STRUCT = -> Array[1] That has the following fields : INDEX INT Array[1790] LAT_MIN FLOAT Array[1790] LAT_MAX FLOAT Array[1790] LON_MIN FLOAT Array[1790] LON_MAX FLOAT Array[1790] ACCEPTED FLOAT Array[1790, 512] ACC_ERR FLOAT Array[1790, 512] REJECTED FLOAT Array[1790, 512] REJ_ERR FLOAT Array[1790, 512] NSPECTRA FLOAT Array[1790] ENERGY FLOAT Array[512] The meaning of these fieds is : INDEX Equal-area pixel number LAT_MIN Minimum latitude of the pixel LAT_MAX Maximum latitude of the pixel LON_MIN Minimum longitude of the pixel LON_MAX Maximum longitude of the pixel ACCEPTED Accepted spectrum (512 channels) mean ACC_ERR Accepted spectrum (512 channels) sdt deviation REJECTED Rejected spectrum (512 channels) mean REJ_ERR Rejected spectrum (512 channels) sdt deviation NSPECTRA Number of 32s-spectra in pixel ENERGY Channel to Energy correspondance CREATED Date when the file was created AUTHORS Authors names Example : How to extract a rectangular array of the accepted 2.2MeV line ? restore,'GRS_spectra_V001.idl' % We see that : grspds_spectra.energy[123] = 2.20271 data=grspds_spectra.accepted[*,123] pixels={lat_min : grspds_spectra.lat_min , $ lat_max : grspds_spectra.lat_max , $ lon_min : grspds_spectra.lon_min , $ lon_max : grspds_spectra.lon_max} make_grid,data,pixels,array % array is a 720 by 360 array that correspond to % 720 longitudes from -180 to +180 % 360 latitudes from -90 t +90 % You may display that array : tvscl,array