OCIAM一维反演项目专注于通过先进的数学模型与算法,对地球物理数据进行精确的一维反演分析,旨在深入理解地壳结构和组成。
```fortran
subroutine startup(itform, descr, modfil, datfil, maxitr,
& tolreq, iruf, nit, pmu, rlast, tlast,
& ifftol)
C-----------------------------------------------------------------------
C OCCAM 2.0: Steven Constable IGPP/SIO La Jolla CA 92093-0225
C Subroutine Revision 1.04, 18 May 1992
include imp.inc
character*80 itform, descr, modfil, datfil
real pmu, rlast, tlast
integer maxitr, tolreq, iruf, nit
logical ifftol
C on input:
C nothing except the parameters in the calling sequence.
C on output:
C sets up OCCAM for a new run or continues from an old one.
call itrin(itform, descr, modfil, datfil,
& maxitr, tolreq, iruf, nit,
& pmu, rlast, tlast,
& ifftol)
return
end
C-----------------------------------------------------------------------
subroutine run(nitmax)
C OCCAM 2.0: Steven Constable IGPP/SIO La Jolla CA 92093-0225
C Subroutine Revision 1.04, 18 May 1992
include imp.inc
integer nitmax, i, ifftol
real pmu, rlast, tlast
character*80 itform, descr, modfil, datfil
C on input:
C nothing except the parameters in the calling sequence.
C on output:
C performs a new run of OCCAM.
do 1 i = 1,nitmax
call startup(itform, descr,
& modfil, datfil,
& maxitr, tolreq, iruf,
& nit, pmu, rlast, tlast,
& ifftol)
if (ifftol) goto 20
1 continue
call itin(itform, descr)
20 return
end
C-----------------------------------------------------------------------
subroutine stopit(pmnew, tolnw)
C OCCAM 2.0: Steven Constable IGPP/SIO La Jolla CA 92093-0225
C Subroutine Revision 1.04, 18 May 1992
include imp.inc
real pmnew(npp), tolnw
character*80 itform, descr, modfil, datfil
C on input:
C nothing except the parameters in the calling sequence.
C on output:
C stops an OCCAM run.
call itout(itform, descr,
& modfil, datfil,
& maxitr, tolreq, iruf,
& nit, pmu, rlast, tlast,
& ifftol)
return
end
C-----------------------------------------------------------------------
subroutine startup2(pmnew)
C OCCAM 2.0: Steven Constable IGPP/SIO La Jolla CA 92093-0225
C Subroutine Revision 1.04, 18 May 1992
include imp.inc
real pmnew(npp)
C on input:
C nothing except the parameters in the calling sequence.
C on output:
C sets up OCCAM for a new run or continues from an old one.
call itout(itform, descr,
& modfil, datfil,
& maxitr, tolreq, iruf,
& nit, pmu, rlast, tlast,
& ifftol)
return
end
C-----------------------------------------------------------------------
subroutine run2(nitmax)
C OCCAM 2.0: Steven Constable IGPP/SIO La Jolla CA 92093-0225
C Subroutine Revision 1.04, 18 May 1992
include imp.inc
integer nitmax, i
C on input:
C nothing except the parameters in the calling sequence.
C on output:
C performs a new run of OCCAM.
do 1 i = 1,nitmax
call startup(itform, descr,
& modfil, datfil,
& maxitr, tolreq, iruf,
& nit, pmu, rlast, tlast,
& ifftol)
1 continue
return
end
C-----------------------------------------------------------------------
subroutine stopit2(pmnew)
C OCCAM 2.0: Steven Constable IGPP/SIO La Jolla CA