PIOMAS Arctic sea ice thickness daily数据下载和IDL处理

忘是亡心i 2022-11-10 10:52 67阅读 0赞

网上好像没有PIOMAS数据处理的文章，最近项目需要进行了处理，因此给大家提供以下参考，不过自己摸索一下应该很快就会了

PIOMAS数据下载

[http://psc.apl.uw.edu/research/projects/arctic-sea-ice-volume-anomaly/data/model\_grid][http_psc.apl.uw.edu_research_projects_arctic-sea-ice-volume-anomaly_data_model_grid]

在表格中的Sea ice thickness (Volume per unit Area) , daily mean即为我下载的海冰厚度数据

点击右侧的[  hiday.H<yyyy>][hiday.H_yyyy]

跳转到[https://pscfiles.apl.washington.edu/zhang/PIOMAS/data/v2.1/hiday/][hiday.H_yyyy]下载数据

我这里是下载了[hiday.H2021][]

这个数据明显比其他数据小，因为其实只有1和2月的数据

PIOMAS数据处理

其实官网提供了读取数据的IDL程序，就是叫做read\_draw\_hi.pro这个程序，但是这个程序包含了制图，因此我觉得对我来说稍微有点复杂了。

我按照它程序提供的信息和思路，重新用IDL写了代码，因此以下代码就是能够将我上述下载的数据读取出来，存成每天一个波段的envi格式的数据（不带地理信息），我也将mask之外的数据存储成nan值了

除了上述的数据外，还需要下载[grid.dat数据和][grid.dat]io.dat\_360\_120.output数据，[http://psc.apl.uw.edu/research/projects/arctic-sea-ice-volume-anomaly/data/model\_grid][http_psc.apl.uw.edu_research_projects_arctic-sea-ice-volume-anomaly_data_model_grid][ 这个网址都有。][http_psc.apl.uw.edu_research_projects_arctic-sea-ice-volume-anomaly_data_model_grid]

以下程序应该还算是比较简明，值得注意的是一些数据的地址要修改，另外就是如果读取的数据改变导致波段改变，59要改成别的值。因为我处理的是北极的信息，因此nx和ny是这个值，但是如果不是的话，可以从read\_draw\_hi.pro这个程序再找出进行修改

之后我使用envi的Georeference from IGM建立地理信息。这也是以下代码需要保存lon和lat的原因。

之后可能会更新批量处理的代码，代码也会加上建立地理信息的步骤。

pro read_piomas_2021
      nx=360 & ny=120;120 ;Arctic Parallel Ocean and sea Ice Model (POIM) grid
      max_levels=30
      desired_level =1
      mask_value=99999.0
      barwidth=200.0
      lon=fltarr(nx,ny) & lat=fltarr(nx,ny)
      data4=fltarr(nx,ny,59, /NOZERO);注意这里的59其实是指数据为1和2月数据一起59天，最后生成59各波段的数据
      kmt=intarr(nx,ny)
      ;读取经纬度
      openr, lunin, 'H:\mission\PIOMAS\grid.dat',/Get_lun  ;打开数据文件
      readf, lunin, lon, format='(10f8.2)'
      readf, lunin, lat, format='(10f8.2)'
      close,lunin
      free_lun,lunin
      ;read grid mask (ocean levels)
      openr,lunin,'H:\mission\PIOMAS\io.dat_360_120.output',/Get_lun  ;打开数据文件
      readf, lunin, kmt,format='(360i2)'
      close,lunin
      free_lun,lunin
      
      
    ;读取数据
     openr,lunin, 'H:\mission\PIOMAS\hiday.H2021', /swap_if_big_endian,/Get_lun  ;打开数据文件'H:\mission\LNG\PIOMAS\monthly\heff.H2021'
     readu,lunin, data4
     close,lunin
     free_lun,lunin
     location=where(kmt lt 1,count)
     if count gt 0 then begin
       numFiles=size(data4)
       numFiles=numFiles[3]
       FOR fidx=0, numFiles-1 DO BEGIN
        temp=data4[*,*,fidx]
        temp[location]=0.0/0
        data4[*,*,fidx]=temp[*,*]
       ENDFOR
     endif
    
     print,'count',count
     
     ;储存数据
     openw,lunou,'H:\mission\PIOMAS\lon',/get_lun
     writeu,lunou,lon;
     close,lunou
     free_lun,lunou
     
     openw,lunou,'H:\mission\PIOMAS\lat',/get_lun
     writeu,lunou,lat;
     close,lunou
     free_lun,lunou
     
     openw,lunou,'H:\mission\PIOMAS\piomas2021',/get_lun
     writeu,lunou,data4;
     close,lunou
     free_lun,lunou
    print,'finished'
    
    
    
    end

PIOMAS的一些详细信息，包括格网可以参考：[http://psc.apl.uw.edu/research/projects/projections-of-an-ice-diminished-arctic-ocean/model/][http_psc.apl.uw.edu_research_projects_projections-of-an-ice-diminished-arctic-ocean_model]

[http_psc.apl.uw.edu_research_projects_arctic-sea-ice-volume-anomaly_data_model_grid]: http://psc.apl.uw.edu/research/projects/arctic-sea-ice-volume-anomaly/data/model_grid
[hiday.H_yyyy]: https://pscfiles.apl.washington.edu/zhang/PIOMAS/data/v2.1/hiday/
[hiday.H2021]: https://pscfiles.apl.washington.edu/zhang/PIOMAS/data/v2.1/hiday/hiday.H2021
[grid.dat]: https://pscfiles.apl.washington.edu/zhang/PIOMAS/utilities/grid.dat
[http_psc.apl.uw.edu_research_projects_projections-of-an-ice-diminished-arctic-ocean_model]: http://psc.apl.uw.edu/research/projects/projections-of-an-ice-diminished-arctic-ocean/model/