Changeset 183 for palm/trunk/SOURCE

Timestamp:

Aug 4, 2008 3:39:12 PM (16 years ago)

Author:

letzel

Message:

• mrun: (optionally) link local input
• diffusion_s: bugfix for calculation of horizontal fluxes at vertical walls
• NCL scripts in trunk/SCRIPTS/NCL updated

Location:

palm/trunk/SOURCE

Files:

• CURRENT_MODIFICATIONS (modified) (2 diffs)
• diffusion_s.f90 (modified) (3 diffs)

palm/trunk/SOURCE/CURRENT_MODIFICATIONS

@@ -81 +81 @@
 Errors:
 ------
+Bugfix: calculation of horizontal fluxes at vertical walls (diffusion_s)
+
 Bugfix: zero assignments to tendency arrays in case of restarts (init_3d_model)
 
@@ -117 +119 @@
 
 
-check_parameters, flow_statistics, init_dvrp, init_3d_model, local_stop, plant_canopy_model, poismg, pres, read_3d_binary, user_interface, write_3d_binary
+check_parameters, diffusion_s, flow_statistics, init_dvrp, init_3d_model, local_stop, plant_canopy_model, poismg, pres, read_3d_binary, user_interface, write_3d_binary
 

palm/trunk/SOURCE/diffusion_s.f90

@@ -4 +4 @@
 ! Actual revisions:
 ! -----------------
-! 
+! bugfix: calculation of fluxes at vertical surfaces
 !
 ! Former revisions:
@@ -90 +90 @@
 
                    tend(k,j,i) = tend(k,j,i)                                  &
-                                          + 0.5 * ( fwxp(j,i) *               &
+                                                + ( fwxp(j,i) * 0.5 *         &
                         ( kh(k,j,i) + kh(k,j,i+1) ) * ( s(k,j,i+1)-s(k,j,i) ) &
                         + ( 1.0 - fwxp(j,i) ) * wall_s_flux(1)                &
-                                                   -fwxm(j,i) *               &
+                                                   -fwxm(j,i) * 0.5 *         &
                         ( kh(k,j,i) + kh(k,j,i-1) ) * ( s(k,j,i)-s(k,j,i-1) ) &
                         + ( 1.0 - fwxm(j,i) ) * wall_s_flux(2)                &
                                                   ) * ddx2                    &
-                                          + 0.5 * ( fwyp(j,i) *               &
+                                                + ( fwyp(j,i) * 0.5 *         &
                         ( kh(k,j,i) + kh(k,j+1,i) ) * ( s(k,j+1,i)-s(k,j,i) ) &
                         + ( 1.0 - fwyp(j,i) ) * wall_s_flux(3)                &
-                                                   -fwym(j,i) *               &
+                                                   -fwym(j,i) * 0.5 *         &
                         ( kh(k,j,i) + kh(k,j-1,i) ) * ( s(k,j,i)-s(k,j-1,i) ) &
                         + ( 1.0 - fwym(j,i) ) * wall_s_flux(4)                &
@@ -200 +200 @@
 
              tend(k,j,i) = tend(k,j,i)                                        &
-                                          + 0.5 * ( fwxp(j,i) *               &
+                                                + ( fwxp(j,i) * 0.5 *         &
                         ( kh(k,j,i) + kh(k,j,i+1) ) * ( s(k,j,i+1)-s(k,j,i) ) &
                         + ( 1.0 - fwxp(j,i) ) * wall_s_flux(1)                &
-                                                   -fwxm(j,i) *               &
+                                                   -fwxm(j,i) * 0.5 *         &
                         ( kh(k,j,i) + kh(k,j,i-1) ) * ( s(k,j,i)-s(k,j,i-1) ) &
                         + ( 1.0 - fwxm(j,i) ) * wall_s_flux(2)                &
                                                   ) * ddx2                    &
-                                          + 0.5 * ( fwyp(j,i) *               &
+                                                + ( fwyp(j,i) * 0.5 *         &
                         ( kh(k,j,i) + kh(k,j+1,i) ) * ( s(k,j+1,i)-s(k,j,i) ) &
                         + ( 1.0 - fwyp(j,i) ) * wall_s_flux(3)                &
-                                                   -fwym(j,i) *               &
+                                                   -fwym(j,i) * 0.5 *         &
                         ( kh(k,j,i) + kh(k,j-1,i) ) * ( s(k,j,i)-s(k,j-1,i) ) &
                         + ( 1.0 - fwym(j,i) ) * wall_s_flux(4)                &