Skip to content

Commit

Permalink
remove vertical interpolation check
Browse files Browse the repository at this point in the history
  • Loading branch information
mbruckner-work authored Nov 21, 2024
1 parent ee0880a commit 72f1ef2
Showing 1 changed file with 0 additions and 15 deletions.
15 changes: 0 additions & 15 deletions melodies_monet/util/sat_l2_swath_utility.py
Original file line number Diff line number Diff line change
Expand Up @@ -291,8 +291,6 @@ def cal_amf_wrfchem(scatw, wrfpreslayer, tpreslev, troppres, wrfno2layer_molec,
vertical_pres = []
vertical_scatw = []
vertical_wrfp = []

runvertcheck = True

for llb in range(len(lb[0])):
yy = lb[0][llb]
Expand All @@ -303,19 +301,6 @@ def cal_amf_wrfchem(scatw, wrfpreslayer, tpreslev, troppres, wrfno2layer_molec,
f = interpolate.interp1d(np.log10(vertical_pres[:]),vertical_scatw[:], fill_value="extrapolate")# relationship between pressure to avk
wrfavk[yy,xx,:] = f(np.log10(vertical_wrfp[:])) #wrf-chem averaging kernel

# MEB: check added for vertical interp
if runvertcheck and (~np.isnan(vertical_pres)).all():
from matplotlib import pyplot as plt
plt.plot(vertical_pres,label='Obs pressure')
plt.plot(vertical_wrfp,label='Model pressure')
plt.legend()
plt.ytitle("Pressure")
plt.xtitle('z index')
plt.title('Tropomi pair check: Pressure dim')
plt.show()
print(wrfavk)
runvercheck = False

for l in range(nz-1):
# check if it's within tropopause
preminus[:,:] = wrfpreslayer[:,:,l] - troppres[:,:]
Expand Down

0 comments on commit 72f1ef2

Please sign in to comment.