Changeset 405 for palm/trunk
- Timestamp:
- Nov 16, 2009 12:20:59 PM (15 years ago)
- Location:
- palm/trunk
- Files:
-
- 4 edited
-
DOC/app/chapter_3.5.4.html (modified) (4 diffs)
-
SCRIPTS/mrun (modified) (2 diffs)
-
SOURCE/CURRENT_MODIFICATIONS (modified) (2 diffs)
-
SOURCE/modules.f90 (modified) (2 diffs)
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/DOC/app/chapter_3.5.4.html
r343 r405 1 1 <!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 Transitional//EN"> 2 2 <html><head> 3 3 4 <meta http-equiv="CONTENT-TYPE" content="text/html; charset=windows-1252"><title>PALM chapter 3.5.4</title> <meta name="GENERATOR" content="StarOffice 7 (Win32)"> <meta name="AUTHOR" content="Siegfried Raasch"> <meta name="CREATED" content="20040802;14050943"> <meta name="CHANGED" content="20041117;12180008"> <meta name="KEYWORDS" content="parallel LES model"> <style> 4 5 <!-- 5 6 @page { size: 21cm 29.7cm } 6 7 --> 7 </style></head> 8 <body style="direction: ltr;" lang="en-US"><h3 style="line-height: 100%;"><font size="4">3.5.4 8 </style></head><body style="direction: ltr;" lang="en-US"><h3 style="line-height: 100%;"><font size="4">3.5.4 9 9 User-defined output quantities<br> 10 10 </font></h3>A very typical request of users is the … … 172 172 IF ( initializing_actions == 'read_restart_data' ) 173 173 THEN</span><br style="font-family: monospace;"><span style="font-family: monospace;"> 174 READ ( 13 ) field_ch r</span><br style="font-family: monospace;"><span style="font-family: monospace;"> 175 DO WHILE ( TRIM( field_ch r ) /= '*** end174 READ ( 13 ) field_char</span><br style="font-family: monospace;"><span style="font-family: monospace;"> 175 DO WHILE ( TRIM( field_char ) /= '*** end 176 176 user ***' )</span><br style="font-family: monospace;"><br style="font-family: monospace;"><span style="font-family: monospace;"> 177 SELECT CASE ( TRIM( field_ch r ) )</span><br style="font-family: monospace;"><br style="font-family: monospace;"><span style="font-family: monospace;"> 177 SELECT CASE ( TRIM( field_char ) )</span><br style="font-family: monospace;"><br style="font-family: monospace;"><span style="font-family: monospace;"> 178 178 CASE ( 'u2_av' )</span><br style="font-family: monospace;"><span style="font-family: monospace;"> 179 179 ALLOCATE( … … 186 186 187 187 188 TRIM( 189 field_chr ), '" found in'</span><br style="font-family: monospace;"><span style="font-family: monospace;"> 188 TRIM( field_char ), '" found in'</span><br style="font-family: monospace;"><span style="font-family: monospace;"> 190 189 PRINT*, 191 190 ' … … 276 275 grid = 'zu'</span><br> <br>The <span style="font-family: monospace;">ELSE</span> case is 277 276 only needed in case that output of time-averaged data is requested.<br><br> 278 As a special case, xy cross section output can also be defined at one single level at height <span style="font-family: monospace;">k=nzb+1</span> on the u-grid. This features is useful for output of surface data (e.g. heat fluxes). In this case, the corresponding 2d data has to be resorted to the array <span style="font-family: monospace;">local_pf(i,j,nzb+1).</span> In addition to this, the grid in <a href="chapter_3.5.1.html#user_define_netcdf_grid"><span style="font-family: monospace;">user_define_netcdf_grid</span></a> as well as in <a href="chapter_3.5.1.html#user_data_output_2d"><span style="font-family: monospace;">user_data_output_2d</span></a> must be set to <span style="font-family: monospace;">grid = 'zu1'</span>.<br><br> 277 As a special case, xy cross section output can also be defined at one single level at height <span style="font-family: monospace;">k=nzb+1</span> 278 on the u-grid. This features is useful for output of surface data (e.g. 279 heat fluxes). In this case, the corresponding 2d data has to be 280 resorted to the array <span style="font-family: monospace;">local_pf(i,j,nzb+1).</span> In addition to this, the grid in <a href="chapter_3.5.1.html#user_define_netcdf_grid"><span style="font-family: monospace;">user_define_netcdf_grid</span></a> as well as in <a href="chapter_3.5.1.html#user_data_output_2d"><span style="font-family: monospace;">user_data_output_2d</span></a> must be set to <span style="font-family: monospace;">grid = 'zu1'</span>.<br><br> 279 281 <span style="font-family: monospace;"> CASE 280 282 ( 'u2_xy' )</span><br style="font-family: monospace;"><span style="font-family: monospace;"> -
palm/trunk/SCRIPTS/mrun
r400 r405 3162 3162 if [[ $threads_per_task != 1 ]] 3163 3163 then 3164 printf "\n threads per task: $threads_per_task" 3164 # increase stack size to unlimited, because large runs 3165 # may abort otherwise 3166 ulimit -Ss unlimited 3167 printf "\n threads per task: $threads_per_task stacksize: unlimited" 3165 3168 fi 3166 3169 if [[ $run_coupled_model = false ]] … … 3197 3200 export MV2_NUM_PORTS=2 3198 3201 # The default setting of MV2_CPU_MAPPING gives best results 3202 # export MV2_ENABLE_AFFINITY=1 3199 3203 # export MV2_CPU_MAPPING=0,1,4,5,2,3,6,7 3200 mpiexec -np $ii ./a.out $ROPTS < runfile_atmos 3204 if [[ $threads_per_task != 1 ]] 3205 then 3206 mpiexec -npernode 1 ./a.out $ROPTS < runfile_atmos 3207 else 3208 mpiexec -np $ii ./a.out $ROPTS < runfile_atmos 3209 fi 3201 3210 fi 3202 3211 elif [[ $host = lcxt4 ]] -
palm/trunk/SOURCE/CURRENT_MODIFICATIONS
r403 r405 7 7 var_ts is replaced by dots_max (modules,init_3d_model) 8 8 9 modules,init_3d_model 9 init_3d_model, modules 10 10 11 11 … … 18 18 Bugfix: calculation of precipitation_rate (calc_precipitation) 19 19 20 sum_up_3d_data, time_integration, calc_precipitation 20 Bugfix: initial data assignments to some dvrp arrays changed due to error 21 messages from gfortran compiler (modules= 22 23 calc_precipitation, modules, sum_up_3d_data, time_integration 21 24 22 25 -
palm/trunk/SOURCE/modules.f90
r394 r405 5 5 ! Current revisions: 6 6 ! ----------------- 7 ! -var_ts: replaced by dots_max 8 ! 7 ! -var_ts: replaced by dots_max, 8 ! initial data assignments to some dvrp arrays changed due to error messages 9 ! from gfortran compiler 9 10 ! 10 11 ! Former revisions: … … 583 584 topography_color = (/ 0.8, 0.7, 0.6 /) 584 585 585 #if defined( __decalpha )586 REAL, DIMENSION(2,10) :: slicer_range_limits_dvrp = (/&586 !#if defined( __decalpha ) 587 REAL, DIMENSION(2,10) :: slicer_range_limits_dvrp = RESHAPE( (/ & 587 588 -1.0, 1.0, -1.0, 1.0, -1.0, 1.0, -1.0, 1.0, & 588 589 -1.0, 1.0, -1.0, 1.0, -1.0, 1.0, -1.0, 1.0, & 589 -1.0, 1.0, -1.0, 1.0 /) 590 591 REAL, DIMENSION(3,10) :: isosurface_color = (/&590 -1.0, 1.0, -1.0, 1.0 /), (/ 2, 10 /) ) 591 592 REAL, DIMENSION(3,10) :: isosurface_color = RESHAPE( (/ & 592 593 0.9, 0.9, 0.9, 0.8, 0.1, 0.1, 0.1, 0.1, 0.8, & 593 594 0.1, 0.8, 0.1, 0.6, 0.1, 0.1, 0.1, 0.1, 0.6, & 594 0.1, 0.6, 0.1, 0.4, 0.1, 0.1, 0.1, 0.1, 0.4 595 0.1, 0.4, 0.1 /) 595 0.1, 0.6, 0.1, 0.4, 0.1, 0.1, 0.1, 0.1, 0.4, & 596 0.1, 0.4, 0.1 /), (/ 3, 10 /) ) 596 597 597 598 REAL(4), DIMENSION(2,100) :: interval_values_dvrp, interval_h_dvrp = & 598 (/ 270.0, 225.0, 225.0, 180.0, 70.0, 25.0, & 599 25.0, -25.0, ( 0.0, i9 = 1, 192 ) /), & 599 RESHAPE( (/ 270.0, 225.0, 225.0, 180.0, & 600 70.0, 25.0, 25.0, -25.0, & 601 ( 0.0, i9 = 1, 192 ) /), & 602 (/ 2, 100 /) ), & 600 603 interval_l_dvrp = 0.5, interval_s_dvrp = 1.0,& 601 604 interval_a_dvrp = 0.0, & 602 605 interval_values_dvrp_prt, & 603 interval_h_dvrp_prt ,&606 interval_h_dvrp_prt = RESHAPE( & 604 607 (/ 270.0, 225.0, 225.0, 180.0, 70.0, 25.0, & 605 608 25.0, -25.0, ( 0.0, i9 = 1, 192 ) /), & 609 (/ 2, 100 /) ), & 606 610 interval_l_dvrp_prt = 0.5, & 607 611 interval_s_dvrp_prt = 1.0, & 608 612 interval_a_dvrp_prt = 0.0 609 #else610 REAL, DIMENSION(2,10) :: slicer_range_limits_dvrp611 612 REAL, DIMENSION(3,10) :: isosurface_color613 614 REAL(4), DIMENSION(2,100) :: interval_values_dvrp, &615 interval_values_dvrp_prt, interval_h_dvrp, &616 interval_h_dvrp_prt, interval_l_dvrp = 0.5, &617 interval_l_dvrp_prt = 0.5, interval_s_dvrp = 1.0, &618 interval_s_dvrp_prt = 1.0, interval_a_dvrp = 0.0, &619 interval_a_dvrp_prt = 0.0620 621 DATA slicer_range_limits_dvrp / -1.0, 1.0, -1.0, 1.0, -1.0, 1.0, &622 -1.0, 1.0, -1.0, 1.0, -1.0, 1.0, &623 -1.0, 1.0, -1.0, 1.0, -1.0, 1.0, &624 -1.0, 1.0 /625 626 DATA isosurface_color / 0.9, 0.9, 0.9, 0.8, 0.1, 0.1, 0.1, 0.1, 0.8, &627 0.1, 0.8, 0.1, 0.6, 0.1, 0.1, 0.1, 0.1, 0.6, &628 0.1, 0.6, 0.1, 0.4, 0.1, 0.1, 0.1, 0.1, 0.4, &629 0.1, 0.4, 0.1 /630 631 DATA interval_h_dvrp / 270.0, 225.0, 225.0, 180.0, 70.0, 25.0, &632 25.0, -25.0, 192 * 0.0 /633 634 DATA interval_h_dvrp_prt / 270.0, 225.0, 225.0, 180.0, 70.0, 25.0, &635 25.0, -25.0, 192 * 0.0 /636 #endif613 !#else 614 ! REAL, DIMENSION(2,10) :: slicer_range_limits_dvrp 615 616 ! REAL, DIMENSION(3,10) :: isosurface_color 617 618 ! REAL(4), DIMENSION(2,100) :: interval_values_dvrp, & 619 ! interval_values_dvrp_prt, interval_h_dvrp, & 620 ! interval_h_dvrp_prt, interval_l_dvrp = 0.5, & 621 ! interval_l_dvrp_prt = 0.5, interval_s_dvrp = 1.0, & 622 ! interval_s_dvrp_prt = 1.0, interval_a_dvrp = 0.0, & 623 ! interval_a_dvrp_prt = 0.0 624 625 ! DATA slicer_range_limits_dvrp / -1.0, 1.0, -1.0, 1.0, -1.0, 1.0, & 626 ! -1.0, 1.0, -1.0, 1.0, -1.0, 1.0, & 627 ! -1.0, 1.0, -1.0, 1.0, -1.0, 1.0, & 628 ! -1.0, 1.0 / 629 630 ! DATA isosurface_color / 0.9, 0.9, 0.9, 0.8, 0.1, 0.1, 0.1, 0.1, 0.8, & 631 ! 0.1, 0.8, 0.1, 0.6, 0.1, 0.1, 0.1, 0.1, 0.6, & 632 ! 0.1, 0.6, 0.1, 0.4, 0.1, 0.1, 0.1, 0.1, 0.4, & 633 ! 0.1, 0.4, 0.1 / 634 635 ! DATA interval_h_dvrp / 270.0, 225.0, 225.0, 180.0, 70.0, 25.0, & 636 ! 25.0, -25.0, 192 * 0.0 / 637 638 ! DATA interval_h_dvrp_prt / 270.0, 225.0, 225.0, 180.0, 70.0, 25.0, & 639 ! 25.0, -25.0, 192 * 0.0 / 640 !#endif 637 641 638 642 REAL(4), DIMENSION(:), ALLOCATABLE :: xcoor_dvrp, ycoor_dvrp, zcoor_dvrp
Note: See TracChangeset
for help on using the changeset viewer.
