subroutine interpolate_TIPEM_TIGCM(ngrid,nslope,nsoil_PEM,nsoil_GCM,TI_PEM,TI_GCM)

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!
!!! Purpose: Transfer the thermal inertia from the PEM vertical  grid to the GCM vertical grid
!!! 
!!! 
!!! Author: LL
!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!


      implicit none


!======================================================================
!  arguments
!  ---------
!  inputs:
      integer,intent(in) :: ngrid                                ! # of horizontal grid points
      integer,intent(in) :: nslope                               ! # of subslope wihtin the mesh
      integer,intent(in) :: nsoil_PEM                            ! # of soil layers in the PEM
      integer,intent(in) :: nsoil_GCM                            ! # of soil layers in the GCM
      real,intent(in) :: TI_PEM(ngrid,nsoil_PEM,nslope)          !  Thermal inertia in the PEM vertical grid [J/m^2/K/s^{1/2}]
      real,intent(inout) :: TI_GCM(ngrid,nsoil_GCM,nslope)       ! Thermal inertia in the GCM vertical grid [J/m^2/K/s^{1/2}]

!local variable
      integer :: ig,islope,iloop                                 ! loop variables
     do ig = 1,ngrid
       do islope = 1,nslope
         do iloop = 1,nsoil_GCM
           TI_GCM(ig,iloop,islope) = TI_PEM(ig,iloop,islope)
         enddo
       enddo
     enddo


  end subroutine interpolate_TIPEM_TIGCM