使用NCL计算锋生函数（修改：已转为python计算）

mozhiting

zengmin 发表于 2021-6-23 21:20
请问楼主用到的数据是哪些，可以提供一下数据列表和链接吗

我用的是ERA5的温度和比湿数据，下载地址：https://cds.climate.copernicus.e ... levels?tab=overview

zengmin

mozhiting 发表于 2021-6-30 09:47
我用的是ERA5的温度和比湿数据，下载地址：https://cds.climate.copernicus.eu/cdsapp#!/dataset/reanaly ...

感谢楼主的解答{:eb511:}
def getVars(var_name,level_index):
        if var_name == 'theta':
                fvar = nc.Dataset('../era5_theta_theta_e.2020.nc')
                var = fvar.variables['pot_tem'][:,level_index,:,:]
        elif var_name == 'theta_e':
                fvar = nc.Dataset('../era5_theta_theta_e.2020.nc')
                var = fvar.variables['pot_tem_equiv'][:,level_index,:,:]
        else:
                fvar = nc.Dataset('../era5_' + var_name + '_pressures.2020.nc')
                var = fvar.variables[var_name][:,level_index,120:361,400:721]
        var = np.array(var)
        return var
这一部分是θ，Θe 和什么的读入？
比湿是哪个变量呢？

mozhiting

zengmin 发表于 2021-7-2 15:14
感谢楼主的解答
def getVars(var_name,level_index):
        if var_name == 'theta':

theta和theta_e是自己算出来的，算theta_e的时候需要用到比湿，算锋生的时候用theta或者theta_e

zengmin

mozhiting 发表于 2021-7-20 09:44
theta和theta_e是自己算出来的，算theta_e的时候需要用到比湿，算锋生的时候用theta或者theta_e

非常感谢楼主的回答

坨坨and便便

你好我下载了你的代码， 你这个代码是不是只是计算锋生然后生成另一个nc文件呀，没有做图部分

坨坨and便便

还有一个问题是，你下载的数据有几层呀，能不能只下载单层，比如想计算700hpa，我只下载700hpa的数据计算

river

下载学习，感谢分享

mozhiting

zengmin 发表于 2021-6-23 21:20
请问楼主用到的数据是哪些，可以提供一下数据列表和链接吗

数据就是ERA5的比湿和温度

wade

mozhiting 发表于 2021-7-20 09:44
theta和theta_e是自己算出来的，算theta_e的时候需要用到比湿，算锋生的时候用theta或者theta_e

算出来是非地转湿q矢量锋生还是经典q矢量锋生？

blueskyhpc

感谢分享，