init_slt.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 init_slt (SLT, &
                     HPROGRAM  &! Program calling unit
       )  

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

!PASSED VARIABLES
!
TYPE(slt_t), INTENT(INOUT) :: SLT
!
 CHARACTER(LEN=6), INTENT(IN) :: HPROGRAM              !Passing unit
!
!LOCAL VARIABLES
 CHARACTER(LEN=4)    :: CRGUNITS              ! type of log-normal geometric mean radius
INTEGER             :: JMODE                 ! Counter for sea salt modes
INTEGER             :: JMODE_IDX             ! Index for sea salt modes
REAL(KIND=JPRB) :: ZHOOK_HANDLE

!get output listing unit
IF (lhook) CALL dr_hook('INIT_SLT',0,zhook_handle)
!
!Allocate memory for the real values which will be used by the model
ALLOCATE(slt%XEMISRADIUS_SLT(nsltmde))
ALLOCATE(slt%XEMISSIG_SLT   (nsltmde))
!
!Get initial size distributions. This is cut and pasted
!from dead routine dstpsd.F90
!Check for different source parameterizations
IF(cemisparam_slt.eq."Vig01") THEN
  crgunits   = 'NUMB'
  xemisradius_ini_slt(:) = (/ 0.2, 2.0, 12.  /)  ! [um]  Number median radius She84 p. 75 Table 1
  xemissig_ini_slt(:) = (/ 1.9, 2.0, 3.00 /)  ! [frc] Geometric standard deviation She84 p. 75 Table 1
ELSE  ! use default of Schultz et al, 2004
  crgunits   = 'MASS'
  xemisradius_ini_slt(:) = 0.5*(/0.28, 2.25, 15.32/) ! [um] Mass median radius
  xemissig_ini_slt(:) =     (/1.59, 2.00, 2.00 /) ! [frc] Geometric standard deviation
ENDIF
!
DO jmode=1,nsltmde
  !
  jmode_idx = jorder_slt(jmode)
  !
  slt%XEMISSIG_SLT   (jmode) = xemissig_ini_slt(jmode_idx)
  slt%XEMISRADIUS_SLT(jmode) = xemisradius_ini_slt(jmode_idx)
  !
  IF (crgunits=="MASS") &
    slt%XEMISRADIUS_SLT(jmode) = slt%XEMISRADIUS_SLT(jmode) * exp(-3.d0 * (log(slt%XEMISSIG_SLT(jmode)))**2)    
  !
ENDDO
!
IF (lhook) CALL dr_hook('INIT_SLT',1,zhook_handle)
!
END SUBROUTINE init_slt