MVLib - Short Writeup - Version 4.0
Example Number 08 - Retrive and plot wordwide Era data precipitation
This example shows how to retrive and plot wordwide ERA project precipitaion data. The routine eraheader read the header of the file, filling the arrays for latitudes and longitudes (lat and lonin the fortran code). The routine eradata read the precipitation field for the selected hour. The routine worldplot will plot the rain array of data.
The source code to produce these plots is available here and the complete script running the code is available here.
Previous Example List of Examples Next Example
MVLib Home Page
	integer ora,giorno,mese,anno
	character data*60,file*80
	parameter (nlon=480,nlat=241)
	real rain(nlon,nlat),lat(nlon,nlat),lon(nlon,nlat)
	
	ora=10 ; giorno=20 ; mese=6 ; anno=2000
	call eraheader(10,ora,giorno,mese,anno,lat,lon)
	call eradata  (10,ora,giorno,mese,anno,rain)
	call datafromhour (ora,giorno,mese,anno,data)
	call worldplot(rain,nlon,nlat,' ',' ','lm0.1c6p4M10')
	call mvsetflags('Tipo di Carattere',-13.0)
	call displayexample('example08','ERA-Interim project - Precipitation '//
	2	'ReAnalysis','precipitation (mm) '//data)
	end