7 & LBi, UBi, LBj, UBj, &
8 & IminS, ImaxS, JminS, JmaxS, &
49 integer,
intent(in) :: ng, tile
50 integer,
intent(in) :: LBi, UBi, LBj, UBj
51 integer,
intent(in) :: IminS, ImaxS, JminS, JmaxS
53 real(r8),
intent(in) :: Zscale
55 real(r8),
intent(in) :: Z(IminS:ImaxS,JminS:JmaxS)
56 real(r8),
intent(out) :: swdk(IminS:ImaxS,JminS:JmaxS)
60 integer :: Jindex, i, j
62 real(r8),
dimension(IminS:ImaxS) :: fac1, fac2, fac3
64# include "set_bounds.h"
73 jindex=int(
mixing(ng)%Jwtype(i,j))
80 swdk(i,j)=exp(z(i,j)*fac1(i))*fac3(i)+ &
81 & exp(z(i,j)*fac2(i))*(1.0_r8-fac3(i))
subroutine lmd_swfrac_tile(ng, tile, lbi, ubi, lbj, ubj, imins, imaxs, jmins, jmaxs, zscale, z, swdk)