- Timestamp:
- Jul 30, 2007 10:49:46 AM (17 years ago)
- Location:
- palm/trunk
- Files:
-
- 8 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/DOC/app/chapter_4.1.html
r97 r103 1880 1880 condition must be used for the potential temperature (see <a href="chapter_4.1.html#bc_pt_t">bc_pt_t</a>), 1881 1881 because 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> </p>1882 the top flux so that a constant flux value cannot be guaranteed.<span style="font-style: italic;"></span> </p> 1883 1883 <p><span style="font-weight: bold;">Note:</span><br>The 1884 1884 application of a top heat flux additionally requires the setting of … … 1887 1887 Prandtl-layer is available at the top boundary so far.</p><p>See 1888 1888 also <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 1890 two-dimensional u-momentum flux field <span style="font-family: monospace;">uswst</span> is 1891 initialized with the value of <span style="font-weight: bold;">top_momentumflux_u</span> as 1892 top (horizontally homogeneous) boundary condition for the u-momentum equation.</p><p><span style="font-weight: bold;">Notes:</span><br>The 1893 application of a top momentum flux additionally requires the setting of 1894 initial 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 Neumann 1896 condition should be used for the u velocity component (see <a href="chapter_4.1.html#bc_uv_t">bc_uv_t</a>), 1897 because otherwise the resolved scale may contribute to 1898 the top flux so that a constant flux value cannot be guaranteed.<span style="font-style: italic;"></span> </p> 1899 <span style="font-weight: bold;"></span><p>No 1900 Prandtl-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 1901 two-dimensional v-momentum flux field <span style="font-family: monospace;">vswst</span> is 1902 initialized with the value of <span style="font-weight: bold;">top_momentumflux_v</span> as 1903 top (horizontally homogeneous) boundary condition for the v-momentum equation.</p><p><span style="font-weight: bold;">Notes:</span><br>The 1904 application of a top momentum flux additionally requires the setting of 1905 initial 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 Neumann 1907 condition should be used for the v velocity component (see <a href="chapter_4.1.html#bc_uv_t">bc_uv_t</a>), 1908 because otherwise the resolved scale may contribute to 1909 the top flux so that a constant flux value cannot be guaranteed.<span style="font-style: italic;"></span> </p> 1910 <span style="font-weight: bold;"></span><p>No 1911 Prandtl-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> 1890 1912 salinityflux</span></td><td style="vertical-align: top;"><p>Kinematic 1891 1913 salinity flux at the top boundary, i.e. the sea surface (in psu m/s). </p> … … 1896 1918 condition must be used for the salinity (see <a href="chapter_4.1.html#bc_sa_t">bc_sa_t</a>), 1897 1919 because 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> </p>1920 the top flux so that a constant flux value cannot be guaranteed.<span style="font-style: italic;"></span> </p> 1899 1921 <p><span style="font-weight: bold;">Note:</span><br>The 1900 1922 application of a salinity flux at the model top additionally requires the setting of … … 2077 2099 = <i>.TRUE.</i>, the user-assigned top fluxes are used 2078 2100 instead 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 2081 2102 heat flux at the top will be automatically set to zero.</p></td></tr><tr> 2082 2103 <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 1020 1020 average'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 1021 1021 more 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> 1023 1023 </td> <td style="vertical-align: top;">R<br> </td> 1024 1024 <td style="vertical-align: top;"><span style="font-style: italic;">9999999.9</span><br> -
palm/trunk/DOC/app/chapter_4.6.html
r97 r103 367 367 interval of data which are subject to temporal averaging of <font face="Thorndale, serif"><font size="3">vertical 368 368 profiles and/or spectra (</font></font>in <font face="Thorndale, serif"><font size="3">s). </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> 370 370 <td>R</td> <td>R</td> <td><span style="font-style: italic;">9999999.9</span></td> 371 371 <td> <p lang="en-GB"><font face="Thorndale"><font face="Thorndale, serif">Temporal interval</font> … … 1371 1371 <td>I</td> <td>C * 40</td> <td><span style="font-style: italic;">'flat'</span></td> <td>Topography 1372 1372 mode.</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">Kinematic1373 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_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 1374 1374 salinity flux at the top boundary, i.e. the sea surface (in psu m/s).</td></tr><tr> 1375 1375 <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 1083 1083 IF ( surface_heatflux == 9999999.9 ) constant_heatflux = .FALSE. 1084 1084 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 1085 1097 1086 1098 ! -
palm/trunk/SOURCE/diffusion_u.f90
r102 r103 204 204 !-- Vertical diffusion at the first gridpoint below the top boundary, 205 205 !-- if the momentum flux at the top is prescribed by the user 206 IF ( use_top_fluxes ) THEN206 IF ( use_top_fluxes .AND. constant_top_momentumflux ) THEN 207 207 k = nzt 208 208 ! … … 368 368 !-- Vertical diffusion at the first gridpoint below the top boundary, 369 369 !-- if the momentum flux at the top is prescribed by the user 370 IF ( use_top_fluxes ) THEN370 IF ( use_top_fluxes .AND. constant_top_momentumflux ) THEN 371 371 k = nzt 372 372 ! -
palm/trunk/SOURCE/diffusion_v.f90
r102 r103 202 202 !-- Vertical diffusion at the first gridpoint below the top boundary, 203 203 !-- if the momentum flux at the top is prescribed by the user 204 IF ( use_top_fluxes ) THEN204 IF ( use_top_fluxes .AND. constant_top_momentumflux ) THEN 205 205 k = nzt 206 206 ! … … 367 367 !-- Vertical diffusion at the first gridpoint below the top boundary, 368 368 !-- if the momentum flux at the top is prescribed by the user 369 IF ( use_top_fluxes ) THEN369 IF ( use_top_fluxes .AND. constant_top_momentumflux ) THEN 370 370 k = nzt 371 371 ! -
palm/trunk/SOURCE/header.f90
r102 r103 5 5 ! ----------------- 6 6 ! Output of informations for coupled model runs (boundary conditions etc.) 7 ! + output of momentumfluxes at the top boundary 7 8 ! 8 9 ! Former revisions: … … 462 463 WRITE ( io, 304 ) 463 464 IF ( coupling_mode == 'uncoupled' ) THEN 465 WRITE ( io, 319 ) top_momentumflux_u, top_momentumflux_v 464 466 IF ( constant_top_heatflux ) THEN 465 467 WRITE ( io, 306 ) top_heatflux … … 1349 1351 318 FORMAT (/' outflow damping layer width: ',I3,' gridpoints with km_', & 1350 1352 'max =',F5.1,' m**2/s') 1353 319 FORMAT (' Predefined constant momentumflux: u: ',F9.6,' m**2/s**2'/ & 1354 ' v: ',F9.6,' m**2/s**2') 1351 1355 320 FORMAT (//' List output:'/ & 1352 1356 ' -----------'// & -
palm/trunk/SOURCE/modules.f90
r102 r103 5 5 ! Actual revisions: 6 6 ! ----------------- 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* 9 10 ! 10 11 ! Former revisions: … … 311 312 conserve_volume_flow = .FALSE., constant_diffusion = .FALSE., & 312 313 constant_heatflux = .TRUE., constant_top_heatflux = .TRUE., & 314 constant_top_momentumflux = .FALSE., & 313 315 constant_top_salinityflux = .TRUE., & 314 316 constant_waterflux = .TRUE., create_disturbances = .TRUE., & … … 396 398 time_do_sla = 0.0, time_dvrp = 0.0, time_prel = 0.0, & 397 399 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, & 400 402 ug_surface = 0.0, u_gtrans = 0.0, & 401 403 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.