Changeset 2756

Timestamp:

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

Author:

suehring

Message:

Fill values for 3D data output of chemical species introduced

Location:

palm/trunk/SOURCE

Files:

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

palm/trunk/SOURCE/data_output_3d.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Fill values for 3D output of chemical species introduced.
+! 
+! 2746 2018-01-15 12:06:04Z suehring
 ! Move flag plant canopy to modules
 ! 
@@ -680 +683 @@
              IF ( .NOT. found  .AND.  air_chemistry )  THEN
                 CALL chem_data_output_3d( av, do3d(av,if), found,              &
-                                          local_pf )
+                                          local_pf, fill_value )
                 resorted = .TRUE.
              ENDIF