pro append_5minfiles_write_abbrev, 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]

;ir_bands_index_id=ncdf_dimid(cdfid,'ir_bands_index')
iwc_id=ncdf_varid(cdfid,'IWC_best_estimate')

if iwc_id lt 0 then begin


ncdf_control, cdfid, /verbose

; need to add dimensions for ir and solar raprad bands
;ir_bands_index_id=ncdf_dimdef(cdfid, 'ir_bands_index',n_elements(rap_ir_wl))
;solar_bands_index_id=ncdf_dimdef(cdfid, 'solar_bands_index',n_elements(rap_solar_wl))

;ncdf_control, cdfid, /endef
ncdf_control, cdfid, /redef

dim_indices2=[height_index_id,time_index_id]
dim_indices3=[height_index_id,time_index_id]


iwc_id=ncdf_vardef(cdfid,'IWC_best_estimate',dim_indices1,/short)

a_string='Best_Estimate_IWC'
ncdf_attput, cdfid, iwc_id, 'long_name', a_string
a_string='iwc (g/m^3)->10.^(float((*)/1000.)'
ncdf_attput, cdfid, iwc_id, 'units', a_string
a_string='>1.e-6 g/m^3, missing data: 8.e-8, no cloud: 9.e-9'
ncdf_attput, cdfid, iwc_id, 'valid_range', a_string
a_string='See IWC_source variable for information on how this quantity was computed'
ncdf_attput, cdfid, iwc_id, 'Comment1', a_string


re_ice_id=ncdf_vardef(cdfid,'re_ice_best_estimate',dim_indices1,/short)

a_string='Best_Estimate_ice_effective_radius'
ncdf_attput, cdfid, re_ice_id, 'long_name', a_string
a_string='effective_radius (m)->10.^(float((*)/1000.)'
ncdf_attput, cdfid, re_ice_id, 'units', a_string
a_string='>1.e-6 m, missing data: 8.e-8, no cloud: 9.e-9'
ncdf_attput, cdfid, re_ice_id, 'valid_range', a_string
a_string='See re_ice_source variable for information on how this quantity was computed'
ncdf_attput, cdfid, re_ice_id, 'Comment1', a_string

ice_error_id=ncdf_vardef(cdfid,'ice_micro_error_estimate',dim_indices1,/short)

a_string='estimated fractional error in the ice microphysical quantities'
ncdf_attput, cdfid, ice_error_id, 'long_name', a_string
a_string='ice_error (unitless fraction)->10.^(float((*)/1000.)'
ncdf_attput, cdfid, ice_error_id, 'units', a_string
a_string='>0., missing data: 8.e-8, no cloud: 9.e-9'
ncdf_attput, cdfid, ice_error_id, 'valid_range', a_string

iwc_source_id=ncdf_vardef(cdfid,'iwc_source',dim_indices1,/short)

a_string='Source index of the ice water content algorithm'
ncdf_attput, cdfid, iwc_source_id, 'long_name', a_string
a_string='none'
ncdf_attput, cdfid, iwc_source_id, 'units', a_string
a_string='1: mace et al 1998, 2: mace et al 2001, 3: ccm3 Parameterization (Kiehl et al 1998)'
ncdf_attput, cdfid, iwc_source_id, 'Comment1', a_string
a_string='4: Liu and Illingworth, 2000, 5: Dong et al 1998 parameterization
ncdf_attput, cdfid, iwc_source_id, 'Comment2', a_string
a_string='6: Dong and Mace 2001 retrieval, 7: frisch et al 1998 lwc retrieval
ncdf_attput, cdfid, iwc_source_id, 'Comment3', a_string

ice_re_source_id=ncdf_vardef(cdfid,'ice_effective_radius_source',dim_indices1,/short)

a_string='Source index of the effective radius algorithm'
ncdf_attput, cdfid, ice_re_source_id, 'long_name', a_string
a_string='none'
ncdf_attput, cdfid, ice_re_source_id, 'units', a_string
a_string='1: mace et al 1998, 2: mace et al 2001, 3: ccm3 Parameterization (Kiehl et al 1998)'
ncdf_attput, cdfid, ice_re_source_id, 'Comment1', a_string
a_string='4: Liu and Illingworth, 2000, 5: Dong et al 1998 parameterization
ncdf_attput, cdfid, ice_re_source_id, 'Comment2', a_string
a_string='6: Dong and Mace 2001 retrieval, 7: frisch et al 1998 lwc retrieval
ncdf_attput, cdfid, ice_re_source_id, 'Comment3', a_string

ice_tau_solar_id=ncdf_vardef(cdfid,'ice_tau_solar',dim_indices1,/short)

a_string='Spectral solar optical depth at 0.5 microns of the ice phase in this vertical bin'
ncdf_attput, cdfid, ice_tau_solar_id, 'long_name', a_string
a_string='optical depth (unitless)->10.^(float((*)/1000.)'
ncdf_attput, cdfid, ice_tau_solar_id, 'units', a_string
a_string='>0.0, missing data: 8.e-8, no cloud: 9.e-9'
ncdf_attput, cdfid, ice_tau_solar_id, 'valid_range', a_string
a_string='See ice_rad_param_source variable for information on how this quantity was computed'
ncdf_attput, cdfid, ice_tau_solar_id, 'Comment1', a_string

ice_omega_solar_id=ncdf_vardef(cdfid,'ice_omega_solar',dim_indices1,/short)

a_string='Spectral solar single scattering albedo at 0.5 microns  of the ice phase in this vertical bin'
ncdf_attput, cdfid, ice_omega_solar_id, 'long_name', a_string
a_string='omega (unitless)->10.^(float((*)/1000.)'
ncdf_attput, cdfid, ice_omega_solar_id, 'units', a_string
a_string='0.0<omega<1.0, missing data: 8.e-8, no cloud: 9.e-9'
ncdf_attput, cdfid, ice_omega_solar_id, 'valid_range', a_string
a_string='See ice_rad_param_source variable for information on how this quantity was computed'
ncdf_attput, cdfid, ice_omega_solar_id, 'Comment1', a_string

ice_g_solar_id=ncdf_vardef(cdfid,'ice_g_solar',dim_indices1,/short)

a_string='Spectral solar assymetry parameter at 0.5 microns of the ice phase in this vertical bin'
ncdf_attput, cdfid, ice_g_solar_id, 'long_name', a_string
a_string='omega (unitless)->10.^(float((*)/1000.)'
ncdf_attput, cdfid, ice_g_solar_id, 'units', a_string
a_string='-1.<omega<1.0, missing data: 8.e-8, no cloud: 9.e-9'
ncdf_attput, cdfid, ice_g_solar_id, 'valid_range', a_string
a_string='See ice_rad_param_source variable for information on how this quantity was computed'
ncdf_attput, cdfid, ice_g_solar_id, 'Comment1', a_string

ice_tau_IR_id=ncdf_vardef(cdfid,'ice_tau_IR',dim_indices1,/short)

a_string='Spectral IR optical depth at 11 microns of the ice phase in this vertical bin'
ncdf_attput, cdfid, ice_tau_IR_id, 'long_name', a_string
a_string='optical depth (unitless)->10.^(float((*)/1000.)'
ncdf_attput, cdfid, ice_tau_IR_id, 'units', a_string
a_string='>0.0, missing data: 8.e-8, no cloud: 9.e-9'
ncdf_attput, cdfid, ice_tau_IR_id, 'valid_range', a_string
a_string='See ice_rad_param_source variable for information on how this quantity was computed'
ncdf_attput, cdfid, ice_tau_IR_id, 'Comment1', a_string

ice_omega_IR_id=ncdf_vardef(cdfid,'ice_omega_IR',dim_indices1,/short)

a_string='Spectral IR single scattering albedo at 11 microns of the ice phase in this vertical bin'
ncdf_attput, cdfid, ice_omega_IR_id, 'long_name', a_string
a_string='omega (unitless)->10.^(float((*)/1000.)'
ncdf_attput, cdfid, ice_omega_IR_id, 'units', a_string
a_string='0.0<omega<1.0, missing data: 8.e-8, no cloud: 9.e-9'
ncdf_attput, cdfid, ice_omega_IR_id, 'valid_range', a_string
a_string='See ice_rad_param_source variable for information on how this quantity was computed'
ncdf_attput, cdfid, ice_omega_IR_id, 'Comment1', a_string

ice_g_IR_id=ncdf_vardef(cdfid,'ice_g_IR',dim_indices1,/short)

a_string='Spectral IR assymetry parameter at 11 microns of the ice phase in this vertical bin'
ncdf_attput, cdfid, ice_g_IR_id, 'long_name', a_string
a_string='omega (unitless)->10.^(float((*)/1000.)'
ncdf_attput, cdfid, ice_g_IR_id, 'units', a_string
a_string='-1.<omega<1.0, missing data: 8.e-8, no cloud: 9.e-9'
ncdf_attput, cdfid, ice_g_IR_id, 'valid_range', a_string
a_string='See ice_rad_param_source variable for information on how this quantity was computed'
ncdf_attput, cdfid, ice_g_IR_id, 'Comment1', a_string

ice_solar_radparam_source_id=ncdf_vardef(cdfid,'ice_RadiationParameterization_source_solar',dim_indices1,/short)

a_string='Source index of the solar radiative parameterization algorithm'
ncdf_attput, cdfid, ice_solar_radparam_source_id, 'long_name', a_string
a_string='none'
ncdf_attput, cdfid, ice_solar_radparam_source_id, 'units', a_string
a_string='1: Fu 1998, 2: Fu 1996, 3: ccm3 Parameterization (Kiehl et al 1998)'
ncdf_attput, cdfid, ice_solar_radparam_source_id, 'Comment1', a_string
a_string='4: Fu and Liou, 1993, 5: Ebert and Curry, 1992
ncdf_attput, cdfid, ice_solar_radparam_source_id, 'Comment2', a_string
a_string='6: Chyleck, 1995
ncdf_attput, cdfid, ice_solar_radparam_source_id, 'Comment3', a_string

ice_ir_radparam_source_id=ncdf_vardef(cdfid,'ice_RadiationParameterization_source_ir',dim_indices1,/short)

a_string='Source index of the IR radiative parameterization algorithm'
ncdf_attput, cdfid, ice_ir_radparam_source_id, 'long_name', a_string
a_string='none'
ncdf_attput, cdfid, ice_ir_radparam_source_id, 'units', a_string
a_string='1: Fu 1998, 2: Fu 1996, 3: ccm3 Parameterization (Kiehl et al 1998)'
ncdf_attput, cdfid, ice_ir_radparam_source_id, 'Comment1', a_string
a_string='4: Fu and Liou, 1993, 5: Ebert and Curry, 1992
ncdf_attput, cdfid, ice_ir_radparam_source_id, 'Comment2', a_string
a_string='6: Chyleck, 1995
ncdf_attput, cdfid, ice_ir_radparam_source_id, 'Comment3', a_string

lwc_id=ncdf_vardef(cdfid,'lwc_best_estimate',dim_indices1,/short)

a_string='Best_Estimate_lwc'
ncdf_attput, cdfid, lwc_id, 'long_name', a_string
a_string='lwc (g/m^3)->10.^(float((*)/1000.)'
ncdf_attput, cdfid, lwc_id, 'units', a_string
a_string='>1.e-6 g/m^3, missing data: 8.e-8, no cloud: 9.e-9'
ncdf_attput, cdfid, lwc_id, 'valid_range', a_string
a_string='See lwc_source variable for information on how this quantity was computed'
ncdf_attput, cdfid, lwc_id, 'Comment1', a_string

re_liquid_id=ncdf_vardef(cdfid,'re_liquid_best_estimate',dim_indices1,/short)

a_string='Best_Estimate_liquid_effective_radius'
ncdf_attput, cdfid, re_liquid_id, 'long_name', a_string
a_string='effective_radius (m)->10.^(float((*)/1000.)'
ncdf_attput, cdfid, re_liquid_id, 'units', a_string
a_string='>1.e-6 m, missing data: 8.e-8, no cloud: 9.e-9'
ncdf_attput, cdfid, re_liquid_id, 'valid_range', a_string
a_string='See re_liquid_source variable for information on how this quantity was computed'
ncdf_attput, cdfid, re_liquid_id, 'Comment1', a_string

liquid_error_id=ncdf_vardef(cdfid,'liquid_micro_error_estimate',dim_indices1,/short)

a_string='estimated fractional error in the liquid microphysical quantities'
ncdf_attput, cdfid, liquid_error_id, 'long_name', a_string
a_string='liquid_error (unitless fraction)->10.^(float((*)/1000.)'
ncdf_attput, cdfid, liquid_error_id, 'units', a_string
a_string='>0., missing data: 8.e-8, no cloud: 9.e-9'
ncdf_attput, cdfid, liquid_error_id, 'valid_range', a_string

lwc_source_id=ncdf_vardef(cdfid,'lwc_source',dim_indices1,/short)

a_string='Source index of the liquid water content algorithm'
ncdf_attput, cdfid, lwc_source_id, 'long_name', a_string
a_string='none'
ncdf_attput, cdfid, lwc_source_id, 'units', a_string
a_string='1: mace et al 1998, 2: mace et al 2001, 3: ccm3 Parameterization (Kiehl et al 1998)'
ncdf_attput, cdfid, lwc_source_id, 'Comment1', a_string
a_string='4: Liu and Illingworth, 2000, 5: Dong et al 1998 parameterization
ncdf_attput, cdfid, lwc_source_id, 'Comment2', a_string
a_string='6: Dong and Mace 2001 retrieval, 7: frisch et al 1998 lwc retrieval
ncdf_attput, cdfid, lwc_source_id, 'Comment3', a_string

liquid_re_source_id=ncdf_vardef(cdfid,'liquid_effective_radius_source',dim_indices1,/short)

a_string='Source index of the effective radius algorithm'
ncdf_attput, cdfid, liquid_re_source_id, 'long_name', a_string
a_string='none'
ncdf_attput, cdfid, liquid_re_source_id, 'units', a_string
a_string='1: mace et al 1998, 2: mace et al 2001, 3: ccm3 Parameterization (Kiehl et al 1998)'
ncdf_attput, cdfid, liquid_re_source_id, 'Comment1', a_string
a_string='4: Liu and Illingworth, 2000, 5: Dong et al 1998 parameterization
ncdf_attput, cdfid, liquid_re_source_id, 'Comment2', a_string
a_string='6: Dong and Mace 2001 retrieval, 7: frisch et al 1998 lwc retrieval
ncdf_attput, cdfid, liquid_re_source_id, 'Comment3', a_string

liquid_tau_solar_id=ncdf_vardef(cdfid,'liquid_tau_solar',dim_indices1,/short)

a_string='Spectral solar optical depth at 0.5 microns of the liquid phase in this vertical bin'
ncdf_attput, cdfid, liquid_tau_solar_id, 'long_name', a_string
a_string='optical depth (unitless)->10.^(float((*)/1000.)'
ncdf_attput, cdfid, liquid_tau_solar_id, 'units', a_string
a_string='>0.0, missing data: 8.e-8, no cloud: 9.e-9'
ncdf_attput, cdfid, liquid_tau_solar_id, 'valid_range', a_string
a_string='See liquid_rad_param_source variable for information on how this quantity was computed'
ncdf_attput, cdfid, liquid_tau_solar_id, 'Comment1', a_string

liquid_omega_solar_id=ncdf_vardef(cdfid,'liquid_omega_solar',dim_indices1,/short)

a_string='Spectral solar single scattering albedo at 0.5 microns of the liquid phase in this vertical bin'
ncdf_attput, cdfid, liquid_omega_solar_id, 'long_name', a_string
a_string='omega (unitless)->10.^(float((*)/1000.)'
ncdf_attput, cdfid, liquid_omega_solar_id, 'units', a_string
a_string='0.0<omega<1.0, missing data: 8.e-8, no cloud: 9.e-9'
ncdf_attput, cdfid, liquid_omega_solar_id, 'valid_range', a_string
a_string='See liquid_rad_param_source variable for information on how this quantity was computed'
ncdf_attput, cdfid, liquid_omega_solar_id, 'Comment1', a_string

liquid_g_solar_id=ncdf_vardef(cdfid,'liquid_g_solar',dim_indices1,/short)

a_string='Spectral solar assymetry parameter at 0.5 microns of the liquid phase in this vertical bin'
ncdf_attput, cdfid, liquid_g_solar_id, 'long_name', a_string
a_string='omega (unitless)->10.^(float((*)/1000.)'
ncdf_attput, cdfid, liquid_g_solar_id, 'units', a_string
a_string='-1.<omega<1.0, missing data: 8.e-8, no cloud: 9.e-9'
ncdf_attput, cdfid, liquid_g_solar_id, 'valid_range', a_string
a_string='See liquid_rad_param_source variable for information on how this quantity was computed'
ncdf_attput, cdfid, liquid_g_solar_id, 'Comment1', a_string

liquid_tau_IR_id=ncdf_vardef(cdfid,'liquid_tau_IR',dim_indices1,/short)

a_string='Spectral IR optical depth at 11 microns of the liquid phase in this vertical bin'
ncdf_attput, cdfid, liquid_tau_IR_id, 'long_name', a_string
a_string='optical depth (unitless)->10.^(float((*)/1000.)'
ncdf_attput, cdfid, liquid_tau_IR_id, 'units', a_string
a_string='>0.0, missing data: 8.e-8, no cloud: 9.e-9'
ncdf_attput, cdfid, liquid_tau_IR_id, 'valid_range', a_string
a_string='See liquid_rad_param_source variable for information on how this quantity was computed'
ncdf_attput, cdfid, liquid_tau_IR_id, 'Comment1', a_string

liquid_omega_IR_id=ncdf_vardef(cdfid,'liquid_omega_IR',dim_indices1,/short)

a_string='Spectral IR single scattering albedo at 11 microns of the liquid phase in this vertical bin'
ncdf_attput, cdfid, liquid_omega_IR_id, 'long_name', a_string
a_string='omega (unitless)->10.^(float((*)/1000.)'
ncdf_attput, cdfid, liquid_omega_IR_id, 'units', a_string
a_string='0.0<omega<1.0, missing data: 8.e-8, no cloud: 9.e-9'
ncdf_attput, cdfid, liquid_omega_IR_id, 'valid_range', a_string
a_string='See liquid_rad_param_source variable for information on how this quantity was computed'
ncdf_attput, cdfid, liquid_omega_IR_id, 'Comment1', a_string

liquid_g_IR_id=ncdf_vardef(cdfid,'liquid_g_IR',dim_indices1,/short)

a_string='Spectral IR assymetry parameter at 11 microns of the liquid phase in this vertical bin'
ncdf_attput, cdfid, liquid_g_IR_id, 'long_name', a_string
a_string='omega (unitless)->10.^(float((*)/1000.)'
ncdf_attput, cdfid, liquid_g_IR_id, 'units', a_string
a_string='-1.<omega<1.0, missing data: 8.e-8, no cloud: 9.e-9'
ncdf_attput, cdfid, liquid_g_IR_id, 'valid_range', a_string
a_string='See liquid_rad_param_source variable for information on how this quantity was computed'
ncdf_attput, cdfid, liquid_g_IR_id, 'Comment1', a_string

liquid_solar_radparam_source_id=ncdf_vardef(cdfid,'liquid_RadiationParameterization_source_solar',dim_indices1,/short)

a_string='Source index of the solar radiative parameterization algorithm'
ncdf_attput, cdfid, liquid_solar_radparam_source_id, 'long_name', a_string
a_string='none'
ncdf_attput, cdfid, liquid_solar_radparam_source_id, 'units', a_string
a_string='1: Fu 1998, 2: Fu 1996, 3: ccm3 Parameterization (Kiehl et al 1998)'
ncdf_attput, cdfid, liquid_solar_radparam_source_id, 'Comment1', a_string
a_string='4: Fu and Liou, 1993, 5: Ebert and Curry, 1992
ncdf_attput, cdfid, liquid_solar_radparam_source_id, 'Comment2', a_string
a_string='6: Chyleck, 1995
ncdf_attput, cdfid, liquid_solar_radparam_source_id, 'Comment3', a_string

liquid_ir_radparam_source_id=ncdf_vardef(cdfid,'liquid_RadiationParameterization_source_ir',dim_indices1,/short)

a_string='Source index of the ir radiative parameterization algorithm'
ncdf_attput, cdfid, liquid_ir_radparam_source_id, 'long_name', a_string
a_string='none'
ncdf_attput, cdfid, liquid_ir_radparam_source_id, 'units', a_string
a_string='1: Fu 1998, 2: Fu 1996, 3: ccm3 Parameterization (Kiehl et al 1998)'
ncdf_attput, cdfid, liquid_ir_radparam_source_id, 'Comment1', a_string
a_string='4: Fu and Liou, 1993, 5: Ebert and Curry, 1992
ncdf_attput, cdfid, liquid_ir_radparam_source_id, 'Comment2', a_string
a_string='6: Chyleck, 1995
ncdf_attput, cdfid, liquid_ir_radparam_source_id, 'Comment3', a_string

ice_fraction_id=ncdf_vardef(cdfid,'ice_fraction_CCM3',dim_indices1,/short)

a_string='Fraction of the water content that is ice - CCM3 param.  Kiehl et al, 1998'
ncdf_attput, cdfid, ice_fraction_id, 'long_name', a_string
a_string='fraction (unitless)->10.^(float((*)/1000.)'
ncdf_attput, cdfid, ice_fraction_id, 'units', a_string
a_string='0.<omega<1.0, missing data: 8.e-8, no cloud: 9.e-9'
ncdf_attput, cdfid, ice_fraction_id, 'valid_range', a_string

ncdf_control, cdfid, /endef

endif else begin	; variables and dimensions already defined, get id's to write

;ir_bands_index_id=ncdf_dimid(cdfid,'ir_bands_index')
;solar_bands_index_id=ncdf_dimid(cdfid,'solar_bands_index')

;solar_bands_id=ncdf_varid(cdfid,'RAPRAD_solar_bands')
;ir_bands_id=ncdf_varid(cdfid,'RAPRAD_ir_bands')
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')
re_liquid_id=ncdf_varid(cdfid,'re_liquid_best_estimate')
liquid_error_id=ncdf_varid(cdfid,'liquid_micro_error_estimate')
lwc_source_id=ncdf_varid(cdfid,'lwc_source')
liquid_re_source_id=ncdf_varid(cdfid,'liquid_effective_radius_source')
liquid_tau_solar_id=ncdf_varid(cdfid,'liquid_tau_solar')
liquid_omega_solar_id=ncdf_varid(cdfid,'liquid_omega_solar')
liquid_g_solar_id=ncdf_varid(cdfid,'liquid_g_solar')
liquid_tau_IR_id=ncdf_varid(cdfid,'liquid_tau_IR')
liquid_omega_IR_id=ncdf_varid(cdfid,'liquid_omega_IR')
liquid_g_IR_id=ncdf_varid(cdfid,'liquid_g_IR')
liquid_solar_radparam_source_id=ncdf_varid(cdfid,'liquid_RadiationParameterization_source_solar')
liquid_ir_radparam_source_id=ncdf_varid(cdfid,'liquid_RadiationParameterization_source_ir')
ice_fraction_id=ncdf_varid(cdfid,'ice_fraction_CCM3')

endelse

;solar_bands_id,ir_bands_id,iwc_id,re_ice_id,ice_error_id,iwc_source_id,ice_re_source_id,$
;ice_tau_solar_id,ice_omega_solar_id,ice_g_solar_id,ice_tau_IR_id,ice_omega_IR_id, $
;ice_g_IR_id,ice_radparam_source_id,lwc_id,re_liquid_id,liquid_error_id,lwc_source_id, $
;liquid_re_source_id,liquid_tau_solar_id, liquid_omega_solar_id,liquid_g_solar_id,liquid_tau_IR_id, $
;liquid_omega_IR_id,liquid_g_IR_id,liquid_radparam_source_id

;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

;ncdf_varput, cdfid, solar_bands_id,rap_solar_wl
;ncdf_varput, cdfid, ir_bands_id, rap_ir_wl

; find the indexes of the solar and ir spectral bands,

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,*]
if max(n_elements(where(fix(1000.*(alog10(temp_array))) gt 2000.))) gt 0 then begin
	print, temp_array[where(fix(1000.*(alog10(temp_array))) gt 2000.)]
endif
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_varput, cdfid, re_liquid_id, fix(1000.*(alog10(re_liq_out)))
ncdf_varput, cdfid, liquid_error_id, fix(1000.*(alog10(liqmicro_error)))
ncdf_varput, cdfid, lwc_source_id, lwc_source
ncdf_varput, cdfid, liquid_re_source_id,liq_size_source
for j=0,num_heights-1 do temp_array[j,*]=liq_tau_solar[index_solar[0],j,*]
ncdf_varput, cdfid, liquid_tau_solar_id, fix(1000.*(alog10(temp_array)))
for j=0,num_heights-1 do temp_array[j,*]=liq_omega_solar[index_solar[0],j,*]
ncdf_varput, cdfid, liquid_omega_solar_id, fix(1000.*(alog10(temp_array)))
for j=0,num_heights-1 do temp_array[j,*]=liq_g_solar[index_solar[0],j,*]
ncdf_varput, cdfid, liquid_g_solar_id, fix(1000.*(alog10(temp_array)))
for j=0,num_heights-1 do temp_array[j,*]=liq_tau_ir[index_ir[0],j,*]
ncdf_varput, cdfid, liquid_tau_IR_id, fix(1000.*(alog10(temp_array)))
for j=0,num_heights-1 do temp_array[j,*]=liq_omega_ir[index_ir[0],j,*]
ncdf_varput, cdfid, liquid_omega_IR_id, fix(1000.*(alog10(temp_array)))
for j=0,num_heights-1 do temp_array[j,*]=liq_g_ir[index_ir[0],j,*]
ncdf_varput, cdfid, liquid_g_IR_id, fix(1000.*(alog10(temp_array)))
ncdf_varput, cdfid, liquid_solar_radparam_source_id, liq_rad_param_source[*,*,0]
ncdf_varput, cdfid, liquid_solar_radparam_source_id, liq_rad_param_source[*,*,1]
ncdf_varput, cdfid, ice_fraction_id, fix(1000.*(alog10(ice_fraction)))



;RAPRAD_solar_bands,RAPRAD_ir_bands,IWC_best_estimate,re_ice_best_estimate,ice_micro_error_estimate,iwc_source
;ice_effective_radius_source,ice_tau_solar,ice_omega_solar,ice_g_solar,ice_tau_IR,ice_omega_IR,ice_g_IR
;ice_RadiationParameterization_source,lwc_best_estimate,re_liquid_best_estimate,liquid_micro_error_estimate
;lwc_source,liquid_effective_radius_source,liquid_tau_solar,liquid_omega_solar,liquid_g_solar
;liquid_tau_IR,liquid_omega_IR,liquid_g_IR,liquid_RadiationParameterization_source

ncdf_close, cdfid

return
end