Changeset 2663 for palm/trunk

Timestamp:

Dec 4, 2017 5:40:50 PM (7 years ago)

Author:

suehring

Message:

Bugfix, wrong coarse-grid index in init_tkefactor used.

File:

• palm/trunk/SOURCE/pmc_interface_mod.f90 (modified) (6 diffs)

palm/trunk/SOURCE/pmc_interface_mod.f90

@@ -26 +26 @@
 ! -----------------
 ! $Id$
+! Bugfix, wrong coarse-grid index in init_tkefactor used.
+! 
+! 2602 2017-11-03 11:06:41Z hellstea
 ! Index-limit related bug (occurred with nesting_mode='vertical') fixed in
 ! pmci_interp_tril_t. Check for too high nest domain added in pmci_setup_parent.   
@@ -2629 +2632 @@
              DO  k = k_wall + 1, nzt
 
-                kc     = kco(k+1)
+                kc     = kco(k) + 1
                 glsf   = ( dx * dy * dzu(k) )**p13
                 glsc   = ( cg%dx * cg%dy *cg%dzu(kc) )**p13
@@ -2650 +2653 @@
              DO  k = k_wall + 1, nzt
 
-                kc     = kco(k+1)
+                kc     = kco(k) + 1
                 glsf   = ( dx * dy * dzu(k) )**p13
                 glsc   = ( cg%dx * cg%dy * cg%dzu(kc) )**p13
@@ -2671 +2674 @@
              DO  k = k_wall + 1, nzt
                 
-                kc     = kco(k+1)
+                kc     = kco(k) + 1
                 glsf   = ( dx * dy * dzu(k) )**p13
                 glsc   = ( cg%dx * cg%dy * cg%dzu(kc) ) ** p13
@@ -2692 +2695 @@
              DO  k = k_wall + 1, nzt
 
-                kc     = kco(k+1)
+                kc     = kco(k) + 1
                 glsf   = ( dx * dy * dzu(k) )**p13
                 glsc   = ( cg%dx * cg%dy * cg%dzu(kc) )**p13
@@ -2713 +2716 @@
              k_wall = get_topography_top_index( j, i, 's' )
 
-             kc     = kco(k+1)
+             kc     = kco(k) + 1
              glsf   = ( dx * dy * dzu(k) )**p13
              glsc   = ( cg%dx * cg%dy * cg%dzu(kc) )**p13