@@ -2504,6 +2504,28 @@ elemental function refraction_angle_atmos_L2_med_f442_r4(deln0,fc,Nmf,H1,H2,a,z0
2504
2504
L22 = analytic_sol_L22_med_ionosphere_f444_r4 (deln0,fc,Nmf,H1,H2,a,z0)
2505
2505
L2 = L21+ L22
2506
2506
end function refraction_angle_atmos_L2_med_f442_r4
2507
+
2508
+ elemental function refraction_angle_atmos_L2_med_f442_r8 (deln0 ,fc ,Nmf ,H1 ,H2 ,a ,z0 ) result(L2)
2509
+ if defined(__INTEL_COMPILER) && ! defined(__GNUC__)
2510
+ ! dir$ optimize:3
2511
+ ! dir$ attributes code_align : 32 :: refraction_angle_atmos_L2_med_f442_r8
2512
+ ! dir$ attributes forceinline :: refraction_angle_atmos_L2_med_f442_r8
2513
+ #endif
2514
+
2515
+ real (kind= dp), intent (in ) :: deln0
2516
+ real (kind= dp), intent (in ) :: fc
2517
+ real (kind= dp), intent (in ) :: Nmf
2518
+ real (kind= dp), intent (in ) :: H1
2519
+ real (kind= dp), intent (in ) :: H2
2520
+ real (kind= dp), intent (in ) :: a
2521
+ real (kind= dp), intent (in ) :: z0
2522
+ real (kind= dp) :: L2
2523
+ real (kind= dp), automatic :: L21, L22
2524
+ L21 = analytic_sol_L21_med_ionosphere_f443_r8 (fc,Nmf,H1,H2,a,z0)
2525
+ L22 = analytic_sol_L22_med_ionosphere_f444_r8 (deln0,fc,Nmf,H1,H2,a,z0)
2526
+ L2 = L21+ L22
2527
+ end function refraction_angle_atmos_L2_med_f442_r8
2528
+
2507
2529
2508
2530
2509
2531
end module emw_refraction
0 commit comments