Changeset 103 for palm


Ignore:
Timestamp:
Jul 30, 2007 10:49:46 AM (17 years ago)
Author:
raasch
Message:

further preliminary changes concerning coupling

Location:
palm/trunk
Files:
8 edited

Legend:

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

    r97 r103  
    18801880condition must be used for the potential temperature (see <a href="chapter_4.1.html#bc_pt_t">bc_pt_t</a>),
    18811881because otherwise the resolved scale may contribute to
    1882 the top flux so that a constant value cannot be guaranteed.<span style="font-style: italic;"></span>&nbsp;</p>
     1882the top flux so that a constant flux value cannot be guaranteed.<span style="font-style: italic;"></span>&nbsp;</p>
    18831883<p><span style="font-weight: bold;">Note:</span><br>The
    18841884application of a top heat flux additionally requires the setting of
     
    18871887Prandtl-layer is available at the top boundary so far.</p><p>See
    18881888also <a href="#surface_heatflux">surface_heatflux</a>.</p>
    1889 </td></tr><tr><td style="vertical-align: top;"><a name="top_salinityflux"></a><span style="font-weight: bold;">top_salinityflux</span></td><td style="vertical-align: top;">R</td><td style="vertical-align: top;"><span style="font-style: italic;">no prescribed<br>
     1889</td></tr><tr><td style="vertical-align: top;"><a name="top_momentumflux_u"></a><span style="font-weight: bold;">top_momentumflux_u</span></td><td style="vertical-align: top;">R</td><td style="vertical-align: top;"><span style="font-style: italic;">no prescribed momentumflux</span></td><td style="vertical-align: top;">Momentum flux along x at the top boundary (in m2/s2).<br><p>If a value is assigned to this parameter, the internal
     1890two-dimensional u-momentum flux field <span style="font-family: monospace;">uswst</span> is
     1891initialized with the value of <span style="font-weight: bold;">top_momentumflux_u</span> as
     1892top (horizontally homogeneous) boundary condition for the u-momentum equation.</p><p><span style="font-weight: bold;">Notes:</span><br>The
     1893application of a top momentum flux additionally requires the setting of
     1894initial parameter <a href="chapter_4.1.html#use_top_fluxes">use_top_fluxes</a>
     1895= .T.. Setting of <span style="font-weight: bold;">top_momentumflux_u</span> requires setting of <a href="#top_momentumflux_v">top_momentumflux_v</a> also.</p><p>A&nbsp;Neumann
     1896condition should be used for the u velocity component (see <a href="chapter_4.1.html#bc_uv_t">bc_uv_t</a>),
     1897because otherwise the resolved scale may contribute to
     1898the top flux so that a constant flux value cannot be guaranteed.<span style="font-style: italic;"></span>&nbsp;</p>
     1899<span style="font-weight: bold;"></span><p>No
     1900Prandtl-layer is available at the top boundary so far.</p></td></tr><tr><td style="vertical-align: top;"><a name="top_momentumflux_v"></a><span style="font-weight: bold;">top_momentumflux_v</span></td><td style="vertical-align: top;">R</td><td style="vertical-align: top;"><span style="font-style: italic;">no prescribed momentumflux</span></td><td style="vertical-align: top;">Momentum flux along y at the top boundary (in m2/s2).<br><p>If a value is assigned to this parameter, the internal
     1901two-dimensional v-momentum flux field <span style="font-family: monospace;">vswst</span> is
     1902initialized with the value of <span style="font-weight: bold;">top_momentumflux_v</span> as
     1903top (horizontally homogeneous) boundary condition for the v-momentum equation.</p><p><span style="font-weight: bold;">Notes:</span><br>The
     1904application of a top momentum flux additionally requires the setting of
     1905initial parameter <a href="chapter_4.1.html#use_top_fluxes">use_top_fluxes</a>
     1906= .T.. Setting of <span style="font-weight: bold;">top_momentumflux_v</span> requires setting of <a href="chapter_4.1.html#top_momentumflux_u">top_momentumflux_u</a> also.</p><p>A&nbsp;Neumann
     1907condition should be used for the v velocity component (see <a href="chapter_4.1.html#bc_uv_t">bc_uv_t</a>),
     1908because otherwise the resolved scale may contribute to
     1909the top flux so that a constant flux value cannot be guaranteed.<span style="font-style: italic;"></span>&nbsp;</p>
     1910<span style="font-weight: bold;"></span><p>No
     1911Prandtl-layer is available at the top boundary so far.</p></td></tr><tr><td style="vertical-align: top;"><a name="top_salinityflux"></a><span style="font-weight: bold;">top_salinityflux</span></td><td style="vertical-align: top;">R</td><td style="vertical-align: top;"><span style="font-style: italic;">no prescribed<br>
    18901912salinityflux</span></td><td style="vertical-align: top;"><p>Kinematic
    18911913salinity flux at the top boundary, i.e. the sea surface (in psu m/s).&nbsp; </p>
     
    18961918condition must be used for the salinity (see <a href="chapter_4.1.html#bc_sa_t">bc_sa_t</a>),
    18971919because otherwise the resolved scale may contribute to
    1898 the top flux so that a constant value cannot be guaranteed.<span style="font-style: italic;"></span>&nbsp;</p>
     1920the top flux so that a constant flux value cannot be guaranteed.<span style="font-style: italic;"></span>&nbsp;</p>
    18991921<p><span style="font-weight: bold;">Note:</span><br>The
    19001922application of a salinity flux at the model top additionally requires the setting of
     
    20772099= <i>.TRUE.</i>, the user-assigned top fluxes are used
    20782100instead
    2079 (see <a href="chapter_4.1.html#top_heatflux">top_heatflux</a>).</p><p>Currently,
    2080 only a value for the sensible heatflux can be assigned. In case of <span style="font-weight: bold;">use_top_fluxes</span> = <span style="font-style: italic;">.TRUE.</span>, the latent
     2101(see <a href="chapter_4.1.html#top_heatflux">top_heatflux</a>, <a href="#top_momentumflux_u">top_momentumflux_u</a>, <a href="#top_momentumflux_v">top_momentumflux_v</a>, <a href="#top_salinityflux">top_salinityflux</a>).</p><p>Currently, no value for the latent heatflux can be assigned. In case of <span style="font-weight: bold;">use_top_fluxes</span> = <span style="font-style: italic;">.TRUE.</span>, the latent
    20812102heat flux at the top will be automatically set to zero.</p></td></tr><tr>
    20822103<td style="vertical-align: top;"> <p><a name="use_ug_for_galilei_tr"></a><b>use_ug_for_galilei_tr</b></p>
  • palm/trunk/DOC/app/chapter_4.2.html

    r102 r103  
    10201020average's statistics. <span lang="en-GB"><font face="Thorndale, serif"><span style="font-weight: bold;"></span><span style="font-weight: bold;"></span></font></span><a href="chapter_4.1.html#dt"><span lang="en-GB"></span></a><font face="Thorndale, serif"><span lang="en-GB"></span></font><span lang="en-GB"></span><br> </p><p>For
    10211021more explanations see parameter <a href="#dt_averaging_input">dt_averaging_input</a>.<a href="chapter_4.1.html#dt"><span lang="en-GB"></span></a><font face="Thorndale, serif"><span lang="en-GB"></span></font></p></td>
    1022 </tr> <tr> <td style="vertical-align: top;"><a name="dt_data_output"></a><span style="font-weight: bold;">dt_data_output</span><br>
     1022</tr> <tr><td style="vertical-align: top;"><a name="dt_coupling"></a><span style="font-weight: bold;">dt_coupling</span></td><td style="vertical-align: top;">R</td><td style="vertical-align: top;"><span style="font-style: italic;">9999999.9</span></td><td style="vertical-align: top;">Temporal interval for the data exchange in case of runs with coupled models (e.g. atmosphere - ocean) (in s).<br><br>This parameter has an effect only in case of a run with coupled models.<br>It is available starting from version 3.3a. A more detailed explanation will be given soon.</td></tr><tr> <td style="vertical-align: top;"><a name="dt_data_output"></a><span style="font-weight: bold;">dt_data_output</span><br>
    10231023</td> <td style="vertical-align: top;">R<br> </td>
    10241024<td style="vertical-align: top;"><span style="font-style: italic;">9999999.9</span><br>
  • palm/trunk/DOC/app/chapter_4.6.html

    r97 r103  
    367367interval of&nbsp;data which are subject to temporal averaging of <font face="Thorndale, serif"><font size="3">vertical
    368368profiles and/or spectra&nbsp;(</font></font>in <font face="Thorndale, serif"><font size="3">s).&nbsp;</font></font></td>
    369 </tr> <tr> <td><a href="chapter_4.2.html#dt_data_output"><span style="font-weight: bold;">dt_data_output</span></a></td>
     369</tr> <tr><td align="undefined" valign="undefined"><a style="font-weight: bold;" href="chapter_4.2.html#dt_coupling">dt_coupling</a></td><td align="undefined" valign="undefined">R</td><td align="undefined" valign="undefined">R</td><td align="undefined" valign="undefined">9999999.9</td><td align="undefined" valign="undefined">Temporal interval for the data exchange in case of runs with coupled models (e.g. atmosphere - ocean) (in s).</td></tr><tr> <td><a href="chapter_4.2.html#dt_data_output"><span style="font-weight: bold;">dt_data_output</span></a></td>
    370370<td>R</td> <td>R</td> <td><span style="font-style: italic;">9999999.9</span></td>
    371371<td> <p lang="en-GB"><font face="Thorndale"><font face="Thorndale, serif">Temporal interval</font>
     
    13711371<td>I</td> <td>C * 40</td> <td><span style="font-style: italic;">'flat'</span></td> <td>Topography
    13721372mode.</td> </tr> <tr><td><a style="font-weight: bold;" href="chapter_4.1.html#top_heatflux">top_heatflux</a></td><td>I</td><td>R</td><td><span style="font-style: italic;">no prescribed heatflux</span></td><td>Kinematic
    1373 sensible heat flux at the top surface (in K m/s).</td></tr><tr><td align="undefined" valign="undefined"><a href="chapter_4.1.html#top_salinityflux"><span style="font-weight: bold;">top_salinityflux</span></a></td><td align="undefined" valign="undefined">I</td><td align="undefined" valign="undefined">R</td><td align="undefined" valign="undefined"><span style="font-style: italic;">no prescribed</span><br style="font-style: italic;"><span style="font-style: italic;">salinityflux</span></td><td align="undefined" valign="undefined">Kinematic
     1373sensible heat flux at the top surface (in K m/s).</td></tr><tr><td align="undefined" valign="undefined"><a href="chapter_4.1.html#top_momentumflux_u"><span style="font-weight: bold;">top_momentumflux_u</span></a></td><td align="undefined" valign="undefined">I</td><td align="undefined" valign="undefined">R</td><td align="undefined" valign="undefined"><span style="font-style: italic;">no prescribed momentumflux</span></td><td align="undefined" valign="undefined">Momentum flux along x at the top boundary (in m2/s2).</td></tr><tr><td align="undefined" valign="undefined"><a style="font-weight: bold;" href="chapter_4.1.html#top_momentumflux_v">top_momentumflux_v</a></td><td align="undefined" valign="undefined">I</td><td align="undefined" valign="undefined">R</td><td align="undefined" valign="undefined"><span style="font-style: italic;">no prescribed momentumflux</span></td><td align="undefined" valign="undefined">Momentum flux along y at the top boundary (in m2/s2).</td></tr><tr><td align="undefined" valign="undefined"><a href="chapter_4.1.html#top_salinityflux"><span style="font-weight: bold;">top_salinityflux</span></a></td><td align="undefined" valign="undefined">I</td><td align="undefined" valign="undefined">R</td><td align="undefined" valign="undefined"><span style="font-style: italic;">no prescribed</span><br style="font-style: italic;"><span style="font-style: italic;">salinityflux</span></td><td align="undefined" valign="undefined">Kinematic
    13741374salinity flux at the top boundary, i.e. the sea surface (in psu m/s).</td></tr><tr>
    13751375<td style="vertical-align: middle;" width="15%"> <p><a href="chapter_4.1.html#ug_surface"><b>ug_surface</b></a></p>
  • palm/trunk/SOURCE/check_parameters.f90

    r102 r103  
    10831083    IF ( surface_heatflux == 9999999.9 )  constant_heatflux     = .FALSE.
    10841084    IF ( top_heatflux     == 9999999.9 )  constant_top_heatflux = .FALSE.
     1085    IF ( top_momentumflux_u /= 9999999.9  .AND.  &
     1086         top_momentumflux_v /= 9999999.9 )  THEN
     1087       constant_top_momentumflux = .TRUE.
     1088    ELSEIF (  .NOT. ( top_momentumflux_u == 9999999.9  .AND.  &
     1089           top_momentumflux_v == 9999999.9 ) )  THEN   
     1090       IF ( myid == 0 )  THEN
     1091          PRINT*, '+++ check_parameters:'
     1092          PRINT*, '    both, top_momentumflux_u AND top_momentumflux_v'
     1093          PRINT*, '    must be set'
     1094       ENDIF
     1095       CALL local_stop
     1096    ENDIF
    10851097
    10861098!
  • palm/trunk/SOURCE/diffusion_u.f90

    r102 r103  
    204204!--          Vertical diffusion at the first gridpoint below the top boundary,
    205205!--          if the momentum flux at the top is prescribed by the user
    206              IF ( use_top_fluxes )  THEN
     206             IF ( use_top_fluxes  .AND.  constant_top_momentumflux )  THEN
    207207                k = nzt
    208208!
     
    368368!--    Vertical diffusion at the first gridpoint below the top boundary,
    369369!--    if the momentum flux at the top is prescribed by the user
    370        IF ( use_top_fluxes )  THEN
     370       IF ( use_top_fluxes  .AND.  constant_top_momentumflux )  THEN
    371371          k = nzt
    372372!
  • palm/trunk/SOURCE/diffusion_v.f90

    r102 r103  
    202202!--          Vertical diffusion at the first gridpoint below the top boundary,
    203203!--          if the momentum flux at the top is prescribed by the user
    204              IF ( use_top_fluxes )  THEN
     204             IF ( use_top_fluxes  .AND.  constant_top_momentumflux )  THEN
    205205                k = nzt
    206206!
     
    367367!--    Vertical diffusion at the first gridpoint below the top boundary,
    368368!--    if the momentum flux at the top is prescribed by the user
    369        IF ( use_top_fluxes )  THEN
     369       IF ( use_top_fluxes  .AND.  constant_top_momentumflux )  THEN
    370370          k = nzt
    371371!
  • palm/trunk/SOURCE/header.f90

    r102 r103  
    55! -----------------
    66! Output of informations for coupled model runs (boundary conditions etc.)
     7! + output of momentumfluxes at the top boundary
    78!
    89! Former revisions:
     
    462463       WRITE ( io, 304 )
    463464       IF ( coupling_mode == 'uncoupled' )  THEN
     465          WRITE ( io, 319 )  top_momentumflux_u, top_momentumflux_v
    464466          IF ( constant_top_heatflux )  THEN
    465467             WRITE ( io, 306 )  top_heatflux
     
    13491351318 FORMAT (/'       outflow damping layer width: ',I3,' gridpoints with km_', &
    13501352                    'max =',F5.1,' m**2/s')
     1353319 FORMAT ('       Predefined constant momentumflux:  u: ',F9.6,' m**2/s**2'/ &
     1354            '                                          v: ',F9.6,' m**2/s**2')
    13511355320 FORMAT (//' List output:'/ &
    13521356             ' -----------'//  &
  • palm/trunk/SOURCE/modules.f90

    r102 r103  
    55! Actual revisions:
    66! -----------------
    7 ! +comm_inter, coupling_char, coupling_mode, dt_coupling, ngp_xy, port_name,
    8 ! time_coupling, type_xy, uswst*, vswst*
     7! +comm_inter, constant_top_momentumflux, coupling_char, coupling_mode,
     8! dt_coupling, ngp_xy, port_name, time_coupling, top_momentumflux_u|v,
     9! type_xy, uswst*, vswst*
    910!
    1011! Former revisions:
     
    311312                conserve_volume_flow = .FALSE., constant_diffusion = .FALSE., &
    312313                constant_heatflux = .TRUE., constant_top_heatflux = .TRUE., &
     314                constant_top_momentumflux = .FALSE., &
    313315                constant_top_salinityflux = .TRUE., &
    314316                constant_waterflux = .TRUE., create_disturbances = .TRUE., &
     
    396398             time_do_sla = 0.0, time_dvrp = 0.0, time_prel = 0.0, &
    397399             time_restart = 9999999.9, time_run_control = 0.0, &
    398              top_heatflux = 9999999.9, top_momentumflux_u = 0.0, &
    399              top_momentumflux_v = 0.0, top_salinityflux = 9999999.9, &
     400             top_heatflux = 9999999.9, top_momentumflux_u = 9999999.9, &
     401             top_momentumflux_v = 9999999.9, top_salinityflux = 9999999.9, &
    400402             ug_surface = 0.0, u_gtrans = 0.0, &
    401403             ups_limit_e = 0.0, ups_limit_pt = 0.0, ups_limit_u = 0.0, &
Note: See TracChangeset for help on using the changeset viewer.