160站数据转化问题

☆涅垚YTZ☆

最近很多同学都用到这个，本人不才，也来献献丑。

单时次转化：

1. program main         
   
2. real vec(160)
   
3. real lat(160),lon(160)
   
4. character*8 stid(160)
   
5. 
   
6. open(1,file='c:\data\r16006.txt',status='old')
   
7. read(1,*) (vec(i),i=1,160)
   
8. close(1)
   
9. 
   
10. open(2,file='c:\data\china.txt')
    
11. do k=1,160
    
12. read(2,*) lat(k),lon(k)
    
13. enddo
    
14. close(2)
    
15. 
    
16. do 2 i=1,160
    
17. 2 stid(i)=char(i)
    
18. 
    
19. open (3,file='c:\data\r16006.grd',form='binary')
    
20. tim=0.0
    
21. nlev=1
    
22. nflag=1
    
23. do i=1,160
    
24. write(3) stid(i),lat(i),lon(i),tim,nlev,nflag,vec(i)
    
25. enddo
    
26. nlev=0
    
27. write(3) stid(i-1),lat(i-1),lon(i-1),tim,nlev,nflag
    
28. close(3)
    
29. end   
    

program main         <br /> real vec(160)<br /> real lat(160),lon(160)<br /> character*8 stid(160)<br /> <br /> open(1,file='c:\data\r16006.txt',status='old')<br /> read(1,*) (vec(i),i=1,160)<br /> close(1)<br /> <br /> open(2,file='c:\data\china.txt')<br /> do k=1,160<br /> read(2,*) lat(k),lon(k)<br /> enddo<br /> close(2)<br /> <br /> do 2 i=1,160<br /> 2 stid(i)=char(i)<br /> <br /> open (3,file='c:\data\r16006.grd',form='binary')<br /> tim=0.0<br /> nlev=1<br /> nflag=1<br /> do i=1,160<br /> write(3) stid(i),lat(i),lon(i),tim,nlev,nflag,vec(i)<br /> enddo<br /> nlev=0<br /> write(3) stid(i-1),lat(i-1),lon(i-1),tim,nlev,nflag<br /> close(3)<br /> end   <br />

多时次转化：

1. program main         
   
2. real vec(160,50)
   
3. real lat(160),lon(160)
   
4. character*8 stid(160)
   
5. 
   
6. open(1,file='r16006.txt',status='old')
   
7. read(1,*) ((vec(i,j),i=1,160),j=1,50)
   
8. close(1)
   
9. 
   
10. open(2,file='china.txt')
    
11. do k=1,160
    
12. read(2,*) lat(k),lon(k)
    
13. enddo
    
14. close(2)
    
15. 
    
16. do 2 i=1,160
    
17. 2 stid(i)=char(i)
    
18. 
    
19. open (3,file=r16006.grd',form='binary')
    
20. do j=1,50
    
21. tim=0.0
    
22. nlev=1
    
23. nflag=1
    
24. do i=1,160
    
25. write(3) stid(i),lat(i),lon(i),tim,nlev,nflag,vec(i,j)
    
26. enddo
    
27. nlev=0
    
28. write(3) stid(i-1),lat(i-1),lon(i-1),tim,nlev,nflag
    
29. enddo
    
30. close(3)
    
31. end   
    

program main         <br /> real vec(160,50)<br /> real lat(160),lon(160)<br /> character*8 stid(160)<br /> <br /> open(1,file='r16006.txt',status='old')<br /> read(1,*) ((vec(i,j),i=1,160),j=1,50)<br /> close(1)<br /> <br /> open(2,file='china.txt')<br /> do k=1,160<br /> read(2,*) lat(k),lon(k)<br /> enddo<br /> close(2)<br /> <br /> do 2 i=1,160<br /> 2 stid(i)=char(i)<br /> <br /> open (3,file=r16006.grd',form='binary')<br /> do j=1,50<br /> tim=0.0<br /> nlev=1<br /> nflag=1<br /> do i=1,160<br /> write(3) stid(i),lat(i),lon(i),tim,nlev,nflag,vec(i,j)<br /> enddo<br /> nlev=0<br /> write(3) stid(i-1),lat(i-1),lon(i-1),tim,nlev,nflag<br /> enddo<br /> close(3)<br /> end   <br />

大家注意改路径

mofangbao

支持分享啊！

topmad

这样的分享对于提高水平很有帮助的

随缘

抓个原始数据图就好了

言深深

还可以，写成grd主要是和grads接口的吗？
另外，比如do 2这样的结构以后尽量用do enddo结构来代替，前者属于77老版本的，可以与你自己的书写习惯也统一起来（你前文有写do enddo结构）

善人/jw

来顶一个

☆涅垚YTZ☆

言深深 发表于 2012-5-19 16:07

                               
登录/注册后可看大图

还可以，写成grd主要是和grads接口的吗？
另外，比如do 2这样的结构以后尽量用do enddo结构来代替，前者属 ...

嗯，是啊~主要是根据书上的写的，自己的改的时候就会写do enddo 了，谢谢给予指导~

☆涅垚YTZ☆

随缘 发表于 2012-5-19 14:31

                               
登录/注册后可看大图

抓个原始数据图就好了

截图原始数据的存放格式是吗？

随缘

☆涅垚YTZ☆ 发表于 2012-5-19 20:42

                               
登录/注册后可看大图

截图原始数据的存放格式是吗？

嗯，是啊是是

yyb_nuist

小曹，支持一个！