Added refraction_angle_atmos_L2_med_f442_r4

Author: bgin

Radiolocation/GMS_emw_refraction.f90

@@ -2482,4 +2482,28 @@ elemental function analytic_sol_L22_med_ionosphere_f444_r8(deln0,fc,Nmf,H1,H2,a,
             L22     = trm1*trm2 
       end function analytic_sol_L22_med_ionosphere_f444_r8
 
+      ! refraction angle whole atmosphere (medium part).
+      ! formula: 4.42, page: 82
+      elemental function refraction_angle_atmos_L2_med_f442_r4(deln0,fc,Nmf,H1,H2,a,z0) result(L2)
+if defined(__INTEL_COMPILER) && !defined(__GNUC__)           
+            !dir$ optimize:3
+            !dir$ attributes code_align : 32 :: refraction_angle_atmos_L2_med_f442_r4
+            !dir$ attributes forceinline :: refraction_angle_atmos_L2_med_f442_r4
+#endif 
+
+            real(kind=sp),   intent(in) :: deln0
+            real(kind=sp),   intent(in) :: fc 
+            real(kind=sp),   intent(in) :: Nmf 
+            real(kind=sp),   intent(in) :: H1 
+            real(kind=sp),   intent(in) :: H2 
+            real(kind=sp),   intent(in) :: a 
+            real(kind=sp),   intent(in) :: z0 
+            real(kind=sp) :: L2
+            real(kind=sp),   automatic :: L21, L22 
+            L21  = analytic_sol_L21_med_ionosphere_f443_r4(fc,Nmf,H1,H2,a,z0)
+            L22  = analytic_sol_L22_med_ionosphere_f444_r4(deln0,fc,Nmf,H1,H2,a,z0)
+            L2   = L21+L22 
+      end function refraction_angle_atmos_L2_med_f442_r4
+     
+
 end module emw_refraction