source: branches/GRISLIv3/SOURCES/initial2-0.4.f90 @ 483

!> \file initial2-0.4.f90
!!Initialisation du modele
!<

!> SUBROUTINE INITIAL2
!! Cette routine permet d'initialiser le modele. 
!! \author ...
!! \date 19 Novembre 1999
!! @note Elle appelle toutes les routines d'initialisation des differents packages.
!! @note  C'est cette routine qui est appellee dans le main
!! @note  Used modules:
!! @note    - use module3D_phy
!! @note    - use param_phy_mod
!<


subroutine INITIAL2()

  USE module3D_phy, only: cde,abl,bdot,bmelt,epsxx,epsyy,epsxy,hwater,calv,hdot,hdotwater, &
                    slope,taub,ubx,uby,uzk,uxbar,uybar,ibase,tpmp,mk,mk0,h,front, &
                    bsoc,flot,sealevel_2D,flotmx,flotmy,gzmx, &
                    gzmy,flgzmx,flgzmy,sdx,sdy,ux,uy,uzr,t,s
  use geography, only: nx,ny,nz,nzm,dx
  USE param_phy_mod,only: ro,row
  implicit none
  real,dimension(NZ) ::cord_vert 

  integer :: i,j,k

  !     INITIALISATION DES TABLEAUX

  cord_vert(1)=0.
  cord_vert(NZ)=1.
  do K=1,NZ
     if ((K.ne.1).and.(K.ne.NZ)) cord_vert(K)=(K-1.)/(NZ-1.)
     CDE(K)=1.-cord_vert(K)
  end do

  ABL(:,:)=0. 
  BDOT(:,:)=0. 
  BMELT(:,:)=0.
  EPSXX(:,:)=0.
  EPSYY(:,:)=0.
  EPSXY(:,:)=0.
  HWATER(:,:)=0.
  CALV(:,:)=0.    
  HDOT(:,:)=0. 
  HDOTWATER(:,:)=0.
  SLOPE(:,:)=0. 
  TAUB(:,:)=0. 
  UBX(:,:)=0.
  UBY(:,:)=0.
  UZK(:,:)=0.
  UXBAR(:,:)=0.
  UYBAR(:,:)=0.
  IBASE(:,:)=1
  TPMP(:,:,1)=0
  TPMP(:,:,NZ)=-0.00087*H(:,:)

  MK(:,:)=MK0(:,:)


  !-------Initialisation des fronts.
  FRONT(:,:)  =0 

  !-------Initialisation des flot... devrait etre mise dans input_topo.
  ! en tout cas on le fait pour mismip

  do J=1,NY
     do I=1,NX
        if ((BSOC(I,J)+H(I,J)*RO/ROW -SEALEVEL_2D(I,J)).LT.0.) then
           FLOT(I,J)=.TRUE.
        else
           FLOT(I,J)=.FALSE.
        endif
     enddo
  enddo


  FLOTMX(:,:)=FLOT(:,:)
  FLOTMY(:,:)=FLOT(:,:)
  GZMX(:,:)=.FALSE.
  GZMY(:,:)=.FALSE.
  FLGZMX(:,:)=.FALSE.
  FLGZMY(:,:)=.FALSE.


  do J=2,NY
     do I=2,NX
        SDX(I,J)=(S(I,J)-S(I-1,J))/DX
        SDY(I,J)=(S(I,J)-S(I,J-1))/DX
     end do
  end do



  do K=1,NZ
     do I=1,NX  
        do J=1,NY
           UX(I,J,K)=0.
           UY(I,J,K)=0.
           UZR(I,J,K)=0.
           T(I,J,K)=0.
           !CT(I,J,K)=6.62E7
           !CP(I,J,K)=2009.
        end do
     end do
  end do

  do K=2,NZ-1
     do I=1,NX  
        do J=1,NY
           TPMP(I,J,K)=-0.00087*(K-1)/(NZ-1)*H(I,J)
        end do
     end do
  end do

  do K=NZ+1,NZ+NZM
     do I=1,NX  
        do J=1,NY
           T(I,J,K)=0. 
        end do
     end do
  end do
end subroutine INITIAL2