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

Timestamp:

May 9, 2018 8:42:38 AM (6 years ago)

Author:

maronga

Message:

series of bugfixes

File:

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

palm/trunk/SOURCE/chemistry_model_mod.f90

@@ -27 +27 @@
 ! -----------------
 ! $Id$
+! Bugfix: nzb_do and nzt_do were not used for 3d data output
+! 
+! 3004 2018-04-27 12:33:25Z Giersch
 ! Comment concerning averaged data output added 
 ! 
@@ -1152 +1155 @@
 
 
-   SUBROUTINE chem_data_output_3d( av, variable, found, local_pf, fill_value )
+   SUBROUTINE chem_data_output_3d( av, variable, found, local_pf, fill_value, nzb_do, nzt_do )
 
 
@@ -1163 +1166 @@
 
       CHARACTER (LEN=*) ::  variable !<
+
+      INTEGER(iwp) ::  av    !<
+      INTEGER(iwp) ::  nzb_do !< lower limit of the data output (usually 0)
+      INTEGER(iwp) ::  nzt_do !< vertical upper limit of the data output (usually nz_do3d)
+
       LOGICAL      ::  found !<
-      INTEGER(iwp) ::  av    !<
 
       REAL(wp) ::  fill_value !<
-      REAL(sp), DIMENSION(nxl:nxr,nys:nyn,nzb:nzt+1) ::  local_pf 
+      REAL(sp), DIMENSION(nxl:nxr,nys:nyn,nzb_do:nzt_do) ::  local_pf 
 
 
@@ -1189 +1196 @@
                 DO  i = nxl, nxr
                    DO  j = nys, nyn
-                      DO  k = nzb, nzt+1
+                      DO  k = nzb_do, nzt_do
                           local_pf(i,j,k) = MERGE(                             &
                                               chem_species(lsp)%conc(k,j,i),   &
@@ -1201 +1208 @@
                 DO  i = nxl, nxr
                    DO  j = nys, nyn
-                      DO  k = nzb, nzt+1
+                      DO  k = nzb_do, nzt_do
                           local_pf(i,j,k) = MERGE(                             &
                                               chem_species(lsp)%conc_av(k,j,i),&