<div dir="ltr">Dear NCL,<div><br></div><div>I want to plot snow cover extent after multiplying it with snow cover area. Here is the code for it. I am also attaching the figure with this email.</div><div><br></div><div><div>;*******************************************</div><div>;Reading in the weekly snow cover extent </div><div>;*******************************************</div><div><br></div><div>a=addfile(&quot;~/Documents/NCL_files/Sea_ice/<a href="http://nhsce_v01r01_19661004_20140602.nc">nhsce_v01r01_19661004_20140602.nc</a>&quot;,&quot;r&quot;)</div><div>sce=a-&gt;snow_cover_extent ; The values are either 0 or 1</div><div><br></div><div>lat=a-&gt;latitude</div><div>lon=a-&gt;longitude</div><div><br></div><div>sce!1=&quot;lat&quot;</div><div>sce!2=&quot;lon&quot;</div><div><br></div><div>sce@lat2d=lat</div><div>sce@lon2d=lon</div><div><br></div><div>;********************************************</div><div>:Reading in the snow cover are</div><div>;********************************************</div><div><br></div><div>sce_area=a-&gt;area</div><div><br></div><div>sce_area!0=&quot;lat&quot;</div><div>sce_area!1=&quot;lon&quot;</div><div><br></div><div>copy_VarAtts(sce,sce_area)</div><div>copy_VarCoords_1(sce,sce_area)</div><div><br></div><div>;*******************************************</div><div>;Muliplying Snow cover extent and Area</div><div>;*******************************************</div><div><br></div><div>B = new((/2487,88,88/),float)</div><div><br></div><div>     do i = 0,2486</div><div>      B(i,:,:)=sce(i,:,:)#sce_area(:,:)</div><div>      end do</div><div><br></div><div>B!0=&quot;time&quot;      </div><div>B!1=&quot;lat&quot;</div><div>B!2=&quot;lon&quot;      </div><div>   </div><div>B@lat2d=lat</div><div>B@lon2d=lon</div><div><br></div><div>copy_VarAtts(sce,B)</div><div>copy_VarCoords_1(sce,B)</div><div><br></div><div>printVarSummary(B)</div><div><br></div><div>B@_FillValue = -9</div><div><br></div><div> </div><div>;************************************************</div><div>; create plot</div><div>;************************************************</div><div>  wks = gsn_open_wks(&quot;pdf&quot; ,&quot;Snow_cover_area&quot;)             </div><div>  res            = True                          ; plot mods desired</div><div>  res@gsnPolar   = &quot;NH&quot;                          ; specify the hemisphere</div><div>  plot = gsn_csm_contour_map_polar(wks,B(0,:,:),res)    ; create the plot</div><div><br></div>
</div></div>