Changeset 134
- Timestamp:
- Nov 21, 2007 7:28:38 AM (17 years ago)
- Location:
- palm/trunk
- Files:
-
- 4 edited
-
DOC/app/chapter_3.4.html (modified) (18 diffs)
-
DOC/app/chapter_4.1.html (modified) (5 diffs)
-
SOURCE/CURRENT_MODIFICATIONS (modified) (2 diffs)
-
SOURCE/init_grid.f90 (modified) (3 diffs)
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/DOC/app/chapter_3.4.html
r108 r134 1 1 <!DOCTYPE html PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN"> 2 <html> 3 <head> 2 <html><head> 4 3 5 4 … … 8 7 9 8 10 <title>PALM chapter 3.4</title> 11 </head> 12 9 <title>PALM chapter 3.4</title></head> 13 10 14 11 <body> … … 3423 3420 <font color="#000000"><span style="font-family: Cumberland AMT;"> 3424 3421 3425 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 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;"> 3429 0 0 0 0 0 3430 0 0 0 0 3431 0 0 0 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;"> 0 12.5 12.5 12.5 12.5 12.5 12.5 12.5 12.5 12.5 12.5 0 0 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;"> 0 12.5 12.5 12.5 12.5 12.5 12.5 12.5 12.5 12.5 12.5 0 0 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;"> 3434 0 12.5 12.5 50 50 50 50 50 3435 50 12.5 12.5 0 0 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 0 12.5 12.5 50 50 50 50 3438 50 50 12.5 12.5 0 0 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 0 12.5 12.5 50 50 50 50 3441 50 50 12.5 12.5 0 0 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;"> 0 12.5 12.5 50 50 50 50 50 50 12.5 12.5 0 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 0 12.5 12.5 50 50 50 50 3445 50 50 12.5 12.5 0 0 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 0 12.5 12.5 50 50 50 50 3448 50 50 12.5 12.5 0 0 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;"> 0 12.5 12.5 12.5 12.5 12.5 12.5 12.5 12.5 12.5 12.5 0 0 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;"> 0 12.5 12.5 12.5 12.5 12.5 12.5 12.5 12.5 12.5 12.5 0 0 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;"> 3452 0 0 0 0 0 3453 0 0 0 0 3454 0 0 0 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 3457 These 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> 3432 3458 3433 3459 … … 3436 3462 3437 3463 3438 <td><br> 3439 3440 3441 </td> 3464 3442 3465 3443 3466 … … 3482 3505 3483 3506 3484 <td style="vertical-align: top;"><br> 3485 3486 3487 </td> 3507 3488 3508 3489 3509 … … 3528 3548 3529 3549 3530 <td style="vertical-align: top;"><br> 3531 3532 3533 </td> 3550 3534 3551 3535 3552 … … 3574 3591 3575 3592 3576 <td style="vertical-align: top;"><br> 3577 3578 3579 </td> 3593 3580 3594 3581 3595 … … 3619 3633 3620 3634 3621 <td style="vertical-align: top;"><br> 3622 3623 3624 </td> 3635 3625 3636 3626 3637 … … 3664 3675 3665 3676 3666 <td style="vertical-align: top;"><br> 3667 3668 3669 </td> 3677 3670 3678 3671 3679 … … 3709 3717 3710 3718 3711 <td style="vertical-align: top;"><br> 3712 3713 3714 </td> 3719 3715 3720 3716 3721 … … 3754 3759 3755 3760 3756 <td style="vertical-align: top;"><br> 3757 3758 3759 </td> 3761 3760 3762 3761 3763 … … 3799 3801 3800 3802 3801 <td style="vertical-align: top;"><br> 3802 3803 3804 </td> 3803 3805 3804 3806 3805 … … 3861 3860 3862 3861 3863 <td align="left" valign="top"></td>3862 3864 3863 3865 3864 … … 3904 3903 3905 3904 3906 <td></td>3905 3907 3906 3908 3907 … … 3948 3947 3949 3948 3950 <td></td>3949 3951 3950 3952 3951 … … 3992 3991 3993 3992 3994 <td></td>3993 3995 3994 3996 3995 … … 4035 4034 4036 4035 4037 <td></td>4036 4038 4037 4039 4038 … … 4067 4066 4068 4067 4069 </body> 4070 </html> 4068 </body></html> -
palm/trunk/DOC/app/chapter_4.1.html
r132 r134 2024 2024 = <span style="font-style: italic;">( ( <a href="#nx">nx</a> + 2025 2025 1 ) * <a href="#dx">dx</a> - <a href="#building_length_x">building_length_x</a> ) / 2</span> 2026 centers the building in x-direction. </td>2026 centers the building in x-direction. <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> 2027 2027 2028 2028 … … 2090 2090 = <span style="font-style: italic;">( ( <a href="#ny">ny</a> + 2091 2091 1 ) * <a href="#dy">dy</a> - <a href="#building_length_y">building_length_y</a> ) / 2</span> 2092 centers the building in y-direction. </td>2092 centers the building in y-direction. <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> 2093 2093 2094 2094 … … 10537 10537 10538 10538 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>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>.<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> 10540 10540 10541 10541 … … 10574 10574 information</font></font><font color="#000000"> 10575 10575 in m. These data <span style="font-style: italic;"></span>must 10576 exactly match the horizontal grid.< /font> </div>10576 exactly 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> 10577 10577 10578 10578 … … 10586 10586 </span><font color="#000000"> 10587 10587 Alternatively, 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>10588 to 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"> 10589 10589 10590 10590 -
palm/trunk/SOURCE/CURRENT_MODIFICATIONS
r133 r134 21 21 ------- 22 22 23 Redefine initial nzb_local as the actual total size of topography (later the 24 extent of topography in nzb_local is reduced by 1dx at the E topography walls 25 and by 1dy at the N topography walls to form the basis for nzb_s_inner); 26 for consistency redefine 'single_building' case. 27 23 28 Vertical profiles now based on nzb_s_inner; they are divided by 24 29 ngp_2dh_s_inner (scalars, procucts of scalars) and ngp_2dh (staggered velocity … … 40 45 humidity/scalar/salinity in a future release. 41 46 42 buoyancy, check_open, data_output_dvrp, diffusion_s, diffusivities, flow_statistics, header, init_3d_model, init_dvrp, modules, prognostic_equations47 buoyancy, check_open, data_output_dvrp, diffusion_s, diffusivities, flow_statistics, header, init_3d_model, init_dvrp, init_grid, modules, prognostic_equations 43 48 44 49 -
palm/trunk/SOURCE/init_grid.f90
r116 r134 4 4 ! Actual revisions: 5 5 ! ----------------- 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. 6 10 ! Calculation of wall flag arrays 7 11 ! … … 365 369 366 370 ! 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. 374 374 nzb_local = 0 375 nzb_local(bys:byn -1,bxl:bxr-1) = bh375 nzb_local(bys:byn,bxl:bxr) = bh 376 376 377 377 CASE ( 'read_from_file' ) … … 468 468 ENDIF 469 469 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,:) 470 493 471 494 !
Note: See TracChangeset
for help on using the changeset viewer.
