Changeset 531 for palm


Ignore:
Timestamp:
Apr 21, 2010 6:47:21 AM (13 years ago)
Author:
heinze
Message:

call of subsidence added in prognostic equations for humidity/passive scalar, some bugfixes

Location:
palm/trunk
Files:
6 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/DOC/app/chapter_4.1.html

    r411 r531  
    11<!DOCTYPE html PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
    22<html><head>
     3
     4
    35
    46
     
    1416 
    1517 
    16   <meta http-equiv="content-type" content="text/html; charset=ISO-8859-1">
    17 
    18 
    19 
    20 
    21 
    22 
    23  
    24  
    25  
    26  
    27  
    28  
    29   <title>PALM chapter 4.1</title></head>
    30 <body>
     18  <meta http-equiv="content-type" content="text/html; charset=ISO-8859-1"><title>PALM chapter 4.1</title></head><body>
    3119
    3220
     
    34133401      <br>
    34143402      <br>
    3415          Sets the time period a precursor run shall run uncoupled. This parameter is used to set up the precursor run control for atmosphere-ocean-coupled runs. It has to be set individually to the atmospheric / oceanic precursor run. The time in the data output will show negative values during the precursor run. See documentation  for further information.
    3416      
    3417       </td>
     3403Sets the time period a precursor run shall run uncoupled. This
     3404parameter is used to set up the precursor run control for
     3405atmosphere-ocean-coupled runs. It has to be set individually to the
     3406atmospheric / oceanic precursor run. The time in the data output will
     3407show negative values during the precursor run. See documentation for
     3408further information. </td>
    34183409
    34193410    </tr>
    34203411
    3421       <td style="vertical-align: top;"><a name="cthf"></a><span style="font-weight: bold;">cthf</span></td>
     3412      <tr><td style="vertical-align: top;"><a name="cthf"></a><span style="font-weight: bold;">cthf</span></td>
    34223413
    34233414      <td style="vertical-align: top;">R</td>
     
    1564315634
    1564415635A more effective control of
    15645 the &ldquo;overshoots&rdquo; can be achieved with parameter <a href="#cut_spline_overshoot">cut_spline_overshoot</a>.
     15636the &#8220;overshoots&#8221; can be achieved with parameter <a href="#cut_spline_overshoot">cut_spline_overshoot</a>.
    1564615637      </td>
    1564715638
     
    1724917240
    1725017241
    17251  </tr>
     17242 
    1725217243
    1725317244 <tr>
     
    1730317294   <p><span style="font-weight: bold;">Attention:</span><br>
    1730417295
    17305       The large scale vertical motion is only applied to the prognostic equation for the potential
    17306       temperature because it cannot be applied to the momentum equations due to incompressibility.
     17296      The large scale vertical motion is only applied to the prognostic equation for the scalar quantities (potential
     17297      temperature, humidity if  <a href="chapter_4.1.html#humidity">humidity</a> = .T. or passive scalar if <a href="chapter_4.1.html#passive_scalar">passive_scalar</a> = .T.) because it cannot be applied to the momentum equations due to incompressibility.
    1730717298      Thus, the model is not mass consistent.</p>
    1730817299
     
    1731317304
    1731417305
    17315  </tr>
     17306 
    1731617307
    1731717308 <tr>
     
    1734417335default values result in a profile which is zero everywhere regardless of the
    1734517336values of <a href="#ws_vertical_gradient">ws_vertical_gradient</a>.
    17346 For the piecewise construction of temperature profiles see <a href="#ws_vertical_gradient">ws_vertical_gradient</a>.</p>
     17337For the piecewise construction of the subsidence/ascent velocity profile see <a href="#ws_vertical_gradient">ws_vertical_gradient</a>.</p>
    1734717338      </td>
    1734817339
  • palm/trunk/DOC/app/chapter_4.2.html

    r524 r531  
    11<!DOCTYPE html PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
    22<html><head>
     3
     4
    35
    46
     
    48074809
    48084810
    4809  </tr>
     4811 </tr><tr>
     4812            <td style="vertical-align: top;"><font color="#ff6666"><i>hyp</i></font></td>
     4813            <td style="vertical-align: top;">Hydrostatic pressure (in dbar)<br>
     4814            </td>
     4815          </tr>
     4816
    48104817
    48114818
  • palm/trunk/SOURCE/check_parameters.f90

    r525 r531  
    44! Current revisions:
    55! -----------------
    6 !
     6! Bugfix: unit of hyp changed to dbar
    77!
    88! Former revisions:
     
    22052205          CASE ( 'hyp' )
    22062206             dopr_index(i) = 72
    2207              dopr_unit(i)  = 'kPa'
     2207             dopr_unit(i)  = 'dbar'
    22082208             hom(:,2,72,:) = SPREAD( zu, 2, statistic_regions+1 )
    22092209
  • palm/trunk/SOURCE/flow_statistics.f90

    r392 r531  
    880880       hom(:,1,39,sr) = sums(:,38) / ( sums(:,32) + 1E-20 )**1.5    ! Sw
    881881       hom(:,1,40,sr) = sums(:,40)     ! p
    882        hom(:,1,45,sr) = sums(:,45)     ! w"q"
     882       hom(:,1,45,sr) = sums(:,45)     ! w"vpt"
    883883       hom(:,1,46,sr) = sums(:,46)     ! w*vpt*       
    884884       hom(:,1,47,sr) = sums(:,45) + sums(:,46)    ! wvpt
     
    907907       hom(:,1,70,sr) = sums(:,70)     ! q*2
    908908       hom(:,1,71,sr) = sums(:,71)     ! prho
    909        hom(:,1,72,sr) = hyp * 1E-4     ! hyp in kPa
     909       hom(:,1,72,sr) = hyp * 1E-4     ! hyp in dbar
    910910
    911911       hom(:,1,pr_palm-1,sr) = sums(:,pr_palm-1)
  • palm/trunk/SOURCE/modules.f90

    r520 r531  
    55! Current revisions:
    66! -----------------
    7 !
     7! character length of dopr_unit enlarged
    88!
    99! Former revisions:
     
    10191019             ( 'unknown', i9 = 1, dots_max-23 ) /)
    10201020
    1021     CHARACTER (LEN=7), DIMENSION(300) ::  dopr_unit = 'unknown'
     1021    CHARACTER (LEN=9), DIMENSION(300) ::  dopr_unit = 'unknown'
    10221022
    10231023    CHARACTER (LEN=7), DIMENSION(0:1,100) ::  do2d_unit, do3d_unit
  • palm/trunk/SOURCE/prognostic_equations.f90

    r484 r531  
    44! Current revisions:
    55! -----------------
    6 !
     6! add call of subsidence in the equation for humidity / passive scalar
    77!
    88! Former revisions:
     
    658658!--          Sink or source of scalar concentration due to canopy elements
    659659             IF ( plant_canopy ) CALL plant_canopy_model( i, j, 5 )
     660
     661!
     662!--          If required compute influence of large-scale subsidence/ascent
     663             IF ( large_scale_subsidence ) THEN
     664                CALL subsidence ( i, j, tend, q, q_init )
     665             ENDIF
    660666
    661667             CALL user_actions( i, j, 'q-tendency' )
     
    11101116          ENDIF
    11111117
    1112 
    11131118          CALL user_actions( i, j, 'pt-tendency' )
    11141119
     
    12231228!--          Sink or source of scalar concentration due to canopy elements
    12241229             IF ( plant_canopy ) CALL plant_canopy_model( i, j, 5 )
     1230
     1231
     1232!--          If required compute influence of large-scale subsidence/ascent
     1233             IF ( large_scale_subsidence ) THEN
     1234                CALL subsidence ( i, j, tend, q, q_init )
     1235             ENDIF
    12251236
    12261237
     
    19071918       IF ( plant_canopy ) CALL plant_canopy_model( 5 )
    19081919
     1920!
     1921!--    If required compute influence of large-scale subsidence/ascent
     1922       IF ( large_scale_subsidence ) THEN
     1923         CALL subsidence ( tend, q, q_init )
     1924       ENDIF
     1925
    19091926       CALL user_actions( 'q-tendency' )
    19101927
Note: See TracChangeset for help on using the changeset viewer.