parameter (ixm=100,jxm=100,kxm=23)
real tgr(ixm,jxm),u(ixm,jxm,kxm),v(ixm,jxm,kxm),lat(ixm,jxm),
2 lon(ixm,jxm)
character data*40,file*30
open (40,file='tmp/MMOUTP_DOMAIN3',form='unformatted',status='old')
call mvsetflags('Plot confini',3.0)
call mvsetflags('Colore confini',20.0)
call mvsetflags('Colore Titoli',20.0)
n=0
do nn=1,10
call mm5readfld (40,'ground t',n,tgr ,ixm,jxm,1,ix,jx,kx)
call kelvin2centi(tgr,ixm,jxm,1)
call datafrommm5index (mm5mif(1,mm5mif(1,1)),data)
call mm5colormap(tgr,ixm,jxm,ix-1,jx-1,'temperatura')
call mvlibsignature
call makefilename('stemp',nn,'.gif',file)
call mvsetflags('X size',200.0)
call scrivisulplot(data,0.5,0.92)
call displayplot(file,'Example 04 - MM5 ground temperature',' ')
call makefilename('btemp',nn,'.gif',file)
call mvsetflags('X size',600.0)
call displayplot(file,' ',' ')
enddo
call viacolvento('stemp*.gif','example041s.gif')
call viacolvento('btemp*.gif','example041.gif')
end
|