File: mapgen_clip

package info (click to toggle)
magics%2B%2B 4.5.3-1
  • links: PTS, VCS
  • area: main
  • in suites: bullseye
  • size: 94,096 kB
  • sloc: cpp: 164,444; xml: 27,055; perl: 6,352; ansic: 4,097; python: 3,761; sh: 739; javascript: 703; f90: 278; makefile: 90
file content (50 lines) | stat: -rwxr-xr-x 1,398 bytes parent folder | download | duplicates (2)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
#!/usr/bin/python3
#
# Copyright 2007 European Centre for Medium-Range Weather Forecasts (ECMWF)
# 
# Licensed under the Apache License, Version 2.0 (the "License"); 
# you may not use this file except in compliance with the License. 
# You may obtain a copy of the License at 
# 
# 	http://www.apache.org/licenses/LICENSE-2.0
# 
# Unless required by applicable law or agreed to in writing, software 
# distributed under the License is distributed on an "AS IS" BASIS, 
# WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. 
# See the License for the specific language governing permissions and 
# limitations under the License.
#
#
# This script is provided as help to clip / reduce size of mapgen data files.
#
import sys

if len(sys.argv) < 6:
    print ("Not enough parameters!")
    print ("mapgen_clip in.mapgen out.mapgen lon_min lon_max lat_min lat_max")
    sys.exit()

lon_min = float(sys.argv[3])
lon_max = float(sys.argv[4])
lat_min = float(sys.argv[5])
lat_max = float(sys.argv[6])

fin = open(sys.argv[1], 'r')
fou = open(sys.argv[2], 'w')

begin = ""

for line in fin:
  if line == "# -b\n":
    begin=line
  else:
    s = line.split()
    if float(s[0]) > lon_min and float(s[0]) < lon_max:
     if float(s[1]) > lat_min and float(s[1]) < lat_max:
      if begin != "":
        fou.write(begin)
        begin = ""
      fou.write(line)

fin.close()
fou.close()