Changeset 3897 for palm/trunk/SOURCE

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

Location:

palm/trunk/SOURCE

Files:

• chemistry_model_mod.f90 (modified) (4 diffs)
• init_pegrid.f90 (modified) (2 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      

palm/trunk/SOURCE/init_pegrid.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Minor revision of multigrid check; give warning instead of an abort.
+! 
+! 3890 2019-04-12 15:59:20Z suehring
 ! Check if grid coarsening is possible on subdomain, in order to avoid that 
 ! multigrid approach effectively reduces to a Gauss-Seidel scheme.
@@ -955 +958 @@
 !--    This case, the maximum number of grid levels is 1, i.e. effectively 
 !--    a Gauss-Seidel scheme is applied rather than a multigrid approach.
-       IF ( maximum_grid_level == 1 )  THEN
-          message_string = 'grid coarsening on subdomain level cannot be ' //  &
-                           'performed'
-          CALL message( 'poismg', 'PA0236', 1, 2, 0, 6, 0 )
+!--    Give a warning in this case.
+       IF ( maximum_grid_level == 1  .AND.  mg_switch_to_pe0_level == -1 )  THEN
+          message_string = 'No grid coarsening possible, multigrid ' //        &
+                           'approach effectively reduces to a Gauss-Seidel ' //&
+                           'scheme.'
+          
+          CALL message( 'poismg', 'PA0648', 0, 1, 0, 6, 0 )
        ENDIF