Writing to External File

[DaveG]

I wonder if someone could point out why this script is not writing to a file. I can see the correct values populating the terminal but when I go to open the output file it wrote to, it was empty.

file = open(file_ls, "w")
# first line contains 3 integers: the grid dimension, nx and ny, and the total number of snapshots nt
print("%d %d %d" % (np.shape(kxx)[1], np.shape(kxx)[0], collapse_steps), open(file_ls, "a"))
# then it follows by the X coordinates of gridded data (from lines 2 to nx+1), Y coordinates of gridded data (from lines nx+1 to nx+ny+1), and the times at which the snapshots are created (from lines nx+ny+2 to nx+ny+nt+1)
for ii in np.arange(np.shape(kxx)[1]):
    print("%.8f" % (kxx[0, ii]), open(file_ls, "a"))
# for ii in arrange(shape(kyy)[0]):
for ii in range(np.shape(kyy)[0] - 1, -1, -1):
    print("%.8f" % (kyy[ii, 0]), open(file_ls, "a"))
for ii in np.arange(collapse_steps):
    mytime = collapse_times[ii] + collapse_starttime
    print("%.8f" % mytime, open(file_ls, "a"))

# all the nt snapshots of seafloor variation are written into a single column one by one from t=t1 to t=tn; for each snapshot, deltah, the data should be written row by row from the left (i=1) to the right (i=nx) from the bottom (j=1) to the top (j=nx)
# deltaH > 0 corresponds to UPLIFT
# deltaH < 0 corresponds to SUBSIDE

deltaH = elev_t0 - elev_t1
###### time - start of collapse => initial bath => zero difference ??? bit stupid ??? check w/comcot sources
for ii in np.arange(np.shape(elev_t1)[0]):
    for jj in np.arange(np.shape(elev_t1)[1]):
        print("%.8f" % 0.0, open(file_ls, "a"))

for ii in np.arange(np.shape(elev_t1)[0]):
    for jj in np.arange(np.shape(elev_t1)[1]):
        print("%.8f" % (deltaH[ii, jj]), open(file_ls, "a"))

file.close()

buran write Mar-29-2022, 07:36 AM:
Please, use proper tags when post code, traceback, output, etc. This time I have added tags for you.
See BBcode help for more info.