Changeset 135


Ignore:
Timestamp:
Nov 22, 2007 12:24:23 PM (15 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:
4 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/DOC/app/chapter_3.5.4.html

    r130 r135  
    192192output of 3d volume data, the gridpoint data of the quantity has to be
    193193resorted to array <span style="font-family: monospace;">local_pf</span>
    194 in subroutine <a href="chapter_3.5.1.html#user_data_output_3d"><span style="font-family: monospace;">user_data_output_3d</span></a>.
    195 Also the vertical grid, on which the quantity is defined, has to be set
    196 again:<br><br><span style="font-family: monospace;">&nbsp;&nbsp;&nbsp;&nbsp;CASE
     194in 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
    197195( 'u2' )</span><br style="font-family: monospace;"><span style="font-family: monospace;">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
    198196IF ( av == 0 )&nbsp; THEN</span><br style="font-family: monospace;"><span style="font-family: monospace;">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
     
    212210ENDDO</span><br style="font-family: monospace;"><span style="font-family: monospace;">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
    213211ENDDO</span><br style="font-family: monospace;"><span style="font-family: monospace;">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
    214 ENDIF</span><br style="font-family: monospace;"><br style="font-family: monospace;"><span style="font-family: monospace;">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
    215 grid = 'zu'</span><br><br>The <span style="font-family: monospace;">ELSE</span> case is
     212ENDIF</span><br><br>The <span style="font-family: monospace;">ELSE</span> case is
    216213only needed in case that output of time-averaged data is requested.<br><br><br></li><li>In
    217214case of job chains, the sum array has to be written to the (binary)
  • palm/trunk/SOURCE/CURRENT_MODIFICATIONS

    r134 r135  
    6565Bugfix: assignment of fluxes at walls
    6666
    67 advec_particles, check_parameters, combine_plot_fields, diffusion_s, init_ocean, poisfft, production_e
     67Bugfix: absolute value of f must be used when calculating the Blackadar mixing length (init_1d_model)
     68
     69advec_particles, check_parameters, combine_plot_fields, diffusion_s, init_ocean, init_1d_model, poisfft, production_e
  • palm/trunk/SOURCE/init_1d_model.f90

    r83 r135  
    44! Actual revisions:
    55! -----------------
    6 !
     6! Bugfix: absolute value of f must be used when calculating the Blackadar
     7! mixing length!
    78!
    89! Former revisions:
     
    7677!--       Blackadar mixing length
    7778          IF ( f /= 0.0 )  THEN
    78              lambda = 2.7E-4 * SQRT( ug(nzt+1)**2 + vg(nzt+1)**2 ) / f + 1E-10
     79             lambda = 2.7E-4 * SQRT( ug(nzt+1)**2 + vg(nzt+1)**2 ) / &
     80                               ABS( f ) + 1E-10
    7981          ELSE
    8082             lambda = 30.0
  • palm/trunk/SOURCE/user_interface.f90

    r130 r135  
    946946!          ENDIF
    947947!
    948 !          grid = 'zu'
    949948
    950949       CASE DEFAULT
Note: See TracChangeset for help on using the changeset viewer.