Added region_cut_snaphu parameter to intf_tops.csh and batch_tops.config

geofis · Aug 3, 2020 · 3c96b1f · 3c96b1f
1 parent b7247e9
commit 3c96b1f
Show file tree

Hide file tree

Showing 2 changed files with 12 additions and 5 deletions.
diff --git a/gmtsar/csh/batch_tops.config b/gmtsar/csh/batch_tops.config
@@ -87,6 +87,11 @@ near_interp = 1
 # use landmask (1 -- yes; else -- no)
 switch_land = 0
 
+# region to unwrap in radar coordinates of the resulting unwrapped phase grid
+# (leave it blank if to process the whole unwrapped phase grid)
+# example format 500/800/200/400  - OPTIONAL
+region_cut_snaphu = 
+
 #
 # Allow phase discontinuity in unrapped phase. This is needed for interferograms having sharp phase jumps.
 # defo_max = 0 - used for smooth unwrapped phase such as interseismic deformation

diff --git a/gmtsar/csh/intf_tops.csh b/gmtsar/csh/intf_tops.csh
@@ -61,6 +61,7 @@
   set threshold_snaphu = `grep threshold_snaphu $2 | awk '{print $3}'`
   set threshold_geocode = `grep threshold_geocode $2 | awk '{print $3}'`
   set region_cut = `grep region_cut $2 | awk '{print $3}'`
+  set region_cut_snaphu = `grep region_cut_snaphu $2 | awk '{print $3}'`
   set switch_land = `grep switch_land $2 | awk '{print $3}'`
   set defomax = `grep defomax $2 | awk '{print $3}'`
   set range_dec = `grep range_dec $2 | awk '{print $3}'`
@@ -192,12 +193,12 @@ if ($stage <= 2) then
 
     if ($threshold_snaphu != 0 ) then
       if ($switch_land == 1) then
-        if ($region_cut == "") then
-          set region_cut = `gmt grdinfo phase.grd -I- | cut -c3-20`
+        if ($region_cut_snaphu == "") then
+          set region_cut_snaphu = `gmt grdinfo phase.grd -I- | cut -c3-20`
         endif
         cd ../../topo
         if (! -f landmask_ra.grd) then
-          landmask.csh $region_cut
+          landmask.csh $region_cut_snaphu
         endif
         cd ../intf
         cd $ref_id"_"$rep_id
@@ -208,9 +209,9 @@ if ($stage <= 2) then
       echo "SNAPHU.CSH - START"
       echo "threshold_snaphu: $threshold_snaphu"
       if ($near_interp == 1) then
-        snaphu_interp.csh $threshold_snaphu $defomax $region_cut
+        snaphu_interp.csh $threshold_snaphu $defomax $region_cut_snaphu
       else
-        snaphu.csh $threshold_snaphu $defomax $region_cut
+        snaphu.csh $threshold_snaphu $defomax $region_cut_snaphu
       endif
       echo "SNAPHU.CSH - END"
     else
@@ -248,3 +249,4 @@ echo "END STACK OF TOPS INTERFEROGRAMS"
 echo ""
 
 
+