Changeset 14836

Timestamp:

05/01/13 16:21:41 (12 years ago)

Author:

adhikari

Message:

CHG: no more t0=0 restriction. Default numtimes=5, for Ivins benchmarks

File:

• issm/trunk-jpl/src/c/modules/GiaDeflectionCorex/stot.f (modified) (3 diffs)

issm/trunk-jpl/src/c/modules/GiaDeflectionCorex/stot.f

@@ -2 +2 @@
       implicit double precision (a-h,o-z)
       parameter (Nafter = 1)
-      parameter (Ntime = 2)
+      parameter (Ntime = 5)
+      parameter (Ntimm = Ntime-1)
       parameter (nhank = 1024)
       parameter (Ntimp = Ntime + Nafter)
       double precision decay(2)
       double precision pset(7)
-      double precision time(Ntimp),bi(Ntime),dmi(Ntime)
+      double precision time(Ntimp),bi(Ntimm),dmi(Ntimm)
       double precision dekay1(nhank),dekay2(nhank),amp0(nhank),
      1amp1(nhank),amp2(nhank)
+c ::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
       common /blockm/ dekay1,dekay2,amp0,amp1,amp2
       common /blockt/ time,bi,dmi
@@ -42 +44 @@
       t = time(Ntimp)
       sumb = 0.0d0
-      ta = 0.0d0
-      tb = time(1)
+      ta = time(1)
+      tb = time(2)
 c ::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
-      do 97 i = 1,Ntime
+      do 97 i = 1,Ntimm
 c ::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
 c note that this "if" prevents adding load
@@ -64 +66 @@
    38 call qwise(t,ta,qjadon,xi0,xi1,xi2,slope,ycept,decay,bhaq)
    39 sumb = bhaq + sumb
-      ta = time(i)
-      tb = time(i + 1)
+      ta = time(i + 1)
+      tb = time(i + 2)
    97 continue
       fltng = sumb