数据处理，HadISST1_SST_1870-1900.txt转nc

lt5201314

数据处理，HadISST1_SST_1870-1900.txt 转 nc

这个问题我做了很久没有找到合适的方法，用ncl自带的asciiread不能直接读，请问论坛各位前辈怎么处理呢？
ncl能处理吗？不能的话用什么办法或者软件可以呢？

井中月

hadley sst有封装好的nc文件啊，为啥不直接下nc文件？

lt5201314

数据下载https://www.metoffice.gov.uk/hadobs/hadisst/data/download.html

edwardli

Python语言，大概不超过10行就能完成数据准备

lt5201314

edwardli 发表于 2021-1-12 23:01
Python语言，大概不超过10行就能完成数据准备

怎么实现呢，我写的代码读不出来

sc0rp1o1

1. loadscript("$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl")
   
2. loadscript("$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl")
   
3. loadscript("$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl")
   
4. begin  
   
5. 
   
6.         lat=fspan(89.5,-89.5,180)
   
7.         lon=fspan(-179.5,179.5,360)
   
8.        
   
9.         lat!0="lat"
   
10.         lat@long_name="latitude"
    
11.         lat@units="degrees_north"
    
12.         lon!0="lon"
    
13.         lon@long_name="longitude"       
    
14.         lon@units="degrees_east"
    
15. 
    
16.         f0=asciiread("/mnt/e/NCL/SST/HadISST1_SST_2004.txt",(/2172,360/),"integer")
    
17.         sst_jan=f0(1:180,:)
    
18.         sst_jan@_FillValue=-32768       
    
19.         sst_jan2float=tofloat(sst_jan)/100.
    
20. 
    
21.         sst_jan2float!0="latitude"
    
22.         sst_jan2float!1="longitude"
    
23.         sst_jan2float&latitude=lat
    
24.         sst_jan2float&longitude=lon
    
25.         sst_jan2float@units="C"
    
26.         sst_jan2float@long_name="sea surface temperature in january 2004"
    
27.         printVarSummary(sst_jan2float)
    
28. 
    
29.         sst_feb=f0(182:361,:)
    
30.         sst_feb@_FillValue=-32768
    
31.         sst_feb2float=tofloat(sst_feb)/100.       
    
32. 
    
33.         sst_feb2float!0="latitude"
    
34.         sst_feb2float!1="longitude"
    
35.         sst_feb2float&latitude=lat
    
36.         sst_feb2float&longitude=lon
    
37.         sst_feb2float@units="C"
    
38.         sst_feb2float@long_name="sea surface temperature in february 2004"
    
39.         printVarSummary(sst_feb2float)
    
40. 
    
41.         system("rm -f /mnt/e/NCL/SST/HadISST1_SST_2004.nc")
    
42.         fo=addfile("/mnt/e/NCL/SST/HadISST1_SST_2004.nc","c")
    
43.         fo->sst_jan=sst_jan2float
    
44.         fo->sst_feb=sst_feb2float
    
45. 
    
46. end

复制代码

sc0rp1o1

我下载了HadISST1_SST_2004.txt.gz这个数据，测试了一下，没什么问题，你可以参考一下 。这个数据格式稍微麻烦了一点就是每隔一个月就有一行说明，所以读的时候需要跳过那一行说明。剩下的月份你可以计算一下行数，读就行了。

Ifare

matlab应该也可以做

lt5201314

sc0rp1o1 发表于 2021-1-13 10:12
我下载了HadISST1_SST_2004.txt.gz这个数据，测试了一下，没什么问题，你可以参考一下 。这个数据格式稍微 ...

这个一共若是年份太多了挨着算有点多诶，可以循环实现吗

sc0rp1o1

lt5201314 发表于 2021-1-14 08:37
这个一共若是年份太多了挨着算有点多诶，可以循环实现吗

循环应该是可以的，你得想一下 sst = f0(:,:) 这一行应该怎么写