forked from drdaphos/LIDAR
-
Notifications
You must be signed in to change notification settings - Fork 2
/
Copy pathcspot_calc.pro
executable file
·29 lines (20 loc) · 1.32 KB
/
cspot_calc.pro
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
function cspot_calc, spot_diam
;The manufacturer’s calibration of PSAP assumes a spot area A(PSAP) of 17.83mm2,
;equivalent to a circular spot diameter of 4.765mm. This is the diameter of the hole in the filter holder,
;which is smaller than the exposed filter area owing to the use of O-rings in the filter holder to provide
;an airtight seal. A(PSAP) is used internally by the instrument in its derivation of absorption.
;B1999 observed some variation in spot size among instruments and recommend that the actual spot
;area A(TRUE) for each instrument should be measured and used to correct the spot area of the manufacturer’s
;reference instrument A(REF). Cspot = Atrue/Aref
;However, Ogren (2010) notes that the spot size correction should in fact compare the actual spot area A(TRUE)
;for the instrument with the spot area assumed by PSAP in its internal calculations, A(PSAP), rather than
;comparing it to the reference instrument. Thus, the spot size correction factor should really
; be Cspot=Atrue/Apsap
; Typical spot diameters are around about 5.3mm for the FAAM instrument.
if ((spot_diam gt 6.0) or (spot_diam lt 4.5)) then $
print, 'Spot_diam of ', strtrim(spot_diam), 'is significantly different to usual for this instrument.'
Atrue=!dpi*(spot_diam/2.0)^2.0
Apsap = 17.83 ;mm^2
Cspot = Atrue/Apsap
return, Cspot
end