Dec 2013
Copied and modified for Intermediate Scale Work

Copy full directory
/scratch/sft2/lucia/grid_gen/arctic/double_offcet_var

mv output output_var_param (variable param for each poly)
cp mesh_adjoin.py mesh_adjoin_var.py

edit mesh_adjoin.py so every poly has same offset param
      # Set is_params_uniform to 0
      # when the params are different for every plygon
      is_params_uniform = 1

run_second.sh
    rm output/*
    cp output_surf_mesh/* output/
    ./mesh_adjoin.py

CHANGE: Keep each poly set in this directory with new names for re-use 
    move output to out_*new_name
    merge all polys into single mesh at elevation 0.  poly_all_quads.inp
    copy all lgi files into the poly set directory, differ for each set 

CHANGE: move dir to top for use in mesh versions
    keep lidar data in input dir and re-use instead of re-write for each poly set
    Process lidar data if not already done
    Use input/lidar_site3.inp and write lidar_site3_zero (zero elev with z_save) 

Run LaGriT commands (this all goes in mesh work directory pointing to appropriate poly sheet)
    CHANGE: instead of 'process_lidar.lgi', 'merge_polys.lgi', 'refine.lgi', 'interpolate.lgi'
    use processed lidar data
    interpolate elevations to merged poly surface



Problem Polys
aquifer.lanl.gov% cp lobster_var_param/param_08.txt work_single_polys
aquifer.lanl.gov% cp lobster_var_param/param_45.txt work_single_polys
aquifer.lanl.gov% cp lobster_var_param/param_56.txt work_single_polys
aquifer.lanl.gov% cp lobster_var_param/param_21.txt work_single_polys