Changeset 220 for palm/trunk


Ignore:
Timestamp:
Dec 18, 2008 7:00:36 AM (13 years ago)
Author:
raasch
Message:

some small bugfixes in user_module, user_read_restart_data, read_3d_binary, flow_statistics and mrun

Location:
palm/trunk
Files:
6 edited

Legend:

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

    r210 r220  
    4747averaged mean profiles as provided with the standard PALM output are
    4848used. The user has to set parameters <a href="chapter_4.2.html#dt_data_output_pr">dt_data_output_pr</a>, <a href="chapter_4.2.html#averaging_interval">averaging_interval</a>,
    49 etc. for the precursor run appropriately, so that an output is done at
    50 the end of the precursor run. The profile information is then contained
     49etc. for the precursor run appropriately, so that an output is done <span style="font-weight: bold;">at
     50the end</span> of the precursor run. The profile information is then contained
    5151in the restart (binary) file created at the end of the precursor run
    52 and can be used by the main run. <span style="font-weight: bold;">It is very important that the mean profiles at the end of the precursor run are in a stationary or quasi-stationary state</span>, because otherwise it may not be justified to use them as constant profiles at the inflow. <span style="font-weight: bold;">Also, turbulence at the end of the precursor run should be fully developed. </span>Otherwise, the main run would need an additional spinup-time at the beginning to get the turbulence to its final stage.<br></li><li>The
    53 main run has to read the binary data from the precursor run .... &nbsp;
     52and can be used by the main run. <span style="font-weight: bold;">It is very important that the mean profiles at the end of the precursor run are in a stationary or quasi-stationary state</span>, because otherwise it may not be justified to use them as constant profiles at the inflow. <span style="font-weight: bold;">Also, turbulence at the end of the precursor run should be fully developed. </span>Otherwise, the main run would need an additional spinup-time at the beginning to get the turbulence to its final stage.<br></li><li>In order to read the binary data from the precursor run, parameter <a href="chapter_4.1.html#initializing_actions">initializing_actions</a> = <span style="font-style: italic;">'read_data_for_recycling'</span> has to be set in the main run. &nbsp; .... &nbsp;
    5453&nbsp; set bc_lr = 'dirichlet/radiation' ... &nbsp;
    5554conserve_volume_flow = .TRUE. ? ....</li><li>The main run allows to use a larger domain along x and y than in the precursor run ...</li><li>Recycling width should be equal to the domain size of the precursor run (or larger?)</li></ul><pre style="line-height: 100%;"><font style="font-size: 10pt;" size="2">#</font><br><font style="font-size: 10pt;" size="2">RUN_CONTROL<span style="font-weight: bold;">_O</span></font><font style="font-size: 10pt;" size="2"> out:loc:tr d3<span style="font-weight: bold;">o</span># ~/palm/current_version/JOBS/$fname/MONITORING <span style="font-weight: bold;">_</span></font><font style="font-size: 10pt;" size="2"><span style="font-weight: bold;">o</span>_</font><font style="font-size: 10pt;" size="2">rc</font><br><font style="font-size: 10pt;" size="2">HEADER</font><font style="font-size: 10pt;" size="2"><span style="font-weight: bold;">_O</span> out:loc:tr d3<span style="font-weight: bold;">o</span># ~/palm/current_version/JOBS/$fname/MONITORING <span style="font-weight: bold;">_</span></font><font style="font-size: 10pt;" size="2"><span style="font-weight: bold;">o</span>_</font><font style="font-size: 10pt;" size="2">header</font><br><font style="font-size: 10pt;" size="2">PLOT1D_PAR<span style="font-weight: bold;">_O</span></font><font style="font-size: 10pt;" size="2"> out:loc:tr pr<span style="font-weight: bold;">o</span># ~/palm/current_version/JOBS/$fname/OUTPUT <span style="font-weight: bold;">_</span></font><font style="font-size: 10pt;" size="2"><span style="font-weight: bold;">o</span>_</font><font style="font-size: 10pt;" size="2">pr_par</font><br><font style="font-size: 10pt;" size="2">PLOT1D_DATA</font><font style="font-size: 10pt;" size="2"><span style="font-weight: bold;">_O</span> out:loc:tr pr<span style="font-weight: bold;">o</span># ~/palm/current_version/JOBS/$fname/OUTPUT <span style="font-weight: bold;">_</span></font><font style="font-size: 10pt;" size="2"><span style="font-weight: bold;">o</span>_</font><font style="font-size: 10pt;" size="2">pr_in</font><br></pre><br><hr>
  • palm/trunk/SCRIPTS/mrun

    r215 r220  
    18961896 then
    18971897
    1898     printf "\n  +++ no compiler specified for \"$localhost $cond1 $cond2\""
     1898    printf "\n  +++ no compiler specified for \"$host $cond1 $cond2\""
    18991899    locat=compiler_name; exit
    19001900
     
    23752375       elif [[ $localhost = lcsgih ]]
    23762376       then
    2377           printf "      compiler is called via ssh on \"hicegate\" \n"
     2377          printf "      compiler is called via ssh on \"hicegate2\" \n"
    23782378          ssh  130.75.4.103  -l $usern  ". /usr/share/modules/init/bash; module load mpt; cd $TEMPDIR; make -f Makefile PROG=a.out  F90=$compiler_name  COPT=\"$cpp_options\"  F90FLAGS=\"$fopts\"  LDFLAGS=\"$lopts\" 2>&1 "
    23792379#          ssh  130.75.4.102  -l $usern  ". /usr/share/modules/init/bash; module load mvapich2; cd $TEMPDIR; make -f Makefile PROG=a.out  F90=$compiler_name  COPT=\"$cpp_options\"  F90FLAGS=\"$fopts\"  LDFLAGS=\"$lopts\" 2>&1 "
  • palm/trunk/SOURCE/flow_statistics.f90

    r198 r220  
    10081008       ts_value(19,sr) = hom(nzb+1,1,pr_palm,sr)    ! u'w'    at k=0
    10091009       ts_value(20,sr) = hom(nzb+2,1,pr_palm,sr)    ! v'w'    at k=0
    1010        ts_value(21,sr) = hom(nzb+12,1,pr_palm,sr)   ! w'q'    at k=0
     1010       ts_value(21,sr) = hom(nzb+12,1,pr_palm,sr)   ! q*
    10111011
    10121012       IF ( ts_value(5,sr) /= 0.0 )  THEN
  • palm/trunk/SOURCE/read_3d_binary.f90

    r216 r220  
    974974!
    975975!--    Read user-defined restart data
    976        CALL user_read_restart_data( nxlfa, nxl_on_file, nxrfa, nxr_on_file, &
    977                                     nynfa, nyn_on_file, nysfa, nys_on_file, &
    978                                     offset_xa, offset_ya, overlap_count(i), &
     976       CALL user_read_restart_data( i, nxlfa, nxl_on_file, nxrfa, nxr_on_file, &
     977                                    nynfa, nyn_on_file, nysfa, nys_on_file,    &
     978                                    offset_xa, offset_ya, overlap_count(i),    &
    979979                                    tmp_2d, tmp_3d )
    980980
  • palm/trunk/SOURCE/user_module.f90

    r211 r220  
    3737    INTEGER ::  dots_num_palm, user_idummy
    3838    LOGICAL ::  user_defined_namelist_found = .FALSE.
    39     REAL    ::  user_dummy
     39    REAL    ::  user_rdummy
    4040
    4141!
  • palm/trunk/SOURCE/user_read_restart_data.f90

    r216 r220  
    1  SUBROUTINE user_read_restart_data( nxlfa, nxl_on_file, nxrfa, nxr_on_file, &
    2                                     nynfa, nyn_on_file, nysfa, nys_on_file, &
    3                                     offset_xa, offset_ya, overlap_count,    &
     1 SUBROUTINE user_read_restart_data( i, nxlfa, nxl_on_file, nxrfa, nxr_on_file, &
     2                                    nynfa, nyn_on_file, nysfa, nys_on_file,    &
     3                                    offset_xa, offset_ya, overlap_count,       &
    44                                    tmp_2d, tmp_3d )
    55
     
    3333    CHARACTER (LEN=20) :: field_char
    3434
    35     INTEGER ::  k, nxlc, nxlf, nxl_on_file, nxrc, nxrf, nxr_on_file, nync, &
     35    INTEGER ::  i, k, nxlc, nxlf, nxl_on_file, nxrc, nxrf, nxr_on_file, nync, &
    3636                nynf, nyn_on_file, nysc, nysf, nys_on_file, overlap_count
    3737
Note: See TracChangeset for help on using the changeset viewer.