[CIG-SEISMO] a problem about Mineos
潘磊
panlei7 at gmail.com
Wed Oct 27 19:16:59 PDT 2010
Dear Sim/Madam,
I'm a student from China. when studying codes of Mineos, I met some
that I couldn't understand. In minos_bran.f, we need call the
subroutine detqn which returns the count of zero crossings. I don't
understand how the count of zero crossings is computed and what
relation between it and the branch number. Is there any document to
explain this part? Could you explain it to me with the simpler
toroidal modes?
Thank you! Look forward to your reply.
Yours,
Pan Lei
some codes including this problem
detqn:
if(knsw.eq.1) then
knt=kount-2
if(jcom.ne.4) then
if(l.ne.1) then
knt=knt+1
irem=mod(knt,2)
if(irem.eq.0.and.det.gt.0.d0) knt=knt+1
if(irem.ne.0.and.det.lt.0.d0) knt=knt+1
end if
end if
end if
tprop:
if(knsw.eq.1) then
if(s(2)*f(2).le.0.d0) then
if(f(2).eq.0.d0) then
tes=-s(2)*s(1)
else
tes=f(2)*s(1)-s(2)*f(1)
endif
if(tes.lt.0.d0) kount=kount+1
if(tes.gt.0.d0) kount=kount-1
end if
end if
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://geodynamics.org/pipermail/cig-seismo/attachments/20101028/a4aceb5f/attachment.htm
More information about the CIG-SEISMO
mailing list