; ; This program reads in the averaged data from the ; .comb file. It is mostly radiation and retrieval data ; ; ; Read in the average data ; pro get_average_data_rad,site,searchdate,avg_int ; ; common containing the data ; common monthly_averaged_data_rad, base_time,time_offset,height,numtimes,numheights,$ 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_solar,ice_RadiationParameterization_source_ir, $ 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_solar,$ liquid_RadiationParameterization_source_ir,$ ice_fraction_CCM3,solar_flux_down_clear,solar_flux_up_clear, $ IR_flux_up_clear,IR_flux_down_clear,downwelling_solar_diffuse_clear, $ downwelling_solar_direct_clear,TOA_solar_up_clear,TOA_solar_down_clear, $ TOA_IR_up_clear,solar_flux_down,solar_flux_up,IR_flux_up,IR_flux_down, $ downwelling_solar_diffuse,downwelling_solar_direct,TOA_solar_up,TOA_solar_down, $ TOA_IR_up ; ; Date to be analyzed ; curdate=strcompress(string(searchdate),/remove_all) month=strmid(curdate,4,2) site=strcompress(site,/remove_all) ; ; Open the netcdf file ; if site eq 'sgp' then begin armdatadir='/data/mace3/arm_data/sgp/sgp.average.'+avg_int+'/' if avg_int eq '1hr' then armdatadir=armdatadir+'full_ret/' arm_file='sgp.average.'+avg_int+'.'+curdate+'*comb.cdf*' endif else if site eq 'nsa' then begin armdatadir='/data/mace2/arm_data/nsa/nsa.average.'+avg_int+'/' arm_file='nsa.average.'+avg_int+'.'+curdate+'*comb.cdf*' endif else if site eq 'twpc1' then begin armdatadir='/data/mace/arm_data/twp/Manus_C1/twpc1.average.'+avg_int+'/' arm_file='twpc1.average.'+avg_int+'.'+curdate+'*comb.cdf*' endif else if site eq 'twpc2' then begin armdatadir='/data/mace/arm_data/twp/Nauru_C2/twpc2.average.'+avg_int+'/' arm_file='twpc2.average.'+avg_int+'.'+curdate+'*comb.cdf*' endif ; ; Unzip the data files ; unix_command='gunzip -fq '+armdatadir+arm_file spawn,unix_command ; ; Get a list of the matching file names ; print, 'filename ',arm_file files=findfile(armdatadir+arm_file,count=num_files) print, 'Found ',num_files,' files' ; ; Loop through all the matching files in the directory ; j=-1 ;counts number of files containing the variable for i=0,num_files-1 do begin ; ; Don't try to process any gif files that may be in the directory ; if strpos(files[i],'gif') eq -1 and $ strpos(files[i],'png') eq -1 then begin j=j+1 print,i,j,files[i] ; ; Open the netcdf file ; cdf_id=ncdf_open(files[i]) ; ; Get the dimension id's ; time_did=ncdf_dimid(cdf_id,'time') height_did=ncdf_dimid(cdf_id,'height') ; ; Get the dimensions ; if time_did ge 0 then ncdf_diminq,cdf_id,time_did,charstring,numtimes if height_did ge 0 then ncdf_diminq,cdf_id,height_did,charstring,numheights ; ; Create some dummy variables ; dummy1d=make_array(numtimes,/float,value=-8888) dummy2d=make_array(numheights,numtimes,/float,value=-8888) ; ; Get the variable id's ; base_time_id=ncdf_varid(cdf_id,'base_time') time_offset_id=ncdf_varid(cdf_id,'time_offset') height_id=ncdf_varid(cdf_id,'height') IWC_best_estimate_id=ncdf_varid(cdf_id,'IWC_best_estimate') re_ice_best_estimate_id=ncdf_varid(cdf_id,'re_ice_best_estimate') ice_micro_error_estimate_id=ncdf_varid(cdf_id,'ice_micro_error_estimate') iwc_source_id=ncdf_varid(cdf_id,'iwc_source') ice_effective_radius_source_id=ncdf_varid(cdf_id,'ice_effective_radius_source') ice_tau_solar_id=ncdf_varid(cdf_id,'ice_tau_solar') ice_omega_solar_id=ncdf_varid(cdf_id, 'ice_omega_solar') ice_g_solar_id=ncdf_varid(cdf_id, 'ice_g_solar') ice_tau_IR_id=ncdf_varid(cdf_id, 'ice_tau_IR') ice_omega_IR_id=ncdf_varid(cdf_id, 'ice_omega_IR') ice_g_IR_id=ncdf_varid(cdf_id, 'ice_g_IR') ice_RadiationParameterization_source_solar_id= $ ncdf_varid(cdf_id,'ice_RadiationParameterization_source_solar') ice_RadiationParameterization_source_ir_id= $ ncdf_varid(cdf_id,'ice_RadiationParameterization_source_ir') lwc_best_estimate_id=ncdf_varid(cdf_id, 'lwc_best_estimate') re_liquid_best_estimate_id=ncdf_varid(cdf_id, 're_liquid_best_estimate') liquid_micro_error_estimate_id=ncdf_varid(cdf_id,'liquid_micro_error_estimate') lwc_source_id=ncdf_varid(cdf_id,'lwc_source') liquid_effective_radius_source_id= $ ncdf_varid(cdf_id,'liquid_effective_radius_source') liquid_tau_solar_id=ncdf_varid(cdf_id,'liquid_tau_solar') liquid_omega_solar_id=ncdf_varid(cdf_id,'liquid_omega_solar') liquid_g_solar_id=ncdf_varid(cdf_id,'liquid_g_solar') liquid_tau_IR_id=ncdf_varid(cdf_id,'liquid_tau_IR') liquid_omega_IR_id=ncdf_varid(cdf_id,'liquid_omega_IR') liquid_g_IR_id=ncdf_varid(cdf_id,'liquid_g_IR') liquid_RadiationParameterization_source_solar_id= $ ncdf_varid(cdf_id,'liquid_RadiationParameterization_source_solar') liquid_RadiationParameterization_source_ir_id= $ ncdf_varid(cdf_id,'liquid_RadiationParameterization_source_ir') ice_fraction_CCM3_id=ncdf_varid(cdf_id,'ice_fraction_CCM3') solar_flux_down_clear_id=ncdf_varid(cdf_id,'solar_flux_down_clear') solar_flux_up_clear_id=ncdf_varid(cdf_id,'solar_flux_up_clear') IR_flux_up_clear_id=ncdf_varid(cdf_id,'IR_flux_up_clear') IR_flux_down_clear_id=ncdf_varid(cdf_id,'IR_flux_down_clear') downwelling_solar_diffuse_clear_id=ncdf_varid(cdf_id,'downwelling_solar_diffuse_clear') downwelling_solar_direct_clear_id=ncdf_varid(cdf_id,'downwelling_solar_direct_clear') TOA_solar_up_clear_id=ncdf_varid(cdf_id,'TOA_solar_up_clear') TOA_solar_down_clear_id=ncdf_varid(cdf_id,'TOA_solar_down_clear') TOA_IR_up_clear_id=ncdf_varid(cdf_id,'TOA_IR_up_clear') solar_flux_down_id=ncdf_varid(cdf_id,'solar_flux_down') solar_flux_up_id=ncdf_varid(cdf_id,'solar_flux_up') IR_flux_up_id=ncdf_varid(cdf_id,'IR_flux_up') IR_flux_down_id=ncdf_varid(cdf_id,'IR_flux_down') downwelling_solar_diffuse_id=ncdf_varid(cdf_id,'downwelling_solar_diffuse') downwelling_solar_direct_id=ncdf_varid(cdf_id,'downwelling_solar_direct') TOA_solar_up_id=ncdf_varid(cdf_id,'TOA_solar_up') TOA_solar_down_id=ncdf_varid(cdf_id,'TOA_solar_down') TOA_IR_up_id=ncdf_varid(cdf_id,'TOA_IR_up') ; ; Get the variables ; if base_time_id ne -1 then ncdf_varget,cdf_id,base_time_id,base_time if time_offset_id ne -1 then ncdf_varget,cdf_id,time_offset_id,time_offset if height_id ne -1 then ncdf_varget,cdf_id,height_id,height if IWC_best_estimate_id ne -1 then $ ncdf_varget,cdf_id,IWC_best_estimate_id,IWC_best_estimate $ else IWC_best_estimate=dummy2d IWC_best_estimate=transpose(IWC_best_estimate) if re_ice_best_estimate_id ne -1 then $ ncdf_varget,cdf_id,re_ice_best_estimate_id,re_ice_best_estimate $ else re_ice_best_estimate=dummy2d re_ice_best_estimate=transpose(re_ice_best_estimate) if ice_micro_error_estimate_id ne -1 then $ ncdf_varget,cdf_id,ice_micro_error_estimate_id,ice_micro_error_estimate $ else ice_micro_error_estimate=dummy2d ice_micro_error_estimate=transpose(ice_micro_error_estimate) if iwc_source_id ne -1 then $ ncdf_varget,cdf_id,iwc_source_id,iwc_source $ else iwc_source=dummy2d iwc_source=transpose(iwc_source) if ice_effective_radius_source_id ne -1 then $ ncdf_varget,cdf_id,ice_effective_radius_source_id,ice_effective_radius_source $ else ice_effective_radius_source=dummy2d ice_effective_radius_source=transpose(ice_effective_radius_source) if ice_tau_solar_id ne -1 then $ ncdf_varget,cdf_id,ice_tau_solar_id,ice_tau_solar $ else ice_tau_solar=dummy2d ice_tau_solar=transpose(ice_tau_solar) if ice_omega_solar_id ne -1 then $ ncdf_varget,cdf_id,ice_omega_solar_id,ice_omega_solar $ else ice_omega_solar=dummy2d ice_omega_solar=transpose(ice_omega_solar) if ice_g_solar_id ne -1 then $ ncdf_varget,cdf_id,ice_g_solar_id,ice_g_solar $ else ice_g_solar=dummy2d ice_g_solar=transpose(ice_g_solar) if ice_tau_IR_id ne -1 then $ ncdf_varget,cdf_id,ice_tau_IR_id,ice_tau_IR $ else ice_tau_IR=dummy2d ice_tau_IR=transpose(ice_tau_IR) if ice_omega_IR_id ne -1 then $ ncdf_varget,cdf_id,ice_omega_IR_id,ice_omega_IR $ else ice_omega_IR=dummy2d ice_omega_IR=transpose(ice_omega_IR) if ice_g_IR_id ne -1 then $ ncdf_varget,cdf_id,ice_g_IR_id,ice_g_IR $ else ice_g_IR=dummy2d ice_g_IR=transpose(ice_g_IR) if ice_RadiationParameterization_source_solar_id ne -1 then $ ncdf_varget,cdf_id,ice_RadiationParameterization_source_solar_id,ice_RadiationParameterization_source_solar $ else ice_RadiationParameterization_source_solar=dummy2d ice_RadiationParameterization_source_solar=transpose(ice_RadiationParameterization_source_solar) if ice_RadiationParameterization_source_ir_id ne -1 then $ ncdf_varget,cdf_id,ice_RadiationParameterization_source_ir_id,ice_RadiationParameterization_source_ir $ else ice_RadiationParameterization_source_ir=dummy2d ice_RadiationParameterization_source_ir=transpose(ice_RadiationParameterization_source_ir) if lwc_best_estimate_id ne -1 then $ ncdf_varget,cdf_id,lwc_best_estimate_id,lwc_best_estimate $ else lwc_best_estimate=dummy2d lwc_best_estimate=transpose(lwc_best_estimate) if re_liquid_best_estimate_id ne -1 then $ ncdf_varget,cdf_id,re_liquid_best_estimate_id,re_liquid_best_estimate $ else re_liquid_best_estimate=dummy2d re_liquid_best_estimate=transpose(re_liquid_best_estimate) if liquid_micro_error_estimate_id ne -1 then $ ncdf_varget,cdf_id,liquid_micro_error_estimate_id,liquid_micro_error_estimate $ else liquid_micro_error_estimate=dummy2d liquid_micro_error_estimate=transpose(liquid_micro_error_estimate) if lwc_source_id ne -1 then $ ncdf_varget,cdf_id,lwc_source_id,lwc_source $ else lwc_source=dummy2d lwc_source=transpose(lwc_source) if liquid_effective_radius_source_id ne -1 then $ ncdf_varget,cdf_id,liquid_effective_radius_source_id,liquid_effective_radius_source $ else liquid_effective_radius_source=dummy2d liquid_effective_radius_source=transpose(liquid_effective_radius_source) if liquid_tau_solar_id ne -1 then $ ncdf_varget,cdf_id,liquid_tau_solar_id,liquid_tau_solar $ else liquid_tau_solar=dummy2d liquid_tau_solar=transpose(liquid_tau_solar) if liquid_omega_solar_id ne -1 then $ ncdf_varget,cdf_id,liquid_omega_solar_id,liquid_omega_solar $ else liquid_omega_solar=dummy2d liquid_omega_solar=transpose(liquid_omega_solar) if liquid_g_solar_id ne -1 then $ ncdf_varget,cdf_id,liquid_g_solar_id,liquid_g_solar $ else liquid_g_solar=dummy2d liquid_g_solar=transpose(liquid_g_solar) if liquid_tau_IR_id ne -1 then $ ncdf_varget,cdf_id,liquid_tau_IR_id,liquid_tau_IR $ else liquid_tau_IR=dummy2d liquid_tau_IR=transpose(liquid_tau_IR) if liquid_omega_IR_id ne -1 then $ ncdf_varget,cdf_id,liquid_omega_IR_id,liquid_omega_IR $ else liquid_omega_IR=dummy2d liquid_omega_IR=transpose(liquid_omega_IR) if liquid_g_IR_id ne -1 then $ ncdf_varget,cdf_id,liquid_g_IR_id,liquid_g_IR $ else liquid_g_IR=dummy2d liquid_g_IR=transpose(liquid_g_IR) if liquid_RadiationParameterization_source_solar_id ne -1 then $ ncdf_varget,cdf_id,liquid_RadiationParameterization_source_solar_id,liquid_RadiationParameterization_source_solar $ else liquid_RadiationParameterization_source_solar=dummy2d liquid_RadiationParameterization_source_solar=transpose(liquid_RadiationParameterization_source_solar) if liquid_RadiationParameterization_source_ir_id ne -1 then $ ncdf_varget,cdf_id,liquid_RadiationParameterization_source_ir_id,liquid_RadiationParameterization_source_ir $ else liquid_RadiationParameterization_source_ir=dummy2d liquid_RadiationParameterization_source_ir=transpose(liquid_RadiationParameterization_source_ir) if ice_fraction_CCM3_id ne -1 then $ ncdf_varget,cdf_id,ice_fraction_CCM3_id,ice_fraction_CCM3 $ else ice_fraction_CCM3=dummy2d ice_fraction_CCM3=transpose(ice_fraction_CCM3) if solar_flux_down_clear_id ne -1 then $ ncdf_varget,cdf_id,solar_flux_down_clear_id,solar_flux_down_clear $ else solar_flux_down_clear=dummy2d solar_flux_down_clear=transpose(solar_flux_down_clear) if solar_flux_up_clear_id ne -1 then $ ncdf_varget,cdf_id,solar_flux_up_clear_id,solar_flux_up_clear $ else solar_flux_up_clear=dummy2d solar_flux_up_clear=transpose(solar_flux_up_clear) if IR_flux_up_clear_id ne -1 then $ ncdf_varget,cdf_id,IR_flux_up_clear_id,IR_flux_up_clear $ else IR_flux_up_clear=dummy2d IR_flux_up_clear=transpose(IR_flux_up_clear) if IR_flux_down_clear_id ne -1 then $ ncdf_varget,cdf_id,IR_flux_down_clear_id,IR_flux_down_clear $ else IR_flux_down_clear=dummy2d IR_flux_down_clear=transpose(IR_flux_down_clear) if downwelling_solar_diffuse_clear_id ne -1 then $ ncdf_varget,cdf_id,downwelling_solar_diffuse_clear_id,downwelling_solar_diffuse_clear $ else downwelling_solar_diffuse_clear=dummy2d downwelling_solar_diffuse_clear=transpose(downwelling_solar_diffuse_clear) if downwelling_solar_direct_clear_id ne -1 then $ ncdf_varget,cdf_id,downwelling_solar_direct_clear_id,downwelling_solar_direct_clear $ else downwelling_solar_direct_clear=dummy2d downwelling_solar_direct_clear=transpose(downwelling_solar_direct_clear) if TOA_solar_up_clear_id ne -1 then $ ncdf_varget,cdf_id,TOA_solar_up_clear_id,TOA_solar_up_clear $ else TOA_solar_up_clear=dummy1d TOA_solar_up_clear=transpose(TOA_solar_up_clear) if TOA_solar_down_clear_id ne -1 then $ ncdf_varget,cdf_id,TOA_solar_down_clear_id,TOA_solar_down_clear $ else TOA_solar_down_clear=dummy1d TOA_solar_down_clear=transpose(TOA_solar_down_clear) if TOA_IR_up_clear_id ne -1 then $ ncdf_varget,cdf_id,TOA_IR_up_clear_id,TOA_IR_up_clear $ else TOA_IR_up_clear=dummy1d TOA_IR_up_clear=transpose(TOA_IR_up_clear) if solar_flux_down_id ne -1 then $ ncdf_varget,cdf_id,solar_flux_down_id,solar_flux_down $ else solar_flux_down=dummy2d solar_flux_down=transpose(solar_flux_down) if solar_flux_up_id ne -1 then $ ncdf_varget,cdf_id,solar_flux_up_id,solar_flux_up $ else solar_flux_up=dummy2d solar_flux_up=transpose(solar_flux_up) if IR_flux_up_id ne -1 then $ ncdf_varget,cdf_id,IR_flux_up_id,IR_flux_up $ else IR_flux_up=dummy2d IR_flux_up=transpose(IR_flux_up) if IR_flux_down_id ne -1 then $ ncdf_varget,cdf_id,IR_flux_down_id,IR_flux_down $ else IR_flux_up=dummy2d IR_flux_down=transpose(IR_flux_down) if downwelling_solar_diffuse_id ne -1 then $ ncdf_varget,cdf_id,downwelling_solar_diffuse_id,downwelling_solar_diffuse $ else downwelling_solar_diffuse=dummy2d downwelling_solar_diffuse=transpose(downwelling_solar_diffuse) if downwelling_solar_direct_id ne -1 then $ ncdf_varget,cdf_id,downwelling_solar_direct_id,downwelling_solar_direct $ else downwelling_solar_direct=dummy2d downwelling_solar_direct=transpose(downwelling_solar_direct) if TOA_solar_up_id ne -1 then $ ncdf_varget,cdf_id,TOA_solar_up_id,TOA_solar_up $ else TOA_solar_up=dummy2d TOA_solar_up=transpose(TOA_solar_up) if TOA_solar_down_id ne -1 then $ ncdf_varget,cdf_id,TOA_solar_down_id,TOA_solar_down $ else TOA_solar_down=dummy2d TOA_solar_down=transpose(TOA_solar_down) if TOA_IR_up_id ne -1 then $ ncdf_varget,cdf_id,TOA_IR_up_id,TOA_IR_up $ else TOA_IR_up=dummy2d TOA_IR_up=transpose(TOA_IR_up) ; ; Close the netcdf file ; ncdf_close,cdf_id ; ; If it is the first file ; if j eq 0 then begin com_btto=long(base_time)+long(time_offset) com_IWC_best_estimate=IWC_best_estimate com_re_ice_best_estimate=re_ice_best_estimate com_ice_micro_error_estimate=ice_micro_error_estimate com_iwc_source=iwc_source com_ice_effective_radius_source=ice_effective_radius_source com_ice_tau_solar=ice_tau_solar com_ice_omega_solar=ice_omega_solar com_ice_g_solar=ice_g_solar com_ice_tau_IR=ice_tau_IR com_ice_omega_IR=ice_omega_IR com_ice_g_IR=ice_g_IR com_ice_RadiationParameterization_source_solar=ice_RadiationParameterization_source_solar com_ice_RadiationParameterization_source_ir=ice_RadiationParameterization_source_ir com_lwc_best_estimate=lwc_best_estimate com_re_liquid_best_estimate=re_liquid_best_estimate com_liquid_micro_error_estimate=liquid_micro_error_estimate com_lwc_source=lwc_source com_liquid_effective_radius_source=liquid_effective_radius_source com_liquid_tau_solar=liquid_tau_solar com_liquid_omega_solar=liquid_omega_solar com_liquid_g_solar=liquid_g_solar com_liquid_tau_IR=liquid_tau_IR com_liquid_omega_IR=liquid_omega_IR com_liquid_g_IR=liquid_g_IR com_liquid_RadiationParameterization_source_solar=liquid_RadiationParameterization_source_solar com_liquid_RadiationParameterization_source_ir=liquid_RadiationParameterization_source_ir com_ice_fraction_CCM3=ice_fraction_CCM3 com_solar_flux_down_clear=solar_flux_down_clear com_solar_flux_up_clear=solar_flux_up_clear com_IR_flux_up_clear=IR_flux_up_clear com_IR_flux_down_clear=IR_flux_down_clear com_downwelling_solar_diffuse_clear=downwelling_solar_diffuse_clear com_downwelling_solar_direct_clear=downwelling_solar_direct_clear com_TOA_solar_up_clear=TOA_solar_up_clear com_TOA_solar_down_clear=TOA_solar_down_clear com_TOA_IR_up_clear=TOA_IR_up_clear com_solar_flux_down=solar_flux_down com_solar_flux_up=solar_flux_up com_IR_flux_up=IR_flux_up com_IR_flux_down=IR_flux_down com_downwelling_solar_diffuse=downwelling_solar_diffuse com_downwelling_solar_direct=downwelling_solar_direct com_TOA_solar_up=TOA_solar_up com_TOA_solar_down=TOA_solar_down com_TOA_IR_up=TOA_IR_up ; ; It if is not the first file ; endif else begin com_btto=[long(com_btto),long(base_time)+long(time_offset)] com_IWC_best_estimate=[com_IWC_best_estimate,IWC_best_estimate] com_re_ice_best_estimate=[com_re_ice_best_estimate,re_ice_best_estimate] com_ice_micro_error_estimate=[com_ice_micro_error_estimate,ice_micro_error_estimate] com_iwc_source=[com_iwc_source,iwc_source] com_ice_effective_radius_source=[com_ice_effective_radius_source,ice_effective_radius_source] com_ice_tau_solar=[com_ice_tau_solar,ice_tau_solar] com_ice_omega_solar=[com_ice_omega_solar,ice_omega_solar] com_ice_g_solar=[com_ice_g_solar,ice_g_solar] com_ice_tau_IR=[com_ice_tau_IR,ice_tau_IR] com_ice_omega_IR=[com_ice_omega_IR,ice_omega_IR] com_ice_g_IR=[com_ice_g_IR,ice_g_IR] com_ice_RadiationParameterization_source_solar= $ [com_ice_RadiationParameterization_source_solar,ice_RadiationParameterization_source_solar] com_ice_RadiationParameterization_source_ir= $ [com_ice_RadiationParameterization_source_ir,ice_RadiationParameterization_source_ir] com_lwc_best_estimate=[com_lwc_best_estimate,lwc_best_estimate] com_re_liquid_best_estimate=[com_re_liquid_best_estimate,re_liquid_best_estimate] com_liquid_micro_error_estimate=[com_liquid_micro_error_estimate,liquid_micro_error_estimate] com_lwc_source=[com_lwc_source,lwc_source] com_liquid_effective_radius_source=[com_liquid_effective_radius_source,liquid_effective_radius_source] com_liquid_tau_solar=[com_liquid_tau_solar,liquid_tau_solar] com_liquid_omega_solar=[com_liquid_omega_solar,liquid_omega_solar] com_liquid_g_solar=[com_liquid_g_solar,liquid_g_solar] com_liquid_tau_IR=[com_liquid_tau_IR,liquid_tau_IR] com_liquid_omega_IR=[com_liquid_omega_IR,liquid_omega_IR] com_liquid_g_IR=[com_liquid_g_IR,liquid_g_IR] com_liquid_RadiationParameterization_source_solar= $ [com_liquid_RadiationParameterization_source_solar,liquid_RadiationParameterization_source_solar] com_liquid_RadiationParameterization_source_ir= $ [com_liquid_RadiationParameterization_source_ir,liquid_RadiationParameterization_source_ir] com_ice_fraction_CCM3=[com_ice_fraction_CCM3,ice_fraction_CCM3] com_solar_flux_down_clear=[com_solar_flux_down_clear,solar_flux_down_clear] com_solar_flux_up_clear=[com_solar_flux_up_clear,solar_flux_up_clear] com_IR_flux_up_clear=[com_IR_flux_up_clear,IR_flux_up_clear] com_IR_flux_down_clear=[com_IR_flux_down_clear,IR_flux_down_clear] com_downwelling_solar_diffuse_clear=[com_downwelling_solar_diffuse_clear,downwelling_solar_diffuse_clear] com_downwelling_solar_direct_clear=[com_downwelling_solar_direct_clear,downwelling_solar_direct_clear] com_TOA_solar_up_clear=[com_TOA_solar_up_clear,TOA_solar_up_clear] com_TOA_solar_down_clear=[com_TOA_solar_down_clear,TOA_solar_down_clear] com_TOA_IR_up_clear=[com_TOA_IR_up_clear,TOA_IR_up_clear] com_solar_flux_down=[com_solar_flux_down,solar_flux_down] com_solar_flux_up=[com_solar_flux_up,solar_flux_up] com_IR_flux_up=[com_IR_flux_up,IR_flux_up] com_IR_flux_down=[com_IR_flux_down,IR_flux_down] com_downwelling_solar_diffuse=[com_downwelling_solar_diffuse,downwelling_solar_diffuse] com_downwelling_solar_direct=[com_downwelling_solar_direct,downwelling_solar_direct] com_TOA_solar_up=[com_TOA_solar_up,TOA_solar_up] com_TOA_solar_down=[com_TOA_solar_down,TOA_solar_down] com_TOA_IR_up=[com_TOA_IR_up,TOA_IR_up] endelse endif ;end of not a gif file endfor ;end of loop through matching files ; ; Create the new arrays ; base_time=long(com_btto[0]) ;base time time_offset=long(com_btto)-long(base_time) IWC_best_estimate=com_IWC_best_estimate re_ice_best_estimate=com_re_ice_best_estimate ice_micro_error_estimate=com_ice_micro_error_estimate iwc_source=com_iwc_source ice_effective_radius_source=com_ice_effective_radius_source ice_tau_solar=com_ice_tau_solar ice_omega_solar=com_ice_omega_solar ice_g_solar=com_ice_g_solar ice_tau_IR=com_ice_tau_IR ice_omega_IR=com_ice_omega_IR ice_g_IR=com_ice_g_IR ice_RadiationParameterization_source_solar=com_ice_RadiationParameterization_source_solar ice_RadiationParameterization_source_ir=com_ice_RadiationParameterization_source_ir lwc_best_estimate=com_lwc_best_estimate re_liquid_best_estimate=com_re_liquid_best_estimate liquid_micro_error_estimate=com_liquid_micro_error_estimate lwc_source=com_lwc_source liquid_effective_radius_source=com_liquid_effective_radius_source liquid_tau_solar=com_liquid_tau_solar liquid_omega_solar=com_liquid_omega_solar liquid_g_solar=com_liquid_g_solar liquid_tau_IR=com_liquid_tau_IR liquid_omega_IR=com_liquid_omega_IR liquid_g_IR=com_liquid_g_IR liquid_RadiationParameterization_source_solar=com_liquid_RadiationParameterization_source_solar liquid_RadiationParameterization_source_ir=com_liquid_RadiationParameterization_source_ir ice_fraction_CCM3=com_ice_fraction_CCM3 solar_flux_down_clear=com_solar_flux_down_clear solar_flux_up_clear=com_solar_flux_up_clear IR_flux_up_clear=com_IR_flux_up_clear IR_flux_down_clear=com_IR_flux_down_clear downwelling_solar_diffuse_clear=com_downwelling_solar_diffuse_clear downwelling_solar_direct_clear=com_downwelling_solar_direct_clear TOA_solar_up_clear=com_TOA_solar_up_clear TOA_solar_down_clear=com_TOA_solar_down_clear TOA_IR_up_clear=com_TOA_IR_up_clear solar_flux_down=com_solar_flux_down solar_flux_up=com_solar_flux_up IR_flux_up=com_IR_flux_up IR_flux_down=com_IR_flux_down downwelling_solar_diffuse=com_downwelling_solar_diffuse downwelling_solar_direct=com_downwelling_solar_direct TOA_solar_up=com_TOA_solar_up TOA_solar_down=com_TOA_solar_down TOA_IR_up=com_TOA_IR_up ; ; Do some conversions ; ;*** IWC_best_estimate IWC_best_estimate=10.0^(float(IWC_best_estimate)/1000.0) result_no_cloud=where(IWC_best_estimate gt 1e-9 and $ IWC_best_estimate lt 9.99e-9,count_no_cloud) result_missing=where(IWC_best_estimate gt 1e-8 and $ IWC_best_estimate lt 9.99e-8,count_missing) if count_no_cloud gt 0 then IWC_best_estimate[result_no_cloud]=-9999 if count_missing gt 0 then IWC_best_estimate[result_missing]=-8888 ;*** lwc_best_estimate lwc_best_estimate=10.0^(float(lwc_best_estimate)/1000.0) result_no_cloud=where(lwc_best_estimate gt 1e-9 and $ lwc_best_estimate lt 9.99e-9,count_no_cloud) result_missing=where(lwc_best_estimate gt 1e-8 and $ lwc_best_estimate lt 9.99e-8,count_missing) if count_no_cloud gt 0 then lwc_best_estimate[result_no_cloud]=-9999 if count_missing gt 0 then lwc_best_estimate[result_missing]=-8888 ;*** re_liquid_best_estimate re_liquid_best_estimate=10.0^(float(re_liquid_best_estimate)/1000.0) result_no_cloud=where(re_liquid_best_estimate gt 1e-9 and $ re_liquid_best_estimate lt 9.99e-9,count_no_cloud) result_missing=where(re_liquid_best_estimate gt 1e-8 and $ re_liquid_best_estimate lt 9.99e-8,count_missing) if count_no_cloud gt 0 then re_liquid_best_estimate[result_no_cloud]=-9999 if count_missing gt 0 then re_liquid_best_estimate[result_missing]=-8888 ;*** re_liquid_best_estimate re_ice_best_estimate=10.0^(float(re_ice_best_estimate)/1000.0) result_no_cloud=where(re_ice_best_estimate gt 1e-9 and $ re_ice_best_estimate lt 9.99e-9,count_no_cloud) result_missing=where(re_ice_best_estimate gt 1e-8 and $ re_ice_best_estimate lt 9.99e-8,count_missing) if count_no_cloud gt 0 then re_ice_best_estimate[result_no_cloud]=-9999 if count_missing gt 0 then re_ice_best_estimate[result_missing]=-8888 ;*** ice_tau_solar ice_tau_solar=10.^(float(ice_tau_solar)/1000.) result_no_cloud=where(ice_tau_solar gt 1e-9 and $ ice_tau_solar lt 9.99e-9,count_no_cloud) result_missing=where(ice_tau_solar gt 1e-8 and $ ice_tau_solar lt 9.99e-8,count_missing) if count_no_cloud gt 0 then ice_tau_solar[result_no_cloud]=-9999 if count_missing gt 0 then ice_tau_solar[result_missing]=-8888 ;*** Convert liquid_tau_solar liquid_tau_solar=10.^(float(liquid_tau_solar)/1000.) result_no_cloud=where(liquid_tau_solar gt 1e-9 and $ liquid_tau_solar lt 9.99e-9,count_no_cloud) result_missing=where(liquid_tau_solar gt 1e-8 and $ liquid_tau_solar lt 9.99e-8,count_missing) if count_no_cloud gt 0 then liquid_tau_solar[result_no_cloud]=-9999 if count_missing gt 0 then liquid_tau_solar[result_missing]=-8888 ; ; Get parameters of new arrays ; numtimes=n_elements(time_offset) numheights=n_elements(height) return end