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

Timestamp:

Jan 16, 2018 6:11:14 PM (6 years ago)

Author:

suehring

Message:

Fill values for 3D data output of chemical species introduced

File:

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

palm/trunk/SOURCE/chemistry_model_mod.f90

@@ -27 +27 @@
 ! -----------------
 ! $Id$
+! Fill values in 3D output introduced.
+! 
+! 2718 2018-01-02 08:49:38Z maronga
 ! Initial revision
 !
@@ -1200 +1203 @@
 
 
-   SUBROUTINE chem_data_output_3d( av, variable, found, local_pf )
+   SUBROUTINE chem_data_output_3d( av, variable, found, local_pf, fill_value )
 
 
@@ -1213 +1216 @@
       LOGICAL      ::  found !<
       INTEGER(iwp) ::  av    !<
-      REAL(sp), DIMENSION(nxl:nxr,nys:nyn,nzb:nzt+1) ::  local_pf !<
+
+      REAL(wp) ::  fill_value !<
+      REAL(sp), DIMENSION(nxl:nxr,nys:nyn,nzb:nzt+1) ::  local_pf 
+
 
       !-- local variables
@@ -1234 +1240 @@
                    DO  j = nys, nyn
                       DO  k = nzb, nzt+1
-                          local_pf(i,j,k) = chem_species(lsp)%conc(k,j,i)
+                          local_pf(i,j,k) = MERGE(                             &
+                                              chem_species(lsp)%conc(k,j,i),   &
+                                              REAL( fill_value, KIND = wp ),   &
+                                              BTEST( wall_flags_0(k,j,i), 0 ) )
                       ENDDO
                    ENDDO
@@ -1243 +1252 @@
                    DO  j = nys, nyn
                       DO  k = nzb, nzt+1
-                          local_pf(i,j,k) = chem_species(lsp)%conc_av(k,j,i)
+                          local_pf(i,j,k) = MERGE(                             &
+                                              chem_species(lsp)%conc_av(k,j,i),&
+                                              REAL( fill_value, KIND = wp ),   &
+                                              BTEST( wall_flags_0(k,j,i), 0 ) )
                       ENDDO
                    ENDDO