Changeset 3897 for palm/trunk/SOURCE/chemistry_model_mod.f90

Timestamp:

Apr 15, 2019 11:51:14 AM (5 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:

• palm/trunk/SOURCE/chemistry_model_mod.f90 (modified) (4 diffs)

palm/trunk/SOURCE/chemistry_model_mod.f90

@@ -27 +27 @@
 ! -----------------
 ! $Id: chemistry_model_mod.f90 3784 2019-03-05 14:16:20Z banzhafs
+! 2D output of emission fluxes
+! 
+! 3784 2019-03-05 14:16:20Z banzhafs
 ! Bugfix, uncomment erroneous commented variable used for dry deposition. 
 ! Bugfix in 3D emission output. 
@@ -1315 +1318 @@
     INTEGER(iwp) ::  j               !< grid index along y-direction
     INTEGER(iwp) ::  k               !< grid index along z-direction
+    INTEGER(iwp) ::  l               !< running index for orientation of surface elements
+    INTEGER(iwp) ::  m               !< running indices for surfaces
     INTEGER(iwp) ::  char_len        !< length of a character string
 !
@@ -1324 +1329 @@
 
     spec_name = TRIM( variable(4:char_len-3) )
+!
+!-- Output of emission values, i.e. surface fluxes cssws. 
+    IF ( variable(1:3) == 'em_' )  THEN
+
+       local_pf = 0.0_wp
+
+       DO  lsp = 1, nvar
+          IF ( TRIM( spec_name ) == TRIM( chem_species(lsp)%name) )  THEN
+!
+!--          No average output for now. 
+             DO  m = 1, surf_lsm_h%ns
+                local_pf(surf_lsm_h%i(m),surf_lsm_h%j(m),nzb+1) =              &
+                              local_pf(surf_lsm_h%i(m),surf_lsm_h%j(m),nzb+1)  &
+                                            + surf_lsm_h%cssws(lsp,m)
+             ENDDO
+             DO  m = 1, surf_usm_h%ns
+                   local_pf(surf_usm_h%i(m),surf_usm_h%j(m),nzb+1) =           &
+                              local_pf(surf_usm_h%i(m),surf_usm_h%j(m),nzb+1)  &
+                                            + surf_usm_h%cssws(lsp,m)
+             ENDDO
+             grid = 'zu'
+             found = .TRUE.
+          ENDIF
+       ENDDO
+
+    ELSE
 
        DO  lsp=1,nspec
@@ -1362 +1393 @@
           ENDIF
        ENDDO
-
-       RETURN
+    ENDIF
+
+    RETURN
 
  END SUBROUTINE chem_data_output_2d