Changeset 134


Ignore:
Timestamp:
Nov 21, 2007 7:28:38 AM (14 years ago)
Author:
letzel
Message:

Redefine initial nzb_local as the actual total size of topography (later the
extent of topography in nzb_local is reduced by 1dx at the E topography walls
and by 1dy at the N topography walls to form the basis for nzb_s_inner);
for consistency redefine 'single_building' case.

Location:
palm/trunk
Files:
4 edited

Legend:

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

    r108 r134  
    11<!DOCTYPE html PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
    2 <html>
    3 <head>
     2<html><head>
    43
    54
     
    87
    98 
    10   <title>PALM chapter 3.4</title>
    11 </head>
    12 
     9  <title>PALM chapter 3.4</title></head>
    1310
    1411<body>
     
    34233420      <font color="#000000"><span style="font-family: Cumberland AMT;">&nbsp; &nbsp;
    34243421&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;
    3425 &nbsp; &nbsp; &nbsp;&nbsp; &nbsp;S</span></font><br>
    3426 
    3427 
    3428       <span style="font-family: Arial;"></span><br>
    3429 
    3430 
    3431 These data must exactly match the horizontal grid.<br>
     3422&nbsp; &nbsp; &nbsp;&nbsp; &nbsp;S<br></span></font>
     3423
     3424
     3425      <br><span style="font-family: Arial;"></span><font color="#000000">Example for a 50m building surrounded by a 12.5m podium on flat ground:<br>
     3426
     3427
     3428      <span style="font-family: Cumberland AMT;"></span><font color="#000000"><font color="#000000"><span style="font-family: Cumberland AMT;">&nbsp;
     3429&nbsp;0 &nbsp; &nbsp;0 &nbsp; &nbsp;0 &nbsp; &nbsp;0 &nbsp; &nbsp;0
     3430&nbsp; &nbsp;0 &nbsp; &nbsp;0 &nbsp; &nbsp;0 &nbsp; &nbsp;0 &nbsp;
     3431&nbsp;0 &nbsp; &nbsp;0 &nbsp; &nbsp;0 &nbsp; &nbsp;0<br></span></font></font></font><font color="#000000"><font color="#000000"><font color="#000000"><font color="#000000"><font color="#000000"><span style="font-family: Cumberland AMT;"> &nbsp; &nbsp;0 12.5 12.5 12.5 12.5 12.5 12.5 12.5 12.5 12.5 12.5 &nbsp; &nbsp;0 &nbsp; &nbsp;0&nbsp;
     3432</span></font></font></font></font><font color="#000000"><font color="#000000"><font color="#000000"><span style="font-family: Cumberland AMT;"><br></span></font></font></font><font color="#000000"><font color="#000000"><font color="#000000"><font color="#000000"><span style="font-family: Cumberland AMT;"></span></font></font><font color="#000000"><font color="#000000"><span style="font-family: Cumberland AMT;"> </span></font></font></font></font></font><font color="#000000"><font color="#000000"><font color="#000000"><font color="#000000"><font color="#000000"><span style="font-family: Cumberland AMT;"> &nbsp; &nbsp;0 12.5 12.5 12.5 12.5 12.5 12.5 12.5 12.5 12.5 12.5 &nbsp; &nbsp;0 &nbsp; &nbsp;0&nbsp;
     3433</span></font></font></font></font><font color="#000000"><font color="#000000"><font color="#000000"><span style="font-family: Cumberland AMT;"><br></span></font></font></font><font color="#000000"><font color="#000000"><font color="#000000"><font color="#000000"><span style="font-family: Cumberland AMT;"></span></font></font><font color="#000000"><font color="#000000"><span style="font-family: Cumberland AMT;"> </span></font></font></font></font></font><font color="#000000"><font color="#000000"><font color="#000000"><font color="#000000"><span style="font-family: Cumberland AMT;">&nbsp;
     3434&nbsp;0 12.5 12.5 &nbsp; 50 &nbsp; 50 &nbsp; 50 &nbsp; 50 &nbsp; 50
     3435&nbsp; 50 12.5 12.5 &nbsp; &nbsp;0 &nbsp; &nbsp;0&nbsp;
     3436<br></span></font></font></font></font><font color="#000000"><font color="#000000"><font color="#000000"><font color="#000000"><span style="font-family: Cumberland AMT;"></span></font></font><font color="#000000"><font color="#000000"><span style="font-family: Cumberland AMT;">
     3437&nbsp; &nbsp;0 12.5 12.5 &nbsp; 50 &nbsp; 50 &nbsp; 50 &nbsp; 50 &nbsp;
     343850 &nbsp; 50 12.5 12.5 &nbsp; &nbsp;0 &nbsp; &nbsp;0&nbsp;
     3439<br></span></font></font></font></font><font color="#000000"><font color="#000000"><font color="#000000"><font color="#000000"><span style="font-family: Cumberland AMT;"></span></font></font><font color="#000000"><font color="#000000"><span style="font-family: Cumberland AMT;">
     3440&nbsp; &nbsp;0 12.5 12.5 &nbsp; 50 &nbsp; 50 &nbsp; 50 &nbsp; 50 &nbsp;
     344150 &nbsp; 50 12.5 12.5 &nbsp; &nbsp;0 &nbsp; &nbsp;0&nbsp;
     3442</span></font></font></font></font><br><font color="#000000"><font color="#000000"><font color="#000000"><span style="font-family: Cumberland AMT;"></span></font></font></font><font color="#000000"><font color="#000000"><font color="#000000"><font color="#000000"><span style="font-family: Cumberland AMT;">&nbsp; &nbsp;0 12.5 12.5 &nbsp; 50 &nbsp; 50 &nbsp; 50 &nbsp; 50 &nbsp; 50 &nbsp; 50 12.5 12.5 &nbsp; &nbsp;0 &nbsp; &nbsp;0 <br>
     3443</span></font></font></font><font color="#000000"><font color="#000000"><font color="#000000"><font color="#000000"><span style="font-family: Cumberland AMT;"></span></font></font><font color="#000000"><font color="#000000"><span style="font-family: Cumberland AMT;">
     3444&nbsp; &nbsp;0 12.5 12.5 &nbsp; 50 &nbsp; 50 &nbsp; 50 &nbsp; 50 &nbsp;
     344550 &nbsp; 50 12.5 12.5 &nbsp; &nbsp;0 &nbsp; &nbsp;0&nbsp;
     3446</span></font></font></font></font><font color="#000000"><font color="#000000"><font color="#000000"><span style="font-family: Cumberland AMT;"><br></span></font></font></font><font color="#000000"><font color="#000000"><font color="#000000"><font color="#000000"><span style="font-family: Cumberland AMT;"></span></font></font><font color="#000000"><font color="#000000"><span style="font-family: Cumberland AMT;">
     3447&nbsp; &nbsp;0 12.5 12.5 &nbsp; 50 &nbsp; 50 &nbsp; 50 &nbsp; 50 &nbsp;
     344850 &nbsp; 50 12.5 12.5 &nbsp; &nbsp;0 &nbsp; &nbsp;0&nbsp;
     3449<br></span></font></font></font></font><font color="#000000"><font color="#000000"><font color="#000000"><font color="#000000"><span style="font-family: Cumberland AMT;"></span></font></font><font color="#000000"><font color="#000000"><span style="font-family: Cumberland AMT;"></span></font></font></font></font></font><font color="#000000"><font color="#000000"><font color="#000000"><font color="#000000"><font color="#000000"><font color="#000000"><span style="font-family: Cumberland AMT;"> &nbsp; &nbsp;0 12.5 12.5 12.5 12.5 12.5 12.5 12.5 12.5 12.5 12.5 &nbsp; &nbsp;0 &nbsp; &nbsp;0&nbsp;
     3450</span></font></font></font></font><font color="#000000"><font color="#000000"><font color="#000000"><span style="font-family: Cumberland AMT;"><br></span></font></font></font><font color="#000000"><font color="#000000"><font color="#000000"><font color="#000000"><span style="font-family: Cumberland AMT;"></span></font></font><font color="#000000"><font color="#000000"><span style="font-family: Cumberland AMT;"> </span></font></font></font></font></font><font color="#000000"><font color="#000000"><font color="#000000"><font color="#000000"><font color="#000000"><span style="font-family: Cumberland AMT;"> &nbsp; &nbsp;0 12.5 12.5 12.5 12.5 12.5 12.5 12.5 12.5 12.5 12.5 &nbsp; &nbsp;0 &nbsp; &nbsp;0&nbsp;
     3451</span></font></font></font></font><font color="#000000"><font color="#000000"><font color="#000000"><span style="font-family: Cumberland AMT;"><br></span></font></font></font><font color="#000000"><font color="#000000"><font color="#000000"><font color="#000000"><span style="font-family: Cumberland AMT;"></span></font></font><font color="#000000"><font color="#000000"><span style="font-family: Cumberland AMT;"> </span></font></font></font></font></font><font color="#000000"><font color="#000000"><font color="#000000"><font color="#000000"><span style="font-family: Cumberland AMT;"></span></font></font></font></font></font><font color="#000000"><font color="#000000"><font color="#000000"><font color="#000000"><span style="font-family: Cumberland AMT;">&nbsp;
     3452&nbsp;0 &nbsp; &nbsp;0 &nbsp; &nbsp;0 &nbsp; &nbsp;0 &nbsp; &nbsp;0
     3453&nbsp; &nbsp;0 &nbsp; &nbsp;0 &nbsp; &nbsp;0 &nbsp; &nbsp;0 &nbsp;
     3454&nbsp;0 &nbsp; &nbsp;0 &nbsp; &nbsp;0 &nbsp; &nbsp;0<br></span></font></font></font><font color="#000000"><font color="#000000"><font color="#000000"><font color="#000000"><font color="#000000"><span style="font-family: Cumberland AMT;"></span></font></font></font></font></font></font><font color="#000000"><font color="#000000"><font color="#000000"><font color="#000000"><font color="#000000"><span style="font-family: Cumberland AMT;"></span></font></font></font></font></font><font color="#000000"><font color="#000000"><font color="#000000"><font color="#000000"><font color="#000000"><span style="font-family: Cumberland AMT;"></span></font></font></font></font><font color="#000000"><font color="#000000"><font color="#000000"><span style="font-family: Cumberland AMT;"></span></font></font></font></font><br>
     3455
     3456
     3457These data must exactly match the horizontal grid. <font color="#000000">Due to the staggered grid the topography will be displaced by -0.5 <a href="chapter_4.1.html#dx">dx</a> in x-direction and -0.5 <a href="chapter_4.1.html#dy">dy</a> in y-direction.</font><br>
    34323458
    34333459
     
    34363462
    34373463
    3438       <td><br>
    3439 
    3440 
    3441       </td>
     3464     
    34423465
    34433466
     
    34823505
    34833506
    3484       <td style="vertical-align: top;"><br>
    3485 
    3486 
    3487       </td>
     3507     
    34883508
    34893509
     
    35283548
    35293549
    3530       <td style="vertical-align: top;"><br>
    3531 
    3532 
    3533       </td>
     3550     
    35343551
    35353552
     
    35743591
    35753592
    3576       <td style="vertical-align: top;"><br>
    3577 
    3578 
    3579       </td>
     3593     
    35803594
    35813595
     
    36193633
    36203634
    3621       <td style="vertical-align: top;"><br>
    3622 
    3623 
    3624       </td>
     3635     
    36253636
    36263637
     
    36643675
    36653676
    3666       <td style="vertical-align: top;"><br>
    3667 
    3668 
    3669       </td>
     3677     
    36703678
    36713679
     
    37093717
    37103718
    3711       <td style="vertical-align: top;"><br>
    3712 
    3713 
    3714       </td>
     3719     
    37153720
    37163721
     
    37543759
    37553760
    3756       <td style="vertical-align: top;"><br>
    3757 
    3758 
    3759       </td>
     3761     
    37603762
    37613763
     
    37993801
    38003802
    3801       <td style="vertical-align: top;"><br>
    3802 
    3803 
    3804       </td>
     3803     
    38053804
    38063805
     
    38613860
    38623861
    3863       <td align="left" valign="top"></td>
     3862     
    38643863
    38653864
     
    39043903
    39053904
    3906       <td></td>
     3905     
    39073906
    39083907
     
    39483947
    39493948
    3950       <td></td>
     3949     
    39513950
    39523951
     
    39923991
    39933992
    3994       <td></td>
     3993     
    39953994
    39963995
     
    40354034
    40364035
    4037       <td></td>
     4036     
    40384037
    40394038
     
    40674066
    40684067
    4069 </body>
    4070 </html>
     4068</body></html>
  • palm/trunk/DOC/app/chapter_4.1.html

    r132 r134  
    20242024= <span style="font-style: italic;">( ( <a href="#nx">nx</a>&nbsp;+
    202520251 ) * <a href="#dx">dx</a> -&nbsp; <a href="#building_length_x">building_length_x</a> ) / 2</span>
    2026 centers the building in x-direction. </td>
     2026centers the building in x-direction.&nbsp;<font color="#000000">Due to the staggered grid the building will be displaced by -0.5 <a href="chapter_4.1.html#dx">dx</a> in x-direction and -0.5 <a href="chapter_4.1.html#dy">dy</a> in y-direction.</font> </td>
    20272027
    20282028
     
    20902090= <span style="font-style: italic;">( ( <a href="#ny">ny</a>&nbsp;+
    209120911 ) * <a href="#dy">dy</a> -&nbsp; <a href="#building_length_y">building_length_y</a> ) / 2</span>
    2092 centers the building in y-direction. </td>
     2092centers the building in y-direction.&nbsp;<font color="#000000">Due to the staggered grid the building will be displaced by -0.5 <a href="chapter_4.1.html#dx">dx</a> in x-direction and -0.5 <a href="chapter_4.1.html#dy">dy</a> in y-direction.</font> </td>
    20932093
    20942094
     
    1053710537
    1053810538
    10539 The building size and location can be specified with the parameters <a href="#building_height">building_height</a>, <a href="#building_length_x">building_length_x</a>, <a href="#building_length_y">building_length_y</a>, <a href="#building_wall_left">building_wall_left</a> and <a href="#building_wall_south">building_wall_south</a>.</div>
     10539The building size and location can be specified with the parameters <a href="#building_height">building_height</a>, <a href="#building_length_x">building_length_x</a>, <a href="#building_length_y">building_length_y</a>, <a href="#building_wall_left">building_wall_left</a> and <a href="#building_wall_south">building_wall_south</a>.<br><font color="#000000">Due to the staggered grid the building will be displaced by -0.5 <a href="chapter_4.1.html#dx">dx</a> in x-direction and -0.5 <a href="chapter_4.1.html#dy">dy</a> in y-direction.</font></div>
    1054010540
    1054110541
     
    1057410574information</font></font><font color="#000000">
    1057510575in m. These data&nbsp;<span style="font-style: italic;"></span>must
    10576 exactly match the horizontal grid.</font> </div>
     10576exactly match the horizontal grid.<br>Due to the staggered grid the topography will be displaced by -0.5 <a href="#dx">dx</a> in x-direction and -0.5 <a href="#dy">dy</a> in y-direction.</font> </div>
    1057710577
    1057810578
     
    1058610586 </span><font color="#000000">
    1058710587Alternatively, the user may add code to the user interface subroutine <a href="chapter_3.5.1.html#user_init_grid">user_init_grid</a>
    10588 to allow further topography modes.<br>
     10588to allow further topography modes. Again, </font><font color="#000000">due to the staggered grid the topography will be displaced by -0.5 <a href="chapter_4.1.html#dx">dx</a> in x-direction and -0.5 <a href="chapter_4.1.html#dy">dy</a> in y-direction.</font><br><font color="#000000">
    1058910589
    1059010590
  • palm/trunk/SOURCE/CURRENT_MODIFICATIONS

    r133 r134  
    2121-------
    2222
     23Redefine initial nzb_local as the actual total size of topography (later the
     24extent of topography in nzb_local is reduced by 1dx at the E topography walls
     25and by 1dy at the N topography walls to form the basis for nzb_s_inner);
     26for consistency redefine 'single_building' case.
     27
    2328Vertical profiles now based on nzb_s_inner; they are divided by
    2429ngp_2dh_s_inner (scalars, procucts of scalars) and ngp_2dh (staggered velocity
     
    4045humidity/scalar/salinity in a future release.
    4146
    42 buoyancy, check_open, data_output_dvrp, diffusion_s, diffusivities, flow_statistics, header, init_3d_model, init_dvrp, modules, prognostic_equations
     47buoyancy, check_open, data_output_dvrp, diffusion_s, diffusivities, flow_statistics, header, init_3d_model, init_dvrp, init_grid, modules, prognostic_equations
    4348
    4449
  • palm/trunk/SOURCE/init_grid.f90

    r116 r134  
    44! Actual revisions:
    55! -----------------
     6! Redefine initial nzb_local as the actual total size of topography (later the
     7! extent of topography in nzb_local is reduced by 1dx at the E topography walls
     8! and by 1dy at the N topography walls to form the basis for nzb_s_inner);
     9! for consistency redefine 'single_building' case.
    610! Calculation of wall flag arrays
    711!
     
    365369
    366370!
    367 !--       Set the individual index arrays for all velocity components and
    368 !--       scalars, taking into account the staggered grid. The horizontal
    369 !--       wind component normal to a wall defines the position of the wall, and
    370 !--       in the respective direction the building is as long as specified in
    371 !--       building_length_?, but in the other horizontal direction (for w and s
    372 !--       in both horizontal directions) the building appears shortened by one
    373 !--       grid length due to the staggered grid.
     371!--       Set the actual total size of the building. Due to the staggered grid,
     372!--       the building will be displaced by -0.5dx in x-direction and by -0.5dy
     373!--       in y-direction compared to the scalar grid.
    374374          nzb_local = 0
    375           nzb_local(bys:byn-1,bxl:bxr-1) = bh
     375          nzb_local(bys:byn,bxl:bxr) = bh
    376376
    377377       CASE ( 'read_from_file' )
     
    468468          ENDIF
    469469       ENDIF
     470
     471!
     472!--    The array nzb_local as defined above describes the actual total size of
     473!--    topography which is defined by u=0 on the topography walls in x-direction
     474!--    and by v=0 on the topography walls in y-direction. However, PALM uses
     475!--    individual arrays nzb_u|v|w|s_inner|outer that are based on nzb_s_inner.
     476!--    Therefore, the extent of topography in nzb_local is now reduced by 1dx
     477!--    at the E topography walls and by 1dy at the N topography walls to form
     478!--    the basis for nzb_s_inner.
     479       DO  j = -gls, ny + gls
     480          DO  i = -gls, nx
     481             nzb_local(j,i) = MIN( nzb_local(j,i), nzb_local(j,i+1) )
     482          ENDDO
     483       ENDDO
     484!--    apply cyclic boundary conditions in x-direction
     485       nzb_local(:,nx+1:nx+gls) = nzb_local(:,0:gls-1)
     486       DO  i = -gls, nx + gls
     487          DO  j = -gls, ny
     488             nzb_local(j,i) = MIN( nzb_local(j,i), nzb_local(j+1,i) )
     489          ENDDO
     490       ENDDO
     491!--    apply cyclic boundary conditions in y-direction
     492       nzb_local(ny+1:ny+gls,:) = nzb_local(0:gls-1,:)
    470493
    471494!
Note: See TracChangeset for help on using the changeset viewer.