Added refraction_angle_whole_atmos_vw5cm3m_f543_r4

Author: bgin

Radiolocation/GMS_emw_refraction.f90

@@ -4998,6 +4998,39 @@ elemental function analytic_sol_L2_whole_atmosphere_f552_r8(fc,Nmf,beta,R0,delnA
                                                                     z0,H0,H1,H2,Hc)
             L2    = L21+L22 
        end function analytic_sol_L2_whole_atmosphere_f552_r8
+
+       ! Formula: 5.43, page: 104
+       elemental function refraction_angle_whole_atmos_vw5cm3m_f543_r4(na,nc,fc,Nmf,beta,R0,delnA,        &
+                                                              z0,H0,H1,H2,Hc) result(alpha)
+if defined(__INTEL_COMPILER) && !defined(__GNUC__)           
+            !dir$ optimize:3
+            !dir$ attributes code_align : 32 :: analytic_sol_L2_whole_atmosphere_f552_r4
+            !dir$ attributes forceinline :: analytic_sol_L2_whole_atmosphere_f552_r4
+#endif 
+            real(kind=sp),    intent(in) :: na 
+            real(kind=sp),    intent(in) :: nc 
+            real(kind=sp),    intent(in) :: fc
+            real(kind=sp),    intent(in) :: Nmf 
+            real(kind=sp),    intent(in) :: beta 
+            real(kind=sp),    intent(in) :: R0 
+            real(kind=sp),    intent(in) :: delnA 
+            real(kind=sp),    intent(in) :: z0 
+            real(kind=sp),    intent(in) :: H0 
+            real(kind=sp),    intent(in) :: H1 
+            real(kind=sp),    intent(in) :: H2 
+            real(kind=sp),    intent(in) :: Hc                     
+            real(kind=sp)                :: alpha 
+            real(kind=sp),   automatic   :: H10, ctgz0 
+            real(kind=sp),   automatic   :: L1, L2
+            real(kind=sp),   automatic   :: L  
+            H10   = H1-H0 
+            ctgz0 = 1.0_sp/tan(z0) 
+            L1    = analytic_sol_L1_whole_atmosphere_f548_r4(beta,R0,delnA,z0,H10)
+            L2    = analytic_sol_L2_whole_atmosphere_f552_r4(fc,Nmf,beta,R0,delnA,        &
+                                                                    z0,H0,H1,H2,Hc)
+            L     = L1+L2 
+            alpha = -log(na/nc)*ctgz0+L                                                        
+       end function refraction_angle_whole_atmos_vw5cm3m_f543_r4