<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("~/Documents/NCL_files/Sea_ice/<a href="http://nhsce_v01r01_19661004_20140602.nc">nhsce_v01r01_19661004_20140602.nc</a>","r")</div><div>sce=a->snow_cover_extent ; The values are either 0 or 1</div><div><br></div><div>lat=a->latitude</div><div>lon=a->longitude</div><div><br></div><div>sce!1="lat"</div><div>sce!2="lon"</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->area</div><div><br></div><div>sce_area!0="lat"</div><div>sce_area!1="lon"</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="time" </div><div>B!1="lat"</div><div>B!2="lon" </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("pdf" ,"Snow_cover_area") </div><div> res = True ; plot mods desired</div><div> res@gsnPolar = "NH" ; specify the hemisphere</div><div> plot = gsn_csm_contour_map_polar(wks,B(0,:,:),res) ; create the plot</div><div><br></div>
</div></div>