#!/bin/csh # map of the Takayama site # Hodaka: since 1968, 1150m; # Tochio: since 1978.11.20 set OUT=site_map.ps set WEST=136.5 set EAST=139 set NORTH=37.5 set SOUTH=34.0 set MARKSIZE1=0.2 set MARKSIZE2=0.2 set MARKSIZE3=0.4 set FONTSIZE=15 set GMTHOME=/usr/local/GMT4.0 gmtset LABEL_FONT_SIZE 22 gmtset ANNOT_FONT_SIZE_PRIMARY 22 gmtset ANNOT_FONT_SIZE_SECONDARY 22 ### site data #### cat site_JMA.txt | awk '{print $8+$9/60.0,$6+$7/60.0}' > site_rain.txt cat site_JMA.txt | grep »Í | awk '{print $8+$9/60.0,$6+$7/60.0}' > site.txt cat site_JMA.txt | grep ´± | awk '{print $8+$9/60.0,$6+$7/60.0}' >> site.txt echo "137.57416 36.25750" > site_Hodaka.txt echo "137.25333 36.1550" > site_TakayamaJMA.txt echo "138.23167 35.3283" > site_Ikawa.txt echo "138.50333 35.915" > site_Kawakami.txt echo "138.34778 36.5236" > site_Sugadaira.txt echo "137.42311 36.1462" > site_Takayama.txt ### color pallette ### cat < height.cpt 0 156 203 219 500 183 234 142 500 183 234 142 1000 230 224 152 1000 230 224 152 1500 248 234 183 1500 248 234 183 2000 255 242 212 2000 255 242 212 2500 255 251 240 2500 255 251 240 3000 255 255 255 3000 255 255 255 3500 255 255 255 3500 255 255 255 4000 255 255 255 B 236 140 255 F 255 255 255 N 116 163 179 EOF ### SCALE ### set SCALEX = `echo "$EAST * 0.4 + $WEST * 0.6" | bc -l` set SCALEY = `echo "$SOUTH * 0.90 + $NORTH * 0.10" | bc -l` set SCALELAT = `echo "( $NORTH + $SOUTH ) / 2" | bc -l` ### main routine ### psbasemap -Jm4.0 -R${WEST}/${EAST}/${SOUTH}/${NORTH} -Ba1f1 -X5.0 -Y5.0 -P -K -V >! $OUT grdimage -Jm $GMTHOME/share/dbase/EastAsia.grd -I$GMTHOME/share/dbase/EastAsia_I.grd -R${WEST}/${EAST}/${SOUTH}/${NORTH} -P -Cheight.cpt -K -O >> $OUT pscoast -Jm -R -Df -P -V -K -W2/0/0/0 -S200/200/255 -Lf${SCALEX}/${SCALEY}/${SCALELAT}/100 -O >> $OUT psxy -Jm -St${MARKSIZE2} -G30/30/255 -R -V -K -O < site_rain.txt >> $OUT psxy -Jm -Sc${MARKSIZE2} -G255/120/0 -R -V -K -O < site.txt >> $OUT psxy -Jm -Sc${MARKSIZE3} -G255/0/200 -R -V -K -O < site_Hodaka.txt >> $OUT #psxy -Jm -Sc${MARKSIZE3} -G250/0/0 -R -V -K -O < site_TakayamaJMA.txt >> $OUT #psxy -Jm -Sc${MARKSIZE2} -G250/250/50 -R -V -K -O < site_TakayamaJMA.txt >> $OUT psxy -Jm -Sc${MARKSIZE3} -G255/0/200 -R -V -K -O < site_Ikawa.txt >> $OUT psxy -Jm -Sc${MARKSIZE3} -G255/0/200 -R -V -K -O < site_Kawakami.txt >> $OUT psxy -Jm -Sc${MARKSIZE3} -G255/0/200 -R -V -K -O < site_Sugadaira.txt >> $OUT psxy -Jm -Sc${MARKSIZE3} -G255/0/200 -R -V -K -O < site_Takayama.txt >> $OUT psxy -Jm $GMTHOME/share/map.pref -R -M -W3 -O -W1p/50/50/50 >> $OUT ### main routine ### set OUT2=site_map_big.ps psbasemap -Jm2.0 -R134/141/33/39 -Ba1f1 -X5.0 -Y5.0 -P -K -V >! $OUT2 grdimage -Jm $GMTHOME/share/dbase/EastAsia.grd -I$GMTHOME/share/dbase/EastAsia_I.grd -R -P -Cheight.cpt -K -O >> $OUT2 pscoast -Jm -R -Df -P -V -K -W2/0/0/0 -S200/200/255 -O >> $OUT2 psxy -Jm -Sc${MARKSIZE3} -G250/0/0 -R -V -K -O < site_Hodaka.txt >> $OUT2 psxy -Jm -Sc${MARKSIZE2} -G250/250/50 -R -V -K -O < site_Hodaka.txt >> $OUT2 psxy -Jm $GMTHOME/share/map.pref -R -M -W3 -O >> $OUT2 #### cat site_JMA.txt | awk '136.5<=$8 && $8<139{print $3,$6+$7/60.0,$10}' > site_lat_alt_rain.txt cat site_JMA.txt | grep »Í | awk '136.5<=$8 && $8<139{print $3,$6+$7/60.0,$10}' > site_lat_alt.txt cat site_JMA.txt | grep ´± | awk '136.5<=$8 && $8<139{print $3,$6+$7/60.0,$10}' >> site_lat_alt.txt echo "Hodaka 36.25750 1150" > site_lat_alt_Hodaka.txt echo "Takayama(JMA) 36.1550 560" > site_lat_alt_TakayamaJMA.txt echo "Ikawa 35.3283 1175" > site_lat_alt_Ikawa.txt echo "Kawakami 35.915 1500" > site_lat_alt_Kawakami.txt echo "Sugadaira 36.5236 1310" > site_lat_alt_Sugadaira.txt echo "Takayama 36.1462 1420" > site_lat_alt_Takayama.txt cat site_lat_alt_Hodaka.txt site_lat_alt_Ikawa.txt site_lat_alt_Kawakami.txt site_lat_alt_Sugadaira.txt site_lat_alt_Takayama.txt > site_lat_alt_univ.txt gnuplot<