计算涡度、涡度平流、温度平流、散度、水汽通量、水汽通量散度、比湿等的FORTRAN程序

river

      今天看到有人发帖问有没有计算涡度、温度平流的程序，搜索了一下论坛，grads因为有自带函数，所以类似的gs很多啊，但是fortran的好像几乎没有。      突然想到大学实习的时候有过一个程序可以一次计算出很多物理量，包括几个平流的计算。只是那会儿没好好学，不知道当时用的什么资料。看了一下，是格点的，可能是micaps的某一类资料吧。需要的人自己研究吧。贴其中的一段出来吧

1. !---------------------------- vorticity of observed wind-------------------------------
   
2. 
   
3. voro=1.e+36
   
4. 
   
5. DO it=1,NT
   
6. DO iz=1,NZ   
   
7.      do i=1+1,NX-1   
   
8.      do j=1+1,NY-1   
   
9. voro(i,j,iz,it)=1./(R*2)*( (v(i+1,j,iz,it)-v(i-1,j,iz,it))/(Delta*pi/180*cos(lat(j)))-(u(i,j+1,iz,it)-u(i,j-1,iz,it))/(Delta*pi/180)+2*u(i,j,iz,it)*tan(lat(j)) )   
   
10.      enddo   
    
11.      enddo
    
12. ENDDO
    
13. ENDDO
    
14. !-------------------------------divergence--------------------------------------
    
15. div=1.e+36

!---------------------------- vorticity of observed wind-------------------------------<br /> <br /> voro=1.e+36<br /> <br /> DO it=1,NT<br /> DO iz=1,NZ   <br />      do i=1+1,NX-1   <br />      do j=1+1,NY-1   <br /> voro(i,j,iz,it)=1./(R*2)*( (v(i+1,j,iz,it)-v(i-1,j,iz,it))/(Delta*pi/180*cos(lat(j)))-(u(i,j+1,iz,it)-u(i,j-1,iz,it))/(Delta*pi/180)+2*u(i,j,iz,it)*tan(lat(j)) )    <br />      enddo   <br />      enddo<br /> ENDDO<br /> ENDDO<br /> !-------------------------------divergence--------------------------------------<br /> div=1.e+36

     废话不多说，上附件。 all.f90 (6.35 KB, 下载次数: 2415)

2012-10-25 16:41 上传

点击文件名下载附件
阅读权限: 20




来自群组: 龙王山皇家气象学院

Aires

river 发表于 2012-10-25 17:15

                               
登录/注册后可看大图

我发了两个贴了，每次都是你第一个来哦。好感动哦

哈哈，因为我经常会来论坛，希望以后多多分享好东西哈

Aires

谢谢river的分享~

mlittlej

谢谢谢谢谢谢啦

cold_wq

谢谢分享你的程序，程序收到

善人/jw

{:eb502:}{:eb502:}{:eb502:}
顶起顶起

奕紘

感謝分享~

freekiller

顶lz               

Ivy

谢谢！你是我们的榜样

river

Aires 发表于 2012-10-25 17:22

                               
登录/注册后可看大图

哈哈，因为我经常会来论坛，希望以后多多分享好东西哈

嗯，有好东西一定会分享的

river

Aires 发表于 2012-10-25 16:52

                               
登录/注册后可看大图

谢谢river的分享~

我发了两个贴了，每次都是你第一个来哦。好感动哦