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

Timestamp:

Dec 14, 2017 5:12:51 PM (6 years ago)

Author:

kanani

Message:

Merge of branch palm4u into trunk

Location:

palm/trunk

Files:

• . (modified) (1 prop)
• SOURCE (modified) (1 prop)
• SOURCE/pmc_interface_mod.f90 (modified) (8 diffs, 1 prop)

palm/trunk/SOURCE/pmc_interface_mod.f90

@@ -1 @@
-MODULE pmc_interface
-
+!> @file pmc_interface_mod.f90
 !------------------------------------------------------------------------------!
-! This file is part of PALM.
+! This file is part of the PALM model system.
 !
 ! PALM is free software: you can redistribute it and/or modify it under the
@@ -26 +25 @@
 ! -----------------
 ! $Id$
+! - Bugfix in init_tke_factor (MS)
+! 
+! 2669 2017-12-06 16:03:27Z raasch
 ! file extension for nested domains changed to "_N##",
 ! created flag file in order to enable combine_plot_fields to process nest data
@@ -194 +196 @@
 ! @todo Data transfer of qc and nc is prepared but not activated
 !-------------------------------------------------------------------------------!
+ MODULE pmc_interface
 
 #if defined( __nopointer )
@@ -2685 +2688 @@
              
              DO  k = k_wall + 1, nzt
-                
+ 
                 kc     = kco(k) + 1
                 glsf   = ( dx * dy * dzu(k) )**p13
@@ -3066 +3069 @@
                    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 ) )
+!                    e(k,j,i)   = MERGE( e(k,j,i), 0.0_wp,                       &
+!                                        BTEST( wall_flags_0(k,j,i), 0 ) )
                    u_p(k,j,i) = MERGE( u_p(k,j,i), 0.0_wp,                     &
                                        BTEST( wall_flags_0(k,j,i), 1 ) )
@@ -3074 +3077 @@
                    w_p(k,j,i) = MERGE( w_p(k,j,i), 0.0_wp,                     &
                                        BTEST( wall_flags_0(k,j,i), 3 ) )
-                   e_p(k,j,i) = MERGE( e_p(k,j,i), 0.0_wp,                     &
-                                       BTEST( wall_flags_0(k,j,i), 0 ) )
+!                    e_p(k,j,i) = MERGE( e_p(k,j,i), 0.0_wp,                     &
+!                                        BTEST( wall_flags_0(k,j,i), 0 ) )
                 ENDDO
              ENDDO
@@ -4492 +4495 @@
       ENDIF
 !
-!--   Store the boundary values also into the other redundant ghost node layers
+!--   Store the boundary values also into the other redundant ghost node layers.
+!--   Please note, in case of only one ghost node layer, e.g. for the PW 
+!--   scheme, the following loops will not be entered.
       IF ( edge == 'l' )  THEN
          DO  ibgp = -nbgp, ib
@@ -4713 +4718 @@
       ENDIF
 !
-!--   Store the boundary values also into the other redundant ghost node layers
+!--   Store the boundary values also into the other redundant ghost node layers.
+!--   Please note, in case of only one ghost node layer, e.g. for the PW 
+!--   scheme, the following loops will not be entered.
       IF ( edge == 's' )  THEN
          DO  jbgp = -nbgp, jb