<div dir="ltr">Hello all,<div><br></div><div>I tried to plot a simple map but I got warnings. Here is my code and warnings:</div><div><br></div><div><div><br></div><div>;*********** From Observation **********</div><div><br></div><div>ana=addfile("NOaerosol_allsky_3DQ.nc","r")</div><div>lat2=ana->lat</div><div>lon2=ana->lon</div><div>lev2=ana->lev/10</div><div>Absp3D=ana->Absp3D ; Atmospheric heating (K/sec)</div><div>Absp=ana->Absp ; Solar radiation absorption in the atmosphere (W/m^2)</div><div>Sfc=ana->Sfc ; Surface net flux (W/m^2)</div><div><br></div><div>delete(ana)</div><div><br></div><div>;********* T's data *************</div><div>;z1 = asciiread("Absp.dat",(/12,128,64,26/),"float")</div><div>z1 = asciiread("Absp.dat",(/128,64,26,12/),"float")</div><div><br></div><div>z1!1="lat"</div><div>z1!0="lon"</div><div>;Coordinate</div><div>z1&lat=lat2</div><div>z1&lon=lon2</div><div><br></div><div><br></div><div>wks = gsn_open_wks("eps","T_Atmospheric_heating_26")</div><div><br></div><div>;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;</div><div> res = True ; plot mods desired</div><div> res@cnFillOn = True ; turn on color </div><div> res@gsnSpreadColors = True ; use full color map</div><div> res@cnLinesOn = False ; no contour lines</div><div> res@cnLineLabelsOn = False ; no line labels</div><div><br></div><div> res@mpMaxLatF = 87.8638 ; specify the plot domain</div><div> res@mpMinLatF = -87.8638 ; </div><div> res@mpMinLonF = 0 ;</div><div> res@mpMaxLonF = 357.1875 </div><div> </div><div> res@gsnAddCyclic = False </div><div><br></div><div> res@cnLevelSelectionMode = "AutomaticLevels" ; AutomaticLevels</div><div><br></div><div><br></div><div> ;res@cnLevelSelectionMode = "ManualLevels" ; AutomaticLevels</div><div> ;res@cnMinLevelValF = -15. ; min level</div><div> ;res@cnMaxLevelValF = 15. ; max level</div><div> ;res@cnLevelSpacingF = 1. ; interval</div><div><br></div><div>;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;</div><div>res@tiMainString = "Atmospheric Heating in Jan at 26th level"</div><div><br></div><div>plot= gsn_csm_contour_map(wks,z1(:,:,25,0),res)</div><div>delete(wks)</div><div><br></div></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div>
<p class="gmail-p1"><span class="gmail-s1">0) check_for_y_lat_coord: Warning: Data either does not contain a valid latitude coordinate array or doesn't contain one at all.</span></p>
<p class="gmail-p1"><span class="gmail-s1">(0) A valid latitude coordinate array should have a 'units' attribute equal to one of the following values: </span></p>
<p class="gmail-p1"><span class="gmail-s1">(0) 'degrees_north' 'degrees-north' 'degree_north' 'degrees north' 'degrees_N' 'Degrees_north' 'degree_N' 'degreeN' 'degreesN' 'deg north'</span></p>
<p class="gmail-p1"><span class="gmail-s1">(0) check_for_lon_coord: Warning: Data either does not contain a valid longitude coordinate array or doesn't contain one at all.</span></p>
<p class="gmail-p1"><span class="gmail-s1">(0) A valid longitude coordinate array should have a 'units' attribute equal to one of the following values: </span></p>
<p class="gmail-p1"><span class="gmail-s1">(0) 'degrees_east' 'degrees-east' 'degree_east' 'degrees east' 'degrees_E' 'Degrees_east' 'degree_E' 'degreeE' 'degreesE' 'deg east'</span></p>
<p class="gmail-p1"><span class="gmail-s1">(0) gsn_csm_map_ce: Warning: you set mpMaxLonF to a value > 180, but</span></p>
<p class="gmail-p1"><span class="gmail-s1">(0) didn't set mpCenterLonF. Setting mpCenterLonF to 178.594</span></p><p class="gmail-p1"><span class="gmail-s1"><br></span></p><p class="gmail-p1"><span class="gmail-s1"><br></span></p><p class="gmail-p1"><span class="gmail-s1"><br></span></p><p class="gmail-p1">Thanks,</p><p class="gmail-p1">Ana</p></div></div>