[cig-commits] r8460 - in seismo/2D/SPECFEM2D/trunk/DATA: . util

[walter at geodynamics.org]

Author: walter
Date: 2007-12-07 15:47:52 -0800 (Fri, 07 Dec 2007)
New Revision: 8460

Added:
   seismo/2D/SPECFEM2D/trunk/DATA/util/
   seismo/2D/SPECFEM2D/trunk/DATA/util/su_create_Paul.python
Log:
added su_create_Paul.python to 2D code


Added: seismo/2D/SPECFEM2D/trunk/DATA/util/su_create_Paul.python
===================================================================
--- seismo/2D/SPECFEM2D/trunk/DATA/util/su_create_Paul.python	2005-02-08 15:11:05 UTC (rev 8459)
+++ seismo/2D/SPECFEM2D/trunk/DATA/util/su_create_Paul.python	2007-12-07 23:47:52 UTC (rev 8460)
@@ -0,0 +1,53 @@
+#!/usr/bin/env python
+#
+#  Python code to generate Seismic Unix files for Ux and Uz and plot them
+#
+from os import *
+import os.path, string
+#
+def createSU():
+	SEM=getcwd()
+	filename=SEM+'/DATA/Par_file'
+	#
+	# Variables to be put in SU header
+	# 
+	if path.exists(filename):
+		variables=['nt','dt']
+	else:
+		print 'No Par_file found !'
+		return 
+	#     
+	# Open the file and get the lines
+	# 
+	f = file(filename,'r')
+	lignes= f.readlines()
+	f.close()
+	# Get the title
+	for ligne in lignes:
+		lsplit=string.split(ligne)
+		if lsplit[0]=='title':
+		   title=' '.join(lsplit[2:])
+		   break 
+	print '#'*50
+	print '#  SU file creation for '+title
+	print  '#'*50
+	#  Get the variables
+	for var in variables:
+		for ligne in lignes:
+			lsplit=string.split(ligne)
+			if lsplit[0]==var:
+			   exec var+'='+string.replace(string.split(''.join(ligne))[2],'d','e') 
+			   break
+	# 
+	chdir(SEM+'/OUTPUT_FILES')
+	labels='label1="Time" label2="Receivers"'
+	# Create headers ans Su file
+	ordres=['suaddhead < Ux_file.bin ns='+str(nt)+' | sushw key=dt a='+str(int(dt*1e6))+' > Ux_file.su']
+	ordres.append('suaddhead < Uz_file.bin ns='+str(nt)+' | sushw key=dt a='+str(int(dt*1e6))+' > Uz_file.su')
+	ordres.append('suxwigb < Ux_file.su perc=96 '+labels+' title=" Ux : '+title+'"&')
+	ordres.append('suxwigb < Uz_file.su xbox=600 perc=96 '+labels+' title=" Uz : '+title+'"&')
+	# 
+	for i in range(len(ordres)):
+		 system(ordres[i]) 
+if __name__=='__main__':
+	createSU()
\ No newline at end of file