- 积分
- 3601
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2012-9-8
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
本帖最后由 hillside 于 2013-6-23 13:08 编辑
现上传下载于“南京大气资料服务中心”官网(“ 国家自然科学基金委员会地球科学部南京信息工程大气资料服务中心 ” ,简称 “ 南京大气资料服务中心 ”)的《南亚高压环流指数》计算程序。简单检索之下,未发现标题博文,如发现此程序已有上传,版主或可建议版主删除此帖。
程序内容摘录:
!用(王盘兴老师定义的)闭合气压系统的方法计算1948-2007年6-8月的南亚高压的面积、强度和中心位置
program main
parameter(ny=60,nm=12,m=143,n=72) !!!!! ny: 年份,nm: 月份,
dimension ys(0:m,0:n,nm,ny),sj(nm,ny),pj(nm,ny),staj(nm,ny),fi(nm,ny)
real lamj(nm,ny)
integer f0(nm)
data f0/1652,1654,1654,1656,1662,1672,1676,1676,1669,1660,1654,1650/ !1~12月特征等值线
!********读取原始资料(100hPa位势高度场)************
ys=0.0 !!!!! 先把所有格点的高度场赋值为零
open(1,file='e:\ng\hgt100.grd',form='binary')
read(1)((((ys(i,j,im,iy),i=0,143),j=0,72),im=1,nm),iy=1,ny-1)
!!!!! 2008年9月在资料中心拷的hgt.mean.mon.nc
!!!!! hgt.mean.mon.nc只到2007年8月,所以必须分开赋值,2007年9-12月赋值为零
read(1)(((ys(i,j,im,ny),i=0,143),j=0,72),im=1,8)
close(1)
!!!!资料存放方式:(纬向:纬度从南极-90到赤道0再到北极90,格点是0——36——72)
注:该官网似乎仅此一个程序可供下载,镇网之宝?
|
-
-
ng.f90
10.96 KB, 下载次数: 259, 下载积分: 金钱 -5
评分
-
查看全部评分
|