Changeset 135 for palm

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

Files:

• DOC/app/chapter_3.5.4.html (modified) (2 diffs)
• SOURCE/CURRENT_MODIFICATIONS (modified) (1 diff)
• SOURCE/init_1d_model.f90 (modified) (2 diffs)
• SOURCE/user_interface.f90 (modified) (1 diff)

palm/trunk/DOC/app/chapter_3.5.4.html

@@ -192 +192 @@
 output of 3d volume data, the gridpoint data of the quantity has to be
 resorted to array <span style="font-family: monospace;">local_pf</span>
-in subroutine <a href="chapter_3.5.1.html#user_data_output_3d"><span style="font-family: monospace;">user_data_output_3d</span></a>.
-Also the vertical grid, on which the quantity is defined, has to be set
-again:<br><br><span style="font-family: monospace;">&nbsp;&nbsp;&nbsp;&nbsp;CASE
+in subroutine <a href="chapter_3.5.1.html#user_data_output_3d"><span style="font-family: monospace;">user_data_output_3d</span></a>.:<br><br><span style="font-family: monospace;">&nbsp;&nbsp;&nbsp;&nbsp;CASE
 ( 'u2' )</span><br style="font-family: monospace;"><span style="font-family: monospace;">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
 IF ( av == 0 )&nbsp; THEN</span><br style="font-family: monospace;"><span style="font-family: monospace;">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
@@ -212 +210 @@
 ENDDO</span><br style="font-family: monospace;"><span style="font-family: monospace;">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
 ENDDO</span><br style="font-family: monospace;"><span style="font-family: monospace;">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
-ENDIF</span><br style="font-family: monospace;"><br style="font-family: monospace;"><span style="font-family: monospace;">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
-grid = 'zu'</span><br><br>The <span style="font-family: monospace;">ELSE</span> case is
+ENDIF</span><br><br>The <span style="font-family: monospace;">ELSE</span> case is
 only needed in case that output of time-averaged data is requested.<br><br><br></li><li>In
 case of job chains, the sum array has to be written to the (binary)

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