Changeset 405


Ignore:
Timestamp:
Nov 16, 2009 12:20:59 PM (15 years ago)
Author:
raasch
Message:

some initial data assignments changed due to problems with gfortran compiler; bugfix in docu chapter 3.5.4; mrun test modification for hybrid mpi-openmp

Location:
palm/trunk
Files:
4 edited

Legend:

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

    r343 r405  
    11<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 Transitional//EN">
    22<html><head>
     3
    34<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>
    45<!--
    56@page { size: 21cm 29.7cm }
    67-->
    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
    99User-defined output quantities<br>
    1010</font></h3>A very typical request of users is the
     
    172172&nbsp;IF ( initializing_actions == 'read_restart_data' )&nbsp;
    173173THEN</span><br style="font-family: monospace;"><span style="font-family: monospace;">&nbsp;&nbsp;&nbsp;&nbsp;
    174 &nbsp; READ ( 13 )&nbsp; field_chr</span><br style="font-family: monospace;"><span style="font-family: monospace;">&nbsp;&nbsp;&nbsp;
    175 &nbsp; &nbsp;DO&nbsp; WHILE ( TRIM( field_chr ) /= '*** end
     174&nbsp; READ ( 13 )&nbsp; field_char</span><br style="font-family: monospace;"><span style="font-family: monospace;">&nbsp;&nbsp;&nbsp;
     175&nbsp; &nbsp;DO&nbsp; WHILE ( TRIM( field_char ) /= '*** end
    176176user ***' )</span><br style="font-family: monospace;"><br style="font-family: monospace;"><span style="font-family: monospace;">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
    177 &nbsp; &nbsp; SELECT CASE ( TRIM( field_chr ) )</span><br style="font-family: monospace;"><br style="font-family: monospace;"><span style="font-family: monospace;">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
     177&nbsp; &nbsp; SELECT CASE ( TRIM( field_char ) )</span><br style="font-family: monospace;"><br style="font-family: monospace;"><span style="font-family: monospace;">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
    178178&nbsp;&nbsp; &nbsp;&nbsp; CASE ( 'u2_av' )</span><br style="font-family: monospace;"><span style="font-family: monospace;">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
    179179&nbsp;&nbsp; &nbsp; ALLOCATE(
     
    186186&nbsp; &nbsp;&nbsp;
    187187&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
    188 TRIM(
    189 field_chr ), '" found in'</span><br style="font-family: monospace;"><span style="font-family: monospace;">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
     188TRIM( field_char ), '" found in'</span><br style="font-family: monospace;"><span style="font-family: monospace;">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
    190189&nbsp; &nbsp;&nbsp; &nbsp; &nbsp; PRINT*,
    191190'&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
     
    276275grid = 'zu'</span><br>&nbsp;<br>The <span style="font-family: monospace;">ELSE</span> case is
    277276only 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>
     277As 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>
     278on the u-grid. This features is useful for output of surface data (e.g.
     279heat fluxes). In this case, the corresponding 2d data has to be
     280resorted 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>
    279281<span style="font-family: monospace;">&nbsp;&nbsp;&nbsp;&nbsp;CASE
    280282( 'u2_xy' )</span><br style="font-family: monospace;"><span style="font-family: monospace;">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
  • palm/trunk/SCRIPTS/mrun

    r400 r405  
    31623162             if [[ $threads_per_task != 1 ]]
    31633163             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"
    31653168             fi
    31663169             if [[ $run_coupled_model = false ]]
     
    31973200                      export MV2_NUM_PORTS=2
    31983201                      #  The default setting of MV2_CPU_MAPPING gives best results
     3202                      # export MV2_ENABLE_AFFINITY=1
    31993203                      #  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
    32013210                   fi
    32023211                elif [[ $host = lcxt4 ]]
  • palm/trunk/SOURCE/CURRENT_MODIFICATIONS

    r403 r405  
    77var_ts is replaced by dots_max (modules,init_3d_model)
    88
    9 modules,init_3d_model
     9init_3d_model, modules
    1010
    1111
     
    1818Bugfix: calculation of precipitation_rate (calc_precipitation)
    1919
    20 sum_up_3d_data, time_integration, calc_precipitation
     20Bugfix: initial data assignments to some dvrp arrays changed due to error
     21messages from gfortran compiler (modules=
     22
     23calc_precipitation, modules, sum_up_3d_data, time_integration
    2124
    2225
  • palm/trunk/SOURCE/modules.f90

    r394 r405  
    55! Current revisions:
    66! -----------------
    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
    910!
    1011! Former revisions:
     
    583584                           topography_color = (/ 0.8, 0.7, 0.6 /)
    584585
    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( (/      &
    587588                                -1.0, 1.0, -1.0, 1.0, -1.0, 1.0, -1.0, 1.0, &
    588589                                -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( (/                  &
    592593                                0.9, 0.9, 0.9,  0.8, 0.1, 0.1,  0.1, 0.1, 0.8, &
    593594                                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 /) )
    596597
    597598    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 /) ),                     &
    600603                                  interval_l_dvrp = 0.5, interval_s_dvrp = 1.0,&
    601604                                  interval_a_dvrp = 0.0,                       &
    602605                                  interval_values_dvrp_prt,                    &
    603                                   interval_h_dvrp_prt,                         &
     606                                  interval_h_dvrp_prt = RESHAPE(               &
    604607                                  (/ 270.0, 225.0, 225.0, 180.0, 70.0, 25.0,   &
    605608                                     25.0, -25.0, ( 0.0, i9 = 1, 192 ) /),     &
     609                                                  (/ 2, 100 /) ),              &
    606610                                  interval_l_dvrp_prt = 0.5,                   &
    607611                                  interval_s_dvrp_prt = 1.0,                   &
    608612                                  interval_a_dvrp_prt = 0.0
    609 #else
    610     REAL, DIMENSION(2,10)     ::  slicer_range_limits_dvrp
    611 
    612     REAL, DIMENSION(3,10)     ::  isosurface_color
    613 
    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.0
    620 
    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 #endif
     613!#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
    637641
    638642    REAL(4), DIMENSION(:), ALLOCATABLE ::  xcoor_dvrp, ycoor_dvrp, zcoor_dvrp
Note: See TracChangeset for help on using the changeset viewer.