ch_init_names.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_init_names (KLUOUT,HSV, SV, OVARSIGI, OVARSIGJ)  
!!    ###########################################
!!
!!*** *CH_INIT_NAMES*
!!
!!    PURPOSE
!!    -------
!!      Read and filter all chemical species into the CSV array
!!     initialize NSV_CHSBEG and  NSV_CHSEND index for the begin and the ending chemical index
!!     
!!
!!    REFERENCE
!!    ---------
!!    
!!    AUTHOR
!!    ------
!!    P. Tulet    *Meteo France*
!!
!!    MODIFICATIONS
!!    -------------
!!    Original 16/10/01
!!    01/12/03    (D.Gazen) change emissions handling for surf. externalization
!!    01/06/05    (P.Tulet) add aerosols list
!!
!!    EXTERNAL
!!    --------
!!
!!    IMPLICIT ARGUMENTS
!!    ------------------
USE modd_sv_n, ONLY : sv_t
!
USE modd_chs_aerosol
!-------------------------------------------------------------------------------
!
!*       0.     DECLARATIONS
!               ------------
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE
!
!*      0.1    declarations of arguments
!

INTEGER,                         INTENT(IN)  :: KLUOUT ! output listing channel
 CHARACTER(LEN=*), DIMENSION(:),  INTENT(IN)  :: HSV    ! name of chemical species
                                                       ! with character # (gas chemistry )
                                                       ! and  character @ (aerosols)
TYPE(sv_t), INTENT(INOUT) :: SV
!
LOGICAL,                         INTENT(OUT) :: OVARSIGI, OVARSIGJ ! type of standard deviation
!
!*      0.2    declarations of local variables
INTEGER :: JSV  !! loop  NBEQ
 CHARACTER        :: YRC1
 CHARACTER(LEN=5) :: YRC2
REAL(KIND=JPRB) :: ZHOOK_HANDLE
!
!-------------------------------------------------------------------------------

IF (lhook) CALL dr_hook('CH_INIT_NAMES',0,zhook_handle)
sv%NBEQ = 0
sv%NAEREQ = 0
sv%NSV_CHSBEG = 0
sv%NSV_AERBEG = 0
sv%NSV_CHSEND = 0
sv%NSV_AEREND = 0
ovarsigi = .false.
ovarsigj = .false.
nsoa = 0


DO jsv=1, SIZE(hsv)

  sv%CSV(jsv) = hsv(jsv)
  yrc1= hsv(jsv)(1:1)
  yrc2 = hsv(jsv)(2:)

  IF (yrc1 == '#') THEN
    sv%CSV(jsv) = trim(yrc2)
    sv%NBEQ = sv%NBEQ + 1
    IF (sv%NBEQ == 1) sv%NSV_CHSBEG=jsv
  ELSE IF (yrc1 == '@') THEN
    sv%CSV(jsv) = trim(yrc2)
    sv%NAEREQ = sv%NAEREQ + 1
    IF (sv%NAEREQ == 1) sv%NSV_AERBEG=jsv
    IF (sv%CSV(jsv) == "M6I") ovarsigi = .true.
    IF (sv%CSV(jsv) == "M6J") ovarsigj = .true.
    IF (sv%CSV(jsv) == "SOA1I") nsoa = 10
  ENDIF

ENDDO

sv%NSV_CHSEND = sv%NSV_CHSBEG + sv%NBEQ -1
sv%NSV_AEREND = sv%NSV_AERBEG + sv%NAEREQ -1

IF (sv%NAEREQ .GT. 0) THEN
DO jsv=1, size(sv%CSV)
   IF (trim(sv%CSV(jsv)) == "M0I") jp_ch_m0i=jsv-sv%NSV_CHSEND
   IF (trim(sv%CSV(jsv)) == "M0J") jp_ch_m0j=jsv-sv%NSV_CHSEND
   IF (trim(sv%CSV(jsv)) == "M6I") jp_ch_m6i=jsv-sv%NSV_CHSEND
   IF (trim(sv%CSV(jsv)) == "M6J") jp_ch_m6j=jsv-sv%NSV_CHSEND
   IF (trim(sv%CSV(jsv)) == "H2OI") jp_ch_h2oi=jsv-sv%NSV_CHSEND
   IF (trim(sv%CSV(jsv)) == "H2OJ") jp_ch_h2oj=jsv-sv%NSV_CHSEND
   IF (trim(sv%CSV(jsv)) == "SO4I") jp_ch_so4i=jsv-sv%NSV_CHSEND
   IF (trim(sv%CSV(jsv)) == "SO4J") jp_ch_so4j=jsv-sv%NSV_CHSEND
   IF (trim(sv%CSV(jsv)) == "NO3I") jp_ch_no3i=jsv-sv%NSV_CHSEND
   IF (trim(sv%CSV(jsv)) == "NO3J") jp_ch_no3j=jsv-sv%NSV_CHSEND
   IF (trim(sv%CSV(jsv)) == "NH3I") jp_ch_nh3i=jsv-sv%NSV_CHSEND
   IF (trim(sv%CSV(jsv)) == "NH3J") jp_ch_nh3j=jsv-sv%NSV_CHSEND
   IF (trim(sv%CSV(jsv)) == "OCI") jp_ch_oci=jsv-sv%NSV_CHSEND
   IF (trim(sv%CSV(jsv)) == "OCJ") jp_ch_ocj=jsv-sv%NSV_CHSEND
   IF (trim(sv%CSV(jsv)) == "BCI") jp_ch_bci=jsv-sv%NSV_CHSEND
   IF (trim(sv%CSV(jsv)) == "BCJ") jp_ch_bcj=jsv-sv%NSV_CHSEND
   IF (trim(sv%CSV(jsv)) == "DSTI") jp_ch_dsti=jsv-sv%NSV_CHSEND
   IF (trim(sv%CSV(jsv)) == "DSTJ") jp_ch_dstj=jsv-sv%NSV_CHSEND
END DO

END IF
IF (lhook) CALL dr_hook('CH_INIT_NAMES',1,zhook_handle)

!
END SUBROUTINE ch_init_names