source: trunk/SOURCES/prop-thermiques_mod.f90 @ 10

!> \file prop-thermiques_mod.f90
!! Module pour definir les proprietes thermiques de la galce
!<

!> \namespace prop_thermiques_ice
!! Module pour definir les proprietes thermiques de la galce
!! \author ...
!! \date ...
!! @note ce block était avant dans icetemp 
!! @note Used module
!! @note   - use module3D_phy
!<



module prop_thermiques_ice

  ! défini les proprietes thermiques de la glace
  ! bloc qui était avant dans icetemp
  ! valeurs standard

  use module3D_phy

  real,dimension(nx,ny,nz) :: CPP     !< CP specific heat capacity (J/(m-3)/K)=ro cp
  real,dimension(nx,ny,nz) :: CTT     !< CT thermal conductivity (J/m/K/a)


contains
  !> SUBROUTINE: thermal_prop
  !! Subroutine pour calculer les proprietes thermiques 
  !>
  subroutine thermal_prop


    do k=1,nz  
       do j=1,ny
          do i=1,nx

             !  temperature de fusion:TPMP (en deg C)
             !  capacite calorifique : CP en J/Kg/K puis multiplie par ro en J/m3/K
             !  conductivite : CT  en J/m/K/an
             !   Attention pour la conductivite c'est la formule de Yin-Chao Yen


             tpmp(i,j,k)=-0.00087*(k-1)*de*h(i,j)

             cpp(i,j,k)=(2115.3+7.79293*t(i,j,k)) 
             cpp(i,j,k)=cpp(i,j,k)*ro              ! en J/m3/K
             ctt(i,j,k)=3.1014e8*exp(-0.0057*(t(i,j,k)+273.15)) !en J/m/K/an

          end do
       end do
    end do

  end subroutine thermal_prop
end module prop_thermiques_ice