Skip to content

Commit

Permalink
Fixed interpolation for different angles.
Browse files Browse the repository at this point in the history
  • Loading branch information
@fnemina
fnemina committed May 6, 2020
1 parent daaa4a5 commit 4da3f38
Showing 1 changed file with 6 additions and 6 deletions.
12 changes: 6 additions & 6 deletions pyOSOAA/osoaahelpers.py
View file
Original file line number Diff line number Diff line change
Expand Up @@ -91,7 +91,7 @@ def RunWavelengths(s, wavelengths=[0.550], angle=0, output="I", taur=False, sun=
# We interpolte the values and add it to a numpy array
#f = interp1d(results['vza'], results[output])
#values = np.append(values, f(angle[idx[0]]))
values = np.append(values, np.interp(angle[idx], results['vza'], results[output]))
values = np.append(values, np.interp(angle[idx[0]], results['vza'], results[output]))
tauv = np.append(tauv, s.outputs.profileatm.tau[-1])
tau[wl] = s.outputs.profileatm
elif taur is True:
Expand All @@ -106,7 +106,7 @@ def RunWavelengths(s, wavelengths=[0.550], angle=0, output="I", taur=False, sun=
# We interpolte the values and add it to a numpy array
#f = interp1d(results['vza'], results[output])
#values = np.append(values, f(angle[idx[0]]))
values = np.append(values, np.interp(angle[idx], results['vza'], results[output]))
values = np.append(values, np.interp(angle[idx[0]], results['vza'], results[output]))
tauv = np.append(tauv, s.outputs.profileatm.tau[-1])
tau[wl] = s.outputs.profileatm
elif taur is False:
Expand All @@ -121,7 +121,7 @@ def RunWavelengths(s, wavelengths=[0.550], angle=0, output="I", taur=False, sun=
# We interpolte the values and add it to a numpy array
#f = interp1d(results['vza'], results[output])
#values = np.append(values, f(angle[idx[0]]))
values = np.append(values, np.interp(angle[idx], results['vza'], results[output]))
values = np.append(values, np.interp(angle[idx[0]], results['vza'], results[output]))
tauv = np.append(tauv, s.outputs.profileatm.tau[-1])
tau[wl] = s.outputs.profileatm
if getTau:
Expand Down Expand Up @@ -152,7 +152,7 @@ def RunWavelengths(s, wavelengths=[0.550], angle=0, output="I", taur=False, sun=
# We interpolte the values and add it to a numpy array
#f = interp1d(results['vza'], results[output])
#values = np.append(values, f(angle[idx]))
values = np.append(values, np.interp(angle[idx], results['vza'], results[output]))
values = np.append(values, np.interp(angle[idx[0]], results['vza'], results[output]))
tauv = np.append(tauv, s.outputs.profileatm.tau[-1])
tau[wl] = s.outputs.profileatm
elif taur is True:
Expand All @@ -169,7 +169,7 @@ def RunWavelengths(s, wavelengths=[0.550], angle=0, output="I", taur=False, sun=
# We interpolte the values and add it to a numpy array
#f = interp1d(results['vza'], results[output])
#values = np.append(values, f(angle[idx]))
values = np.append(values, np.interp(angle[idx], results['vza'], results[output]))
values = np.append(values, np.interp(angle[idx[0]], results['vza'], results[output]))
tauv = np.append(tauv, s.outputs.profileatm.tau[-1])
tau[wl] = s.outputs.profileatm
elif taur is False:
Expand All @@ -186,7 +186,7 @@ def RunWavelengths(s, wavelengths=[0.550], angle=0, output="I", taur=False, sun=
# We interpolte the values and add it to a numpy array
#f = interp1d(results['vza'], results[output])
#values = np.append(values, f(angle[idx]))
values = np.append(values, np.interp(angle[idx], results['vza'], results[output]))
values = np.append(values, np.interp(angle[idx[0]], results['vza'], results[output]))
tauv = np.append(tauv, s.outputs.profileatm.tau[-1])
tau[wl] = s.outputs.profileatm
if getTau:
Expand Down

0 comments on commit 4da3f38

Please sign in to comment.