The program generates the run script and .INN files required to iterate maplets using lithosP.
make_script.in - Text file containing a list of the landmarks to be processed, appended with the end-of-file identifier, 'END'.
INIT_LITHOS - make_scriptP reads the USRMX (maximum number of core processors) value.
- run.sh - Controlling run script for batch maplet iteration using parallel processing.
- run script for each core processor:
- etc ...
- run script for each core processor:
.INN files - lithosP command seed file for each maplet.
- rem_script.b - Script for removal of temporary files (working directory clean-up).
.OOT files - standard output from lithosP for each landmark.]
processed files outputs:
You must set the correct number of core processors for USRMX, a variable contained in INIT_LITHOS. Refer to that entry for further information.
We need to develop a standard way to put in these cross references.
For more information on iteration procedures, refer to the 'How-To': Standard_iteration.
- Create Input Files
The make_script.in file lists the names for the landmarks that should be processed.
Here is a sample make_script.in file:
EE0001 EE0002 EE0003 etc... END
Be sure that the list is ordered by:
- highest-to-lowest maplet ground sample distance
name (region identifier (e.g., FF) + unique numeric identifier (e.g., 0001))
Here is a sample make_scriptP.seed file (see that entry for for further information):
# This script performs two maplet iterations in parallel mode. # It reinitializes slope and albedo. # It recomputes maplet overlaps. # It recomputes limb fits. # It recomputes central maplet vector. # Topography determination from slopes is conditioned by: # Initial topography, # Overlapping topography, # Shadowing, # Differential stereo, # Limb heghts. # Dataless area slopes are set from the shape model. # # The script hides most of the screen output. In order to obtain the # full display replace H* by #H*. # # Modified 2011_12_31. Header and comments added. # Modified 2014_02_26. Hide screen output added. # H* x .025 # # Initialize slope and albedo. # a y y # # Begin first iteration. # 0 40 1 0 1 n 0 y u 0 40 # # Begin topo determination. # 2 8 fill no data with shape slope 2 condition with existing topo .01 7 condition with shadows .025 1 condition with overlapping maps 4 condition with limbs 1, 2.5, 3 6 condition with differential stereo y 0 begin integration .025 15 0 u 1 0 2 n 0 y 1 0 1 n 0 y u # # Begin second iteration. # 0 40 1 0 1 n 0 y u 0 40 # # Begin topo determination. # 2 8 2 .01 7 .025 1 4 1, 2.5, 3 6 y 0 .025 15 0 u 1 0 2 n 0 y H* 1 0 1 n 0 y H* u # # Begin overlaps and limb fits. # o RECENT y 1 o RECENT n 3 y 1, 3, 5 # # Central vector determination. # v RECENT q END
- Run make_scriptP
make_scriptP generates a .INN file for each maplet.
Here is a sample .INN file for maplet EE0001.INN:
01 i EE0001 H* x .025 0 40 etc ...
The first three lines of the .INN file supply a processor identifier (e.g., 01), and load the landmark in lithosP. The remainder of the .INN file lists a copy of the lithosP commands contained in make_scriptP.seed.
make_scriptP also generates the controller run script run.sh, which looks like this:
chmod +x run_script01.b chmod +x run_script02.b etc ... cd SHAPEFILES chmod +w SHAPE.TXT cp SHAPE.TXT TSHP01.TXT cp SHAPE.TXT TSHP02.TXT etc ... chmod -w SHAPE.TXT cd .. ./run_script01.b & ./run_script02.b & etc ...
make_scriptP also generates the run_script01.b, run_script02.b etc. files, which look like this:
rm -f EE0001.OOT /usr/local/bin/LITHOSP < EE0001.INN > EE0001.OOT cp LMRK_DISPLAY01.pgm TESTFILES1/EE0001.pgm rm -f EE0007.OOT /usr/local/bin/LITHOSP < EE0007.INN > EE0007.OOT cp LMRK_DISPLAY01.pgm TESTFILES1/EE0007.pgm rm -f EE0013.OOT etc ...
- Conduct Standard Iteration
Here is a sample command line for running run.sh:
You can also run this in the background using this command line:
nohup sh run.sh &
- Quality Check Iterated Maplets
You must review the .OOT files to ascertain the success of the iteration process. You can detect problem landmarks using find_nofitP. Refer to that entry for more information.
- Clean Working Directory
You can use rem_script.b to clean out the directory of make_scriptP working files once iteration is complete and you have quality checked the iterated maplets.
Here is a sample rem_script.b file:
rm -f *.INN rm -f *.OOT rm -f run_script*
(Compiled by DL)