[CIG-SEISMO] mineos 1.0.0 - syndat and displacements

[daniel peter]

hello,

your mineos-package is really nice. the code is very clean and it  
comes with a well written documentation. really well done,  
congratulations!

i just have a little question, maybe stupid, but would be of great  
help: how is it possible, that the syndat-program produces  
displacements instead of accelerations?

actually, you describe it in the documentation (page 30, chapter 3.4):

"...
Line 5: datatype
This defines the type of synthetic data. If datatype = 0 the output  
synthetic waveforms are accelero-
grams in nm/s2 . This is recommended native Mineos output. If  
datatype = 1 the output is velocity
waveform in nm/s, and if datatype = 2 the output is displacement in  
nm. This is done by additional
conversion of accelerograms to velocity or displacement.
...
"


looking at the syndat.f-file (with linenumbers), it seems, that this  
datatype-option is missing:

...
      47 c
      48 c    read input parameters
      49 c
      50 c....read in file name within event and moment tensor info
      51       write(*,*) '============== Program syndat  
=================='
      52       write(*,*) 'enter input CMT file name:'
      53       read(*,'(a256)') fname1
      54       write(*,*) fname1(1:lnblnk(fname1))
      55 c.... setup tensor type: moment/Mo,nodal2/Mo,nodal2
      56       write(*,*) 'enter tensor type: 0 - moment, ',
      57      *     '1 - nodal plane 1, 2 - nodal plane 2'
      58       read(*,*) itptens
      59       write(*,*) itptens
      60 c.... i/o database name
      61       write(*,*) 'enter input dbname'
      62       read(*,'(a256)') dbin
      63       write(*,*) dbin(1:lnblnk(dbin))
      64       write(*,*) 'enter output dbname'
      65       read(*,'(a256)') dbout
      66       write(*,*) dbout(1:lnblnk(dbout))
      67 c
      68 c    read source and moment
      69 c
      70       open(12,file=fname1,status='old')
...


maybe adding something like this would do it:
...
      64       write(*,*) 'enter output dbname'
      65       read(*,'(a256)') dbout
      66       write(*,*) dbout(1:lnblnk(dbout))

       write(*,*) 'enter output datatype'
       read(*,*) idout
       write(*,*) idout


      67 c
      68 c    read source and moment
      69 c
      70       open(12,file=fname1,status='old')
...


unfortunately, this seems to produce wrong results for displacements.  
you can check it with the same setup as in the DEMO2 example (see the  
attached trace-figure for LHZ component).

i would be very glad about some hints how to fix this :)

best regards,
daniel

-------------- next part --------------
A non-text attachment was scrubbed...
Name: TLY.displacement.LHZ.ASC.pdf
Type: application/pdf
Size: 21007 bytes
Desc: not available
Url : http://geodynamics.org/pipermail/cig-seismo/attachments/20070619/21f4476a/TLY.displacement.LHZ.ASC-0001.pdf
-------------- next part --------------

--
ETH Zurich
Daniel Peter
Institute of Geophysics
HPP O 13.2
Schafmattstr. 30
CH-8093 Zurich

Tel.: +41 44 633 38 61
Email: dpeter at erdw.ethz.ch
url: http://www.sg.geophys.ethz.ch/geodynamics/dpeter