<div dir="ltr">Dear NCL,<div><br></div><div>I am trying to plot SWE data and the lat lon dimensions don&#39;t show up.</div><div>I would be grateful if you could guide me on this</div><div>*********************************************************</div><div><div>;********************************************</div><div>;This is reading in Snow water equivalent which goes from 0 to 250 mm</div><div>; from 1982 to 2012</div><div>;********************************************</div><div><br></div><div>dira = &quot;/Users/ipshita/Desktop/GlobsnowV2.0/version02/&quot;</div><div>dirb = dira + &quot;2011/&quot;</div><div>;dira = &quot;./&quot;</div><div>;dirb = &quot;./&quot;</div><div>a = addfile(dira+&quot;<a href="http://monthly_82_13.nc">monthly_82_13.nc</a>&quot;,&quot;r&quot;)</div><div>b = addfile(dirb+&quot;<a href="http://GlobSnow_SWE_L3A_20110101_v2.0.nc">GlobSnow_SWE_L3A_20110101_v2.0.nc</a>&quot;,&quot;r&quot;)</div><div><br></div><div>swe  = a-&gt;swe</div><div>time = a-&gt;time</div><div>lat  = b-&gt;lat</div><div>lon  = b-&gt;lon</div><div>swe@lat2d = lat</div><div>swe@lon2d = lon</div><div><br></div><div>print(dimsizes(lat))</div><div>print(dimsizes(lon))</div><div><br></div><div>swe_FillValue=-9.99e+08</div><div><br></div><div>utc_date = cd_calendar(time, 0)</div><div>year=utc_date(:,0)</div><div>month=utc_date(:,1)</div><div><br></div><div>;********************************************</div><div>;This is to extract 1982-2012 of swe data</div><div>;********************************************</div><div>swe_82_12=swe(0:371,:,:)</div><div><br></div><div>;********************************************</div><div><br></div><div>;********************************************</div><div>;This is to read in All India Rainfall</div><div>;For March , March , June</div><div>;********************************************</div><div><br></div><div>MJ_1982_2014=asciiread(&quot;~/Documents/NCL_files/SST/Monthly/MJ.txt&quot;,(/34,1/),&quot;float&quot;)</div><div>MJ_1982_2012=MJ_1982_2014(0:30,0)</div><div><br></div><div>print(dimsizes(MJ_1982_2012))</div><div><br></div><div>;*******************************************</div><div>;Extracting monthly data for sce</div><div>;*******************************************</div><div><br></div><div>March_swe = new((/31,721,721/),float)</div><div><br></div><div>March_swe@_FillValue = -9.99e+08</div><div><br></div><div>do nyr=0,371,12</div><div><br></div><div>March_swe(nyr/12,:,:) =swe_82_12(nyr+2,:,:)</div><div><br></div><div>end do</div><div>;=================================================</div><div>March_swe@_FillValue = -9.99e+08</div><div>;**********************************************</div><div><br></div><div>March_dt = dtrend_msg(March_swe&amp;time,March_swe(lat|:,lon|:,time|:),True,False)</div><div><br></div><div>printMinMax(March_dt,False)</div><div><br></div><div>;****************************************</div><div>;Calculating Correlation</div><div>;***************************************</div><div>x1d=ndtooned(MJ_1982_2012)</div><div><br></div><div>corr_March=escorc(March_dt,x1d)</div><div>corr_March!0=&quot;lat&quot; </div><div>corr_March!1=&quot;lon&quot; </div><div>corr_March&amp;lat=lat</div><div>corr_March&amp;lon=lon</div><div><br></div><div>copy_VarCoords(swe,corr_March)<br></div><div><br></div><div><br></div><div>;**********************************************</div><div><br></div><div>;************************************************</div><div>; create plot</div><div>;************************************************</div><div>  wks = gsn_open_wks(&quot;x11&quot;,&quot;SWE_March_MJ_correlation_significance&quot;)               ; </div><div>  </div><div>    res                     = True              ; Plot modes desired.</div><div>    res@gsnMaximize         = True              ; Maximize plot</div><div>    res@cnFillOn            = True              ; color plot desired</div><div>    res@cnFillPalette       = &quot;BrownBlue12&quot;</div><div>    res@cnLinesOn           = False             ; turn off contour lines</div><div>    res@pmLabelBarWidthF    = 0.9               ; make wider</div><div>    res@pmLabelBarHeightF   = 0.1               ; default is taller</div><div>    res@lbLabelFontHeightF  = .018              ; default is HUGE</div><div>    </div><div> </div><div>;*******************************************</div><div>; georeferencing: plot on polar projection</div><div>;*******************************************</div><div>    res@trGridType   = &quot;TriangularMesh&quot;        ; allow missing coordinates</div><div>    res@gsnAddCyclic = False</div><div>    res@cnFillMode   = &quot;RasterFill&quot;</div><div><br></div><div>    res@gsnPolar   = &quot;NH&quot;                          ; specify the hemisphere</div><div>    res@mpMinLatF  = 35</div><div>    </div><div>    res@tiMainString = &quot;March SWE and MJ correlation&quot;</div><div><br></div><div>  plot = gsn_csm_contour(wks,corr_March,res)</div><div>  </div><div>;*******************************************************</div><div><br></div><div><div>warning:gsnPolar is not a valid resource in SWE_March_MJ_correlation_significance_contour at this time</div><div>warning:mpMinLatF is not a valid resource in SWE_March_MJ_correlation_significance_contour at this time</div><div><br></div></div><div><br></div><div class="gmail_signature"><div dir="ltr"><br></div></div>
</div></div>