Fortran批量为nc写ctl

黎大页子

有的NetCDF文件不符合COARES约定，不能用sdfopen直接打开，需要自己配写ctl文件。详见兰溪的帖子【免费派送】“SDF file has no discernable X coordinate” 问题详解http://bbs.06climate.com/forum.php?mod=viewthread&tid=6008&fromuid=1368

但是网上下载数据的多半是一年一个nc文件，如果是日资料，需要判断闰年，这种情况我不会用一个ctl描述全部文件。一个一个写又太麻烦。所以可以用Fortran批量生产。
1.先在excel里面写好要生产的ctl路径+文件名，斜杠一定要用“\”，用"/"Fortran读的时候会出错，然后Ctrl+下拉.....
2.复制放进namelist.txt，然后用fortran读namelist，读了路径+文件名放进数组name()里面。
3.加入闰年判断子程序
4.依次把ctl文件内容write进去，运行就可以看见可爱的ctl们一个一个自己蹦出来啦，真是神清气爽。

程序很简单，唯一让我百度了一下的地方是不换行的问题，有一个现成的advance=“no”，很好用。advance是write（）里面的控制语句，默认为yes,也就是换行，改成no可以接着往下写，但是一定要指明写的内容的格式。

这个问题不难，经验总结比较弱，自己mark下来，也给同为小白的同志们共勉。

PS：后来发现，其实GrADS考虑了闰年的情况，所以无论是日资料，还是月资料，都可以一个ctl搞定，详见
http://bbs.06climate.com/forum.php?mod=viewthread&tid=6663
但是因为我的数据处理是每年分开做的，所以ctl分开比较方便，不然每次还要掐着指头算第X年在几千天当中究竟是哪一段.....大家根据自己数据处理的需要选择合适的方法吧。

1. !This program is to write ctl for monthly mean documents:''sh/lh_oaflux_xxxx.nc''
   
2. Program writectl
   
3. implicit none
   
4. integer::i,j,y,var,days
   
5. character(len=100)::name(54)
   
6. 
   
7. open(10,file='e:/thesis/lh/ctl_name.txt')!读取批处理文件名
   
8. rewind(10)
   
9.   do i=1,26,1
   
10.   y=1985+i-1
    
11.   call year(y,days)
    
12. 
    
13.   read(10,*)name(i) !把文件名放入数组
    
14. !----写文件--------------------------------------------
    
15.   write(*,*)name(i)
    
16.   open(37,file=name(i),form='formatted',iostat=var)
    
17.   write(37,'(a)',advance="no")'dset e:/thesis/lh/lh_oaflux_'
    
18.   write(37,'(i4)',advance="no")y
    
19.   write(37,'(a)',advance="no")'.nc'
    
20.   write(37,*)'title monthly mean surface latent heat flux'
    
21.   write(37,*)'undef 32766'
    
22.   write(37,*)'xdef lon 360 linear 0.5 1.0'
    
23.   write(37,*)'ydef lat 180 linear -89.5 1.0'
    
24.   write(37,'(a)',advance="no")'tdef time '
    
25.   write(37,'(i3)',advance="no")days
    
26.   write(37,'(a)',advance="no")' linear 00Z01JAN'
    
27.   write(37,'(i4)',advance="no")y
    
28.   write(37,*)'1dy'
    
29.   write(37,*)'vars 2'
    
30.   write(37,*)'lhtfl=>lhtfl 0 99 daily mean surface latent heat'
    
31.   write(37,*)'err=>err 0 99 daily mean estimated error of latent heat'
    
32.   write(37,*)'endvars'
    
33. !---------------------------------------------------------
    
34.   close(37)
    
35.   end do
    
36. close(10)
    
37. 
    
38. stop
    
39. end
    
40. 
    
41. !闰年判断-------------------------------
    
42. subroutine year(m,dys)
    
43. implicit none
    
44. integer::m,dys
    
45. dys=365
    
46. if(mod(m,100)==0)then
    
47.    if(mod(m,400)==0)then
    
48.    dys=dys+1
    
49.    end if
    
50. else
    
51.    if(mod(m,4)==0)then
    
52.    dys=dys+1
    
53.    end if
    
54. end if
    
55. 
    
56. return
    
57. stop
    
58. end subroutine
    
59. end

!This program is to write ctl for monthly mean documents:''sh/lh_oaflux_xxxx.nc'' <br /> Program writectl<br /> implicit none<br /> integer::i,j,y,var,days<br /> character(len=100)::name(54)<br /> <br /> open(10,file='e:/thesis/lh/ctl_name.txt')!读取批处理文件名<br /> rewind(10)<br />   do i=1,26,1<br />   y=1985+i-1<br />   call year(y,days)<br /> <br />   read(10,*)name(i) !把文件名放入数组<br /> !----写文件--------------------------------------------<br />   write(*,*)name(i)<br />   open(37,file=name(i),form='formatted',iostat=var)<br />   write(37,'(a)',advance="no")'dset e:/thesis/lh/lh_oaflux_'<br />   write(37,'(i4)',advance="no")y<br />   write(37,'(a)',advance="no")'.nc'<br />   write(37,*)'title monthly mean surface latent heat flux'<br />   write(37,*)'undef 32766'<br />   write(37,*)'xdef lon 360 linear 0.5 1.0'<br />   write(37,*)'ydef lat 180 linear -89.5 1.0'<br />   write(37,'(a)',advance="no")'tdef time '<br />   write(37,'(i3)',advance="no")days<br />   write(37,'(a)',advance="no")' linear 00Z01JAN'<br />   write(37,'(i4)',advance="no")y<br />   write(37,*)'1dy'<br />   write(37,*)'vars 2'<br />   write(37,*)'lhtfl=>lhtfl 0 99 daily mean surface latent heat'<br />   write(37,*)'err=>err 0 99 daily mean estimated error of latent heat'<br />   write(37,*)'endvars'<br /> !---------------------------------------------------------<br />   close(37)<br />   end do<br /> close(10)<br /> <br /> stop<br /> end<br /> <br /> !闰年判断-------------------------------<br /> subroutine year(m,dys)<br /> implicit none<br /> integer::m,dys<br /> dys=365<br /> if(mod(m,100)==0)then<br />    if(mod(m,400)==0)then<br />    dys=dys+1<br />    end if<br /> else<br />    if(mod(m,4)==0)then<br />    dys=dys+1<br />    end if<br /> end if<br /> <br /> return<br /> stop<br /> end subroutine<br /> end

topmad

artmunich 发表于 2012-5-4 16:30

                               
登录/注册后可看大图

深深你还我沙发。。。。。

你们这....人民内部矛盾....结果被升级到阶级矛盾了

言深深

哎呀，沙发啊

黎大页子

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

                               
登录/注册后可看大图

哎呀，沙发啊

你真的好快！，一看就是天天泡在网上的技术宅男呀，嘿嘿。

言深深

黎大页子 发表于 2012-5-4 16:21

                               
登录/注册后可看大图

你真的好快！，一看就是天天泡在网上的技术宅男呀，嘿嘿。

吓···哪有，碰巧一下子登录进来就看到了呗

artmunich

深深你还我沙发。。。。。

mofangbao

真好，很省力啊

言深深

artmunich 发表于 2012-5-4 16:30

                               
登录/注册后可看大图

深深你还我沙发。。。。。

哈····

传说中的谁

围观一下，好久木有原创帖子了

黎大页子

传说中的谁 发表于 2012-5-4 16:50

                               
登录/注册后可看大图

围观一下，好久木有原创帖子了

什么叫做大战A？