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

Timestamp:

Jul 26, 2018 12:55:23 PM (6 years ago)

Author:

suehring

Message:

Further revision of 2D surface output for radiation and chemistry quantities; bugfix for commit 3170

File:

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

palm/trunk/SOURCE/chemistry_model_mod.f90

@@ -27 +27 @@
 ! -----------------
 ! $Id$
+! Revise output of surface quantities in case of overhanging structures
+! 
+! 3045 2018-05-28 07:55:41Z Giersch
 ! error messages revised
 ! 
@@ -1257 +1260 @@
     CHARACTER (LEN=*) :: variable !< 
  
+    LOGICAL      ::  match_def !< flag indicating natural-type surface
+    LOGICAL      ::  match_lsm !< flag indicating natural-type surface
+    LOGICAL      ::  match_usm !< flag indicating urban-type surface
 
     INTEGER(iwp) ::  i                  !< grid index x direction
@@ -1288 +1294 @@
              ENDDO
           ELSEIF ( TRIM(variable(4:)) == TRIM('cssws*') )        THEN
-             DO  m = 1, surf_def_h(0)%ns
-                 i = surf_def_h(0)%i(m)
-                 j = surf_def_h(0)%j(m)
-                 chem_species(lsp)%cssws_av(j,i) = chem_species(lsp)%cssws_av(j,i) + surf_def_h(0)%cssws(lsp,m)
-             ENDDO 
-             DO  m = 1, surf_lsm_h%ns
-                 i = surf_lsm_h%i(m)
-                 j = surf_lsm_h%j(m)
-                 chem_species(lsp)%cssws_av(j,i) = chem_species(lsp)%cssws_av(j,i) + surf_lsm_h%cssws(lsp,m)
+             DO  i = nxl, nxr
+                DO  j = nys, nyn
+                   match_def = surf_def_h(0)%start_index(j,i) <=               &
+                               surf_def_h(0)%end_index(j,i)
+                   match_lsm = surf_lsm_h%start_index(j,i) <=                  &
+                               surf_lsm_h%end_index(j,i)
+                   match_usm = surf_usm_h%start_index(j,i) <=                  &
+                               surf_usm_h%end_index(j,i)
+
+                   IF ( match_def )  THEN
+                      m = surf_def_h(0)%end_index(j,i)
+                      chem_species(lsp)%cssws_av(j,i) =                        &
+                                             chem_species(lsp)%cssws_av(j,i) + &
+                                             surf_def_h(0)%cssws(lsp,m)
+                   ELSEIF ( match_lsm  .AND.  .NOT. match_usm )  THEN
+                      m = surf_lsm_h%end_index(j,i)
+                      chem_species(lsp)%cssws_av(j,i) =                        &
+                                             chem_species(lsp)%cssws_av(j,i) + &
+                                             surf_lsm_h%cssws(lsp,m)
+                   ELSEIF ( match_usm )  THEN
+                      m = surf_usm_h%end_index(j,i)
+                      chem_species(lsp)%cssws_av(j,i) =                        &
+                                             chem_species(lsp)%cssws_av(j,i) + &
+                                             surf_usm_h%cssws(lsp,m)
+                   ENDIF
+                ENDDO
              ENDDO
-             DO  m = 1, surf_usm_h%ns
-                 i = surf_usm_h%i(m)
-                 j = surf_usm_h%j(m)
-                 chem_species(lsp)%cssws_av(j,i) = chem_species(lsp)%cssws_av(j,i) + surf_usm_h%cssws(lsp,m)
-             ENDDO
-
           ENDIF
        ENDDO