source: branches/GRISLIv3/SOURCES/initial-0.3.f90 @ 399

!> \file initial-0.3.f90
!!Initialisation du modele
!<

!> SUBROUTINE initial
!! 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 module_choix
!! @note    - use icetempmod
!! @note    - use diagno_mod 
!! @note    - use resolmeca_SIA_L1
!! @note    - use sorties_ncdf_grisli
!! @note    - use util_recovery
!<


subroutine initial

  use module3D_phy, only:time,H,hmx,hmy,T,tpmp,bmelt,geoplace,icompteur,num_rep_42
  use runparam, only: tbegin
  use module_choix, only:input_topo,init_forclim,input_clim,init_ablation,init_tracer,&
       input_rsl,init_eaubasale,init_bmelt,init_outshort,init_iso,init_deformation,initial_heino,&
       init_furst_schoof,init_sliding,init_spinup,init_dragging,init_icethick,initial_matrice,&
       init_calving
  use icetempmod, only:init_icetemp
  use diagno_mod, only:init_diagno
  use resolmeca_SIA_L1, only:init_resol_meca
  use flottab_mod, only:
  use sorties_ncdf_grisli, only:sortie_ncdf_cat,init_out_ncdf
  use util_recovery, only:init_recovery
  use bilan_flux_mod, only:init_bilan_flux

  !-------------------------------------------------------------------------------------
  ! nouvelles entrees modulaires  janvier 2006
  call initial_phy          ! La physique
  call initial2             ! Initialisation des tableaux
  !
  call  input_topo          ! La geographie

  !------------------------------------------------------
  ! initialisation du climat  (reference et forcage)
!cdc  call input_climat_ref    ! Le climat de reference
  call init_forclim        ! parametres du forcage
  call input_clim          ! lecture fichiers de forcage
  call init_ablation       ! parametres du calcul de l'ablation
  !-------------------------------------------------------------------------------------

  ! aurel neem -> pour initialisation traceurs
  call init_tracer  

  !initialisation du niveau marin local
  call input_rsl


  ! call  input_profile()     ! Les fichiers profils
  !

  if (icompteur.ne.0) then ! reprise  d un fichier cptr
     call read_recovery(icompteur)
     time = tbegin
     call sortie_ncdf_cat
!  else if (icompteur.eq.2) then ! reprise d'un fichier cptr mais sans la topo
!     call read_recovery_temp

!  else if (icompteur.eq.3) then ! reprise d'un fichier cptr mais sans la topo ni l'eau
!     call read_recovery_temp
!     hwater(:,:)=0.
  else
     call read_no_recovery
  endif

  call sortie_ncdf_cat

  ! calcul de Hmx et Hmy -> shift=-1, dim=1 -> H(i-1,j)

  H(:,:)   = max(H(:,:),0.)
  hmx(:,:) = 0.5*(H(:,:)+eoshift(H(:,:),shift=-1,boundary=0d0,dim=1))
  hmy(:,:) = 0.5*(H(:,:)+eoshift(H(:,:),shift=-1,boundary=0d0,dim=2))
  hmx(:,:) = max(hmx(:,:),0.)
  hmy(:,:) = max(hmy(:,:),0.)


  call init_eaubasale                  ! initialise le calcul des temperatures

  if (geoplace(1:5).ne.'mism3') then
     call init_icetemp(num_rep_42)        ! initialise le calcul de temperature
     call thermal_prop_icetemp            ! appel aux proprietes thermiques
     call init_bmelt                      ! Appel aux tableaux d'initialisation des ice shelves

  else ! initialisation dans le cas mismip
     T(:,:,:)    = 0.
     TPMP(:,:,:) = 0.
     bmelt(:,:)  = 0.
  end if

  !-------------------------------------------------------------------------------------

  !-------------------------------------------------------------------------------------
  call init_outshort        ! initialisation sorties temporelles
  call init_out_ncdf        ! initialisation des sorties netcdf
  call init_recovery        ! initialisation des sorties cptr ou nc

  !-------------------------------------------------------------------------------------
  !Appels a l'origine (2005) dans initial
  !
  call init_iso ! initialisation de l'isostasie cette routine est dans le module
  !               isostasie_mod et dans noisostasie_mod



  call init_deformation ! initialisation de la deformation. cette routine est 
  !                       dans le module deformation_mod

  call init_resol_meca  ! type d'association SIA L1

  call initial_heino   ! a mettre avant les init sliding et deformation

  call init_furst_schoof ! initialisation furst schoof gr line

  call init_diagno     ! initialisation de la resolution equation elliptique vitesses 

  call init_sliding    ! initialisation du glissement

  call init_spinup     ! initialisation du spinup

  call init_dragging   ! initialisation du frottement

  call init_icethick   ! interverti avec init_dragging

  call initial_matrice !   cette routine est dans le module eq_elliptique_mod

  call init_calving

  call init_bilan_flux
!
!-------------------------------------------------------------------------------------
! ecriture netcdf apres initialisation

!!$  if (iter_beta.eq.0) then
!!$
!!$     if (itracebug.eq.1)  call tracebug(' Avant appel routine icethick3')
!!$     call icethick3
!!$     debug_3D(:,:,88) = S(:,:)
!!$     if (itracebug.eq.1)  call tracebug(' Apres appel routine icethick3')
!!$  end if
  return

end subroutine initial