Re: FW: problem with plotting ECMWF data

Hi Basit,
I replied to this yesterday; I am guessing you did not see the reply. I
copied it below.. Adam

On 03/11/2013 10:28 AM, Adam Phillips wrote:
> Hi Basit,
> If the latitudes and longitudes of your ECMWF data are 2-dimensional,
> then you should attach them as lat2d/lon2d attributes to your msl
> array. If g0_lat_1/g0_lon_2 are 2-dimensional, change this:
> msl&g0_lat_1 = g0_lat_1
> msl&g0_lon_2 = g0_lon_2
>
> to this:
> msl@lat2d = g0_lat_1
> msl@lon2d = g0_lon_2
>
> If the latitude and longitude arrays are not 2-dimensional: Are you
> getting any warning messages? Is the data on a special map projection?
>
> If the above does not help let ncl-talk know.. Adam
>
> On 03/11/2013 06:57 AM, BasitAli Khan wrote:
>> 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 attached) 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 below:
>>
>>
>> Thanks and best regards,
>> Basit
>>
>> ;*************************************************
>> 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
>>
>>
>> ----
>>
>> Basit A. Khan, Ph.D.
>>
>> Postdoctoral Research Fellow
>>
>> Division of Physical Sciences & Engineering
>>
>> Office# 3204, Level 3, Building 1,
>>
>> King Abdullah University of Science & Technology
>>
>> 4700 King Abdullah Blvd, Box 2753, Thuwal 23955 --6900,
>>
>> Kingdom of Saudi Arabia.
>>
>> Office: +966(0)2 808 0276, Mobile: +966(0)5 0860 3617
>>
>> E-mail: basitali.khan@kaust.edu.sa <mailto:basitali.khan@kaust.edu.sa>
>>
>> Skype name: basit.a.khan
>>
>> ------------------------------------------------------------------------
>>
>> 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
>
> --
>

-- 
______________________________________________________________
Adam Phillips                                asphilli@ucar.edu
NCAR/Climate and Global Dynamics Division       (303) 497-1726
P.O. Box 3000				
Boulder, CO 80307-3000    http://www.cgd.ucar.edu/cas/asphilli

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Tue Mar 12 07:26:54 2013