Changeset 220 for palm/trunk
- Timestamp:
- Dec 18, 2008 7:00:36 AM (16 years ago)
- Location:
- palm/trunk
- Files:
-
- 6 edited
-
DOC/app/chapter_3.9.html (modified) (1 diff)
-
SCRIPTS/mrun (modified) (2 diffs)
-
SOURCE/flow_statistics.f90 (modified) (1 diff)
-
SOURCE/read_3d_binary.f90 (modified) (1 diff)
-
SOURCE/user_module.f90 (modified) (1 diff)
-
SOURCE/user_read_restart_data.f90 (modified) (2 diffs)
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/DOC/app/chapter_3.9.html
r210 r220 47 47 averaged mean profiles as provided with the standard PALM output are 48 48 used. 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 at50 the end of the precursor run. The profile information is then contained49 etc. for the precursor run appropriately, so that an output is done <span style="font-weight: bold;">at 50 the end</span> of the precursor run. The profile information is then contained 51 51 in 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 .... 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>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. .... 54 53 set bc_lr = 'dirichlet/radiation' ... 55 54 conserve_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 1896 1896 then 1897 1897 1898 printf "\n +++ no compiler specified for \"$ localhost $cond1 $cond2\""1898 printf "\n +++ no compiler specified for \"$host $cond1 $cond2\"" 1899 1899 locat=compiler_name; exit 1900 1900 … … 2375 2375 elif [[ $localhost = lcsgih ]] 2376 2376 then 2377 printf " compiler is called via ssh on \"hicegate \" \n"2377 printf " compiler is called via ssh on \"hicegate2\" \n" 2378 2378 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 " 2379 2379 # 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 1008 1008 ts_value(19,sr) = hom(nzb+1,1,pr_palm,sr) ! u'w' at k=0 1009 1009 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=01010 ts_value(21,sr) = hom(nzb+12,1,pr_palm,sr) ! q* 1011 1011 1012 1012 IF ( ts_value(5,sr) /= 0.0 ) THEN -
palm/trunk/SOURCE/read_3d_binary.f90
r216 r220 974 974 ! 975 975 !-- 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), & 979 979 tmp_2d, tmp_3d ) 980 980 -
palm/trunk/SOURCE/user_module.f90
r211 r220 37 37 INTEGER :: dots_num_palm, user_idummy 38 38 LOGICAL :: user_defined_namelist_found = .FALSE. 39 REAL :: user_ dummy39 REAL :: user_rdummy 40 40 41 41 ! -
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, & 4 4 tmp_2d, tmp_3d ) 5 5 … … 33 33 CHARACTER (LEN=20) :: field_char 34 34 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, & 36 36 nynf, nyn_on_file, nysc, nysf, nys_on_file, overlap_count 37 37
Note: See TracChangeset
for help on using the changeset viewer.
