提取nc資料的逐月資料 grads

2021-07-06 10:35:29 字數 1723 閱讀 5404

**原帖位址

**我還在嘗試搞清楚出來的grd在matlab裡怎麼排布。。。想不明白0 0!

'reinit'

*先設定各月的名稱

mon.1=jan;mon.2=feb;mon.3=mar;mon.4=apr;mon.5=may;mon.6=jun;

mon.7=jul;mon.8=aug;mon.9=sep;mon.10=oct;mon.11=nov;mon.12=dec

*檔案讀入,準備輸出

'sdfopen e:\f\essay\data\data\ncep\pres.mon.mean.nc'

'set gxout fwrite'

imon=1

while(imon<=12)

'set fwrite e:\f\essay\data\data\ncep\div\pres\'mon.imon'_p.grd'

y=1948*資料起始年

y1=1960

y2=2013

while(y1<=y2)

ta=(y1-y)*12+imon

'set t 'ta''

nz=1

*地面氣壓僅為一層資料,故不迴圈,而shum,u,vwnd都是8層

*while(nz<=8)

'set z 'nz''

'set x 29 57'

'set y 43 60'

*這裡注意一下,當x,y輸出的時候設定了其實,ctl檔案描述相應的要進行修改

*(xdef 28 liner 70 2.5

ydef 17 liner 15 2.5)

'd pres'

*nz=nz+1

*endwhile

y1=y1+1

endwhile

'disable fwrite'

imon=imon+1

endwhile ;

********************fortran處理成txt***************=

program main

implicit none

integer it,iz,n,m,i,j

real::shum(29,17,8,792)

n=29

m=17

open (10,file='e:/llx/shum.grd',form='binary')

do it=1,792

do iz=1,8

read(10 )((shum(i,j,iz,it),i=1,n),j=1,m)

end do

end do

close(10)

open(12,file='e:/llx/shuiqi/shum1.txt',status='new')

do it=139,792 !我所需要擷取的時間

do iz=1,8

write(12,*) ((shum(i,j,iz,it),i=1,n),j=1,m)

end do

end do

close(12)

end其中①it,iz是grd檔案中時間和高度的格點數;

②n,m分別是x方向和y方向上格點數,為「大數-小數+1」

i.e.:'set x 29 57'則i=57-29+1

③write()函式在寫入grd時寫「write(12)」、寫入txt時為「write(12,*)」

④寫成grd的時候必須時間,高度,緯度,經度的順序,怎麼寫進去怎麼讀出來

資料提取 常用的手機資料提取方法

閒暇時間,隨手記錄,願與諸位朋友分享。學識有限,不當之處,懇請各位大神不吝賜教,也是對我自己的學習提高過程!言歸正文 手機取證,從字面理解,可以分為取和證兩個過程。取,把資料原原本本的從手機中提取出來。證,通過資料檢索 挖掘 分析,尋找與案件有關的線索和證據。對於我們來說,二者缺一不可。但萬事開頭難...

資料提取 資料分析 第七步提取相應的資料指標

sr research中國代表處官網 通過前面的5或6步操作後,基本操作的預處理過程已經完成,對於某些領域的實驗,例如閱讀中,還會存在步驟為注視過濾等操作剔除小於或大於一定範圍值的資料等的操作 在後面的文章中會依次介紹到 完成上述操作後,就可以進行眼動資料指標的提取了,對於ia興趣區指標的提取可參考...

NC14661 簡單的資料結構

因為要從隊頭隊尾插入刪除元素,可以用stl中的deque 雙端佇列 便於操作.這些是deque的一些基本操作 1.把x壓入後 前端 push back x push front x 2.訪問 不刪除 後 前端元素 back front 3.刪除後 前端元素 pop back pop front 4....