@@ -5204,4 +5204,27 @@ elemental function analytic_sol_del12_whole_atmos_f560_r8(fc,Nmf,z0,H10,Hc0,beta
5204
5204
del12 = trm1* trm2- trm3+ trm4
5205
5205
end function analytic_sol_del12_whole_atmos_f560_r8
5206
5206
5207
+ ! Formula: 5.55, page: 106
5208
+ elemental function analytic_sol_del1_whole_atmos_f555_r4 (fc ,Nmf ,delnA ,z0 ,H10 ,Hc0 ,beta ,d ) result(del1)
5209
+ if defined(__INTEL_COMPILER) && ! defined(__GNUC__)
5210
+ ! dir$ optimize:3
5211
+ ! dir$ attributes code_align : 32 :: analytic_sol_del1_whole_atmos_f555_r4
5212
+ ! dir$ attributes forceinline :: analytic_sol_del1_whole_atmos_f555_r4
5213
+ #endif
5214
+
5215
+ real (kind= sp), intent (in ) :: fc
5216
+ real (kind= sp), intent (in ) :: Nmf
5217
+ real (kind= sp), intent (in ) :: delnA
5218
+ real (kind= sp), intent (in ) :: z0
5219
+ real (kind= sp), intent (in ) :: H10
5220
+ real (kind= sp), intent (in ) :: Hc0
5221
+ real (kind= sp), intent (in ) :: beta
5222
+ real (kind= sp), intent (in ) :: d
5223
+ real (kind= sp) :: del1
5224
+ real (kind= sp), automatic :: del11, del12
5225
+ del11 = analytic_sol_del11_whole_atmos_f559_r4 (delnA,z0,H10,Hc0,beta)
5226
+ del12 = analytic_sol_del12_whole_atmos_f560_r4 (fc,Nmf,z0,H10,Hc0,beta,d)
5227
+ del1 = del11+ del12
5228
+ end function analytic_sol_del1_whole_atmos_f555_r4
5229
+
5207
5230
end module emw_refraction
0 commit comments