JAG2025-cloud-top-properties

This repository provides necessary code and datasets for transparency and reproducibility of our work -- using MODIS/Aqua and CALIOP/CALIPSO to estimate cloud top height, tempearture, and pressure.CTH2017065.0650.npy, CTP2017065.0650.npy, CTT2017065.0650.npy is the CTH, CTP, CTT produced using the proposed algorithm, saved as numpy npy format. time stamp: UTC 06:50, 2017 day 65.LUT2010.txt, LUT2012.txt, LUT2014.txt are used to form the total LUT. And they are generated using the code "BuildBUT.py"BuildBUT.py is the code used for generating the LUT. The meaning of each colomn of the LUTs can be found in this file. And the key code is: Xdata = Xdata[[0,1,5,6,19,21,22,23,28,30,32,33,36],:,:] and LUT[line] = np.hstack((data,doy,dem,sza, CTH,CTT,CTP))ApplyLUT.py provides a way to apply the propsed LUT algorithm.Validate.py is used to validate the proposed algorithm and MODIS operational product(MYD06) against CALIOP/CALIPSO measurements. And the key code is Validate[line] = [calipso_cth,calipso_ctt,calipso_ctp, lut_cth,lut_ctt,lut_ctp, cb_cth,cb_ctt,cb_ctp, modis_cth,modis_ctt,modis_ctp, DoY,sza,type0,cot,layer]The corresponding result of this code is Validate2024.txt.CloudPhase1.py further provides the way to discuss the accuracy of the algirothm in different conditions. And the key code is Validate[line] = [calipso_cth,calipso_ctt,calipso_ctp, lut_cth,lut_ctt,lut_ctp, modis_cth,modis_ctt,modis_ctp, DoY,sza,type0,cot,layer,cph,cfr,cgt]The corresponding result of this code is Validate2025.txt.