subroutine xTi_profile use modglobal, only : ib, ie, jb, je, kb, ke use modfields, only : thl0, thlm, thlprof integer j, k ! set ghost cell ! do j = jb - 1, je + 1 ! do k = kb, ke + 1 ! thl0(ib - 1, j, k) = 2*thlprof(k) - thl0(ib, j, k) ! thlm(ib - 1, j, k) = 2*thlprof(k) - thlm(ib, j, k) ! end do ! end do do j = jb - 1, je + 1 do k = kb, ke + 1 thl0(ib - 1, j, k) = thlprof(k) thlm(ib - 1, j, k) = thlprof(k) end do end do ! set first internal cell as well do j = jb - 1, je + 1 do k = kb, ke thl0(ib, j, k) = thlprof(k) thlm(ib, j, k) = thlprof(k) end do end do end subroutine xTi_profile