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

Timestamp:

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

Author:

maronga

Message:

series of bugfixes

File:

• palm/trunk/SOURCE/turbulence_closure_mod.f90 (modified) (14 diffs)

palm/trunk/SOURCE/turbulence_closure_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Bugfix: nzb_do and nzt_do were not used for 3d data output
+! 
+! 3004 2018-04-27 12:33:25Z Giersch
 ! Further allocation checks implemented
 ! 
@@ -626 +629 @@
     REAL(wp) ::  fill_value = -999.0_wp  !< value for the _FillValue attribute
 
-    REAL(wp), DIMENSION(nxl:nxr,nys:nyn,nzb:nzt+1) ::  local_pf !< local
+    REAL(wp), DIMENSION(nxl:nxr,nys:nyn,nzb_do:nzt_do) ::  local_pf !< local
        !< array to which output data is resorted to
 
@@ -723 +726 @@
 !> Define 3D output variables.
 !------------------------------------------------------------------------------!
- SUBROUTINE tcm_data_output_3d( av, variable, found, local_pf )
+ SUBROUTINE tcm_data_output_3d( av, variable, found, local_pf, nzb_do, nzt_do )
  
 
@@ -733 +736 @@
     CHARACTER (LEN=*) ::  variable   !< 
 
-    INTEGER(iwp) ::  av   !< 
-    INTEGER(iwp) ::  i    !< 
-    INTEGER(iwp) ::  j    !< 
-    INTEGER(iwp) ::  k    !< 
+    INTEGER(iwp) ::  av     !< 
+    INTEGER(iwp) ::  i      !< 
+    INTEGER(iwp) ::  j      !< 
+    INTEGER(iwp) ::  k      !< 
+    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   !< 
@@ -742 +747 @@
     REAL(wp) ::  fill_value = -999.0_wp  !< value for the _FillValue attribute
 
-    REAL(sp), DIMENSION(nxl:nxr,nys:nyn,nzb:nzt+1) ::  local_pf   !< local
+    REAL(sp), DIMENSION(nxl:nxr,nys:nyn,nzb_do:nzt_do) ::  local_pf   !< local
        !< array to which output data is resorted to
 
@@ -756 +761 @@
              DO  i = nxl, nxr
                 DO  j = nys, nyn
-                   DO  k = nzb, nzt+1
+                   DO  k = nzb_do, nzt_do
                       local_pf(i,j,k) = diss(k,j,i)
                    ENDDO
@@ -768 +773 @@
              DO  i = nxl, nxr
                 DO  j = nys, nyn
-                   DO  k = nzb, nzt+1
+                   DO  k = nzb_do, nzt_do
                       local_pf(i,j,k) = diss_av(k,j,i)
                    ENDDO
@@ -779 +784 @@
              DO  i = nxl, nxr
                 DO  j = nys, nyn
-                   DO  k = nzb, nzt+1
+                   DO  k = nzb_do, nzt_do
                       local_pf(i,j,k) = kh(k,j,i)
                    ENDDO
@@ -791 +796 @@
              DO  i = nxl, nxr
                 DO  j = nys, nyn
-                   DO  k = nzb, nzt+1
+                   DO  k = nzb_do, nzt_do
                       local_pf(i,j,k) = kh_av(k,j,i)
                    ENDDO
@@ -802 +807 @@
              DO  i = nxl, nxr
                 DO  j = nys, nyn
-                   DO  k = nzb, nzt+1
+                   DO  k = nzb_do, nzt_do
                       local_pf(i,j,k) = km(k,j,i)
                    ENDDO
@@ -814 +819 @@
              DO  i = nxl, nxr
                 DO  j = nys, nyn
-                   DO  k = nzb, nzt+1
+                   DO  k = nzb_do, nzt_do
                       local_pf(i,j,k) = km_av(k,j,i)
                    ENDDO
@@ -825 +830 @@
              DO  i = nxl, nxr
                 DO  j = nys, nyn
-                   DO  k = nzb, nzt+1
+                   DO  k = nzb_do, nzt_do
                       local_pf(i,j,k) = dummy1(k,j,i)
                    ENDDO
@@ -836 +841 @@
              DO  i = nxl, nxr
                 DO  j = nys, nyn
-                   DO  k = nzb, nzt+1
+                   DO  k = nzb_do, nzt_do
                       local_pf(i,j,k) = dummy2(k,j,i)
                    ENDDO
@@ -847 +852 @@
              DO  i = nxl, nxr
                 DO  j = nys, nyn
-                   DO  k = nzb, nzt+1
+                   DO  k = nzb_do, nzt_do
                       local_pf(i,j,k) = dummy3(k,j,i)
                    ENDDO