da_reflectivity.inc

1 subroutine da_reflectivity(rf, rho, qrn)
2 
3    !-----------------------------------------------------------------------
4    ! Purpose: TBD
5    !-----------------------------------------------------------------------
6 
7    implicit none
8 
9    real, intent(in)  :: rho, qrn
10    real, intent(out) :: rf
11 
12    real, parameter   :: leh1=43.1
13    real, parameter   :: leh2=17.5
14 
15    rf = leh1 + leh2 * ALOG10(rho * qrn * 1.E+3)
16 
17 end subroutine da_reflectivity
18 
19