Changeset 3876 for palm/trunk/SOURCE/chemistry_model_mod.f90

Timestamp:

Apr 8, 2019 6:41:49 PM (5 years ago)

Author:

knoop

Message:

Moved "photolysis_scheme", "chem_species" and "phot_frequen" to chem_modules

File:

• palm/trunk/SOURCE/chemistry_model_mod.f90 (modified) (3 diffs)

palm/trunk/SOURCE/chemistry_model_mod.f90

@@ -318 +318 @@
 
     LOGICAL ::  nest_chemistry = .TRUE.  !< flag for nesting mode of chemical species, independent on parent or not
-!
-!-  Define chemical variables within chem_species
-    TYPE species_def
-       CHARACTER(LEN=15)                            ::  name         !< name of chemical species
-       CHARACTER(LEN=15)                            ::  unit         !< unit (ppm for gases, kg m^-3 for aerosol tracers)
-       REAL(kind=wp), POINTER, DIMENSION(:,:,:)     ::  conc         !< concentrations of trace gases
-       REAL(kind=wp), POINTER, DIMENSION(:,:,:)     ::  conc_av      !< averaged concentrations
-       REAL(kind=wp), POINTER, DIMENSION(:,:,:)     ::  conc_p       !< conc at prognostic time level
-       REAL(kind=wp), POINTER, DIMENSION(:,:,:)     ::  tconc_m      !< weighted tendency of conc for previous sub-timestep (Runge-Kutta)
-       REAL(kind=wp), ALLOCATABLE, DIMENSION(:,:)   ::  cssws_av     !< averaged fluxes of trace gases at surface
-       REAL(kind=wp), ALLOCATABLE, DIMENSION(:,:)   ::  flux_s_cs    !< 6th-order advective flux at south face of grid box of chemical species (='cs')
-       REAL(kind=wp), ALLOCATABLE, DIMENSION(:,:)   ::  diss_s_cs    !< artificial numerical dissipation flux at south face of grid box of chemical species
-       REAL(kind=wp), ALLOCATABLE, DIMENSION(:,:,:) ::  flux_l_cs    !< 6th-order advective flux at left face of grid box of chemical species (='cs')
-       REAL(kind=wp), ALLOCATABLE, DIMENSION(:,:,:) ::  diss_l_cs    !< artificial numerical dissipation flux at left face of grid box of chemical species
-       REAL(kind=wp), ALLOCATABLE, DIMENSION(:)     ::  conc_pr_init !< initial profile of chemical species
-    END TYPE species_def
-!
-!-- Define photolysis frequencies in phot_frequen
-    TYPE photols_def                                                            
-       CHARACTER(LEN=15)                            :: name          !< name of pgotolysis frequency
-       CHARACTER(LEN=15)                            :: unit          !< unit (1/s)
-       REAL(kind=wp), POINTER, DIMENSION(:,:,:)     :: freq          !< photolysis frequency
-    END TYPE photols_def
-
-
-    TYPE(species_def), ALLOCATABLE, DIMENSION(:), TARGET, PUBLIC ::  chem_species
-    TYPE(photols_def), ALLOCATABLE, DIMENSION(:), TARGET, PUBLIC ::  phot_frequen  
 
     REAL(kind=wp), ALLOCATABLE, DIMENSION(:,:,:,:), TARGET ::  spec_conc_1  !< pointer for swapping of timelevels for conc
@@ -374 +347 @@
 
     REAL(kind=wp), PUBLIC ::  cs_time_step = 0._wp
-    CHARACTER(10), PUBLIC ::  photolysis_scheme
-                              !< 'constant',
-                              !< 'simple' (Simple parameterisation from MCM, Saunders et al., 2003, Atmos. Chem. Phys., 3, 161-180
-                              !< 'fastj'  (Wild et al., 2000, J. Atmos. Chem., 37, 245-282) STILL NOT IMPLEMENTED
-
 
 !
@@ -445 +413 @@
     PUBLIC nspec               !< number of gas phase chemical species including constant compound (e.g. N2)
     PUBLIC nvar                !< number of variable gas phase chemical species (nvar <= nspec)
-    PUBLIC photols_def         !< type defining phot_frequen
-    PUBLIC species_def         !< type defining chem_species
     PUBLIC spc_names           !< names of gas phase chemical species (come from KPP) (come from KPP)
     PUBLIC spec_conc_2