7 & LBi, UBi, LBj, UBj, &
8 & IminS, ImaxS, JminS, JmaxS, &
9 & Zscale, Z, tl_Z, tl_swdk)
51 integer,
intent(in) :: ng, tile
52 integer,
intent(in) :: LBi, UBi, LBj, UBj
53 integer,
intent(in) :: IminS, ImaxS, JminS, JmaxS
55 real(r8),
intent(in) :: Zscale
57 real(r8),
intent(in) :: Z(IminS:ImaxS,JminS:JmaxS)
58 real(r8),
intent(in) :: tl_Z(IminS:ImaxS,JminS:JmaxS)
60 real(r8),
intent(out) :: tl_swdk(IminS:ImaxS,JminS:JmaxS)
64 integer :: Jindex, i, j
66 real(r8) :: cff1, cff2
67 real(r8) :: tl_cff1, tl_cff2
69 real(r8),
dimension(IminS:ImaxS) :: fac1, fac2, fac3
71# include "set_bounds.h"
80 jindex=int(
mixing(ng)%Jwtype(i,j))
87 cff1=exp(z(i,j)*fac1(i))
88 tl_cff1=fac1(i)*tl_z(i,j)*cff1
89 cff2=exp(z(i,j)*fac2(i))
90 tl_cff2=fac2(i)*tl_z(i,j)*cff2
94 tl_swdk(i,j)=tl_cff1*fac3(i)+ &
95 & tl_cff2*(1.0_r8-fac3(i))
subroutine tl_lmd_swfrac_tile(ng, tile, lbi, ubi, lbj, ubj, imins, imaxs, jmins, jmaxs, zscale, z, tl_z, tl_swdk)