use chymdata , only : chymreadrec parameter (nlon=330,nlat=280) parameter (slon=6.50,elon=9.46,slat=44.10,elat=46.61) character vsource*50 real dem(nlon,nlat) open(10,file='tmp/Po.chym',status='old',form='unformatted') call chymheader(10) call chymreadrec(10,'dem',vsource,dem,idate) call mvsetflags('Palette di Colori',5.0) call mvsetflags('Scrivi Max-Min',0.0) call settaiconfini(slon,elon,slat,elat) call mvsetflags('Curva di Livello', 200.0) call mvsetflags('Curva di Livello',1000.0) call mvsetflags('Curva di Livello', 200.0) call mvsetflags('colore isolinee',4.0) call contorni(dem,nlon,nlat,nlon,nlat) call mvsetflags('Size delle linee',1.5) call mvsetflags('Colore confini',1.0) call boundaries(slon,elon,slat,elat,1) call drawlatlon(slon,elon,slat,elat) call displayexample('example29','How to use the subroutine contorni',' ') end