Changeset 116 for palm


Ignore:
Timestamp:
Oct 11, 2007 2:30:27 AM (17 years ago)
Author:
raasch
Message:

further preliminary updates concerning particle sorting and documentation

Location:
palm/trunk
Files:
12 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/DOC/app/chapter_4.1.html

    r115 r116  
    11<!DOCTYPE html PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
    2 <html>
    3 <head>
     2<html><head>
    43
    54
     
    2019 
    2120 
    22   <title>PALM chapter 4.1</title>
    23 </head>
     21  <title>PALM chapter 4.1</title></head>
    2422
    2523<body>
     
    21532151
    21542152
    2155 Cloud droplets require to use the particle package (<span style="font-weight: bold;">mrun</span>-option <span style="font-family: monospace;">-p particles</span>),
    2156 so in this case a particle corresponds to a droplet. The droplet
     2153<span style="font-weight: bold;"></span><span style="font-family: monospace;"></span>
     2154
     2155
     2156
     2157
     2158Cloud droplets require to use&nbsp;particles (i.e. the NAMELIST group <span style="font-family: Courier New,Courier,monospace;">particles_par</span> has to be included in the parameter file<span style="font-family: monospace;"></span>). Then each particle is a representative for a certain number of droplets. The droplet
    21572159features (number of droplets, initial radius, etc.) can be steered with
    21582160the&nbsp; respective particle parameters (see e.g. <a href="#chapter_4.2.html#radius">radius</a>).
     
    1260212604
    1260312605
    12604 </body>
    12605 </html>
     12606</body></html>
  • palm/trunk/DOC/app/chapter_4.2.html

    r108 r116  
    11<!DOCTYPE html PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
    2 <html>
    3 <head>
     2<html><head>
    43
    54
     
    1211 
    1312 
    14   <title>PALM chapter 4.2</title>
    15 </head>
    16 
     13  <title>PALM chapter 4.2</title></head>
    1714
    1815<body>
     
    86868683
    86878684
    8688     <tr>
     8685    <tr><td style="vertical-align: top;"><span style="font-weight: bold;"><a name="dt_sort_particles"></a>dt_sort_particles</span></td><td style="vertical-align: top;">R</td><td style="vertical-align: top;"><span style="font-style: italic;">0.0</span></td><td style="vertical-align: top;">Temporal interval for sorting particles (in s).<br><br>By
     8686default, particles are sorted in memory in a way that their order
     8687follows the order in which the gridpoint values are stored. This may
     8688improve cache coherence in case of larger numbers of particles and
     8689gridpoints. However, since the sorting itself is time consuming and
     8690since the requirement of sorting depends on the strength of mixing in
     8691the flow, performance can be improved if the sorting is applied only
     8692after certain time intervals. The proper length of this interval <span style="font-weight: bold;">dt_sort_particles</span> must be determined empirically by carrying out test runs with different intervals. Check file <a href="chapter_3.4.html#CPU_MEASURES">CPU_MEASURES</a> to find the value of <span style="font-weight: bold;">dt_sort_particles</span> which gives the best performance.<br><br><span style="font-weight: bold;">Note:</span><br>In case of <a href="chapter_4.1.html#cloud_droplets">cloud_droplets</a> = <span style="font-style: italic;">.T.</span>, any given non-zero value of <span style="font-weight: bold;">dt_sort_particles</span> will be reset to zero and a corresponding warning message will appear in the job protocol.</td></tr><tr>
    86898693
    86908694
     
    1159011594
    1159111595
    11592 </body>
    11593 </html>
     11596</body></html>
  • palm/trunk/DOC/app/chapter_4.6.html

    r115 r116  
    11<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 Transitional//EN">
    2 <html>
    3 <head>
     2<html><head>
    43
    54
     
    1514 
    1615 
    17   <title>PALM chapter 4.6</title>
    18   <meta name="GENERATOR" content="StarOffice 7 (Win32)">
     16 
     17  <title>PALM chapter 4.6</title><meta name="GENERATOR" content="StarOffice 7 (Win32)">
    1918
    2019
     
    5049@page { size: 21cm 29.7cm }
    5150-->
    52   </style>
    53 </head>
    54 
     51  </style></head>
    5552
    5653<body style="direction: ltr;" lang="en-US">
     
    37323729
    37333730
    3734  <tr>
     3731 <tr><td><a href="chapter_4.2.html#dt_sort_particles"><span style="font-weight: bold;">dt_sort_particles</span></a></td><td>P</td><td>R</td><td><span style="font-style: italic;">0.0</span></td><td>Temporal interval for sorting particles (in s).</td></tr><tr>
    37353732
    37363733
     
    1182011817
    1182111818
    11822 </body>
    11823 </html>
     11819</body></html>
  • palm/trunk/SCRIPTS/mbuild

    r83 r116  
    9090     #                    from configuration file
    9191     #                    host identifier (local_host) is read from config file
     92     # 10/10/07 - Siggi - bugfix: handling of comment lines
    9293
    9394
     
    185186 while read line
    186187 do
    187     if [[ "$line" != ""  ||  $(echo $line | cut -c1) != "#" ]]
     188    if [[ "$line" != ""  &&  $(echo $line | cut -c1) != "#" ]]
    188189    then
    189190       HOSTNAME=`echo $line | cut -d" " -s -f2`
     
    228229 while read line
    229230 do
    230     if [[ "$line" != ""  ||  $(echo $line | cut -c1) != "#" ]]
     231    if [[ "$line" != ""  &&  $(echo $line | cut -c1) != "#" ]]
    231232    then
    232233       local_username=`echo $line | cut -d" " -s -f2`
     
    252253 while read line
    253254 do
    254     if [[ "$line" != ""  &&  "$(echo $line | cut -d" " -s -f3)" = "" ]]
     255    if [[ "$line" != ""  &&  "$(echo $line | cut -d" " -s -f3)" = ""  &&  $(echo $line | cut -c1) != "#" ]]
    255256    then
    256257       var=`echo $line | cut -d" " -s -f1 | cut -c2-`
     
    721722       # IP-ADRESSE DES REMOTE-RECHNERS BESTIMMEN
    722723    case  $remote_host  in
    723         (lcmuk)          remote_addres=130.75.105.3;;
     724        (lcmuk)          remote_addres=130.75.105.2;;
    724725        (lctit)          remote_addres=172.17.75.161;;
    725726        (decalpha)       remote_addres=165.132.26.56;;
     
    741742       # REMOTE-USERNAMEN ERMITTELN
    742743    line=""
     744    found=false
    743745    grep  "$remote_host_string" $config_file | grep "%remote_username"  >  tmp_mbuild
    744746    while read line1
    745747    do
     748
    746749       if [[ $(echo $line1 | cut -d" " -s -f3-) = "$remote_host_string" ]]
    747750       then
    748751          line="$line1"
    749752       fi
     753
     754       if [[ "$line" != ""  &&  $(echo $line | cut -c1) != "#" ]]
     755       then
     756          remote_username=`echo $line | cut -d" " -s -f2`
     757          found=true
     758       fi
     759
    750760    done < tmp_mbuild
    751761
    752     if [[ "$line" = ""  ||  $(echo $line | cut -c1) = "#" ]]
     762    if [[ $found = false ]]
    753763    then
    754764       printf "\n  +++ no remote username found in configuration file"
    755765       printf "\n      for \"$remote_host_string\" "
    756766       locat=config_file; exit
    757     else
    758        remote_username=`echo $line | cut -d" " -s -f2`
    759767    fi
    760768
     
    766774    while read line1
    767775    do
     776
    768777       if [[ $(echo $line1 | cut -d" " -s -f3-) = "$remote_host_string" ]]
    769778       then
    770779          line="$line1"
    771780       fi
     781
     782       if [[ "$line" != ""  &&  $(echo $line | cut -c1) != "#" ]]
     783       then
     784          remote_source_path=`echo $line | cut -d" " -s -f2`
     785       fi
     786
    772787    done < tmp_mbuild
    773788
    774     if [[ "$line" != ""  &&  $(echo $line | cut -c1) != "#" ]]
    775     then
    776        remote_source_path=`echo $line | cut -d" " -s -f2`
    777     fi
    778789    if [[ "$remote_source_path" = "" ]]
    779790    then
     
    798809    while read line1
    799810    do
     811
    800812       if [[ $(echo $line1 | cut -d" " -s -f3-) = "$remote_host_string" ]]
    801813       then
    802814          line="$line1"
    803815       fi
     816
     817       if [[ "$line" != ""  &&  $(echo $line | cut -c1) != "#" ]]
     818       then
     819          remote_md=`echo $line | cut -d" " -s -f2`
     820       fi
     821
    804822    done < tmp_mbuild
    805823
    806     if [[ "$line" != ""  &&  $(echo $line | cut -c1) != "#" ]]
    807     then
    808        remote_md=`echo $line | cut -d" " -s -f2`
    809     fi
    810824    if [[ "$remote_md" = "" ]]
    811825    then
     
    826840       # COMPILERNAMEN ERMITTELN
    827841    line=""
     842    found=false
    828843    grep  "$remote_host_string" $config_file | grep "%compiler_name "  >  tmp_mbuild
    829844    while read line1
    830845    do
     846
    831847       if [[ $(echo $line1 | cut -d" " -s -f3-) = "$remote_host_string" ]]
    832848       then
    833849          line="$line1"
    834850       fi
     851
     852       if [[ "$line" != ""  &&  $(echo $line | cut -c1) != "#" ]]
     853       then
     854          compiler_name=`echo $line | cut -d" " -s -f2`
     855          found=true
     856       fi
     857
    835858    done < tmp_mbuild
    836859
    837     if [[ "$line" = ""  ||  $(echo $line | cut -c1) = "#" ]]
     860    if [[ $found = false ]]
    838861    then
    839862       printf "\n  +++ no compiler name found in configuration file"
    840863       printf "\n      for \"$remote_host_string\" "
    841864       locat=config_file; exit
    842     else
    843        compiler_name=`echo $line | cut -d" " -s -f2`
    844865    fi
    845866
     
    850871    then
    851872       line=""
     873       found=false
    852874       grep  "$remote_host_string" $config_file | grep "%compiler_name_ser"  >  tmp_mbuild
    853875       while read line1
    854876       do
     877
    855878          if [[ $(echo $line1 | cut -d" " -s -f3-) = "$remote_host_string" ]]
    856879          then
    857880             line="$line1"
    858881          fi
     882
     883          if [[ "$line" != ""  &&  $(echo $line | cut -c1) != "#" ]]
     884          then
     885             compiler_name_ser=`echo $line | cut -d" " -s -f2`
     886             found=true
     887          fi
     888
    859889       done < tmp_mbuild
    860890
    861        if [[ "$line" = ""  ||  $(echo $line | cut -c1) = "#" ]]
     891       if [[ $found = false ]]
    862892       then
    863893          printf "\n  +++ no serial compiler name found in configuration file"
    864894          printf "\n      for \"$remote_host_string\" "
    865895          locat=config_file; exit
    866        else
    867           compiler_name_ser=`echo $line | cut -d" " -s -f2`
    868896       fi
    869897    else
     
    875903       # PRAEPROZESSOR-OPTIONEN/DIREKTIVEN ERMITTELN
    876904    line=""
     905    found=false
    877906    grep  "$remote_host_string" $config_file | grep "%cpp_options"  >  tmp_mbuild
    878907    while read line1
    879908    do
     909
    880910       if [[ $(echo $line1 | cut -d" " -s -f3-) = "$remote_host_string" ]]
    881911       then
    882912          line="$line1"
    883913       fi
     914
     915       if [[ "$line" != ""  &&  $(echo $line | cut -c1) != "#" ]]
     916       then
     917             # EVENTUELLE DOPPELPUNKTE AUS OPTIONSSTRING ENTFERNEN
     918          cpp_options=`echo $line | cut -d" " -s -f2 | sed 's/::/%DUM%/g' | sed 's/:/ /g' | sed 's/%DUM%/:/g'`
     919          found=true
     920       fi
     921
    884922    done < tmp_mbuild
    885     if [[ "$line" = ""  ||  $(echo $line | cut -c1) = "#" ]]
     923
     924    if [[ $found = false ]]
    886925    then
    887926       printf "\n  +++ no preprocessor options found in configuration file"
    888927       printf "\n      for \"$remote_host_string\" "
    889928       locat=config_file; exit
    890     else
    891           # EVENTUELLE DOPPELPUNKTE AUS OPTIONSSTRING ENTFERNEN
    892        cpp_options=`echo $line | cut -d" " -s -f2 | sed 's/::/%DUM%/g' | sed 's/:/ /g' | sed 's/%DUM%/:/g'`
    893     fi
     929    fi
     930
    894931
    895932       # RECHNERSPEZIFISCHE CPP-DIREKTIVEN HINZUFUEGEN
     
    922959    while read line1
    923960    do
     961
    924962       if [[ $(echo $line1 | cut -d" " -s -f3-) = "$remote_host_string" ]]
    925963       then
    926964          line="$line1"
    927965       fi
     966
     967       if [[ "$line" != ""  &&  $(echo $line | cut -c1) != "#" ]]
     968       then
     969             # EVENTUELLE DOPPELPUNKTE AUS OPTIONSSTRING ENTFERNEN
     970          netcdf_inc=`echo $line | cut -d" " -s -f2 | sed 's/::/%DUM%/g' | sed 's/:/ /g' | sed 's/%DUM%/:/g'`
     971       fi
     972
    928973    done < tmp_mbuild
    929     if [[ "$line" != ""  &&  $(echo $line | cut -c1) != "#" ]]
    930     then
    931           # EVENTUELLE DOPPELPUNKTE AUS OPTIONSSTRING ENTFERNEN
    932        netcdf_inc=`echo $line | cut -d" " -s -f2 | sed 's/::/%DUM%/g' | sed 's/:/ /g' | sed 's/%DUM%/:/g'`
    933     fi
    934974
    935975    line=""
     
    937977    while read line1
    938978    do
     979
    939980       if [[ $(echo $line1 | cut -d" " -s -f3-) = "$remote_host_string" ]]
    940981       then
    941982          line="$line1"
    942983       fi
     984
     985       if [[ "$line" != ""  &&  $(echo $line | cut -c1) != "#" ]]
     986       then
     987             # EVENTUELLE DOPPELPUNKTE AUS OPTIONSSTRING ENTFERNEN
     988          netcdf_lib=`echo $line | cut -d" " -s -f2 | sed 's/::/%DUM%/g' | sed 's/:/ /g' | sed 's/%DUM%/:/g'`
     989       fi
     990
    943991    done < tmp_mbuild
    944     if [[ "$line" != ""  &&  $(echo $line | cut -c1) != "#" ]]
    945     then
    946           # EVENTUELLE DOPPELPUNKTE AUS OPTIONSSTRING ENTFERNEN
    947        netcdf_lib=`echo $line | cut -d" " -s -f2 | sed 's/::/%DUM%/g' | sed 's/:/ /g' | sed 's/%DUM%/:/g'`
    948     fi
    949992
    950993
     
    952995       # COMPILEROPTIONEN ERMITTELN
    953996    line=""
     997    found=false
    954998    grep  "$remote_host_string" $config_file | grep "%fopts"  >  tmp_mbuild
    955999    while read line1
    9561000    do
     1001
    9571002       if [[ $(echo $line1 | cut -d" " -s -f3-) = "$remote_host_string" ]]
    9581003       then
    9591004          line="$line1"
    9601005       fi
     1006
     1007       if [[ "$line" != ""  &&  $(echo $line | cut -c1) != "#" ]]
     1008       then
     1009             # EVENTUELLE DOPPELPUNKTE AUS DIREKTIVENSTRING ENTFERNEN
     1010          compiler_options=`echo $line | cut -d" " -s -f2 | sed 's/::/%DUM%/g' | sed 's/:/ /g' | sed 's/%DUM%/:/g'`
     1011          found=true
     1012
     1013             # NETCDF-INCLUDEVERZEICHNIS HINZUFUEGEN
     1014          compiler_options="$compiler_options $netcdf_inc"
     1015       fi
     1016
    9611017    done < tmp_mbuild
    962     if [[ "$line" = ""  ||  $(echo $line | cut -c1) = "#" ]]
     1018
     1019    if [[ $found = false ]]
    9631020    then
    9641021       printf "\n  +++ no compiler options found in configuration file"
    9651022       printf "\n      for \"$remote_host_string\" "
    9661023       locat=config_file; exit
    967     else
    968           # EVENTUELLE DOPPELPUNKTE AUS DIREKTIVENSTRING ENTFERNEN
    969        compiler_options=`echo $line | cut -d" " -s -f2 | sed 's/::/%DUM%/g' | sed 's/:/ /g' | sed 's/%DUM%/:/g'`
    970 
    971           # NETCDF-INCLUDEVERZEICHNIS HINZUFUEGEN
    972        compiler_options="$compiler_options $netcdf_inc"
    9731024    fi
    9741025
     
    9771028       # LADER-OPTIONEN ERMITTELN
    9781029    line=""
     1030    found=false
    9791031    grep  "$remote_host_string" $config_file | grep "%lopts"  >  tmp_mbuild
    9801032    while read line1
    9811033    do
     1034
    9821035       if [[ $(echo $line1 | cut -d" " -s -f3-) = "$remote_host_string" ]]
    9831036       then
    9841037          line="$line1"
    9851038       fi
     1039
     1040       if [[ "$line" != ""  &&  $(echo $line | cut -c1) != "#" ]]
     1041       then
     1042             # EVENTUELLE DOPPELPUNKTE AUS DIREKTIVENSTRING ENTFERNEN
     1043          loader_options=`echo $line | cut -d" " -s -f2 | sed 's/::/%DUM%/g' | sed 's/:/ /g' | sed 's/%DUM%/:/g'`
     1044          found=true
     1045
     1046             # NETCDF-LIBRARY HINZUFUEGEN
     1047          loader_options="$loader_options $netcdf_lib"
     1048       fi
     1049
    9861050    done < tmp_mbuild
    9871051
    988     if [[ "$line" = ""  ||  $(echo $line | cut -c1) = "#" ]]
     1052    if [[ $found = false ]]
    9891053    then
    9901054       printf "\n  +++ no loader options found in configuration file"
    9911055       printf "\n      for \"$remote_host_string\" "
    9921056       locat=config_file; exit
    993     else
    994           # EVENTUELLE DOPPELPUNKTE AUS DIREKTIVENSTRING ENTFERNEN
    995        loader_options=`echo $line | cut -d" " -s -f2 | sed 's/::/%DUM%/g' | sed 's/:/ /g' | sed 's/%DUM%/:/g'`
    996 
    997           # NETCDF-LIBRARY HINZUFUEGEN
    998        loader_options="$loader_options $netcdf_lib"
    9991057    fi
    10001058
  • palm/trunk/SOURCE/CURRENT_MODIFICATIONS

    r114 r116  
    77(user_interface).
    88
    9 init_grid, init_pegrid, modules, user_interface
     9Frequence of sorting particles can be controlled with new particles_par
     10parameter dt_sort_particles. Sorting is moved from the SGS timestep loop in
     11advec_particles after the end of this loop.
     12
     13advec_particles, check_parameters, init_grid, init_pegrid, modules, package_parin, read_var_list, user_interface, write_var_list
    1014
    1115
  • palm/trunk/SOURCE/advec_particles.f90

    r114 r116  
    44! Actual revisions:
    55! -----------------
     6! Sorting of particles is controlled by dt_sort_particles and moved from
     7! the SGS timestep loop after the end of this loop.
    68! Bugfix: pleft/pright changed to pnorth/psouth in sendrecv of particle tail
    79! numbers along y
     
    32543256!    ENDDO
    32553257
    3256 !
    3257 !--    Sort particles in the sequence the gridboxes are stored in the memory
    3258        CALL sort_particles
    3259 
    32603258!    WRITE ( 9, * ) '*** advec_particles: ##9'
    32613259!    CALL local_flush( 9 )
     
    32863284
    32873285    ENDDO   ! timestep loop
     3286
     3287
     3288!
     3289!-- Sort particles in the sequence the gridboxes are stored in the memory
     3290    time_sort_particles = time_sort_particles + dt_3d
     3291    IF ( time_sort_particles >= dt_sort_particles )  THEN
     3292       CALL sort_particles
     3293       time_sort_particles = MOD( time_sort_particles, &
     3294                                  MAX( dt_sort_particles, dt_3d ) )
     3295    ENDIF
    32883296
    32893297
  • palm/trunk/SOURCE/check_parameters.f90

    r114 r116  
    44! Actual revisions:
    55! -----------------
    6 ! Multigrid solver allows topography
     6! Multigrid solver allows topography, checking of dt_sort_particles
    77!
    88! Former revisions:
     
    15001500
    15011501!
     1502!-- Check the interval for sorting particles.
     1503!-- Using particles as cloud droplets requires sorting after each timestep.
     1504    IF ( dt_sort_particles /= 0.0  .AND.  cloud_droplets )  THEN
     1505       dt_sort_particles = 0.0
     1506       IF ( myid == 0 )  THEN
     1507          PRINT*, '+++ WARNING: check_parameters:'
     1508          PRINT*, '    dt_sort_particles is reset to 0.0 because ', &
     1509                       'of cloud_droplets = .TRUE.'
     1510       ENDIF
     1511    ENDIF
     1512
     1513!
    15021514!-- Set the default intervals for data output, if necessary
    15031515!-- NOTE: dt_dosp has already been set in package_parin
  • palm/trunk/SOURCE/init_grid.f90

    r114 r116  
    789789!--       Bit 4:  wall to the left
    790790!--       Bit 5:  wall to the right
    791 !--       Bit 6:  inside / outside building (1/0)
     791!--       Bit 6:  inside building
    792792
    793793          flags = 0
  • palm/trunk/SOURCE/modules.f90

    r114 r116  
    55! Actual revisions:
    66! -----------------
    7 ! +flags, wall_flags_1..10
     7! +dt_sort_particles, time_sort_particles, flags, wall_flags_1..10
    88!
    99! Former revisions:
     
    876876    LOGICAL, DIMENSION(:), ALLOCATABLE ::  particle_mask, tail_mask
    877877
    878     REAL    ::  c_0 = 3.0, dt_min_part = 0.0002,                               &
     878    REAL    ::  c_0 = 3.0, dt_min_part = 0.0002, dt_sort_particles = 0.0,      &
    879879                dt_write_particle_data = 9999999.9, dvrp_psize = 9999999.9,    &
    880880                end_time_prel = 9999999.9, initial_weighting_factor = 1.0,     &
     
    883883                particle_advection_start = 0.0, sgs_wfu_part = 0.3333333,      &
    884884                sgs_wfv_part = 0.3333333, sgs_wfw_part = 0.3333333,            &
    885                 time_write_particle_data = 0.0
     885                time_sort_particles = 0.0, time_write_particle_data = 0.0
    886886
    887887    REAL, DIMENSION(max_number_of_particle_groups) ::  &
  • palm/trunk/SOURCE/package_parin.f90

    r77 r116  
    4545    NAMELIST /particles_par/      bc_par_b, bc_par_lr, bc_par_ns, bc_par_t,    &
    4646                                  density_ratio, radius, dt_dopts,             &
    47                                   dt_min_part, dt_prel,                        &
     47                                  dt_min_part, dt_prel, dt_sort_particles,     &
    4848                                  dt_write_particle_data, dvrp_psize,          &
    4949                                  end_time_prel, initial_weighting_factor,     &
  • palm/trunk/SOURCE/read_var_list.f90

    r110 r116  
    44! Actual revisions:
    55! -----------------
    6 !
     6! +time_sort_particles
    77!
    88! Former revisions:
     
    4646    USE indices
    4747    USE model_1d
     48    USE particle_attributes
    4849    USE pegrid
    4950    USE profil_parameter
     
    395396          CASE ( 'time_run_control' )
    396397             READ ( 13 )  time_run_control
     398          CASE ( 'time_sort_particles' )
     399             READ ( 13 )  time_sort_particles
    397400          CASE ( 'timestep_scheme' )
    398401             READ ( 13 )  timestep_scheme
  • palm/trunk/SOURCE/write_var_list.f90

    r110 r116  
    44! Actual revisions:
    55! -----------------
    6 !
     6! +time_sort_particles
    77!
    88! Former revisions:
     
    4646    USE indices
    4747    USE model_1d
     48    USE particle_attributes
    4849    USE pegrid
    4950    USE profil_parameter
     
    341342    WRITE ( 14 )  'time_run_control              '
    342343    WRITE ( 14 )  time_run_control
     344    WRITE ( 14 )  'time_sort_particles           '
     345    WRITE ( 14 )  time_sort_particles
    343346    WRITE ( 14 )  'timestep_scheme               '
    344347    WRITE ( 14 )  timestep_scheme
Note: See TracChangeset for help on using the changeset viewer.