用python循环计算sst.mon.mean.nc几十年来每个月这个区域内大于28摄氏度的温度格点数以及平均值

要计算sst.mon.mean.nc文件中某个区域内大于28摄氏度的温度格点数以及平均值，可以按照以下步骤进行：

1. 导入需要的库和模块：

import numpy as np
from netCDF4 import Dataset

2. 打开sst.mon.mean.nc文件并读取需要的变量：

nc_file = Dataset('sst.mon.mean.nc')
sst = nc_file.variables['sst'][:]
lat = nc_file.variables['lat'][:]
lon = nc_file.variables['lon'][:]

3. 定义需要计算的区域：

lat_min = -10
lat_max = 10
lon_min = 120
lon_max = 160

lat_idx = np.where((lat >= lat_min) & (lat <= lat_max))[0]
lon_idx = np.where((lon >= lon_min) & (lon <= lon_max))[0]

sst_region = sst[:, lat_idx, :][:, :, lon_idx]

4. 循环计算每个月大于28摄氏度的温度格点数以及平均值：

n_points = []
mean_temp = []

for i in range(sst_region.shape[0]):
    sst_month = sst_region[i, :, :]
    n_points_month = np.sum(sst_month > 28)
    mean_temp_month = np.mean(sst_month[sst_month > 28])
    n_points.append(n_points_month)
    mean_temp.append(mean_temp_month)

其中，n_points列表保存了每个月大于28摄氏度的温度格点数，mean_temp列表保存了每个月大于28摄氏度的温度格点的平均值。

完整代码如下：

import numpy as np
from netCDF4 import Dataset

nc_file = Dataset('sst.mon.mean.nc')
sst = nc_file.variables['sst'][:]
lat = nc_file.variables['lat'][:]
lon = nc_file.variables['lon'][:]

lat_min = -10
lat_max = 10
lon_min = 120
lon_max = 160

lat_idx = np.where((lat >= lat_min) & (lat <= lat_max))[0]
lon_idx = np.where((lon >= lon_min) & (lon <= lon_max))[0]

sst_region = sst[:, lat_idx, :][:, :, lon_idx]

n_points = []
mean_temp = []

for i in range(sst_region.shape[0]):
    sst_month = sst_region[i, :, :]
    n_points_month = np.sum(sst_month > 28)
    mean_temp_month = np.mean(sst_month[sst_month > 28])
    n_points.append(n_points_month)
    mean_temp.append(mean_temp_month)