日志
- 分享
Python处理era中的总云量，出现问题
-
2024-3-7 11:13
-
一开始用的是import xarray as xr 这个读取nc数据，别的数据都能正常读取，但是在读取总云量时，出现了问题，输出的数据都是inf，无穷大，用ncl重新输出了这个变量的详细信息，发现存在scale_factor和add_offset 也就是缩放和偏移，所以 from netCDF4 import Dataset需要导入这个函数，能够将正确的数据读取出来
-
11 次阅读|0 个评论
- 分享
出中山大学气象学874初试资料
-
2023-10-6 09:46
-
出中山大学气象学874初试资料，复试资料也出，可以加qq：1186999490，备注好信息
-
150 次阅读|0 个评论
- 分享
NCL插值，格点变多，精细化
-
2023-5-23 23:54
-
begin ;这边定义的glon和glat是你想要插值的区域的经纬度以及格点数 glon=fspan(73,106,111);定义一个浮点型数据，从90-112，分为221个格点 glat=fspan(24,43,61) glon!0="lon" glon@long_name="lon" glon@units="degrees_east" glonlon=glon glat!0="lat" glat@lon ...
-
138 次阅读|0 个评论
- 分享
提取变量 编写成一个nc文档
-
2023-5-23 23:42
-
begin ; f1=addfile("evspsbl_Amon_E3SM-1-1-ECA_historical_r1i1p1f1_gr_197001-197912.nc","r") ; f2=addfile("evspsbl_Amon_E3SM-1-1-ECA_historical_r1i1p1f1_gr_198001-198912.nc","r") ; f3=addfile("evspsbl_Amon_E3SM-1-1-ECA_historical_r1i1p1f1_gr_199001-199912.nc","r") ; f4=addfile("evspsbl_Amon_E3SM ...
-
87 次阅读|0 个评论
- 分享
将多个逐年数据编写成一个nc数据 ncl
-
2023-5-23 23:35
-
begin Files=systemfunc("ls evspsbl_Amon_EC-Earth3-Veg-LR_historical_r1i1p1f1_gr_*01-*12.nc") f=addfiles(Files,"r") ;printVarSummary(f) evp0=f -evspsbl(:,{24:43},{73:107}) evp1=f -evspsbl(:,{24:43},{73:107}) evp2=f -evspsbl(:,{24:43},{73:107}) evp3=f -evspsbl(:,{24:43},{73:107}) ...
-
72 次阅读|0 个评论