Changeset 135 for palm/trunk/SOURCE

Timestamp:

Nov 22, 2007 12:24:23 PM (17 years ago)

Author:

raasch

Message:

bug concerning mixing length calculation in case of negative f removed; inconsistencies concerning variable grid removed

Location:

palm/trunk/SOURCE

Files:

• CURRENT_MODIFICATIONS (modified) (1 diff)
• init_1d_model.f90 (modified) (2 diffs)
• user_interface.f90 (modified) (1 diff)

palm/trunk/SOURCE/CURRENT_MODIFICATIONS

@@ -65 +65 @@
 Bugfix: assignment of fluxes at walls
 
-advec_particles, check_parameters, combine_plot_fields, diffusion_s, init_ocean, poisfft, production_e
+Bugfix: absolute value of f must be used when calculating the Blackadar mixing length (init_1d_model)
+
+advec_particles, check_parameters, combine_plot_fields, diffusion_s, init_ocean, init_1d_model, poisfft, production_e

palm/trunk/SOURCE/init_1d_model.f90

@@ -4 +4 @@
 ! Actual revisions:
 ! -----------------
-! 
+! Bugfix: absolute value of f must be used when calculating the Blackadar
+! mixing length!
 !
 ! Former revisions:
@@ -76 +77 @@
 !--       Blackadar mixing length
           IF ( f /= 0.0 )  THEN
-             lambda = 2.7E-4 * SQRT( ug(nzt+1)**2 + vg(nzt+1)**2 ) / f + 1E-10
+             lambda = 2.7E-4 * SQRT( ug(nzt+1)**2 + vg(nzt+1)**2 ) / &
+                               ABS( f ) + 1E-10
           ELSE
              lambda = 30.0

palm/trunk/SOURCE/user_interface.f90

@@ -946 +946 @@
 !          ENDIF
 !
-!          grid = 'zu'
 
        CASE DEFAULT