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

File:
1 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     
Note: See TracChangeset for help on using the changeset viewer.