python循环读取NCEP提供的日平均地表辐射通量nc文件并存入文本文件中

a1s2d3

先把要读取的nc文件放在同一文件夹中，再用python循环读取，读取文件为从NCEP下载的向下地表平均日长波辐射通量，并计算其月平均值，存放至一个文本文件中
NCEP日平均地表向下长波辐射通量下载网址：https://psl.noaa.gov/cgi-bin/db_search/DBListFiles.pl?did=195&tid=98666&vid=1234

1. import netCDF4 as nc
   
2. import numpy as np
   
3. from os import listdir
   
4. from os.path import join
   
5.                
   
6. data = np.zeros((5,12,94,192)) #建立年份*月份*纬向格点数*经向格点数的数列
   
7. path = r'E:' #nc文件存放文件夹
   
8. 
   
9. a = 0
   
10. for filename in listdir(path):
    
11.     #读取nc文件
    
12.     file = join(path,filename)
    
13.     nc_obj = nc.Dataset(file)
    
14.     x = nc_obj.variables['dlwrf'][:] #dlwrf为nc文件中的变量名
    
15.     nc_obj.close()
    
16.     #计算月平均值
    
17.     num = 0
    
18.     for m in range(12):
    
19.         if m in[1,3,5,7,8,10,12]:
    
20.             mnum = 31
    
21.         elif m == 2:
    
22.             if (a+1979) % 4 == 0 and (a+1979) % 100 != 0:
    
23.                 mnum = 29
    
24.             else:
    
25.                 mnum = 28
    
26.         else:
    
27.             mnum = 30
    
28.         for i in range(94):
    
29.             for j in range(192):
    
30.                 data[a,m,i,j] = sum(x[num:num+mnum,i,j]) / mnum
    
31.         num = num + mnum
    
32.     a = a + 1
    
33. 
    
34. #将月平均值按年份、月份、纬度、经度的顺序存放至dat文件中
    
35. with open(r'ulwrf.dat','w') as fp:
    
36.     for n in range(39):
    
37.         for m in range(12):
    
38.             for i in range(94):
    
39.                 for j in range(192):
    
40.                     fp.write(str(data[n,m,i,j])+'  ')

复制代码

详情也可查看https://blog.csdn.net/qq_47672397/article/details/124534856?spm=1001.2014.3001.5501