NASA官网下载6小时频次，42个气压水平下全中国的空气温度数据，空间精度50km*60km，时间范围2006至2021年
（1）nc文件数据提取为csv文件（2）将NASA下载下的6000多站点和地级市经纬度矢量比较，得到每个地级市范围内的站点，范围内无站点，取离其最近站点数据
（3）取每一个地级市求范围内站点不同气压下大气温度的均值
（4）根据数据的特点(较高海拔处，低气压水平下本身无温度数据，因此NASA里将其的值设置的很高）参考有关文献，选有实际数据的前两层温度进行比较，第二层比第一层温度高，计为逆温
（5）每天有四个样本，四种方法判定是否为逆温天：四个样本只要有一个出现逆温即判定为逆温天；四个样本出现两天及以上逆温判定为逆温天，依次类推。之后汇总为年度数据
（6）数据里每年有354个地级市逆温数据
（7）示例为2021年部分城市数据
（8）参考文献：
Shihe Fu, V Brian Viard, Peng Zhang. Air Pollution and Manufacturing Firm Productivity: Nationwide Estimates for China. The Economic Journal, ueab033, https://doi.org/10.1093/ej/ueab033.


陈帅, 张丹丹. 空气污染与劳动生产率——基于监狱工厂数据的实证分析. 经济学季刊，2020.