Hi guys,
I am trying to plot a simple 2D variable from ECMWF data on selected area. NCL plots the SLP on the map but the background map is larger than the plotted SLP (png image copied below) although lat/lon coordinates for map and slp are same. I guess it is because of some projection issue, i tried but could not fix it. I will appreciate if anyone could help me fixing this. The NCL script is copied after the image:
Thanks and best regards,
Basit
[cid:E76F4C2D-C33F-4A58-8416-88538D4B176A]
;*************************************************
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
;************************************************
begin
;************************************************
; create pointer to file and read in data
;************************************************
in = addfile("/project/k14/basit/wrf-sims/wrf_166/ei_oper_an_sfc_075x075_20100526_20100719.nc","r")
minLat = 19.0
maxLat = 25.0
minLon = 36.5
maxLon = 44.3
g0_lat_1 = in->g0_lat_1 ; for _SFC
g0_lon_2 = in->g0_lon_2 ; for _SFC
initial_time0_hours = in->initial_time0_hours ;time in julian date format
yymmddhh = cd_calendar(initial_time0_hours,3)
msl = in->MSL_GDS0_SFC
dims_msl = dimsizes(msl)
msl&g0_lat_1 = g0_lat_1
msl&g0_lon_2 = g0_lon_2
do it=136, 136, 1
;************************************************
; create default plot
;************************************************
wks = gsn_open_wks("png","ecmwf_plot") ; open a ps file
gsn_define_colormap(wks,"BlAqGrYeOrRe") ; choose colormap
res = True ; plot mods desired
res@cnFillOn = True ; turn on color fill
res@cnLinesOn = False ; turn off contour lines
res@pmTickMarkDisplayMode = "Always"; use NCL default lat/lon labels
res@gsnAddCyclic = False ; data already has cyclic point
res@mpGeophysicalLineColor = "Black"
res@mpLimbLineColor = "Blue"
res@mpGeophysicalLineThicknessF = 1.5
res@mpDataBaseVersion = "HighRes" ; Use the high-res database
res@mpDataResolution = "FinestResolution"
res@pmTickMarkDisplayMode = "Always" ; Turn on map tickmarks.
res@mpOutlineOn = True
res@mpFillOn = False
res@mpGridLineThicknessF = 0.5
res@mpGridSpacingF = 3
res@mpLimitMode = "LatLon"
res@mpMinLatF = 19. ; range to zoom in on
res@mpMaxLatF = 25.
res@mpMinLonF = 36.5
res@mpMaxLonF = 44.3
res@gsnSpreadColors = True
res@gsnSpreadColorEnd = -3
res@gsnSpreadColors = True ; use full range of color map
plot = gsn_csm_contour_map(wks,msl(it,{19.:25.},{36.5:44.3}),res)
frame(wks)
end do
end
________________________________
This message and its contents including attachments are intended solely for the original recipient. If you are not the intended recipient or have received this message in error, please notify me immediately and delete this message from your computer system. Any unauthorized use or distribution is prohibited. Please consider the environment before printing this email.
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
This archive was generated by hypermail 2.1.8 : Wed Mar 13 2013 - 14:19:38 MDT