IN1='3scm_mf7_1.5polsub_2iterations_x_z.xyz' OUT1='thebault_fig8_comparison_w_MF7_thickness_model.ps' # cat < mag.cpt # # cpt file created by: grd2cpt -Cseis -I thick1.grd -L2/60 #COLOR_MODEL = RGB # 0 0 0 205 6 0 0 205 6 0 80 255 9 0 80 255 9 0 240 110 12 0 240 110 12 90 255 30 16 90 255 30 16 255 255 0 20 255 255 0 20 255 255 0 27 255 255 0 27 255 170 0 32 255 170 0 32 255 85 0 38 255 85 0 38 255 0 0 45 255 0 0 45 170 0 0 70 170 0 0 B 255 255 255 F 0 0 0 N 128 128 128 END # gmtset MEASURE_UNIT cm # /usr/bin/awk '{print $2,$1,$7}' $IN1 >thick.xyz surface thick.xyz -Gthick.grd -R-180/180/-60/78 -I1/1 -fx,y grdlandmask -R -Dl -I1 -N0/1/1/1/1 -Gland_mask.grd grdmath thick.grd land_mask.grd MUL = thick1.grd #grd2cpt -Cseis -I thick1.grd -L2/60 >mag.cpt grdsample thick1.grd -Gthick2.grd -I0.25 -R -Lg grdimage thick2.grd -JM22.3 -R -K \ -Cmag.cpt -Y4 > $OUT1 # pscoast -R -JM -A20000 -W -O -K -S210 -N1 -Bf10a90/f10a30 >> $OUT1 # psscale -D12/-1/15/1h -Cmag.cpt -O -L -B::/:"km": >> $OUT1 # rm thick.xyz rm thick.grd rm thick1.grd rm thick2.grd rm mag.cpt rm land_mask.grd