source: anuga_work/production/pt_hedland_2006/convert2latlon.py @ 5498

Last change on this file since 5498 was 4949, checked in by sexton, 16 years ago

(a) minor updates for Cairns example in manual (b) scripts for stratification investigation for other model areas

File size: 1.4 KB
Line 
1"""
2Convert Easting Northing to lon/lat for ultimately get kml file
3"""
4
5def alter_file(filename,filenameout):
6
7    from anuga.coordinate_transforms.lat_long_UTM_conversion import UTMtoLL
8   
9    fid = open(filename)
10    lines = fid.readlines()
11    fid.close()
12   
13    fid_out = open(filenameout, 'w')
14    s = '%s,%s,%s\n' %('lon','lat','stage')
15    fid_out.write(s)
16    for line in lines[1:]:
17       fields = line.split(',')
18       easting = float(fields[0])
19       northing = float(fields[1])
20       stage = float(fields[3])
21       lat, lon = UTMtoLL(northing,easting, 50)
22       s = '%f,%f,%f\n' %(lon,lat,stage)
23       fid_out.write(s)
24
25    fid_out.close() 
26
27    return
28
29import project
30from os import sep
31
32timedir='20070613_044204_run_final_3.6_exmouth_nbartzis'
33#timedir='20070613_061053_run_final_3.6_onslow_nbartzis'
34
35#filename = project.outputdir + timedir + sep + 'stage0_interp.csv'
36filename = project.outputdir + timedir + sep + 'stage5_interp.csv'
37#filename = project.outputdir+ timedir + sep + 'stage20_interp.csv'
38#filename = project.outputdir+  timedir + sep + 'stage50_interp.csv'
39#filename_out = project.outputdir + timedir + sep + 'stage0_convert.csv'
40filename_out = project.outputdir + timedir + sep + 'stage5_convert.csv'
41#filename_out = project.outputdir + timedir + sep + 'stage20_convert.csv'
42#filename_out = project.outputdir + timedir + sep + 'stage50_convert.csv'
43alter_file(filename,filename_out)
Note: See TracBrowser for help on using the repository browser.