[ncl-talk] Plot a contour

[saeed bayat]

Dear NCL Users
I am pretty new in NCL. I am using output of RegCM4 to plot topo in NCL. I
want to add mask to the map and country borders and then change the numbers
of X and Y axis to the real lat and lon.
 I attached the plot I got using one template codes in NCL website (and the
script). Also attached the header information of file I am using.
Please help me to solve the problem.

-- 
Abdolhalim Bayat
Ph.D Candidate of Synoptic Climatology
The Tehran University of Iran
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20150611/b96bbf41/attachment-0001.html 
-------------- next part --------------
netcdf lorestan-3_DOMAIN000 {
dimensions:
	iy = 60 ;
	jx = 120 ;
	kz = 19 ;
variables:
	float sigma(kz) ;
		sigma:standard_name = "atmosphere_sigma_coordinate" ;
		sigma:long_name = "Sigma at model layer midpoints" ;
		sigma:units = "1" ;
		sigma:axis = "Z" ;
		sigma:positive = "down" ;
		sigma:formula_terms = "sigma: sigma ps: ps ptop: ptop" ;
	float ptop ;
		ptop:standard_name = "air_pressure" ;
		ptop:long_name = "Pressure at model top" ;
		ptop:units = "hPa" ;
	float iy(iy) ;
		iy:standard_name = "projection_y_coordinate" ;
		iy:long_name = "y-coordinate in Cartesian system" ;
		iy:units = "km" ;
	float jx(jx) ;
		jx:standard_name = "projection_x_coordinate" ;
		jx:long_name = "x-coordinate in Cartesian system" ;
		jx:units = "km" ;
	float xlat(iy, jx) ;
		xlat:standard_name = "latitude" ;
		xlat:long_name = "Latitude at cross points" ;
		xlat:units = "degrees_north" ;
	float xlon(iy, jx) ;
		xlon:standard_name = "longitude" ;
		xlon:long_name = "Longitude at cross points" ;
		xlon:units = "degrees_east" ;
	float dlat(iy, jx) ;
		dlat:standard_name = "latitude" ;
		dlat:long_name = "Latitude at dot points" ;
		dlat:units = "degrees_north" ;
	float dlon(iy, jx) ;
		dlon:standard_name = "longitude" ;
		dlon:long_name = "Longitude at dot points" ;
		dlon:units = "degrees_east" ;
	float topo(iy, jx) ;
		topo:standard_name = "surface_altitude" ;
		topo:long_name = "Domain surface elevation" ;
		topo:units = "m" ;
		topo:coordinates = "xlon xlat" ;
	float landuse(iy, jx) ;
		landuse:legend = "1  => Crop/mixed farming\n",
			"2  => Short grass\n",
			"3  => Evergreen needleleaf tree\n",
			"4  => Deciduous needleleaf tree\n",
			"5  => Deciduous broadleaf tree\n",
			"6  => Evergreen broadleaf tree\n",
			"7  => Tall grass\n",
			"8  => Desert\n",
			"9  => Tundra\n",
			"10 => Irrigated Crop\n",
			"11 => Semi-desert\n",
			"12 => Ice cap/glacier\n",
			"13 => Bog or marsh\n",
			"14 => Inland water\n",
			"15 => Ocean\n",
			"16 => Evergreen shrub\n",
			"17 => Deciduous shrub\n",
			"18 => Mixed Woodland\n",
			"19 => Forest/Field mosaic\n",
			"20 => Water and Land mixture\n",
			"21 => Urban\n",
			"22 => Sub-Urban" ;
		landuse:standard_name = "land_type" ;
		landuse:long_name = "Landuse category as defined in BATS1E" ;
		landuse:units = "1" ;
		landuse:coordinates = "xlon xlat" ;
	float xmap(iy, jx) ;
		xmap:standard_name = "map_factor" ;
		xmap:long_name = "Map factor in domain cross points" ;
		xmap:units = "1" ;
		xmap:coordinates = "xlon xlat" ;
	float dmap(iy, jx) ;
		dmap:standard_name = "map_factor" ;
		dmap:long_name = "Map factor in domain dot points" ;
		dmap:units = "1" ;
		dmap:coordinates = "dlon dlat" ;
	float coriol(iy, jx) ;
		coriol:standard_name = "coriolis_parameter" ;
		coriol:long_name = "Coriolis force parameter" ;
		coriol:units = "s-1" ;
		coriol:coordinates = "xlon xlat" ;
	float snowam(iy, jx) ;
		snowam:standard_name = "snowfall_amount" ;
		snowam:long_name = "Snow initial amount in mm" ;
		snowam:units = "kg m-2" ;
		snowam:coordinates = "xlon xlat" ;
	float mask(iy, jx) ;
		mask:standard_name = "land_binary_mask" ;
		mask:long_name = "Land Sea mask" ;
		mask:units = "1" ;
		mask:coordinates = "xlon xlat" ;

// global attributes:
		:title = "ICTP Regional Climatic model V4 domain" ;
		:institution = "ICTP" ;
		:source = "RegCM Model simulation Terrain output" ;
		:Conventions = "CF-1.4" ;
		:history = "2015-03-21 09:16:29 : Created by RegCM terrain" ;
		:references = "http://eforge.escience-lab.org/gf/project/regcm" ;
		:experiment = "lorestan-3" ;
		:projection = "LAMCON" ;
		:grid_size_in_meters = 60000. ;
		:latitude_of_projection_origin = 36. ;
		:longitude_of_projection_origin = 31. ;
		:standard_parallel = 25.f, 46.f ;
		:boundary_smoothing = "No" ;
		:minimum_h2o_pct_for_water = 75. ;
		:input_dataset_resolution_in_minutes = 2 ;
		:landuse_fudging = "No" ;
		:grid_factor = 0.584011211247294 ;
}
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Screenshot from 2015-06-11 13:26:53.png
Type: image/png
Size: 195131 bytes
Desc: not available
Url : http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20150611/b96bbf41/attachment-0001.png 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: topo.ncl
Type: application/octet-stream
Size: 2428 bytes
Desc: not available
Url : http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20150611/b96bbf41/attachment-0001.obj