From 72f1ef2cb21ebb64cbd0281a90a6d97d077b2b74 Mon Sep 17 00:00:00 2001 From: Maggie Bruckner Date: Thu, 21 Nov 2024 14:59:46 -0700 Subject: [PATCH] remove vertical interpolation check --- melodies_monet/util/sat_l2_swath_utility.py | 15 --------------- 1 file changed, 15 deletions(-) diff --git a/melodies_monet/util/sat_l2_swath_utility.py b/melodies_monet/util/sat_l2_swath_utility.py index 0fbda70b..8458a788 100644 --- a/melodies_monet/util/sat_l2_swath_utility.py +++ b/melodies_monet/util/sat_l2_swath_utility.py @@ -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] @@ -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[:,:]