AtHeartEngineer/kaiju
1
0
Fork 0
mirror of https://github.com/JHUAPL/kaiju.git synced 2026-01-09 15:17:56 -05:00
Code Issues Packages Projects Releases Wiki Activity

Changed to use firsst and last times.

Browse Source
This commit is contained in:
Eric Winter
2022-05-11 15:43:15 -04:00
parent fc326a3e2e
commit ac542a123e
1 changed files with 15 additions and 10 deletions
Download Patch File Download Diff File Expand all files Collapse all files

25
scripts/datamodel/helioSatComp.py
View File

@@ -157,28 +157,34 @@ if __name__ == "__main__":
if debug:
print("gamUT = %s" % gamUT)
# Use the first non zero time as the inital MJD.
# Use the first (non zero?) time as the inital MJD.
# N.B. THIS SKIPS THE FIRST SIMULATION STEP SINCE IT TYPICALLY HAS
# gamT[0] = 0.
loc = np.argwhere(gamT > 0.0)[0][0]
t0 = gamUT[loc] # First non-0 time
t1 = gamUT[-1] # Last time
# loc = np.argwhere(gamT > 0.0)[0][0]
# t0 = gamUT[loc] # First non-0 time
t0 = gamUT[0]
t1 = gamUT[-1]
if debug:
print("t0 = %s" % t0)
print("t1 = %s" % t1)
# Compute the time interval (in integer code units) between step 1 and step 2.
deltaT = np.round(gamT[loc + 1] - gamT[loc])
if debug:
print("deltaT = %s" % deltaT)
# deltaT = np.round(gamT[loc + 1] - gamT[loc])
# if debug:
# print("deltaT = %s" % deltaT)
# <HACK>
deltaT = 3600.0 # 1 hour
# </HACK>
# Save the (float) MJD of the first step.
mjdFileStart = gamMJD[loc]
# mjdFileStart = gamMJD[loc]
mjdFileStart = gamMJD[0]
if debug:
print("mjdFileStart = %s" % mjdFileStart)
# Save the elapsed simulation time in code units of the first step.
secFileStart = gamT[loc]
# secFileStart = gamT[loc]
secFileStart = gamT[0]
if debug:
print("secFileStart = %s" % secFileStart)
@@ -221,7 +227,6 @@ if __name__ == "__main__":
data, scIds[scId], scId, mjdFileStart, secFileStart, fdir, ftag,
cmd, numSegments, keep
)
# scutils.matchUnits(data)
cdfname = os.path.join(fdir, scId + ".comp.cdf")
if os.path.exists(cdfname):
if verbose: