pro write_forcing_histos, fname, log_iwp_bins, forcing_bins


common total_forcing_arrays, log_iwp_toa_net_forcing_M_histo_total, log_iwp_toa_solar_forcing_M_histo_total, $
log_iwp_toa_IR_forcing_M_histo_total, log_iwp_atm_net_forcing_histo_total, log_iwp_atm_solar_forcing_histo_total, $
log_iwp_atm_IR_forcing_histo_total, log_iwp_sfc_net_forcing_obs_histo_total, log_iwp_sfc_solar_forcing_obs_histo_total, $
log_iwp_sfc_IR_forcing_obs_histo_total, log_iwp_sfc_net_forcing_histo_total, log_iwp_sfc_IR_forcing_histo_total, $
log_iwp_sfc_solar_forcing_histo_total, log_iwp_toa_net_forcing_histo_total, log_iwp_toa_solar_forcing_histo_total, $
log_iwp_toa_IR_forcing_histo_total

; 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
cdfid=ncdf_open(fname,/write)
print, 'writing to ',fname

log_iwp_toa_net_forcing_M_histo_total_id=ncdf_varid(cdfid, 'log_iwp_toa_net_forcing_M_histo_total')
log_iwp_toa_solar_forcing_M_histo_total_id=ncdf_varid(cdfid, 'log_iwp_toa_solar_forcing_M_histo_total')
log_iwp_toa_IR_forcing_M_histo_total_id=ncdf_varid(cdfid, 'log_iwp_toa_IR_forcing_M_histo_total')
log_iwp_atm_net_forcing_histo_total_id=ncdf_varid(cdfid, 'log_iwp_atm_net_forcing_histo_total')
log_iwp_atm_solar_forcing_histo_total_id=ncdf_varid(cdfid, 'log_iwp_atm_solar_forcing_histo_total')
log_iwp_atm_IR_forcing_histo_total_id=ncdf_varid(cdfid, 'log_iwp_atm_IR_forcing_histo_total')
log_iwp_sfc_net_forcing_obs_histo_total_id=ncdf_varid(cdfid, 'log_iwp_sfc_net_forcing_obs_histo_total')
log_iwp_sfc_solar_forcing_obs_histo_total_id=ncdf_varid(cdfid, 'log_iwp_sfc_solar_forcing_obs_histo_total')
log_iwp_sfc_IR_forcing_obs_histo_total_id=ncdf_varid(cdfid, 'log_iwp_sfc_IR_forcing_obs_histo_total')
log_iwp_sfc_net_forcing_histo_total_id=ncdf_varid(cdfid, 'log_iwp_sfc_net_forcing_histo_total')
log_iwp_sfc_IR_forcing_histo_total_id=ncdf_varid(cdfid, 'log_iwp_sfc_IR_forcing_histo_total')
log_iwp_sfc_solar_forcing_histo_total_id=ncdf_varid(cdfid, 'log_iwp_sfc_solar_forcing_histo_total')
log_iwp_toa_net_forcing_histo_total_id=ncdf_varid(cdfid, 'log_iwp_toa_net_forcing_histo_total')
log_iwp_toa_solar_forcing_histo_total_id=ncdf_varid(cdfid, 'log_iwp_toa_solar_forcing_histo_total')
log_iwp_toa_IR_forcing_histo_total_id=ncdf_varid(cdfid, 'log_iwp_toa_IR_forcing_histo_total')

ncdf_varput, cdfid, log_iwp_toa_net_forcing_M_histo_total_id, log_iwp_toa_net_forcing_M_histo_total
ncdf_varput, cdfid, log_iwp_toa_solar_forcing_M_histo_total_id, log_iwp_toa_solar_forcing_M_histo_total
ncdf_varput, cdfid, log_iwp_toa_IR_forcing_M_histo_total_id,log_iwp_toa_IR_forcing_M_histo_total
ncdf_varput, cdfid, log_iwp_atm_net_forcing_histo_total_id, log_iwp_atm_net_forcing_histo_total
ncdf_varput, cdfid, log_iwp_atm_solar_forcing_histo_total_id, log_iwp_atm_solar_forcing_histo_total
ncdf_varput, cdfid, log_iwp_atm_IR_forcing_histo_total_id,log_iwp_atm_IR_forcing_histo_total
ncdf_varput, cdfid, log_iwp_sfc_net_forcing_obs_histo_total_id,log_iwp_sfc_net_forcing_obs_histo_total
ncdf_varput, cdfid, log_iwp_sfc_solar_forcing_obs_histo_total_id, log_iwp_sfc_solar_forcing_obs_histo_total
ncdf_varput, cdfid, log_iwp_sfc_IR_forcing_obs_histo_total_id, log_iwp_sfc_IR_forcing_obs_histo_total
ncdf_varput, cdfid, log_iwp_sfc_net_forcing_histo_total_id, log_iwp_sfc_net_forcing_histo_total
ncdf_varput, cdfid, log_iwp_sfc_IR_forcing_histo_total_id, log_iwp_sfc_IR_forcing_histo_total
ncdf_varput, cdfid, log_iwp_sfc_solar_forcing_histo_total_id, log_iwp_sfc_solar_forcing_histo_total
ncdf_varput, cdfid, log_iwp_toa_net_forcing_histo_total_id, log_iwp_toa_net_forcing_histo_total
ncdf_varput, cdfid, log_iwp_toa_solar_forcing_histo_total_id, log_iwp_toa_solar_forcing_histo_total
ncdf_varput, cdfid, log_iwp_toa_IR_forcing_histo_total_id, log_iwp_toa_IR_forcing_histo_total

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

; file does not exist, create it
cdfid=ncdf_create(fname,/clobber)
print, 'created ',fname
forcing_bins_did=ncdf_dimdef(cdfid,'forcing_bins',n_elements(forcing_bins))
log_iwp_bins_did=ncdf_dimdef(cdfid,'log_iwp_bins',n_elements(log_iwp_bins))

forcing_bins_id=ncdf_vardef(cdfid,'forcing_bins',$
	[forcing_bins_did],/float)
log_iwp_bins_id=ncdf_vardef(cdfid,'log_iwp_bins',$
	[log_iwp_bins_did],/float)
log_iwp_toa_net_forcing_M_histo_total_id=ncdf_vardef(cdfid,'log_iwp_toa_net_forcing_M_histo_total',$
	[log_iwp_bins_did,forcing_bins_did],/long)
log_iwp_toa_solar_forcing_M_histo_total_id=ncdf_vardef(cdfid,'log_iwp_toa_solar_forcing_M_histo_total',$
	[log_iwp_bins_did,forcing_bins_did],/long)
log_iwp_toa_IR_forcing_M_histo_total_id=ncdf_vardef(cdfid,'log_iwp_toa_IR_forcing_M_histo_total',$
	[log_iwp_bins_did,forcing_bins_did],/long)
log_iwp_atm_net_forcing_histo_total_id=ncdf_vardef(cdfid,'log_iwp_atm_net_forcing_histo_total',$
	[log_iwp_bins_did,forcing_bins_did],/long)
log_iwp_atm_solar_forcing_histo_total_id=ncdf_vardef(cdfid,'log_iwp_atm_solar_forcing_histo_total',$
	[log_iwp_bins_did,forcing_bins_did],/long)
log_iwp_atm_IR_forcing_histo_total_id=ncdf_vardef(cdfid,'log_iwp_atm_IR_forcing_histo_total',$
	[log_iwp_bins_did,forcing_bins_did],/long)
log_iwp_sfc_net_forcing_obs_histo_total_id=ncdf_vardef(cdfid,'log_iwp_sfc_net_forcing_obs_histo_total',$
	[log_iwp_bins_did,forcing_bins_did],/long)
log_iwp_sfc_solar_forcing_obs_histo_total_id=ncdf_vardef(cdfid,'log_iwp_sfc_solar_forcing_obs_histo_total',$
	[log_iwp_bins_did,forcing_bins_did],/long)
log_iwp_sfc_IR_forcing_obs_histo_total_id=ncdf_vardef(cdfid,'log_iwp_sfc_IR_forcing_obs_histo_total',$
	[log_iwp_bins_did,forcing_bins_did],/long)
log_iwp_sfc_net_forcing_histo_total_id=ncdf_vardef(cdfid,'log_iwp_sfc_net_forcing_histo_total',$
	[log_iwp_bins_did,forcing_bins_did],/long)
log_iwp_sfc_IR_forcing_histo_total_id=ncdf_vardef(cdfid,'log_iwp_sfc_IR_forcing_histo_total',$
	[log_iwp_bins_did,forcing_bins_did],/long)
log_iwp_sfc_solar_forcing_histo_total_id=ncdf_vardef(cdfid,'log_iwp_sfc_solar_forcing_histo_total',$
	[log_iwp_bins_did,forcing_bins_did],/long)
log_iwp_toa_net_forcing_histo_total_id=ncdf_vardef(cdfid,'log_iwp_toa_net_forcing_histo_total',$
	[log_iwp_bins_did,forcing_bins_did],/long)
log_iwp_toa_solar_forcing_histo_total_id=ncdf_vardef(cdfid,'log_iwp_toa_solar_forcing_histo_total',$
	[log_iwp_bins_did,forcing_bins_did],/long)
log_iwp_toa_IR_forcing_histo_total_id=ncdf_vardef(cdfid,'log_iwp_toa_IR_forcing_histo_total',$
	[log_iwp_bins_did,forcing_bins_did],/long)

ncdf_control, cdfid, /endef
ncdf_varput, cdfid, forcing_bins_id, forcing_bins
ncdf_varput, cdfid, log_iwp_bins_id, log_iwp_bins

ncdf_varput, cdfid, log_iwp_toa_net_forcing_M_histo_total_id, log_iwp_toa_net_forcing_M_histo_total
ncdf_varput, cdfid, log_iwp_toa_solar_forcing_M_histo_total_id, log_iwp_toa_solar_forcing_M_histo_total
ncdf_varput, cdfid, log_iwp_toa_IR_forcing_M_histo_total_id,log_iwp_toa_IR_forcing_M_histo_total
ncdf_varput, cdfid, log_iwp_atm_net_forcing_histo_total_id, log_iwp_atm_net_forcing_histo_total
ncdf_varput, cdfid, log_iwp_atm_solar_forcing_histo_total_id, log_iwp_atm_solar_forcing_histo_total
ncdf_varput, cdfid, log_iwp_atm_IR_forcing_histo_total_id,log_iwp_atm_IR_forcing_histo_total
ncdf_varput, cdfid, log_iwp_sfc_net_forcing_obs_histo_total_id,log_iwp_sfc_net_forcing_obs_histo_total
ncdf_varput, cdfid, log_iwp_sfc_solar_forcing_obs_histo_total_id, log_iwp_sfc_solar_forcing_obs_histo_total
ncdf_varput, cdfid, log_iwp_sfc_IR_forcing_obs_histo_total_id, log_iwp_sfc_IR_forcing_obs_histo_total
ncdf_varput, cdfid, log_iwp_sfc_net_forcing_histo_total_id, log_iwp_sfc_net_forcing_histo_total
ncdf_varput, cdfid, log_iwp_sfc_IR_forcing_histo_total_id, log_iwp_sfc_IR_forcing_histo_total
ncdf_varput, cdfid, log_iwp_sfc_solar_forcing_histo_total_id, log_iwp_sfc_solar_forcing_histo_total
ncdf_varput, cdfid, log_iwp_toa_net_forcing_histo_total_id, log_iwp_toa_net_forcing_histo_total
ncdf_varput, cdfid, log_iwp_toa_solar_forcing_histo_total_id, log_iwp_toa_solar_forcing_histo_total
ncdf_varput, cdfid, log_iwp_toa_IR_forcing_histo_total_id, log_iwp_toa_IR_forcing_histo_total

endelse		;OPENR, 1, fname, ERROR = err


NCDF_CLOSE, cdfid

return
end