Changeset 2582 for palm/trunk/SOURCE/pmc_interface_mod.f90

Timestamp:

Oct 26, 2017 1:19:46 PM (7 years ago)

Author:

hellstea

Message:

pmc bugfix and gfortran workaround

File:

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

palm/trunk/SOURCE/pmc_interface_mod.f90

@@ -26 +26 @@
 ! -----------------
 ! $Id$
+! Resetting of e within buildings / topography in pmci_parent_datatrans removed
+! as unnecessary since e is not anterpolated, and as incorrect since it overran
+! the default Neumann condition (bc_e_b). 
+! 
+! 2359 2017-08-21 07:50:30Z hellstea
 ! A minor indexing error in pmci_init_loglaw_correction is corrected. 
 ! 
@@ -3640 +3645 @@
 
 !
-!--          Inside buildings/topography reset velocities and TKE back to zero.
-!--          Other scalars (pt, q, s, km, kh, p, sa, ...) are ignored at
+!--          Inside buildings/topography reset velocities back to zero.
+!--          Scalars (pt, q, s, km, kh, p, sa, ...) are ignored at
 !--          present, maybe revise later.
+!--          Resetting of e is removed as unnecessary since e is not 
+!--          anterpolated, and as incorrect since it overran the default 
+!--          Neumann condition (bc_e_b). 
              DO   i = nxlg, nxrg
                 DO   j = nysg, nyng
@@ -3652 +3660 @@
                       w(k,j,i)  = MERGE( w(k,j,i), 0.0_wp,                     &
                                          BTEST( wall_flags_0(k,j,i), 3 ) )
-                      e(k,j,i)  = MERGE( e(k,j,i), 0.0_wp,                     &
-                                         BTEST( wall_flags_0(k,j,i), 0 ) )
 !
 !--                TO_DO: zero setting of temperature within topography creates