|
3 | 3 |
|
4 | 4 | def vectoa(Xg,Yg,X,Y,U,V,corrlenx,corrleny,err,b=0): |
5 | 5 | ''' |
| 6 | + (Adapted from Filipe Fernandes function) |
| 7 | + |
6 | 8 | Vectoa is a vectorial objective analysis function. |
7 | 9 |
|
8 | 10 | It interpolates a velocity field (U and V, east and north velocity components) |
@@ -37,7 +39,6 @@ def vectoa(Xg,Yg,X,Y,U,V,corrlenx,corrleny,err,b=0): |
37 | 39 | PYTHON VERSION by: |
38 | 40 | Iury Sousa and Hélio Almeida - 30 May 2016 |
39 | 41 | Laboratório de Dinâmica Oceânica - IOUSP |
40 | | -
|
41 | 42 | ======================================================================''' |
42 | 43 | # making sure that the input variables aren't changed |
43 | 44 | xc,yc,x,y,u,v=Xg.copy(),Yg.copy(),X.copy(),Y.copy(),U.copy(),V.copy() |
@@ -112,6 +113,7 @@ def vectoa(Xg,Yg,X,Y,U,V,corrlenx,corrleny,err,b=0): |
112 | 113 |
|
113 | 114 | def scaloa(xc, yc, x, y, t=[], corrlenx=None,corrleny=None, err=None, zc=None): |
114 | 115 | """ |
| 116 | + (Adapted from Filipe Fernandes function) |
115 | 117 | Scalar objective analysis. Interpolates t(x, y) into tp(xc, yc) |
116 | 118 | Assumes spatial correlation function to be isotropic and Gaussian in the |
117 | 119 | form of: C = (1 - err) * np.exp(-d**2 / corrlen**2) where: |
|
0 commit comments