fortran版diamond4转grads的源程序

明-春

谢谢分享……

浅水

感谢楼主 分享  学习了哈   

慕容玄悟

试试看，学习中

云淡ぐ风轻

学习学习！

zgzggzgz

这么牛逼？

红河岸

谢谢楼主，学习下

cold_wq

我自己也写过，我读出来后，直接按照**.000文件里面描述的起始经纬度和经纬向间隔来做ctl文件，结果画出来完全和直接用Meteoinfo画的不同啊，不知道是不是自己这样直接做不对...
以下是自己的code，求解：

1. !Micaps文件头：
   
2. !diamond 4 12年09月15日08时850百帕温度平流场
   
3. !12 09 15 08   0  850
   
4. !5.313 -4.000    .0 170.0  80.0  12.0
   
5. !33  18     2.0  -200.0   200.0  1  0
   
6. 
   
7. !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
   
8. integer ri,IREC,cc
   
9. character*2 yea
   
10. character*2 mon,hou
    
11. character*1 ri0
    
12. character*2 ri1
    
13. character*80 a0,b0,c0,d0
    
14. real::tad(18,33),qad(18,33)
    
15. character cc0(11)*4
    
16. data cc0/'1000','925','850','700','500','400','300','250','200','150','100'/
    
17. yea="12"
    
18. mon="11"
    
19. ri1="14"
    
20. hou="08"
    
21. IREC=1
    
22. ! ****高空:温度平流和湿度平流读写************************
    
23.   OPEN (11,FILE="D:\1113_14_fog\20121114_08_phy.dat",form="UNFORMATTED",ACCESS='DIRECT',RECL=33*18)
    
24. 
    
25. !     温度平流tad
    
26. cc=3
    
27. OPEN (8,FILE="C:\data\MICAPS\physic\tadv"//trim(cc0(cc))//""//yea//mon//ri1//hou//".000",status="old")
    
28. read(8,'(a80)') a0
    
29. read(8,'(a80)') b0
    
30. read(8,'(a80)') c0
    
31. read(8,'(a80)') d0
    
32. do lat=1,18
    
33.    read(8,'(10f8.1)') (tad(lat,lon),lon=1,10)
    
34.    read(8,'(10f8.1)') (tad(lat,lon),lon=11,20)
    
35.    read(8,'(10f8.1)') (tad(lat,lon),lon=21,30)
    
36.    read(8,'(3f8.1)') (tad(lat,lon),lon=31,33)
    
37. enddo
    
38. close(8)
    
39. 
    
40. !     水汽平流qad
    
41. OPEN (8,FILE="C:\data\MICAPS\physic\qadv"//trim(cc0(cc))//""//yea//mon//ri1//hou//".000",status="old")
    
42. read(8,'(a80)') a0
    
43. read(8,'(a80)') b0
    
44. read(8,'(a80)') c0
    
45. read(8,'(a80)') d0
    
46. do lat=1,18
    
47.    read(8,'(10f8.1)') (qad(lat,lon),lon=1,10)
    
48.    read(8,'(10f8.1)') (qad(lat,lon),lon=11,20)
    
49.    read(8,'(10f8.1)') (qad(lat,lon),lon=21,30)
    
50.    read(8,'(3f8.1)') (qad(lat,lon),lon=31,33)
    
51. enddo
    
52. close(8)
    
53. 
    
54. 
    
55. write(11,REC=IREC) ((tad(LAT,LON),LON=1,33),LAT=18,1,-1)
    
56. close(9)
    
57. end

复制代码

还有对应的ctl：
dset D:\20120915_08_phy.dat
undef 9999
title physic
options yrev
xdef 33 linear 0 5.313
ydef 18 linear 12 4.000
zdef 1 levels 850
tdef 1 linear 08z15Sep2012 1dy
vars 1
tadv 1 99 temperature advection
endvars

kongfeng0824

请问清风，fortran中处理完数据后用什么画图会更好一些呢？surfer、ncl还是grads？

笑哈乐呵

看哈宝贝。

怿燃

先保存，在好好学习