Last change
on this file since 54 was
11,
checked in by lahlod, 10 years ago
|
ajout d'un fichier test Laura
|
-
Property svn:executable set to
*
|
File size:
1.2 KB
|
Line | |
---|
1 | #!/usr/bin/env python |
---|
2 | # -*- coding: utf-8 -*- |
---|
3 | import numpy as np |
---|
4 | from numpy import * |
---|
5 | |
---|
6 | def newgridns(xx, yy, zz, m, vnx, nscan, pos, z0, z1, y11): |
---|
7 | |
---|
8 | x=xx |
---|
9 | y=yy |
---|
10 | z=zz |
---|
11 | dx=0.1 |
---|
12 | dy=1 |
---|
13 | dparx=0.05 |
---|
14 | dpary=0.5 |
---|
15 | n0=0 |
---|
16 | # y11= -50 |
---|
17 | x0, x1 = -180, 180 |
---|
18 | y0, y1 = -90, 90 |
---|
19 | dx=0.1 |
---|
20 | dy=1.0 |
---|
21 | n=size(x) |
---|
22 | xvec = arange(x0,x1,dx) |
---|
23 | yvec = arange(y0,y11,dy) |
---|
24 | nx = size(xvec) |
---|
25 | ny = size(yvec) |
---|
26 | ix = zeros(n, float) |
---|
27 | iy = zeros(n, float) |
---|
28 | |
---|
29 | for kk in range(0,n): |
---|
30 | if x[kk] == -180: |
---|
31 | ix[kk] = 3600 |
---|
32 | else: |
---|
33 | ix[kk] = math.ceil((x[kk] - x0)/dx)-1 |
---|
34 | |
---|
35 | for kk in range(0,n): |
---|
36 | if y[kk] == -90: |
---|
37 | iy[kk] = 0 |
---|
38 | else: |
---|
39 | iy[kk] = math.ceil((y[kk] - y0)/dy)-1 |
---|
40 | |
---|
41 | |
---|
42 | if nscan==1: |
---|
43 | ipos = ((pos <=5) | (pos >=26)) |
---|
44 | if nscan==2: |
---|
45 | ipos = ((pos >=6) & (pos <=10))| ((pos >=21) & (pos<=25)) |
---|
46 | if nscan==3: |
---|
47 | ipos = ((pos >=11) & (pos<=20)) |
---|
48 | |
---|
49 | inx = (ix >= 0) & (ix < nx) |
---|
50 | iny = (iy >= 0) & (iy < ny) |
---|
51 | inz = (z >= z0) & (z <= z1) |
---|
52 | inn = inx & iny & inz & ipos |
---|
53 | iix = ix[inn] |
---|
54 | iiy = iy[inn] |
---|
55 | zz = z[inn] |
---|
56 | mm = m[inn] |
---|
57 | N = size(iix)# how many datapoints are left now? |
---|
58 | #calculation of the each arc legth |
---|
59 | |
---|
60 | return zz, mm |
---|
Note: See
TracBrowser
for help on using the repository browser.