Skip to content
Snippets Groups Projects
Commit da7be4dc authored by Jacques Grelet's avatar Jacques Grelet
Browse files

first version of write ascii

parent 5e05bdb7
No related branches found
No related tags found
No related merge requests found
...@@ -3,23 +3,8 @@ ...@@ -3,23 +3,8 @@
import os import os
import tools import tools
def writeAscii(cfg, device, fe, r, variables_1D):
if not os.path.exists(cfg['global']['ascii']):
os.makedirs(cfg['global']['ascii'])
fileName = "{}/{}.{}".format(cfg['global']['ascii'],
cfg['cruise']['cycleMesure'], device.lower())
print('writing header file: {}'.format(fileName))
writeHeader(fileName, cfg, fe, r, variables_1D, device.lower())
fileName = "{}/{}_{}".format(cfg['global']['ascii'],
cfg['cruise']['cycleMesure'], device.lower())
print('writing data file: {}'.format(fileName), end='')
writeData(fileName, cfg, fe, r, variables_1D)
print(' done...')
def writeHeader(hdrFile, cfg, fe, r, variables_1D, device): def writeHeader(hdrFile, cfg, device, fe, r):
f = open(hdrFile, 'w') f = open(hdrFile, 'w')
# first line, header ex: # first line, header ex:
# PIRATA-FR30 THALASSA IRD SBE911+ 09P-1263 BOURLES # PIRATA-FR30 THALASSA IRD SBE911+ 09P-1263 BOURLES
...@@ -46,8 +31,41 @@ def writeHeader(hdrFile, cfg, fe, r, variables_1D, device): ...@@ -46,8 +31,41 @@ def writeHeader(hdrFile, cfg, fe, r, variables_1D, device):
f.close() f.close()
def writeData(dataFile, cfg, fe, r, variables_1D): def writeData(dataFile, cfg, device, fe, r):
f = open(dataFile, 'w') f = open(dataFile, 'w')
f.write("ASCII file OK !")
f.write("{} {} {} {} {} {}\n".format(cfg['cruise']['cycleMesure'],
cfg['cruise']['plateforme'], cfg['cruise']['institute'],
cfg[device]['typeInstrument'], cfg[device]['instrumentNumber'],
cfg['cruise']['pi']))
for i in range(fe.n):
f.write("{:0>5d} {:>-5} {:>8.2} {:>2.5} {:>3.5} {} {:>4.1}\n".format(fe['PROFILE'][i],
-1, fe['TIME'][i], fe['LATITUDE'][i], fe['LONGITUDE'][i],
tools.julian2dt(fe['TIME'][i]).strftime("%Y%m%d%H%M%S"), fe['BATH'][i]))
for m in range(fe.m):
#print(fe['PRES'])
if fe['PRES'][i][m] < 1e35:
f.write("{:0>5} ".format(fe['PROFILE'][i]))
for k in fe.keys:
f.write("{:0>5.2f} ".format( fe[k][i][m]))
f.write("\n")
f.close() f.close()
def writeAscii(cfg, device, fe, r):
if not os.path.exists(cfg['global']['ascii']):
os.makedirs(cfg['global']['ascii'])
fileName = "{}/{}.{}".format(cfg['global']['ascii'],
cfg['cruise']['cycleMesure'], device.lower())
print('writing header file: {}'.format(fileName))
writeHeader(fileName, cfg, device.lower(), fe, r, )
fileName = "{}/{}_{}".format(cfg['global']['ascii'],
cfg['cruise']['cycleMesure'], device.lower())
print('writing data file: {}'.format(fileName), end='')
writeData(fileName, cfg, device.lower(), fe, r)
print(' done...')
\ No newline at end of file
...@@ -19,7 +19,7 @@ EXIT_FAILURE = 1 ...@@ -19,7 +19,7 @@ EXIT_FAILURE = 1
# typeInstrument is a dictionary as key: files extension # typeInstrument is a dictionary as key: files extension
typeInstrument = {'CTD': ('cnv', 'CNV'), 'XBT': ( typeInstrument = {'CTD': ('cnv', 'CNV'), 'XBT': (
'EDF', 'edf'), 'LADCP': ('lad', 'LAD'), 'TSG': 'COLCOR'} 'EDF', 'edf'), 'LADCP': ('lad', 'LAD'), 'TSG': 'COLCOR'}
variables_1D = ['TIME', 'LATITUDE', 'LONGITUDE','BATH'] #variables_1D = ['TIME', 'LATITUDE', 'LONGITUDE','BATH']
ti = typeInstrument # an alias ti = typeInstrument # an alias
filesBrowsePosition_row = 2 filesBrowsePosition_row = 2
filesBrowsePosition_column = 1 filesBrowsePosition_column = 1
...@@ -320,10 +320,10 @@ if __name__ == "__main__": ...@@ -320,10 +320,10 @@ if __name__ == "__main__":
# print(fe.disp()) # print(fe.disp())
# write ASCII hdr and data files # write ASCII hdr and data files
ascii.writeAscii(cfg, device, fe, r, variables_1D) ascii.writeAscii(cfg, device, fe, r)
# write the NetCDF file # write the NetCDF file
netcdf.writeNetCDF(cfg, device, fe, r, variables_1D) netcdf.writeNetCDF(cfg, device, fe, r)
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment