MVLib - Short Writeup - Version 4.0 Example Number 01 - How to read and produce a simple 2-d plot with 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 temperature,
the call to mm5colormap
allow to produce a color map.
The source code to produce these plots is available here, the complete script running the code is available here. |
|||||
| ||||||
parameter (ixm=100,jxm=100) real tgr(ixm,jxm) ; character data*40 n=0 open (40,file='tmp/MMOUTP_DOMAIN3',form='unformatted',status='old') call mm5readfld (40,'ground t',n,tgr ,ixm,jxm,1,ix,jx,kx) call datafrommm5index (mm5mif(1,mm5mif(1,1)),data) call kelvin2centi(tgr,ixm,jxm,1) ! conv. da kelvin a centigradi call mvsetflags('Plot Confini',3.0) call mvsetflags('Colore Titoli',20.0) call mvsetflags('Colore Confini',20.0) call mvsetflags('Primo colore',2.0) call mvsetflags('Ogni quanti colori',2.0) call mm5colormap(tgr,ixm,jxm,ix-1,jx-1,'temperatura') call scrivisulplot(data,0.5,0.92) call displayexample('example01','Example 1 - Plot of MM5 Model 2D fields',' ') end |