#!/bin/bash # IN1='brtps_d1_170_dipole_wake_tail_alt30_v14.xy' OUT1='moon_bscalar_d170_alt30_coestimation_v14.grd' # OUTF='brtps_d1_170_dipole_wake_tail_alt30_v14.ps' # # Br first # /bin/gawk '{print $1,$2,$3}' $IN1 > tmp1.xyz xyz2grd tmp1.xyz -I1/0.5 -Rg -V -Gtmp1a.grd grdsample tmp1a.grd -I0.1 -Lg -Gbr_near.grd cp br_near.grd br_far.grd grdedit br_near.grd -S -Rd -fg # Now Bt /bin/gawk '{print $1,$2,$4}' $IN1 > tmp2.xyz xyz2grd tmp2.xyz -I1/0.5 -Rg -V -Gtmp2a.grd grdsample tmp2a.grd -I0.1 -Lg -Gbt_near.grd cp bt_near.grd bt_far.grd grdedit bt_near.grd -S -Rd -fg # Finally Bscalar /bin/gawk '{print $1,$2,$6}' $IN1 > tmp3.xyz xyz2grd tmp3.xyz -I1/0.5 -Rg -V -Gtmp3a.grd grdsample tmp3a.grd -I0.1 -Lg -Gbscalar_near.grd cp bscalar_near.grd bscalar_far.grd cp bscalar_near.grd $OUT1 grdedit bscalar_near.grd -S -Rd -fg # #fi # gmtset PAPER_MEDIA letter+ LABEL_FONT_SIZE 8 ANOT_FONT_SIZE 8 BASEMAP_TYPE plain OUTPUT_DEGREE_FORMAT -DF TICK_LENGTH 0.2c MEASURE_UNIT cm PAGE_ORIENTATION portrait # cat < mag1.cpt # #COLOR_MODEL = RGB -40. 0 30 255 -3. 0 66 255 -3. 0 66 255 -2. 0 96 255 -2. 0 96 255 -1.5 0 126 255 -1.5 0 126 255 -1.2 0 163 255 -1.2 0 163 255 -1. 0 182 255 -1. 0 182 255 -0.8 0 224 255 -0.8 0 224 255 -0.7 0 255 81 -0.7 0 255 81 -0.6 45 255 15 -0.6 45 255 15 -0.5 53 255 13 -0.5 53 255 13 -0.4 112 255 0 -0.4 112 255 0 -0.3 134 255 0 -0.3 134 255 0 -0.2 171 255 0 -0.2 171 255 0 -0.1 178 255 0 -0.1 178 255 0 0.0 244 255 0 0.0 244 255 0 0.1 255 255 0 L 0.1 255 255 0 0.2 255 248 0 0.2 255 248 0 0.3 255 214 0 0.3 255 214 0 0.4 255 193 0 0.4 255 193 0 0.5 255 140 0 0.5 255 140 0 0.6 255 132 0 0.6 255 132 0 0.7 255 66 0 0.7 255 66 0 0.8 255 0 0 0.8 255 0 0 1 255 0 23 1 255 0 23 1.2 255 0 33 1.2 255 0 33 1.5 255 0 66 1.5 255 0 66 2. 255 0 132 2. 255 0 132 3. 244 0 210 3. 244 0 210 40. 244 0 255 B 0 30 255 F 244 0 255 N 190 190 190 END # cat <mag3.cpt # cpt file created by: grd2cpt -I tmp3a.grd -Cseis -Z #COLOR_MODEL = RGB # 0 0 0 205 0.3 0 72 250 0.3 0 72 250 0.4 0 208 139 0.4 0 208 139 0.5 63 250 54 0.5 63 250 54 0.6 189 255 12 0.6 189 255 12 0.7 255 255 0 0.7 255 255 0 0.9 255 221 0 0.9 255 221 0 1.2 255 144 0 1.2 255 144 0 1.8 255 68 0 1.8 255 68 0 5 246 0 0 5 246 0 0 12 170 0 0 B 0 0 255 F 0 0 0 N 128 128 128 END # cat < reiner_gamma.xy 301.5 7.5 END # psbasemap -R0/18/0/31 -Jx1 -B0 -U -K -X1 -Y2 > $OUTF pstext -R -Jx -O -K << EOF >> $OUTF 15 3.5 16 0 4 LT Br EOF pstext -R -Jx -O -K << EOF >> $OUTF 15 12 16 0 4 LT Btheta EOF pstext -R -Jx -O -K << EOF >> $OUTF 15 20 16 0 4 LT Bmag EOF # Br psbasemap -Rd -JA0/0/7.5 -X1.5 -Y3.5 \ -Ba30f30g30/a30f30g30WSen -O -K -V >> $OUTF grdimage br_near.grd -Cmag1.cpt -JA -Rd -O -K \ -Ba30f30g30/a30f30g30WSen >> $OUTF psxy reiner_gamma.xy -Rd -JA -Sc0.1c -O -K >> $OUTF psbasemap -Rg -JA180/0/7.5 -X8 \ -Ba30f30g30/a30f30g30WSen -O -K -V >> $OUTF grdimage br_far.grd -Cmag1.cpt -JA -Rg -O -K \ -Ba30f30g30/a30f30g30WSen >> $OUTF psxy reiner_gamma.xy -Rg -JA -Sc0.1c -O -K >> $OUTF # Bt psbasemap -Rd -JA0/0/7.5 -X-8 -Y8.25 \ -Ba30f30g30/a30f30g30Wsen -O -K -V >> $OUTF grdimage bt_near.grd -Cmag1.cpt -JA -Rd -O -K \ -Ba30f30g30/a30f30g30WSen >> $OUTF psxy reiner_gamma.xy -Rd -JA -Sc0.1c -O -K >> $OUTF psbasemap -Rg -JA180/0/7.5 -X8 \ -Ba30f30g30/a30f30g30Wsen -O -K -V >> $OUTF grdimage bt_far.grd -Cmag1.cpt -JA -Rg -O -K \ -Ba30f30g30/a30f30g30WSen >> $OUTF psxy reiner_gamma.xy -Rg -JA -Sc0.1c -O -K >> $OUTF # B Scalar psbasemap -Rd -JA0/0/7.5 -X-8 -Y8.25 \ -Ba30f30g30/a30f30g30Wsen -O -K -V >> $OUTF grdimage bscalar_near.grd -Cmag3.cpt -JA -Rd -O -K \ -Ba30f30g30/a30f30g30:."Near side":WSen >> $OUTF psxy reiner_gamma.xy -Rd -JA -Sc0.1c -O -K >> $OUTF psbasemap -Rg -JA180/0/7.5 -X8 \ -Ba30f30g30/a30f30g30Wsen -O -K -V >> $OUTF grdimage bscalar_far.grd -Cmag3.cpt -JA -Rg -O -K \ -Ba30f30g30/a30f30g30:."Far side":WSen >> $OUTF psxy reiner_gamma.xy -Rg -JA -Sc0.1c -O -K >> $OUTF psscale -Cmag1.cpt -D-5/-17.8/6/0.75h -B::/:'nT vec': -E -P -L -S -O -K -V >> $OUTF psscale -Cmag3.cpt -D4/-17.8/6/0.75h -B::/:'nT B': -E -P -L -S -O -V >> $OUTF # rm tmp* mag1.cpt mag3.cpt reiner_gamma.xy b*grd