Changeset 4154


Ignore:
Timestamp:
Aug 13, 2019 1:35:59 PM (2 years ago)
Author:
suehring
Message:

Bugfix in chemistry decycling; Replace global arrays to setup chemical emissions by local ones

Location:
palm/trunk/SOURCE
Files:
2 edited

Legend:

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

    r4144 r4154  
    2727! -----------------
    2828! $Id$
     29! Replace global arrays for emissions by local ones.
     30!
     31! 4144 2019-08-06 09:11:47Z raasch
    2932! relational operators .EQ., .NE., etc. replaced by ==, /=, etc.
    3033!
     
    860863
    861864                IF  (  .NOT.  ALLOCATED( emis_distribution) )  THEN
    862                    ALLOCATE( emis_distribution(1,0:ny,0:nx,n_matched_vars) )
     865                   ALLOCATE( emis_distribution(1,nys:nyn,nxl:nxr,n_matched_vars) )
    863866                ENDIF
    864867
     
    868871
    869872                IF  (  .NOT.  ALLOCATED( emis_distribution) )  THEN
    870                    ALLOCATE( emis_distribution(1,0:ny,0:nx,n_matched_vars) )
     873                   ALLOCATE( emis_distribution(1,nys:nyn,nxl:nxr,n_matched_vars) )
    871874                ENDIF
    872875
     
    876879
    877880                IF  (  .NOT.  ALLOCATED( emis_distribution) )  THEN
    878                    ALLOCATE( emis_distribution(nzb:nzt+1,0:ny,0:nx,n_matched_vars) )
     881!
     882!--                Note, at the moment emissions are considered only by surface fluxes
     883!--                rather than by volume sources. Therefore, no vertical dimension is
     884!--                required and is thus allocated with 1. Later when volume sources
     885!--                are considered, the vertical dimension will increase.                 
     886                   !ALLOCATE( emis_distribution(nzb:nzt+1,nys:nyn,nxl:nxr,n_matched_vars) )
     887                   ALLOCATE( emis_distribution(1,nys:nyn,nxl:nxr,n_matched_vars) )
    879888                ENDIF
    880889 
  • palm/trunk/SOURCE/chemistry_model_mod.f90

    r4115 r4154  
    2727! -----------------
    2828! $Id$
     29! Bugfix in decycling
     30!
     31! 4115 2019-07-24 12:50:49Z suehring
    2932! Fix faulty IF statement in decycling initialization
    3033!
     
    27292732          DO WHILE ( TRIM( cs_name( lsp_usr ) ) /= 'novalue' )
    27302733             IF ( TRIM(chem_species(lsp)%name) == TRIM(cs_name(lsp_usr)) )  THEN
    2731 
    2732                 CALL chem_boundary_conds( chem_species(lsp)%conc_p,                     &
    2733                                                   chem_species(lsp)%conc_pr_init )
     2734!
     2735!--             As chem_exchange_horiz_bounds is called at the beginning
     2736!--             of prognostic_equations, boundary conditions are set on
     2737!--             %conc.
     2738                CALL chem_boundary_conds( chem_species(lsp)%conc,              &
     2739                                          chem_species(lsp)%conc_pr_init )
    27342740             
    27352741             ENDIF
    2736              lsp_usr = lsp_usr +1
     2742             lsp_usr = lsp_usr + 1
    27372743          ENDDO
    27382744
Note: See TracChangeset for help on using the changeset viewer.