Changeset 3876 for palm/trunk/SOURCE/chem_modules.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/chem_modules.f90 (modified) (2 diffs)

palm/trunk/SOURCE/chem_modules.f90

@@ -94 +94 @@
     CHARACTER (LEN=80) ::  mode_emis      = 'PARAMETERIZED'     !< namelist parameter: mode of chemistry emissions - DEFAULT, EXPERT, PARAMETERIZED
     CHARACTER (LEN=80) ::  time_fac_type  = 'MDH'               !< namelist parameter: type of time treatment in the mode_emis DEFAULT - HOUR, MDH
-    
+    CHARACTER (LEN=10) ::  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
+
     CHARACTER (LEN=11), DIMENSION(99) ::  cs_name             = 'novalue'  !< namelist parameter: names of species with given fluxes (see csflux)
     CHARACTER (LEN=11), DIMENSION(99) ::  cs_profile_name     = 'novalue'  !< namelist parameter: tbc...???
@@ -183 +186 @@
     REAL, PARAMETER ::  xm_NO3   = xm_N + xm_O * 3         !< NO3      molecular weight (kg/mol)
     REAL, PARAMETER ::  xm_SO4   = xm_S + xm_O * 4         !< SO4      molecular weight (kg/mol)
+!
+!-  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 ::  chem_species
+    TYPE(photols_def), ALLOCATABLE, DIMENSION(:), TARGET ::  phot_frequen
 
     SAVE