CyberShake output data headers - Revision history

Scottcal: /* Seismogram header */

2017-10-25T17:38:42Z

Seismogram header

2017-10-10T21:45:53Z

Sample C code

2017-10-10T21:42:57Z

Sample C code

2017-10-10T21:42:25Z

Seismogram header

2017-10-10T21:37:15Z

Sample C code

2017-10-09T19:30:23Z

2017-10-09T19:25:18Z

2017-10-09T19:24:09Z

Sample Python code

2017-10-09T19:21:36Z

Sample Python code

2017-10-09T19:19:29Z

Sample Python code

@@ Line 28: / Line 28: @@
 *DT: the timestep size used in the seismogram.
 *NT: the number of timesteps in the seismogram.
-*Comps: This tracks the components in the seismogram.  There are three flags, one for each component (X=1, Y=2, Z=4), and the flags are ANDed together to produce the value here.
+*Comps: This tracks the components in the seismogram.  There are three flags, one for each component (X=1, Y=2, Z=4), and the flags are ORed together to produce the value here.
 *Det_max_freq: the maximum frequency of the deterministic part of the seismogram.  This was 0.5 for studies before 15.4, and 1.0 for 15.4, 15.12, and 17.3.
 *Stoch_max_freq: the maximum frequency of the stochastic part of the seismogram.  For studies with no stochastic component this is -1; for studies 1.4 and 15.12 it is 10.0.

@@ Line 70: / Line 70: @@
      int z_flag = shead.comps & 4;
      printf("Version = %s\n", shead.version);
-     printf("Site = %s\n", shead.site);
+     printf("Site = %s\n", shead.site_name);
      printf("Source ID = %d\n", shead.source_id);
      printf("Rupture ID = %d\n", shead.rupture_id);
      printf("Rupture Variation ID = %d\n", shead.rup_var_id);
-     printf("DT = %f\n", dt);
+     printf("DT = %f\n", shead.dt);
-     printf("NT = %d\n", nt);
+     printf("NT = %d\n", shead.nt);
      printf("X component? %d\n", x_flag!=0);
      printf("Y component? %d\n", y_flag!=0);
      printf("Z component? %d\n", z_flag!=0);
-     printf("Maximum deterministic frequency = %f\n", det_max_freq);
+     printf("Maximum deterministic frequency = %f\n", shead.det_max_freq);
-     printf("Maximum stochastic frequency = %f\n", stoch_max_freq);
+     printf("Maximum stochastic frequency = %f\n", shead.stoch_max_freq);
      return 0;
+  }

@@ Line 40: / Line 40: @@
   #include <stdio.h>
   #include <stdlib.h>
-−
   struct seisheader {
    char version[8];
@@ Line 55: / Line 55: @@
    float stoch_max_freq;
   };
-−
   int main(int argc, char** argv) {
      if (argc<2) {

@@ Line 36: / Line 36: @@
 === Sample C code ===
-Below is sample C code to read and
+Below is sample C code to read in and print the header information.
 === Sample Python code ===

← Older revision		Revision as of 21:37, 10 October 2017
Line 36:		Line 36:
	=== Sample C code ===		=== Sample C code ===

		+	Below is sample C code to read and

	=== Sample Python code ===		=== Sample Python code ===

← Older revision		Revision as of 19:30, 9 October 2017
Line 94:		Line 94:
	== RotD header ==		== RotD header ==

		+	The format for the RotD header is exactly the same as the seismogram and PSA headers, given above. However, after the above header, there is an additional value, defined as:

		+	int num_periods;
		+
		+	This value contains the number of periods which are to follow, making it easier to parse.

	== Duration header ==		== Duration header ==
		+
		+	The format for the duration header is exactly the same as the seismogram, PSA, and RotD headers, given above.

← Older revision		Revision as of 19:25, 9 October 2017
Line 90:		Line 90:
	== PSA header ==		== PSA header ==

		+	The format for the PSA header is exactly the same as the seismogram header, given above.

	== RotD header ==		== RotD header ==

@@ Line 65: / Line 65: @@
      nt = struct.unpack('i', header_str[40:44])
      comps = struct.unpack('i', header_str[44:48])
-     x_flag = comps & 1
+     x_flag = y_flag = z_flag = False
-     y_flag = comps & 2
+    if (comps & 1)==1:
-     z_flag = comps & 4
+        x_flag = True
-     det_max_freq = struct.unpack('f', header_str[48:52])
+     if (comps & 2)==2:
-     stoch_max_freq = struct.unpack('f', header_str[52:56])
+        y_flag = True
      print "Version = %s" % version
      print "Site = %s" % site

@@ Line 55: / Line 55: @@
      header_str = fp_in.read(56)
      version = header_str[0:8]
-     if version[0:6]!="12.10":
+     if version[0:5]!="12.10":
          print "Error: version does not match expected string '12.10', aborting."
          sys.exit(2)

@@ Line 41: / Line 41: @@
 This is a Python script to read in and print the header information.
-  #!/usr/bin/env python
+  #!/usr/bin/env python2
   import sys