modd_dstn.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.
MODULE modd_dst_n

!Purpose: 
!Declare variables and constants necessary to do the dust calculations
!Here are only the variables which depend on the grid!

!Author: Alf Grini <alf.grini@cnrm.meteo.fr>

!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE

TYPE dst_t
  !
  INTEGER, DIMENSION(:), POINTER     :: nvt_dst           !MASK: dust vegetation number to vegetation number
  INTEGER, DIMENSION(:), POINTER   :: nsize_patch_dst   !Number of points for a patch and a vegetation class
  INTEGER, DIMENSION(:,:), POINTER :: nr_patch_dst      !Mask from patch-points to dust-points
  REAL,DIMENSION(:), POINTER         :: z0_erod_dst       !Roughness length momentum over erodible dust emitter sfc
  CHARACTER(LEN=6), DIMENSION(:), POINTER  :: csv_dst     !Name of scalar variables 
  REAL, DIMENSION(:,:),POINTER    :: xsfdst             !Dust variables to be send to output
  REAL, DIMENSION(:,:),POINTER    :: xsfdstm            !Dust variables to be send to output
  REAL,DIMENSION(:), POINTER   :: xemisradius_dst         !Number median radius for each source mode
  REAL,DIMENSION(:), POINTER   :: xemissig_dst            !sigma for each source mode
  REAL,DIMENSION(:), POINTER   :: xmss_frc_src            !Mass fraction of each source mode  
  !
END TYPE dst_t
!
TYPE dst_np_t
!
TYPE(dst_t), DIMENSION(:), POINTER :: al=>null()
!
END TYPE dst_np_t
!
CONTAINS
!
SUBROUTINE dst_init(YDST)
TYPE(dst_t), INTENT(INOUT) :: YDST
REAL(KIND=JPRB) :: ZHOOK_HANDLE
IF (lhook) CALL dr_hook("MODD_DST_N:DST_INIT",0,zhook_handle)
  NULLIFY(ydst%NVT_DST)
  NULLIFY(ydst%NSIZE_PATCH_DST)
  NULLIFY(ydst%NR_PATCH_DST)
  NULLIFY(ydst%Z0_EROD_DST)
  NULLIFY(ydst%CSV_DST)
  NULLIFY(ydst%XSFDST)
  NULLIFY(ydst%XSFDSTM)
  NULLIFY(ydst%XEMISRADIUS_DST)
  NULLIFY(ydst%XEMISSIG_DST)
  NULLIFY(ydst%XMSS_FRC_SRC)
IF (lhook) CALL dr_hook("MODD_DST_N:DST_INIT",1,zhook_handle)
END SUBROUTINE dst_init
!
SUBROUTINE dst_np_init(YDST_NP,KPATCH)
TYPE(dst_np_t), INTENT(INOUT) :: YDST_NP
INTEGER, INTENT(IN) :: KPATCH
INTEGER :: JP
REAL(KIND=JPRB) :: ZHOOK_HANDLE
!
IF (lhook) CALL dr_hook("MODD_DST_N:DST_NP_INIT",0,zhook_handle)
!
IF (ASSOCIATED(ydst_np%AL)) THEN
  DO jp = 1,kpatch
    CALL dst_init(ydst_np%AL(jp))
  ENDDO
  DEALLOCATE(ydst_np%AL)
ELSE
  ALLOCATE(ydst_np%AL(kpatch))
  DO jp = 1,kpatch
    CALL dst_init(ydst_np%AL(jp))
  ENDDO
ENDIF
!
IF (lhook) CALL dr_hook("MODD_DST_N:DST_NP_INIT",1,zhook_handle)
!
END SUBROUTINE dst_np_init
!
END MODULE modd_dst_n