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

is working

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: 1757
SVN UUID:     67feda4a-a26e-11df-9d6e-31afc202ad0c
parent ad011381
c this is <readanigemini.f>
c ----------------------------------------------------------------------------
c ($Id: readanigemini.f,v 1.2 2005-06-16 11:43:30 tforb Exp $)
c ($Id: readanigemini.f,v 1.3 2005-06-16 12:56:23 tforb Exp $)
c
c Copyright (c) 2005 by Thomas Forbriger (BFO Schiltach)
c
......@@ -16,18 +16,19 @@ c
c
c this routines reads a gemini model file
c
character filename*(*), text*72
character filename*(*), text*(*)
integer nlayer, maxlayer, lu
double precision qm(maxlayer), qk(maxlayer), rho(maxlayer, 4)
double precision vpv(maxlayer, 4), vsv(maxlayer, 4)
double precision vph(maxlayer, 4), vsh(maxlayer, 4)
double precision eta(maxlayer, 4)
double precision rb(0:maxlayer)
integer iflso(maxlayer), nco(maxlayer), nlay
integer iflso(maxlayer), nco(maxlayer)
cE
character cnlay*2,form*11
double precision fref, rboc
integer i,j,ifanis
character junk*80
c
open(lu,file=filename,status='old',err=99)
c
......@@ -61,8 +62,10 @@ c------------------------------------------------------
c Read header with comments and print them on stdout
c ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
1 read(lu,'(a72)') text
if (text(1:1).eq.'#') then
text='-'
1 read(lu,'(a72)') junk
if (text.eq.'-') text=junk
if (junk(1:1).eq.'#') then
c print '(a72)', text
goto 1
endif
......@@ -109,10 +112,10 @@ c S-velo-hori, Qmu, Qkappa, eta
enddo
read(lu,*) rb(nlay) ! Get Earth radius
read(lu,*) rb(nlayer) ! Get Earth radius
c Check for ocean.
if(iflso(nlay).eq.1) stop '~~~~ I don`t like oceans. ~~~~'
if(iflso(nlayer).eq.1) stop '~~~~ I don`t like oceans. ~~~~'
rboc = rb(nloc) ! Radius of outer core
close(lu)
......
Supports Markdown
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