#! /bin/csh
gmt begin
set hit = 0
set loc=AA
set start_lat = 35
set start_lon = 103
set end_lat = 35
set end_lon = 124
set cpt = vel_new.cpt
gmt figure "$start_lat" png A1c
gmt set FONT_TITLE 10p
gmt set FONT_LABEL 6p
gmt set FONT_ANNOT_PRIMARY 5p
gmt set MAP_TICK_LENGTH_PRIMARY 3p
gmt set MAP_FRAME_PEN 0.8p
echo plot $loc slice
set dep = 5371
set reg1=0
#gawk ' $2>35.99 && $2 < 36.01 { print $3,$1,$4*100}' model_vper > lat36
gawk '{ print $3, $2, $1, $4*100}' model_vper > model.txt
gmt project -C$start_lon/$start_lat -E$end_lon/$end_lat -W-0.3/0.3 -S model.txt -Lw -Fxz > NCC_data.txt
#################### 插值,网格化 ################################
gawk '{print $1, $2, $3}' NCC_data.txt |gmt surface -Gasc.grd -I0.02/2 -T0 -C0.2 -R$start_lon/$end_lon/$dep/6371
gawk '{print $1, $2, $3}' NCC_data.txt > 1.txt
#gmt blockmean -R$start_lon/$end_lon/$dep/6371 -I5m
gmt grd2xyz asc.grd > wg_data
set roaa = `echo $start_lon $end_lon | gawk '{print ($2-$1)*0.5+$1}'`
echo $roaa "roaa"
gmt grdview asc.grd -C$cpt -JPa2.6i/$roaa'z' -R$start_lon/$end_lon/$dep/6371 -Qs
gmt basemap -JPa2.6i/$roaa'z' -R$start_lon/$end_lon/$dep/6371 -Bxa500f50 -Bya200f50+l"Depth" -BWbtr --MAP_FRAME_PEN=1p --MAP_ANNOT_MIN_SPACING=0
gmt colorbar -C$cpt -Bxa0.5+l"P-wave velocity anomaly (%)" -DJMB+w3c/0.3c -Y-0.3c
gmt end show
Could you please post a quick summary of the problem? This will save us time in having to search through your code/figure for issues.