ch_conversion_factor.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 ch_conversion_factor (CHN, &
                                       hconversion,prhoa)
!     #######################################
!
!!****  *CH_CONVERSION_FACTOR
!!
!!    PURPOSE
!!    -------
!     Determines the correct conversion factor to produce emissions in
!     Molec/m2/s
!
!!**  METHOD
!!    ------
!!    
!!    
!!    AUTHOR
!!    ------
!!      S.QUEGUINER 
!!
!!    MODIFICATIONS
!!    -------------
!!      Original        11/2011
!!      A. Alias        07/2013 add MODI_ABOR1_SFX
!!      M. Leriche      04/2014 correct conversion factor
!!-----------------------------------------------------------------------------
!
!*       0.    DECLARATIONS
!
!
USE modd_ch_snap_n, ONLY : ch_emis_snap_t
!
USE modd_csts,       ONLY : xavogadro, xmd
USE modi_abor1_sfx
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
!
IMPLICIT NONE
!
!*       0.1   declarations of arguments
!
!
TYPE(ch_emis_snap_t), INTENT(INOUT) :: chn
!
 CHARACTER(LEN=3),  INTENT(IN)  :: hconversion ! Unit conversion code
REAL, DIMENSION(:),INTENT(IN)  :: prhoa       ! air density
!
!*       0.2   declarations of local variables
!
REAL(KIND=JPRB) :: zhook_handle
!-------------------------------------------------------------------------------
IF (lhook) CALL dr_hook('CH_CONVERSION_FACTOR',0,zhook_handle)
!
! determine the conversion factor
 chn%XCONVERSION(:) = 1.
SELECT CASE (hconversion)
  CASE ('MIX') ! flux given ppp*m/s,  conversion to molec/m2/s
    chn%XCONVERSION(:) = xavogadro * prhoa(:) / xmd
  CASE ('CON') ! flux given in molecules/cm2/s, conversion to molec/m2/s 
    chn%XCONVERSION(:) =  1e4
  CASE ('MOL') ! flux given in microMol/m2/day, conversion to molec/m2/s  
    chn%XCONVERSION(:) = 1e-6 * xavogadro / 86400.
  CASE default
    CALL abor1_sfx('CH_BUILDEMISSN: UNKNOWN CONVERSION FACTOR')
END SELECT
!
IF (lhook) CALL dr_hook('CH_CONVERSION_FACTOR',1,zhook_handle)
!-------------------------------------------------------------------------------
!
END SUBROUTINE ch_conversion_factor