<html>

<head>

<meta http-equiv="Content-Type" content="text/html; charset=us-ascii">

</head>

<body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; color: rgb(0, 0, 0); font-size: 14px; font-family: Calibri, sans-serif; ">

<div>Check that the dimensions of 't' and the dimensions of 'msk' are the same</div>

<div><br>

</div>

<div>Will</div>

<div><br>

</div>

<span id="OLK_SRC_BODY_SECTION">

<div style="font-family:Calibri; font-size:11pt; text-align:left; color:black; BORDER-BOTTOM: medium none; BORDER-LEFT: medium none; PADDING-BOTTOM: 0in; PADDING-LEFT: 0in; PADDING-RIGHT: 0in; BORDER-TOP: #b5c4df 1pt solid; BORDER-RIGHT: medium none; PADDING-TOP: 3pt">

<span style="font-weight:bold">From: </span>Adv &lt;<a href="mailto:advita6@gmail.com">advita6@gmail.com</a>&gt;<br>

<span style="font-weight:bold">Date: </span>Wednesday, 18 February 2015 4:27 PM<br>

<span style="font-weight:bold">To: </span>&quot;<a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a>&quot; &lt;<a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a>&gt;<br>

<span style="font-weight:bold">Subject: </span>[ncl-talk] Sub: How to mask ocean part in my plot. ..<br>

</div>

<div><br>

</div>

<div>

<div>

<div dir="ltr">

<div>

<div>

<div>Hi,<br>

</div>

I would like to plot correlation only over Northern US. Could someone help me.? My script plots over ocean part too. So i wan't to mask ocean. I tried to use mask file with mask function.

<br>

<br>

</div>

Error looks like<br>

fatal:mask: dimension sizes&nbsp; of parameter 0 and parameter 1 do not match<br>

fatal:Execute: Error occurred at or near line 95 in file testcorr.ncl<br>

&nbsp; <br>

</div>

<div>Help and suggestions are appreciated.<br>

<br>

</div>

<div>Thanks<br>

</div>

<div><br>

<br>

t=t(lat|:,lon|:)<br>

printVarSummary(t)<br>

copy_VarMeta(a,t)<br>

<br>

&nbsp;in = addfile(&quot;<a href="http://lsmask.nc">lsmask.nc</a>&quot;,&quot;r&quot;)<br>

&nbsp;msk=in-&gt;mask<br>

printVarSummary(in)<br>

<br>

; use mask function to mask out land and then ocean data<br>

; ocean=0,land=1,sea_ice=2<br>

;************************************************<br>

&nbsp; land_only&nbsp; = t&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; trick to keep cv's and atts<br>

&nbsp; ocean_only = t<br>

<br>

&nbsp; land_only&nbsp; = mask(t,msk(0,:,:),1)&nbsp;&nbsp; ; Return the values of ts where oro=1<br>

;&nbsp; ocean_only = mask(ts,oro,0)&nbsp;&nbsp; ; Return the values of ts where oro=0<br>

<br>

wks = gsn_open_wks(&quot;ps&quot;,&quot;proj&quot;)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; open a ncgm file<br>

&nbsp; gsn_define_colormap(wks,&quot;gui_default&quot;)<br>

<br>

;**************************************************<br>

&nbsp; res&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = True<br>

&nbsp; res@gsnFrame&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = False&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; don't advance frame yet<br>

&nbsp; res@mpDataBaseVersion&nbsp;&nbsp;&nbsp;&nbsp; = &quot;Ncarg4_1&quot;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; higher res data base<br>

&nbsp; res@vpWidthF&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0.90&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; change aspect ratio of plot<br>

&nbsp; res@vpHeightF&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0.90<br>

&nbsp; res@vpXF&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0.05<br>

&nbsp; res@vpYF&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0.95<br>

&nbsp; res@mpProjection&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = &quot;LambertConformal&quot;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; choose projection<br>

&nbsp; res@mpLambertParallel1F&nbsp;&nbsp;&nbsp; = 33.0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; two parallels<br>

&nbsp; res@mpLambertParallel2F&nbsp;&nbsp;&nbsp; = 45.0<br>

&nbsp; res@mpLambertMeridianF&nbsp;&nbsp;&nbsp;&nbsp; = -95.0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; central meridian<br>

&nbsp; res@mpLimitMode&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = &quot;LatLon&quot;<br>

&nbsp; res@mpMinLatF&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 24.0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; map area<br>

&nbsp; res@mpMaxLatF&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 50.0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; latitudes<br>

&nbsp; res@mpMinLonF&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = -125.0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; and<br>

&nbsp; res@mpMaxLonF&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = -65.0<br>

&nbsp; res@mpFillOn&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = True&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; turn on map fill<br>

<br>

&nbsp;res@mpLandFillColor&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = &quot;Background&quot;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; set to white<br>

&nbsp; res@mpPerimOn&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = False&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; don't draw box around map<br>

&nbsp; res@mpGridAndLimbOn&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = False<br>

&nbsp; res@mpPerimOn&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = False&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; turn off box around plot<br>

<br>

&nbsp; res@cnLinesOn&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = False&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; turn off contour lines<br>

<br>

&nbsp; res@gsnSpreadColors&nbsp;&nbsp;&nbsp;&nbsp; = True&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; use full range of color map<br>

;&nbsp; res@gsnSpreadColorStart = 2<br>

;&nbsp; res@gsnSpreadColorEnd&nbsp;&nbsp; = -3<br>

;&nbsp; res@lbLabelFontHeightF&nbsp; = 0.015&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; label bar font height<br>

<br>

&nbsp; res@tiMainString&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = &quot;Example of Orthographic Projection&quot;&nbsp;&nbsp; ; add a title<br>

;&nbsp; res@tiMainFontHeightF&nbsp; = .018&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; font height<br>

<br>

&nbsp; plot = gsn_csm_contour_map(wks,t,res)&nbsp; ; create the plot<br>

</div>

end<br>

<div><br>

<br>

</div>

</div>

</div>

</div>

</span>

<p style="font-size:10pt; line-height:10pt; calibri"><br>

<br>

University of Tasmania Electronic Communications Policy (December, 2014). <br>

This email is confidential, and is for the intended recipient only. Access, disclosure, copying, distribution, or reliance on any of it by anyone outside the intended recipient organisation is prohibited and may be a criminal offence. Please delete if obtained

 in error and email confirmation to the sender. The views expressed in this email are not necessarily the views of the University of Tasmania, unless clearly intended otherwise.

</p>

</body>

</html>