;pro plot_ro_dimm_mass
;
;This script will go through and plot up the various CN2's and compute the 
;seeing based upon them

pro plot_ro_dimm_mass,hr1=hr1,date1=date1,min1=min1,hr2=hr2,date2=date2,$
                      min2=min2

CN2=return_mass_spread(min1,hr1,date1,min2,hr2,date2)
CN2=.1*CN2
;NOTE CN2's in table are off by factor of 10

ground_layer_CN2=CN2(0,*)
Free_atmosphere_CN2=total(CN2(1:5,*),1)

ro_ground=0.185*(4*!pi^2/((2*!pi/500e-9)^2*(ground_layer_CN2)))^(3./5)
ro_FA=0.185*(4*!pi^2/((2*!pi/500e-9)^2*(Free_atmosphere_CN2)))^(3./5)

rad_deg=(360*3600./(2*!pi))
ground_seeing=(.98*500e-9/ro_ground)*rad_deg
FA_seeing=(.98*500e-9/ro_FA)*rad_deg
stop

end