[cig-commits] [commit] devel, master: compressed the large topo files that are unused by default and added some conversion and display scripts (a5f174d)

Thu Nov 6 08:15:06 PST 2014

Repository : https://github.com/geodynamics/specfem3d_globe

On branches: devel,master
Link       : https://github.com/geodynamics/specfem3d_globe/compare/bc58e579b3b0838a0968725a076f5904845437ca...be63f20cbb6f462104e949894dbe205d2398cd7f

>---------------------------------------------------------------

commit a5f174dbcf94e02d6d1cee7b8b3cfe1393b01d98
Author: Dimitri Komatitsch <komatitsch at lma.cnrs-mrs.fr>
Date:   Fri May 9 01:25:57 2014 +0200

    compressed the large topo files that are unused by default and added some conversion and display scripts


>---------------------------------------------------------------

a5f174dbcf94e02d6d1cee7b8b3cfe1393b01d98
 topo_bathy/GMTglobe.cpt                            | 53 ++++++++++++++++++++++
 .../convert_etopo2_to_specfem_topo_format.csh      |  6 +++
 ...splay_topography_of_SPECFEM3D_GLOBE_with_GMT.sh | 29 ++++++++++++
 ..._the_original_unmodified_data_file_from_the_Web |  0
 topo_bathy/smooth_topo_bathy_PPM_image.f90         |  2 +-
 5 files changed, 89 insertions(+), 1 deletion(-)

diff --git a/topo_bathy/GMTglobe.cpt b/topo_bathy/GMTglobe.cpt
new file mode 100644
index 0000000..615ac94
--- /dev/null
+++ b/topo_bathy/GMTglobe.cpt
@@ -0,0 +1,53 @@
+#       $Id: GMT_globe.cpt,v 1.1 2001/09/23 23:11:20 pwessel Exp $
+#
+# Colormap using in global relief maps
+# Bathymetry colours manually redefined for blue-shade effect and
+# new topography colour scheme for use with Generic Mapping Tools.
+# Designed by Lester M. Anderson (CASP, UK) lester.anderson at casp.cam.ac.uk
+# COLOR_MODEL = RGB
+-10000	153	0	255	-9500	153	0	255
+-9500	153	0	255	-9000	153	0	255
+-9000	153	0	255	-8500	153	0	255
+-8500	136	17	255	-8000	136	17	255  
+-8000	119	34	255	-7500	119	34	255
+-7500	102	51	255	-7000	102	51	255
+-7000	85	68	255	-6500	85	68	255
+-6500	68	85	255	-6000	68	85	255
+-6000	51	102	255	-5500	51	102	255
+-5500	34	119	255	-5000	34	119	255  
+-5000	17	136	255	-4500	17	136	255
+-4500	0	153	255	-4000	0	153	255
+-4000	27	164	255	-3500	27	164	255
+-3500	54	175	255	-3000	54	175	255
+-3000	81	186	255	-2500	81	186	255
+-2500	108	197	255	-2000	108	197	255
+-2000	134	208	255	-1500	134	208	255
+-1500	161	219	255	-1000	161	219	255
+-1000	188	230	255	-500	188	230	255
+-500	215	241	255	-200	215	241	255
+-200	241	252	255	0	241	252	255
+0	51	102	0	100	51	204	102
+100	51	204	102	200	187	228	146
+200	187	228	146	500	255	220	185
+500	255	220	185	1000	243	202	137
+1000	243	202	137	1500	230	184	88
+1500	230	184	88	2000	217	166	39
+2000	217	166	39	2500	168	154	31
+2500	168	154	31	3000	164	144	25
+3000	164	144	25	3500	162	134	19
+3500	162	134	19	4000	159	123	13
+4000	159	123	13	4500	156	113	7
+4500	156	113	7	5000	153	102	0
+5000	153	102	0	5500	162	89	89
+5500	162	89	89	6000	178	118	118
+6000	178	118	118	6500	183	147	147
+6500	183	147	147	7000	194	176	176
+7000	194	176	176	7500	204	204	204
+7500	204	204	204	8000	229	229	229
+8000	229	229	229	8500	242	242	242
+8500	242	242	242	9000	255	255	255
+9000	255	255	255	9500	255	255	255
+9500	255	255	255	10000	255	255	255
+F	255	255	255				
+B	0	0	0
+N	128	128	128
diff --git a/topo_bathy/convert_etopo2_to_specfem_topo_format.csh b/topo_bathy/convert_etopo2_to_specfem_topo_format.csh
new file mode 100755
index 0000000..51e3b39
--- /dev/null
+++ b/topo_bathy/convert_etopo2_to_specfem_topo_format.csh
@@ -0,0 +1,6 @@
+#!/bin/csh
+
+### Dimitri Komatitsch, CNRS Marseille, France, May 2014
+
+grd2xyz -Rg etopo2.grd | awk '{ print $3 }' > topo_bathy_etopo2_original_unmodified_unsmoothed.dat
+
diff --git a/topo_bathy/display_topography_of_SPECFEM3D_GLOBE_with_GMT.sh b/topo_bathy/display_topography_of_SPECFEM3D_GLOBE_with_GMT.sh
new file mode 100755
index 0000000..0e725c3
--- /dev/null
+++ b/topo_bathy/display_topography_of_SPECFEM3D_GLOBE_with_GMT.sh
@@ -0,0 +1,29 @@
+#!/bin/bash
+
+### Dimitri Komatitsch, CNRS Marseille, France, May 2014
+
+# The first plotting command needs and _only_ needs "-K"
+# All following but the last need both "-K" and "-O"
+# The last one needs _only_ "-O"
+
+ps=topography_of_SPECFEM_mesh.ps
+
+central_meridian=0
+
+rm -f $ps
+
+gmtset ANNOT_FONT_SIZE_PRIMARY 10p HEADER_FONT_SIZE 18p PLOT_DEGREE_FORMAT ddd:mm:ssF
+
+# convert the non-regular cubed-sphere grid output of SPECFEM to a regular grid for GMT
+ surface observation_grid_long_lat_topo_for_GMT.txt -Gtopography_of_SPECFEM_mesh.grd -Rd -I4m -f0x,1y -V
+
+## DK DK offset of 4.5cm
+## DK DK -Rd is an alias for -R-180/180/-90/90
+grdimage topography_of_SPECFEM_mesh.grd -CGMTglobe.cpt -Rd -JK$central_meridian/9i -Y4.5c -K -V > $ps
+
+## DK DK offset of -1.5 inch
+pscoast -Rd -JK$central_meridian/9i -B45g30:."Topography of SPECFEM3D_GLOBE mesh": -W -Dc -A1000 -U/-0.75i/-1.5i/"SPECFEM3D_GLOBE calculations by Dimitri Komatitsch" -V -O >> $ps
+
+# Clean up
+rm -f .gmt* topography_of_SPECFEM_mesh.grd
+
diff --git a/topo_bathy/note_etopo2.grd_is_the_original_unmodified_data_file_from_the_Web b/topo_bathy/note_etopo2.grd_is_the_original_unmodified_data_file_from_the_Web
new file mode 100644
index 0000000..e69de29
diff --git a/topo_bathy/smooth_topo_bathy_PPM_image.f90 b/topo_bathy/smooth_topo_bathy_PPM_image.f90
index 90356f2..f0cc700 100644
--- a/topo_bathy/smooth_topo_bathy_PPM_image.f90
+++ b/topo_bathy/smooth_topo_bathy_PPM_image.f90
@@ -32,7 +32,7 @@
 
 ! this program uses a very simple window filter
 ! it works fine because the model has a very high resolution
-! in principle it would be better to use a circular cap to take
+! in principle it would be better to use a circular cap (i.e. a Gaussian filter) to take
 ! the notion of angular distance between points into account in the filter
 ! in practice though this simple filter works perfectly fine
 

