Ignore:
Timestamp:
Jun 29, 2020 8:49:58 AM (7 months ago)
Author:
suehring
Message:

mesoscale nesting: omit explicit pressure forcing via geostrophic wind components; chemistry: enable profile output of vertical fluxes; urban-surface: bugfix in initialization in case of cyclic_fill

File:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/chem_modules.f90

    r4577 r4581  
    2626! -----------------
    2727! $Id$
     28! Enable output of vertical fluxes of chemical species.
     29!
     30! 4577 2020-06-25 09:53:58Z raasch
    2831! further re-formatting concerning Fortran parameter variables
    2932!
     
    163166                                            !< 4: non-cyclic
    164167
    165     INTEGER(iwp) ::  cs_pr_count                           = 0      !< counter for chemical species
    166                                                                     !< profiles
     168    INTEGER(iwp) ::  cs_pr_count_fl_res                    = 0      !< counter for vertical flux profiles of chemical species (resolved-scale)
     169    INTEGER(iwp) ::  cs_pr_count_fl_sgs                    = 0      !< counter for vertical flux profiles of chemical species (SGS)
     170    INTEGER(iwp) ::  cs_pr_count_sp                        = 0      !< counter for chemical species profiles
    167171    INTEGER(iwp) ::  cs_vertical_gradient_level_ind(99,10) = -9999  !< grid index values of
    168172                                                                    !< cs_vertical_gradient_level
     
    174178    INTEGER(iwp) ::  main_street_id                        = 0      !< namelist parameter: lower bound of main street IDs
    175179                                                                    !< (OpenStreetMaps) for PARAMETERIZED mode
    176     INTEGER(iwp) ::  max_pr_cs                             = 0      !<
     180    INTEGER(iwp) ::  max_pr_cs                             = 0      !< number of chemistry profiles in output
    177181    INTEGER(iwp) ::  max_street_id                         = 0      !< namelist parameter: upper bound of main street IDs
    178182                                                                    !< (OpenStreetMaps) for PARAMETERIZED mode
     
    182186                                                                    !< (OpenStreetMaps) for PARAMETERIZED mode
    183187
    184     INTEGER(iwp), DIMENSION(99) ::  cs_pr_index  = 0      !< index for chemical species profiles
    185     INTEGER(iwp), DIMENSION(:)  ::  match_spec_nox(1:2)   !< results of matching the input and
    186                                                           !< model's NOx
    187     INTEGER(iwp), DIMENSION(:)  ::  match_spec_pm(1:3)    !< results of matching the input and
    188                                                           !< model's PMs
    189     INTEGER(iwp), DIMENSION(:)  ::  match_spec_sox(1:2)   !< results of matching the input and
    190                                                           !< model's SOx!
     188    INTEGER(iwp), DIMENSION(99) ::  cs_pr_index_sp      = 0 !< index for chemical species mean profiles
     189    INTEGER(iwp), DIMENSION(99) ::  cs_pr_index_fl_res  = 0 !< index for chemical species sgs-flux profiles
     190    INTEGER(iwp), DIMENSION(99) ::  cs_pr_index_fl_sgs  = 0 !< index for chemical species resolved-scale flux profiles
     191    INTEGER(iwp), DIMENSION(99) ::  hom_index_fl_res    = 0 !< index of the resolved-scale flux profile with respect to the hom array
     192    INTEGER(iwp), DIMENSION(99) ::  hom_index_fl_sgs    = 0 !< index of the SGS flux profile with respect to the hom array
     193    INTEGER(iwp), DIMENSION(99) ::  hom_index_spec      = 0 !< index of the profile with respect to the hom array
     194    INTEGER(iwp), DIMENSION(:)  ::  match_spec_nox(1:2)     !< results of matching the input and model's NOx
     195    INTEGER(iwp), DIMENSION(:)  ::  match_spec_pm(1:3)      !< results of matching the input and model's PMs
     196    INTEGER(iwp), DIMENSION(:)  ::  match_spec_sox(1:2)     !< results of matching the input and model's SOx!
    191197
    192198    INTEGER(iwp), DIMENSION(:), ALLOCATABLE ::  match_spec_input      !< index of input chem species
     
    262268    REAL(wp), DIMENSION(:), ALLOCATABLE ::  css          !< scaling parameter for chem species
    263269
    264     REAL(wp), DIMENSION(:,:,:,:), ALLOCATABLE ::  emis_distribution  !< emissions final values
    265                                                                      !< (main module output) ???
     270    REAL(wp), DIMENSION(:,:,:),   ALLOCATABLE, TARGET ::  sums_ws_l  !< subdomain sum of vertical chemistry flux w'ch'
     271
     272    REAL(wp), DIMENSION(:,:,:,:), ALLOCATABLE ::  emis_distribution  !< emissions final values (main module output) ???
    266273
    267274    REAL(wp), DIMENSION(:,:,:,:), ALLOCATABLE, TARGET ::  cs_1  !< pointer for swapping of
Note: See TracChangeset for help on using the changeset viewer.