pro get_crf_comparison_histos_6hr, fname, forcing_bins


common total_forcing_comparison_arrays, toa_net_forcing_comparison_total, toa_solar_forcing_comparison_total, $
toa_ir_forcing_comparison_total, sfc_net_forcing_comparison_total, $
sfc_solar_forcing_comparison_total,sfc_ir_forcing_comparison_total, atm_net_forcing_comparison_total, $
atm_solar_forcing_comparison_total,atm_ir_forcing_comparison_total

common total_forcing_comparison_arrays2, toa_net_forcing_comparison_total_ncep, toa_solar_forcing_comparison_total_ncep, $
toa_ir_forcing_comparison_total_ncep, sfc_net_forcing_comparison_total_ncep, $
sfc_solar_forcing_comparison_total_ncep,sfc_ir_forcing_comparison_total_ncep, atm_net_forcing_comparison_total_ncep, $
atm_solar_forcing_comparison_total_ncep,atm_ir_forcing_comparison_total_ncep

; check to see if the file exists

OPENR, 1, fname, ERROR = err

if err eq 0 then begin
close, 1

; file exists, read the data
print, 'attempting to open ',fname, ' err is ',err
cdfid=ncdf_open(fname,/nowrite)
print, 'reading ',fname

toa_net_forcing_comparison_total_id=ncdf_varid(cdfid, 'toa_net_forcing_comparison_total')
toa_solar_forcing_comparison_total_id=ncdf_varid(cdfid, 'toa_solar_forcing_comparison_total')
toa_ir_forcing_comparison_total_id=ncdf_varid(cdfid, 'toa_ir_forcing_comparison_total')
sfc_net_forcing_comparison_total_id=ncdf_varid(cdfid, 'sfc_net_forcing_comparison_total')
sfc_solar_forcing_comparison_total_id=ncdf_varid(cdfid, 'sfc_solar_forcing_comparison_total')
sfc_ir_forcing_comparison_total_id=ncdf_varid(cdfid, 'sfc_ir_forcing_comparison_total')
atm_net_forcing_comparison_total_id=ncdf_varid(cdfid, 'atm_net_forcing_comparison_total')
atm_solar_forcing_comparison_total_id=ncdf_varid(cdfid, 'atm_solar_forcing_comparison_total')
atm_ir_forcing_comparison_total_id=ncdf_varid(cdfid, 'atm_ir_forcing_comparison_total')
toa_net_forcing_comparison_total_ncep_id=ncdf_varid(cdfid, 'toa_net_forcing_comparison_total_ncep')
toa_solar_forcing_comparison_total_ncep_id=ncdf_varid(cdfid, 'toa_solar_forcing_comparison_total_ncep')
toa_ir_forcing_comparison_total_ncep_id=ncdf_varid(cdfid, 'toa_ir_forcing_comparison_total_ncep')
sfc_net_forcing_comparison_total_ncep_id=ncdf_varid(cdfid, 'sfc_net_forcing_comparison_total_ncep')
sfc_solar_forcing_comparison_total_ncep_id=ncdf_varid(cdfid, 'sfc_solar_forcing_comparison_total_ncep')
sfc_ir_forcing_comparison_total_ncep_id=ncdf_varid(cdfid, 'sfc_ir_forcing_comparison_total_ncep')
atm_net_forcing_comparison_total_ncep_id=ncdf_varid(cdfid, 'atm_net_forcing_comparison_total_ncep')
atm_solar_forcing_comparison_total_ncep_id=ncdf_varid(cdfid, 'atm_solar_forcing_comparison_total_ncep')
atm_ir_forcing_comparison_total_ncep_id=ncdf_varid(cdfid, 'atm_ir_forcing_comparison_total_ncep')

ncdf_varget, cdfid, toa_net_forcing_comparison_total_id, toa_net_forcing_comparison_total
ncdf_varget, cdfid, toa_solar_forcing_comparison_total_id, toa_solar_forcing_comparison_total
ncdf_varget, cdfid, toa_ir_forcing_comparison_total_id, toa_ir_forcing_comparison_total
ncdf_varget, cdfid, sfc_net_forcing_comparison_total_id, sfc_net_forcing_comparison_total
ncdf_varget, cdfid, sfc_solar_forcing_comparison_total_id, sfc_solar_forcing_comparison_total
ncdf_varget, cdfid, sfc_ir_forcing_comparison_total_id, sfc_ir_forcing_comparison_total
ncdf_varget, cdfid, atm_net_forcing_comparison_total_id, atm_net_forcing_comparison_total
ncdf_varget, cdfid, atm_solar_forcing_comparison_total_id, atm_solar_forcing_comparison_total
ncdf_varget, cdfid, atm_ir_forcing_comparison_total_id, atm_ir_forcing_comparison_total
ncdf_varget, cdfid, toa_net_forcing_comparison_total_ncep_id, toa_net_forcing_comparison_total_ncep
ncdf_varget, cdfid, toa_solar_forcing_comparison_total_ncep_id, toa_solar_forcing_comparison_total_ncep
ncdf_varget, cdfid, toa_ir_forcing_comparison_total_ncep_id, toa_ir_forcing_comparison_total_ncep
ncdf_varget, cdfid, sfc_net_forcing_comparison_total_ncep_id, sfc_net_forcing_comparison_total_ncep
ncdf_varget, cdfid, sfc_solar_forcing_comparison_total_ncep_id, sfc_solar_forcing_comparison_total_ncep
ncdf_varget, cdfid, sfc_ir_forcing_comparison_total_ncep_id, sfc_ir_forcing_comparison_total_ncep
ncdf_varget, cdfid, atm_net_forcing_comparison_total_ncep_id, atm_net_forcing_comparison_total_ncep
ncdf_varget, cdfid, atm_solar_forcing_comparison_total_ncep_id, atm_solar_forcing_comparison_total_ncep
ncdf_varget, cdfid, atm_ir_forcing_comparison_total_ncep_id, atm_ir_forcing_comparison_total_ncep

NCDF_CLOSE, cdfid


endif else begin	;OPENR, 1, fname, ERROR = err

; file does not exist, create it

write_CRF_comparison_histos_6hr, fname, log_iwp_bins, forcing_bins

endelse		;OPENR, 1, fname, ERROR = err



return
end