MVLib - Short Writeup - Version 4.0
Example Number 02 - How to read and produce a simple 2-d plot of circulation from 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, the call to mm5freccette allow to produce a map of ground circlation. 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 u(ixm,jxm,kxm),v(ixm,jxm,kxm),lat(ixm,jxm),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)
	close(40)
        call mvsetflags('Colore Freccette',11.0)
	call mvsetflags('Plot confini',3.0)
        call mvsetflags('Colore Confini',20.0)
        call mvsetflags('Colore Titoli',20.0)
        call mvsetflags('Step Freccette',2.0)
        call mm5freccette(u,v,lat,lon,ixm,jxm,kxm,ix,jx)
        call datafrommm5index (mm5mif(1,mm5mif(1,1)),data)
	call displayexample('example02','Example 02 - How to use mm5freccette'//
	2	' subroutine',data)
        end