pro append_5minfiles_write_abbrev_icemicro, fname

; this routine writes the radiative properties at only 0.5 microns and 11 microns.

common cldmicro_output, iwc_out, re_ice_out, icemicro_error, iwc_source, ice_size_source, ice_tau_solar, $
ice_omega_solar, ice_g_solar,ice_tau_ir, ice_omega_ir,ice_g_ir,ice_rad_param_source, lwc_out, re_liq_out, liqmicro_error, $
lwc_source, liq_size_source, liq_tau_solar, liq_omega_solar, liq_g_solar, liq_tau_ir, liq_omega_ir, $
liq_g_ir,liq_rad_param_source, ice_fraction, rap_solar_wl, rap_ir_wl


if n_elements(where(iwc_out eq -9999.)) gt 1 then iwc_out[where(iwc_out eq -9999.)]=9.e-9
if n_elements(where(iwc_out eq -8888.)) gt 1 then iwc_out[where(iwc_out eq -8888.)]=8.e-8
if n_elements(where(iwc_out le 0.)) gt 1 then iwc_out[where(iwc_out le 0.)]=9.e-9

if n_elements(where(re_ice_out eq -9999.)) gt 1 then re_ice_out[where(re_ice_out eq -9999.)]=9.e-9
if n_elements(where(re_ice_out eq -8888.)) gt 1 then re_ice_out[where(re_ice_out eq -8888.)]=8.e-8
if n_elements(where(re_ice_out le 0.)) gt 1 then re_ice_out[where(re_ice_out le 0.)]=9.e-9

if n_elements(where(icemicro_error eq -9999.)) gt 1 then icemicro_error[where(icemicro_error eq -9999.)]=9.e-9
if n_elements(where(icemicro_error eq -8888.)) gt 1 then icemicro_error[where(icemicro_error eq -8888.)]=8.e-8
if n_elements(where(icemicro_error le 0.)) gt 1 then icemicro_error[where(icemicro_error le 0.)]=9.e-9

if n_elements(where(ice_tau_solar eq -9999.)) gt 1 then ice_tau_solar[where(ice_tau_solar eq -9999.)]=9.e-9
if n_elements(where(ice_tau_solar eq -8888.)) gt 1 then ice_tau_solar[where(ice_tau_solar eq -8888.)]=8.e-8
if n_elements(where(ice_tau_solar le 0.)) gt 1 then ice_tau_solar[where(ice_tau_solar le 0.)]=9.e-9

if n_elements(where(ice_omega_solar eq -9999.)) gt 1 then ice_omega_solar[where(ice_omega_solar eq -9999.)]=9.e-9
if n_elements(where(ice_omega_solar eq -8888.)) gt 1 then ice_omega_solar[where(ice_omega_solar eq -8888.)]=8.e-8
if n_elements(where(ice_omega_solar le 0.)) gt 1 then ice_omega_solar[where(ice_omega_solar le 0.)]=9.e-9

if n_elements(where(ice_g_solar eq -9999.)) gt 1 then ice_g_solar[where(ice_g_solar eq -9999.)]=9.e-9
if n_elements(where(ice_g_solar eq -8888.)) gt 1 then ice_g_solar[where(ice_g_solar eq -8888.)]=8.e-8
if n_elements(where(ice_g_solar le 0.)) gt 1 then ice_g_solar[where(ice_g_solar le 0.)]=9.e-9

if n_elements(where(ice_tau_ir eq -9999.)) gt 1 then ice_tau_ir[where(ice_tau_ir eq -9999.)]=9.e-9
if n_elements(where(ice_tau_ir eq -8888.)) gt 1 then ice_tau_ir[where(ice_tau_ir eq -8888.)]=8.e-8
if n_elements(where(ice_tau_ir le 0.)) gt 1 then ice_tau_ir[where(ice_tau_ir le 0.)]=9.e-9

if n_elements(where(ice_omega_ir eq -9999.)) gt 1 then ice_omega_ir[where(ice_omega_ir le 0.)]=9.e-9
if n_elements(where(ice_omega_ir eq -8888.)) gt 1 then ice_omega_ir[where(ice_omega_ir eq -8888.)]=8.e-8
if n_elements(where(ice_omega_ir le 0.)) gt 1 then ice_omega_ir[where(ice_omega_ir le 0.)]=9.e-9

if n_elements(where(ice_g_ir eq -9999.)) gt 1 then ice_g_ir[where(ice_g_ir eq -9999.)]=9.e-9
if n_elements(where(ice_g_ir eq -8888.)) gt 1 then ice_g_ir[where(ice_g_ir eq -8888.)]=8.e-8
if n_elements(where(ice_g_ir le 0.)) gt 1 then ice_g_ir[where(ice_g_ir le 0.)]=9.e-9

if n_elements(where(lwc_out eq -9999.)) gt 1 then lwc_out[where(lwc_out eq -9999.)]=9.e-9
if n_elements(where(lwc_out eq -8888.)) gt 1 then lwc_out[where(lwc_out eq -8888.)]=8.e-8
if n_elements(where(lwc_out le 0.)) gt 1 then lwc_out[where(lwc_out le 0.)]=9.e-9

if n_elements(where(re_liq_out eq -9999.)) gt 1 then re_liq_out[where(re_liq_out eq -9999.)]=9.e-9
if n_elements(where(re_liq_out eq -8888.)) gt 1 then re_liq_out[where(re_liq_out eq -8888.)]=8.e-8
if n_elements(where(re_liq_out le 0.)) gt 1 then re_liq_out[where(re_liq_out le 0.)]=9.e-9

if n_elements(where(liqmicro_error eq -9999.)) gt 1 then liqmicro_error[where(liqmicro_error eq -9999.)]=9.e-9
if n_elements(where(liqmicro_error eq -8888.)) gt 1 then liqmicro_error[where(liqmicro_error eq -8888.)]=8.e-8
if n_elements(where(liqmicro_error le 0.)) gt 1 then liqmicro_error[where(liqmicro_error le 0.)]=9.e-9

if n_elements(where(liq_tau_solar eq -9999.)) gt 1 then liq_tau_solar[where(liq_tau_solar eq -9999.)]=9.e-9
if n_elements(where(liq_tau_solar eq -8888.)) gt 1 then liq_tau_solar[where(liq_tau_solar eq -8888.)]=8.e-8
if n_elements(where(liq_tau_solar le 0.)) gt 1 then liq_tau_solar[where(liq_tau_solar le 0.)]=9.e-9

if n_elements(where(liq_omega_solar eq -9999.)) gt 1 then liq_omega_solar[where(liq_omega_solar eq -9999.)]=9.e-9
if n_elements(where(liq_omega_solar eq -8888.)) gt 1 then liq_omega_solar[where(liq_omega_solar eq -8888.)]=8.e-8
if n_elements(where(liq_omega_solar le 0.)) gt 1 then liq_omega_solar[where(liq_omega_solar le 0.)]=9.e-9

if n_elements(where(liq_g_solar eq -9999.)) gt 1 then liq_g_solar[where(liq_g_solar eq -9999.)]=9.e-9
if n_elements(where(liq_g_solar eq -8888.)) gt 1 then liq_g_solar[where(liq_g_solar eq -8888.)]=8.e-8
if n_elements(where(liq_g_solar le 0.)) gt 1 then liq_g_solar[where(liq_g_solar le 0.)]=9.e-9

if n_elements(where(liq_tau_ir eq -9999.)) gt 1 then liq_tau_ir[where(liq_tau_ir eq -9999.)]=9.e-9
if n_elements(where(liq_tau_ir eq -8888.)) gt 1 then liq_tau_ir[where(liq_tau_ir eq -8888.)]=8.e-8
if n_elements(where(liq_tau_ir le 0.)) gt 1 then liq_tau_ir[where(liq_tau_ir le 0.)]=9.e-9

if n_elements(where(liq_omega_ir eq -9999.)) gt 1 then liq_omega_ir[where(liq_omega_ir eq -9999.)]=9.e-9
if n_elements(where(liq_omega_ir eq -8888.)) gt 1 then liq_omega_ir[where(liq_omega_ir eq -8888.)]=8.e-8
if n_elements(where(liq_omega_ir le 0.)) gt 1 then liq_omega_ir[where(liq_omega_ir le 0.)]=9.e-9

if n_elements(where(liq_g_ir eq -9999.)) gt 1 then liq_g_ir[where(liq_g_ir eq -9999.)]=9.e-9
if n_elements(where(liq_g_ir eq -8888.)) gt 1 then liq_g_ir[where(liq_g_ir eq -8888.)]=8.e-8
if n_elements(where(liq_g_ir le 0.)) gt 1 then liq_g_ir[where(liq_g_ir le 0.)]=9.e-9

if n_elements(where(ice_fraction eq -9999.)) gt 1 then ice_fraction[where(ice_fraction eq -9999.)]=9.e-9
if n_elements(where(ice_fraction eq -8888.)) gt 1 then ice_fraction[where(ice_fraction eq -8888.)]=8.e-8
if n_elements(where(ice_fraction le 0.)) gt 1 then ice_fraction[where(ice_fraction le 0.)]=9.e-9

cdfid=ncdf_open(fname,/write)

time_index_id=ncdf_dimid(cdfid,'time')
height_index_id=ncdf_dimid(cdfid,'height')

ncdf_diminq, cdfid, height_index_id, char_strng, num_heights
ncdf_diminq, cdfid, time_index_id, char_strng, num_times

dim_indices1=[height_index_id,time_index_id]

iwc_id=ncdf_varid(cdfid,'IWC_best_estimate')
re_ice_id=ncdf_varid(cdfid,'re_ice_best_estimate')
ice_error_id=ncdf_varid(cdfid,'ice_micro_error_estimate')
iwc_source_id=ncdf_varid(cdfid,'iwc_source')
ice_re_source_id=ncdf_varid(cdfid,'ice_effective_radius_source')
ice_tau_solar_id=ncdf_varid(cdfid,'ice_tau_solar')
ice_omega_solar_id=ncdf_varid(cdfid,'ice_omega_solar')
ice_g_solar_id=ncdf_varid(cdfid,'ice_g_solar')
ice_tau_IR_id=ncdf_varid(cdfid,'ice_tau_IR')
ice_omega_IR_id=ncdf_varid(cdfid,'ice_omega_IR')
ice_g_IR_id=ncdf_varid(cdfid,'ice_g_IR')
ice_solar_radparam_source_id=ncdf_varid(cdfid,'ice_RadiationParameterization_source_solar')
ice_ir_radparam_source_id=ncdf_varid(cdfid,'ice_RadiationParameterization_source_ir')
lwc_id=ncdf_varid(cdfid,'lwc_best_estimate')

index_solar=where(rap_solar_wl ge 0.5)
index_ir=where(rap_ir_wl ge 11.)
temp_array=fltarr(num_heights,num_times)

ncdf_varput, cdfid, iwc_id, fix(1000.*(alog10(iwc_out)))
ncdf_varput, cdfid, re_ice_id, fix(1000.*(alog10(re_ice_out)))
ncdf_varput, cdfid, ice_error_id, fix(1000.*(alog10(icemicro_error)))
ncdf_varput, cdfid, iwc_source_id, iwc_source
ncdf_varput, cdfid, ice_re_source_id, ice_size_source
for j=0,num_heights-1 do temp_array[j,*]=ice_tau_solar[index_solar[0],j,*]
ncdf_varput, cdfid, ice_tau_solar_id, fix(1000.*(alog10(temp_array)))
for j=0,num_heights-1 do temp_array[j,*]=ice_omega_solar[index_solar[0],j,*]
ncdf_varput, cdfid, ice_omega_solar_id, fix(1000.*(alog10(temp_array)))
for j=0,num_heights-1 do temp_array[j,*]=ice_g_solar[index_solar[0],j,*]
ncdf_varput, cdfid, ice_g_solar_id,fix(1000.*(alog10(temp_array)))
for j=0,num_heights-1 do temp_array[j,*]=ice_tau_ir[index_ir[0],j,*]
ncdf_varput, cdfid, ice_tau_IR_id, fix(1000.*(alog10(temp_array)))
for j=0,num_heights-1 do temp_array[j,*]=ice_omega_ir[index_ir[0],j,*]
ncdf_varput, cdfid, ice_omega_IR_id, fix(1000.*(alog10(temp_array)))
for j=0,num_heights-1 do temp_array[j,*]=ice_g_ir[index_ir[0],j,*]
ncdf_varput, cdfid, ice_g_IR_id, fix(1000.*(alog10(temp_array)))
ncdf_varput, cdfid, ice_solar_radparam_source_id, ice_rad_param_source[*,*,0]
ncdf_varput, cdfid, ice_ir_radparam_source_id, ice_rad_param_source[*,*,1]
ncdf_varput, cdfid, lwc_id, fix(1000.*(alog10(lwc_out)))

ncdf_close, cdfid

return
end