surfex_alloc.F90

Go to the documentation of this file.

!SFX_LIC Copyright 1994-2014 CNRS, Meteo-France and Universite Paul Sabatier
!SFX_LIC This is part of the SURFEX software governed by the CeCILL-C licence
!SFX_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt  
!SFX_LIC for details. version 1.
SUBROUTINE surfex_alloc(YDSURFEX)
!
USE modd_teb_par, ONLY : nteb_patch_max
USE modd_data_cover_par, ONLY : nvegtype_ecosg, ntilesfc
!
USE modd_surfex_n, ONLY : surfex_t
!
USE modd_ch_emis_field_n, ONLY : ch_emis_field_init
USE modd_ch_snap_n, ONLY : ch_emis_snap_init
USE modd_ch_surf_n, ONLY : ch_surf_init
USE modd_data_cover_n, ONLY : data_cover_init
USE modd_dummy_surf_fields_n, ONLY : dummy_surf_fields_init
USE modd_emis_gr_field_n, ONLY : emis_gr_field_init
USE modd_sfx_grid_n, ONLY : grid_init, grid_np_init
USE modd_canopy_n, ONLY : canopy_init
USE modd_grid_conf_proj_n, ONLY : grid_conf_proj_init
USE modd_surf_atm_grid_n, ONLY : surf_atm_grid_init
USE modd_surf_atm_n, ONLY : surf_atm_init
USE modd_sso_n, ONLY : sso_init, sso_np_init
USE modd_sv_n, ONLY : sv_init
!
USE modd_data_tsz0_n, ONLY : data_tsz0_init
!
USE modd_ideal_n, ONLY : ideal_init
!
USE modd_dst_n, ONLY : dst_np_init
USE modd_slt_n, ONLY : slt_init
!
USE modd_diag_n, ONLY : diag_init, diag_np_init, diag_options_init
!
USE modd_diag_evap_isba_n, ONLY : diag_evap_isba_init, diag_evap_isba_np_init
USE modd_diag_misc_isba_n, ONLY : diag_misc_isba_init, diag_misc_isba_np_init
USE modd_diag_ocean_n, ONLY : diag_ocean_init
USE modd_diag_misc_seaice_n, ONLY : diag_misc_seaice_init
USE modd_diag_misc_flake_n, ONLY : diag_misc_flake_init
USE modd_diag_misc_teb_options_n, ONLY : diag_misc_teb_options_init
USE modd_diag_utci_teb_n, ONLY : diag_utci_teb_init
USE modd_diag_misc_teb_n, ONLY : diag_misc_teb_np_init
!
USE modd_data_bem_n, ONLY : data_bem_init
USE modd_bem_option_n, ONLY : bem_options_init
USE modd_bld_description_n, ONLY : bld_desc_init
USE modd_ch_teb_n, ONLY : ch_teb_init
USE modd_data_teb_n, ONLY : data_teb_init
USE modd_teb_irrig_n, ONLY : teb_irrig_init
USE modd_teb_option_n, ONLY : teb_options_init
USE modd_teb_panel_n, ONLY : teb_panel_init
USE modd_bem_n, ONLY : bem_np_init
USE modd_teb_n, ONLY : teb_np_init
!
USE modd_ch_flake_n, ONLY : ch_flake_init
USE modd_flake_n, ONLY : flake_init
!
USE modd_ch_watflux_n, ONLY : ch_watflux_init
USE modd_watflux_n, ONLY : watflux_init
!
USE modd_ch_isba_n, ONLY : ch_isba_init, ch_isba_np_init
USE modd_agri_n, ONLY : agri_np_init
USE modd_data_isba_n, ONLY : data_isba_init
USE modd_gr_biog_n, ONLY : gr_biog_init, gr_biog_np_init
USE modd_isba_options_n, ONLY : isba_options_init
USE modd_isba_n, ONLY : isba_s_init, isba_k_init, isba_p_init, &
                        isba_nk_init, isba_np_init, isba_npe_init
!
USE modd_ch_seaflux_n, ONLY : ch_seaflux_init
USE modd_data_seaflux_n, ONLY : data_seaflux_init
USE modd_seaflux_n, ONLY : seaflux_init
USE modd_ocean_n, ONLY : ocean_init
USE modd_ocean_rel_n, ONLY : ocean_rel_init
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE
!
type(surfex_t), INTENT (INOUT) :: ydsurfex
REAL(KIND=JPRB) :: ZHOOK_HANDLE
IF (lhook) CALL dr_hook("SURFEX_ALLOC",0,zhook_handle)
  !
  CALL diag_options_init(ydsurfex%FM%DFO)
  CALL diag_init(ydsurfex%FM%DF)
  CALL diag_init(ydsurfex%FM%DFC)
  CALL diag_misc_flake_init(ydsurfex%FM%DMF)
  ! 
  CALL grid_init(ydsurfex%FM%G)  
  CALL canopy_init(ydsurfex%FM%SB)  
  CALL ch_flake_init(ydsurfex%FM%CHF)
  CALL flake_init(ydsurfex%FM%F)
  !  
  !
  CALL diag_options_init(ydsurfex%WM%DWO)
  CALL diag_init(ydsurfex%WM%DW)
  CALL diag_init(ydsurfex%WM%DWC)
  !  
  CALL grid_init(ydsurfex%WM%G)
  CALL canopy_init(ydsurfex%WM%SB)  
  CALL ch_watflux_init(ydsurfex%WM%CHW)
  CALL watflux_init(ydsurfex%WM%W)
  !
  !
  CALL diag_options_init(ydsurfex%SM%SD%O) 
  CALL diag_init(ydsurfex%SM%SD%D)
  CALL diag_init(ydsurfex%SM%SD%DC)
  CALL diag_init(ydsurfex%SM%SD%DI)
  CALL diag_init(ydsurfex%SM%SD%DIC)
  CALL diag_ocean_init(ydsurfex%SM%SD%GO)  
  CALL diag_misc_seaice_init(ydsurfex%SM%SD%DMI)
  !  
  CALL data_seaflux_init(ydsurfex%SM%DTS)
  CALL grid_init(ydsurfex%SM%G)
  CALL canopy_init(ydsurfex%SM%SB)  
  CALL ch_seaflux_init(ydsurfex%SM%CHS)
  CALL seaflux_init(ydsurfex%SM%S)
  CALL seaflux_init(ydsurfex%SM%S)  
  CALL ocean_init(ydsurfex%SM%O)
  CALL ocean_rel_init(ydsurfex%SM%OR)
  !
  !
  CALL diag_options_init(ydsurfex%IM%ID%O)  
  CALL diag_init(ydsurfex%IM%ID%D)
  CALL diag_init(ydsurfex%IM%ID%DC)
  CALL diag_np_init(ydsurfex%IM%ID%ND,nvegtype_ecosg)
  CALL diag_np_init(ydsurfex%IM%ID%NDC,nvegtype_ecosg)
  CALL diag_evap_isba_init(ydsurfex%IM%ID%DE)
  CALL diag_evap_isba_init(ydsurfex%IM%ID%DEC)
  CALL diag_evap_isba_np_init(ydsurfex%IM%ID%NDE,nvegtype_ecosg)
  CALL diag_evap_isba_np_init(ydsurfex%IM%ID%NDEC,nvegtype_ecosg) 
  CALL diag_misc_isba_init(ydsurfex%IM%ID%DM)
  CALL diag_misc_isba_np_init(ydsurfex%IM%ID%NDM,nvegtype_ecosg)
  !
  CALL data_isba_init(ydsurfex%IM%DTV)
  CALL canopy_init(ydsurfex%IM%SB)
  CALL isba_options_init(ydsurfex%IM%O)
  CALL isba_s_init(ydsurfex%IM%S)  
  CALL ch_isba_init(ydsurfex%IM%CHI)
  CALL ch_isba_np_init(ydsurfex%IM%NCHI,nvegtype_ecosg)
  CALL gr_biog_init(ydsurfex%IM%GB)
  CALL gr_biog_np_init(ydsurfex%IM%NGB,nvegtype_ecosg)
  CALL sso_init(ydsurfex%IM%ISS)
  CALL sso_np_init(ydsurfex%IM%NISS,nvegtype_ecosg)
  CALL grid_init(ydsurfex%IM%G)
  CALL grid_np_init(ydsurfex%IM%NG,nvegtype_ecosg)
  CALL isba_k_init(ydsurfex%IM%K)
  CALL isba_nk_init(ydsurfex%IM%NK,nvegtype_ecosg)
  CALL isba_np_init(ydsurfex%IM%NP,nvegtype_ecosg)
  CALL isba_npe_init(ydsurfex%IM%NPE,nvegtype_ecosg)  
  CALL agri_np_init(ydsurfex%IM%NAG,nvegtype_ecosg)
  !
  !
  CALL diag_np_init(ydsurfex%GDM%VD%ND,nteb_patch_max)  
  CALL diag_evap_isba_np_init(ydsurfex%GDM%VD%NDE,nteb_patch_max)
  CALL diag_evap_isba_np_init(ydsurfex%GDM%VD%NDEC,nteb_patch_max)
  CALL diag_misc_isba_np_init(ydsurfex%GDM%VD%NDM,nteb_patch_max)  
  !  
  CALL data_isba_init(ydsurfex%GDM%DTV)
  CALL isba_options_init(ydsurfex%GDM%O)
  CALL isba_s_init(ydsurfex%GDM%S)  
  CALL gr_biog_init(ydsurfex%GDM%GB)
  CALL isba_k_init(ydsurfex%GDM%K)
  CALL isba_p_init(ydsurfex%GDM%P)
  CALL isba_npe_init(ydsurfex%GDM%NPE,nteb_patch_max)
  !
  !
  CALL diag_np_init(ydsurfex%GRM%VD%ND,nteb_patch_max)
  CALL diag_evap_isba_np_init(ydsurfex%GRM%VD%NDE,nteb_patch_max)
  CALL diag_evap_isba_np_init(ydsurfex%GRM%VD%NDEC,nteb_patch_max)
  CALL diag_misc_isba_np_init(ydsurfex%GRM%VD%NDM,nteb_patch_max)
  !
  CALL data_isba_init(ydsurfex%GRM%DTV)
  CALL isba_options_init(ydsurfex%GRM%O)
  CALL isba_s_init(ydsurfex%GRM%S)  
  CALL gr_biog_init(ydsurfex%GRM%GB)
  CALL isba_k_init(ydsurfex%GRM%K)
  CALL isba_p_init(ydsurfex%GRM%P)
  CALL isba_npe_init(ydsurfex%GRM%NPE,nteb_patch_max)
  !
  !
  CALL diag_options_init(ydsurfex%TM%TD%O)
  CALL diag_init(ydsurfex%TM%TD%D)
  CALL diag_misc_teb_options_init(ydsurfex%TM%TD%MTO)
  CALL diag_misc_teb_np_init(ydsurfex%TM%TD%NDMT,nteb_patch_max)   
  CALL diag_misc_teb_np_init(ydsurfex%TM%TD%NDMTC,nteb_patch_max)
  CALL diag_utci_teb_init(ydsurfex%TM%TD%DUT)
  ! 
  CALL data_teb_init(ydsurfex%TM%DTT)
  CALL teb_options_init(ydsurfex%TM%TOP)
  CALL canopy_init(ydsurfex%TM%SB)
  CALL grid_init(ydsurfex%TM%G)  
  CALL ch_teb_init(ydsurfex%TM%CHT)
  CALL teb_panel_init(ydsurfex%TM%TPN)
  CALL teb_irrig_init(ydsurfex%TM%TIR)    
  CALL teb_np_init(ydsurfex%TM%NT,nteb_patch_max)  
  !
  CALL data_bem_init(ydsurfex%TM%DTB)    
  CALL bem_options_init(ydsurfex%TM%BOP)  
  CALL bld_desc_init(ydsurfex%TM%BDD)
  CALL bem_np_init(ydsurfex%TM%NB,nteb_patch_max)  
  !
  !
  CALL data_cover_init(ydsurfex%DTCO)
  CALL data_tsz0_init(ydsurfex%DTZ)
  CALL dummy_surf_fields_init(ydsurfex%DUU)
  !
  CALL grid_conf_proj_init(ydsurfex%GCP)
  CALL surf_atm_grid_init(ydsurfex%UG)
  CALL surf_atm_init(ydsurfex%U)
  CALL diag_options_init(ydsurfex%DUO) 
  CALL diag_init(ydsurfex%DU) 
  CALL diag_init(ydsurfex%DUC)
  CALL diag_np_init(ydsurfex%DUP,ntilesfc) 
  CALL diag_np_init(ydsurfex%DUPC,ntilesfc)
  CALL sso_init(ydsurfex%USS)
  CALL canopy_init(ydsurfex%SB)
  !
  CALL diag_init(ydsurfex%DL)
  CALL diag_init(ydsurfex%DLC)
  CALL ideal_init(ydsurfex%L)
  !
  CALL sv_init(ydsurfex%SV)
  CALL ch_surf_init(ydsurfex%CHU)  
  CALL ch_emis_field_init(ydsurfex%CHE)
  CALL ch_emis_snap_init(ydsurfex%CHN)
  CALL emis_gr_field_init(ydsurfex%EGF)  
  CALL dst_np_init(ydsurfex%NDST,nvegtype_ecosg)
  CALL slt_init(ydsurfex%SLT)
  !
IF (lhook) CALL dr_hook("SURFEX_ALLOC",1,zhook_handle)
!
END SUBROUTINE surfex_alloc