MVLib - Short Writeup - Version 4.0
Example Number 03 - How to read and produce a simple 2-d plot of fields of MM5 model output
This example shows the simplest way to produce a map with MM5 model output. The call to mm5readfields routine allows to retrieve data of of different fields, a combination of call to mm5freccette, mm5freccetelle and mm5colormap allow to produce a map of ground circlation and temperature. The source code to produce these plots is available here, the complete script running the code is available here.
Previous Example MVLib Home Page Next Example
	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
        open (40,file='tmp/MMOUTP_DOMAIN3',form='unformatted',status='old')
        n=0
        call mm5readfld (40,'latitdot',n,lat,ixm,jxm,1,ix,jx,kx)
        call mm5readfld (40,'longidot',n,lon,ixm,jxm,1,ix,jx,kx)
        call mm5readfld (40,'u' ,n,u,ixm,jxm,kxm,ix,jx,kx)
        call mm5readfld (40,'v' ,n,v,ixm,jxm,kxm,ix,jx,kx)
        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 mvsetflags('Plot confini',3.0)
	call mvsetflags('Colore confini',7.0)
	call mvsetflags('Label Bar',1.0)
	call mvsetflags('Label Direction',2.0)
        call mvsetflags('Colore Titoli',20.0)
	call mm5freccette(u,v,lat,lon,ixm,jxm,kxm,ix-1,jx-1)
        call mm5colormap(tgr,ixm,jxm,ix-1,jx-1,'temperatura')
	call displayexample('example03','How to use mm5colormap and mm5freccette',
	2	data)

	call mvsetflags('Colore freccetelle',3.0)
	call mvsetflags('step freccette',2.0)
        call mm5colormap(tgr,ixm,jxm,ix-1,jx-1,'temperatura')
	call mm5freccetelle(u,v,ixm,jxm,kxm,ix-1,jx-1)
	call displayexample('example03','How to use mm5colormap and mm5freccetelle',
	2	data)
        end