- 积分
- 754
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2013-1-13
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
本帖最后由 男紫汗 于 2013-11-11 11:05 编辑
箭头长度好长。。。。求救各位大神
gs:
'reinit'
'set mapdraw off'
'open u3.ctl'
'open v3.ctl'
'set dfile 2'
'set lon 0 357.5'
'set lat -90 90'
'd u.1;v.2'
;
u3.ctl:
DSET ^u3.dat
OPTIONS sequential
UNDEF -2.56E33
TITLE 5 Days of Sample Model Output
XDEF 144 LINEAR 0.0 2.5
YDEF 73 LINEAR -90.0 2.5
ZDEF 1 LEVELS 1000
TDEF 42 LINEAR 02JAN1960 1yr
VARS 1
u 1 99 Temperature
ENDVARS
Fortran程序:program ssimplicit none
integer i,j
real(8)::u1(73,144),u2(73,144)=0,v1(73,144),v2(73,144)=0,u3(73,144)=0,v3(73,144)=0
OPEN(11,FILE='vwind.dat')
do i=1,73
do j=1,144
read(11,"(E21.14)") v1(i,j)
enddo
enddo
close(11)
open(20,file='uwind.dat')
do i=1,73
do j=1,144
read(20,"(E21.14)") u1(i,j)
enddo
enddo
close(20)
print*,u1(1,1)
print*,u1(73,144)
do j=1,144
do i=1,73
if(2.le.j.and.j.le.143.and.2.le.i.and.i.le.72) then
u2(i,j)=u1(i,j)/2.+(u1(i,j-1)+u1(i,j+1)+u1(i-1,j)+u1(i+1,j))/8.
v2(i,j)=v1(i,j)/2.+(v1(i,j-1)+v1(i,j+1)+v1(i-1,j)+v1(i+1,j))/8.
else
u2(i,j)=u1(i,j)
v2(i,j)=v1(i,j)
endif
enddo
enddo
u3=u1-u2
v3=v1-v2
open(12,file='u3.dat',form='unformatted',access='SEQUENTIAL')
open(13,file='v3.dat',form='unformatted',access='SEQUENTIAL')
open(14,file='u2.dat',form='unformatted',access='SEQUENTIAL')
open(15,file='v2.dat',form='unformatted',access='SEQUENTIAL')
open(16,file='u1.dat',form='unformatted',access='SEQUENTIAL')
open(17,file='v1.dat',form='unformatted',access='SEQUENTIAL')
do i=1,144
do j=1,73
write(12) u3(j,i)
write(13) v3(j,i)
write(14) u2(j,i)
write(15) v2(j,i)
write(16) u1(j,i)
write(17) v1(j,i)
enddo
enddo
end
这个有问题吗?
|
-
|