Changeset 3897 for palm


Ignore:
Timestamp:
Apr 15, 2019 11:51:14 AM (2 years ago)
Author:
suehring
Message:

2D output of emission fluxes in chemistry; revise check for multigrid coarsening and give a warning instead of an error

Location:
palm/trunk/SOURCE
Files:
2 edited

Legend:

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

    r3889 r3897  
    2727! -----------------
    2828! $Id: chemistry_model_mod.f90 3784 2019-03-05 14:16:20Z banzhafs
     29! 2D output of emission fluxes
     30!
     31! 3784 2019-03-05 14:16:20Z banzhafs
    2932! Bugfix, uncomment erroneous commented variable used for dry deposition.
    3033! Bugfix in 3D emission output.
     
    13151318    INTEGER(iwp) ::  j               !< grid index along y-direction
    13161319    INTEGER(iwp) ::  k               !< grid index along z-direction
     1320    INTEGER(iwp) ::  l               !< running index for orientation of surface elements
     1321    INTEGER(iwp) ::  m               !< running indices for surfaces
    13171322    INTEGER(iwp) ::  char_len        !< length of a character string
    13181323!
     
    13241329
    13251330    spec_name = TRIM( variable(4:char_len-3) )
     1331!
     1332!-- Output of emission values, i.e. surface fluxes cssws.
     1333    IF ( variable(1:3) == 'em_' )  THEN
     1334
     1335       local_pf = 0.0_wp
     1336
     1337       DO  lsp = 1, nvar
     1338          IF ( TRIM( spec_name ) == TRIM( chem_species(lsp)%name) )  THEN
     1339!
     1340!--          No average output for now.
     1341             DO  m = 1, surf_lsm_h%ns
     1342                local_pf(surf_lsm_h%i(m),surf_lsm_h%j(m),nzb+1) =              &
     1343                              local_pf(surf_lsm_h%i(m),surf_lsm_h%j(m),nzb+1)  &
     1344                                            + surf_lsm_h%cssws(lsp,m)
     1345             ENDDO
     1346             DO  m = 1, surf_usm_h%ns
     1347                   local_pf(surf_usm_h%i(m),surf_usm_h%j(m),nzb+1) =           &
     1348                              local_pf(surf_usm_h%i(m),surf_usm_h%j(m),nzb+1)  &
     1349                                            + surf_usm_h%cssws(lsp,m)
     1350             ENDDO
     1351             grid = 'zu'
     1352             found = .TRUE.
     1353          ENDIF
     1354       ENDDO
     1355
     1356    ELSE
    13261357
    13271358       DO  lsp=1,nspec
     
    13621393          ENDIF
    13631394       ENDDO
    1364 
    1365        RETURN
     1395    ENDIF
     1396
     1397    RETURN
    13661398
    13671399 END SUBROUTINE chem_data_output_2d     
  • palm/trunk/SOURCE/init_pegrid.f90

    r3890 r3897  
    2525! -----------------
    2626! $Id$
     27! Minor revision of multigrid check; give warning instead of an abort.
     28!
     29! 3890 2019-04-12 15:59:20Z suehring
    2730! Check if grid coarsening is possible on subdomain, in order to avoid that
    2831! multigrid approach effectively reduces to a Gauss-Seidel scheme.
     
    955958!--    This case, the maximum number of grid levels is 1, i.e. effectively
    956959!--    a Gauss-Seidel scheme is applied rather than a multigrid approach.
    957        IF ( maximum_grid_level == 1 )  THEN
    958           message_string = 'grid coarsening on subdomain level cannot be ' //  &
    959                            'performed'
    960           CALL message( 'poismg', 'PA0236', 1, 2, 0, 6, 0 )
     960!--    Give a warning in this case.
     961       IF ( maximum_grid_level == 1  .AND.  mg_switch_to_pe0_level == -1 )  THEN
     962          message_string = 'No grid coarsening possible, multigrid ' //        &
     963                           'approach effectively reduces to a Gauss-Seidel ' //&
     964                           'scheme.'
     965         
     966          CALL message( 'poismg', 'PA0648', 0, 1, 0, 6, 0 )
    961967       ENDIF
    962968
Note: See TracChangeset for help on using the changeset viewer.