Commit da0bcadd authored by thomas.forbriger's avatar thomas.forbriger Committed by thomas.forbriger
Browse files

bugfix: sff_cs3 has opposite sign to ZQ

This is a legacy commit from before 2015-03-01.
It may be incomplete as well as inconsistent.
See COPYING.legacy and README.history for details.


SVN Path:     http://gpitrsvn.gpi.uni-karlsruhe.de/repos/TFSoftware/trunk
SVN Revision: 3998
SVN UUID:     67feda4a-a26e-11df-9d6e-31afc202ad0c
parent 505d16d0
......@@ -28,6 +28,7 @@ c 10/02/97 added output component selection
c 27/06/97 changed coordinate units (psource) to meters (as defined by sff)
c 25/02/99 introduced spherical coordinates
c 26/02/99 sensible date setting
c 09/06/11 corrected sign for sff_sc3 in refmet_psource
c
c======================================================================
c
......@@ -393,6 +394,8 @@ c
c 27/06/97 changed coordinate units (psource) to meters (as defined by sff)
c 28/04/00 set y2k save nowdate
c 29/04/05 add 1900 not 1999 to year value
c 09/06/11 ZQ means depth below surface, sff_sc3 means height above
c surface; wrong sign was passed; corrected
c
subroutine refmet_psource(ZQ, The, typ, sff_scs, sff_sdate, nowdate,
& sff_stime, sff_stype, sff_sc1, sff_sc2, sff_sc3, radius)
......@@ -411,12 +414,12 @@ c
sff_scs='S'
sff_sc1=90.
sff_sc2=0.
sff_sc3=sngl(ZQ)
sff_sc3=-sngl(ZQ)
else
sff_scs='C'
sff_sc1=0.
sff_sc2=0.
sff_sc3=sngl(ZQ)*1000.
sff_sc3=-sngl(ZQ)*1000.
endif
call idate(nowdate)
nowdate(3)=nowdate(3)-100*int(nowdate(3)/100)
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment