<div dir="ltr"><div class="gmail_default" style="font-size:small">Dear Sara,</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">The problem is that you are trying to use wrf_xxxx functions on non-WRF data. The wrf_xxxx functions are tailored to work specifically with WRF-ARW output data, which is usually named something like "wrfout_d01_xxxxxxx".</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">I don't know what your data looks like, so I can't tell you how to correctly plot it. However, if you visit this page:</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default"><a href="http://www.ncl.ucar.edu/Applications/plot_data_on_map.shtml">http://www.ncl.ucar.edu/Applications/plot_data_on_map.shtml</a><br></div><div class="gmail_default"><br></div><div class="gmail_default">then my guess is that you'll be able to use dataonmap_1.ncl if your data is rectilinear, or dataonmap_3.ncl if your data is curvilinear. The definition of rectilinear and curvilinear are provided on that page.</div><div class="gmail_default"><br></div><div class="gmail_default">Don't worry about the fact that the data was read off a GRIB file for the dataonmap_3.ncl example. The script will be useful for any curvilinear data.</div><div class="gmail_default"><br></div><div class="gmail_default">If you don't understand the examples or have problems with them, please post back to ncl-talk.</div><div class="gmail_default"><br></div><div class="gmail_default">Good luck,</div><div class="gmail_default"><br></div><div class="gmail_default">--Mary</div><div class="gmail_default"><br></div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small"></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Jul 26, 2018 at 2:48 AM, SARA CORNEJO BUENO <span dir="ltr"><<a href="mailto:sarcorne@ucm.es" target="_blank">sarcorne@ucm.es</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><font color="#0000ff">Hello,</font><div><font color="#0000ff"><br></font></div><div><font color="#0000ff">I'm trying to make a sea level pressure plot from a Regrid; this is part of the code</font></div><div><font color="#0000ff"><br></font></div><div><div><font color="#0000ff">;</font><font color="#000000">;;;;; Set basic resources ;;;;;;</font></div><div><font color="#000000"> res = True</font></div><div><font color="#000000"> res@MainTitle = "Sea Level Pressure_811"</font></div><div><font color="#000000"> time = wrf_user_getvar(a,"time",-1) ; get all time in the file</font></div><div><font color="#000000"> ntime = dimsizes(time) </font></div><div><font color="#000000"> do it = 0,ntime-1,1 ; TIME LOOP</font></div><div><font color="#000000"> res@TimeLabel = time(it) ; Set Valid time to use on plots</font></div><div><font color="#000000">;exit()</font></div><div><font color="#000000"> ; Open workstation</font></div><div><font color="#000000"> wks_type = "png"</font></div><div><font color="#000000"> wks_type@wkWidth = 2000</font></div><div><font color="#000000"> wks_type@wkHeight = 2000</font></div><div><font color="#000000"> wks = gsn_open_wks(wks_type, "/Volumes/SimuPlots/slp_811_" +time(it)+ "")</font></div><div><font color="#000000"><br></font></div><div><font color="#000000"> gsn_define_colormap(wks,"BlAqG<wbr>rYeOrRe")</font></div><div><font color="#000000">;exit()</font></div><div><font color="#000000"><br></font></div><div><font color="#000000">;;;;;; Set plot and map resources ;;;;;;</font></div><div><font color="#000000"> pltres = True</font></div><div><font color="#000000"> mpres = True</font></div><div><font color="#000000"> mpres@mpOutlineBoundarySets = "National"</font></div><div><font color="#000000"> mpres@mpDataBaseVersion = "MediumRes"</font></div><div><font color="#000000"> mpres@mpDataSetName = "Earth..4"</font></div><div><font color="#000000"> mpres@mpGeophysicalLineColor = "Black"</font></div><div><font color="#000000"> mpres@mpGeophysicalLineThickne<wbr>ssF = 4.0</font></div><div><font color="#000000">;exit()</font></div><div><font color="#000000">LON =a->longitude </font></div><div><font color="#000000">;print(lon)</font></div><div><font color="#000000">LAT=a->latitude</font></div><div><font color="#000000">;print(lat)</font></div><div><font color="#000000">LON = LON-360 ;así son longitudes W</font></div><div><font color="#000000">;print(LON)</font></div><div><font color="#000000">LAT = LAT(::-1) ; overwrite: reverse using NCL syntax (para recolocar la latitud (tiene que estar S-N))</font></div><div><font color="#000000">;exit()</font></div><div><font color="#000000"> if (isvar("zoomCoords")) then</font></div><div><font color="#000000"> ; lat2d = wrf_user_Agetvar(a,"latitude",<wbr>0)</font></div><div><font color="#000000"> ;lon2d = wrf_user_getvar(a,"longitude",<wbr>0)</font></div><div><font color="#000000">;exit()</font></div><div><font color="#000000"> ; Coordenadas de la zona del mapa que quiero dibujar</font></div><div><font color="#000000"> mpres@mpLimitMode = "Corners"</font></div><div><font color="#000000"> mpres@mpLeftCornerLatF = zoomCoords(0)</font></div><div><font color="#000000"> mpres@mpLeftCornerLonF = zoomCoords(1)</font></div><div><font color="#000000"> mpres@mpRightCornerLatF = zoomCoords(2)</font></div><div><font color="#000000"> mpres@mpRightCornerLonF = zoomCoords(3)</font></div><div><font color="#000000"> end if</font></div><div><font color="#000000">;exit()</font></div><div><font color="#000000"> ; Get the variables</font></div><div><font color="#000000"><br></font></div><div><font color="#000000"> slp = wrf_user_getvar(a,"SLP",it)</font></div><div><font color="#000000"> slp@units ="hPa"</font></div></div><div><font color="#0000ff"><br></font></div><div><font color="#0000ff"><br></font></div><div><font color="#0000ff">All goes right untill the following part:</font></div><div><br></div><div><font color="#0000ff">the problem appears at the end, on this part of the code:</font></div><div><br></div><div><div>;;;;;; Plotting options for slp ;;;;;</div><div> opts2 = res</div><div> opts2@cnLevelSelectionMode = "ExplicitLevels"</div><div> opts2@gsnContourLineThicknesse<wbr>sScale = 7.0</div><div> opts2@cnLineColor = "Blue"</div><div> opts2@cnLineLabelInterval = 1 ; default = 2</div><div> opts2@cnLineLabelPlacementMode = "constant"</div><div> opts2@cnLineLabelDensityF = 2 ;densidad de info-boxes en una linea de contorno</div><div>; opts2@cnLevels = (/ 980, 984, 988, 990, 994, 998, 1000, 1002, 1004, 1008, 1010, 1012, 1018, 1020, 1022 /)</div><div> opts2@ContourParameters = (/ 980, 1022, 2 /)</div><div> opts2@cnFillOn = False</div><div> contour_slp = wrf_contour(a, wks, slp, opts2)</div><div> delete(opts2)</div><div><br></div><div><font color="#0000ff">It appears these warnings:</font></div><div>
<font color="#0000ff"><span></span>
</font><p class="m_-1402359263384047780m_-4229415893897382798gmail-p1" style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-weight:normal;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;background-color:rgb(255,255,255)"><span class="m_-1402359263384047780m_-4229415893897382798gmail-s1" style="font-variant-ligatures:no-common-ligatures"><font color="#0000ff">warning:Attempt to reference attribute (TITLE) which is undefined</font></span></p>
<p class="m_-1402359263384047780m_-4229415893897382798gmail-p1" style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-weight:normal;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;background-color:rgb(255,255,255)"><span class="m_-1402359263384047780m_-4229415893897382798gmail-s1" style="font-variant-ligatures:no-common-ligatures"><font color="#0000ff">warning:Attempt to reference attribute (DX) which is undefined</font></span></p>
<p class="m_-1402359263384047780m_-4229415893897382798gmail-p1" style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-weight:normal;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;background-color:rgb(255,255,255)"><span class="m_-1402359263384047780m_-4229415893897382798gmail-s1" style="font-variant-ligatures:no-common-ligatures"><font color="#0000ff">warning:Attempt to reference attribute (WEST-EAST_GRID_DIMENSION) which is undefined</font></span></p>
<p class="m_-1402359263384047780m_-4229415893897382798gmail-p1" style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-weight:normal;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;background-color:rgb(255,255,255)"><span class="m_-1402359263384047780m_-4229415893897382798gmail-s1" style="font-variant-ligatures:no-common-ligatures"><font color="#0000ff">warning:Attempt to reference attribute (SOUTH-NORTH_GRID_DIMENSION) which is undefined</font></span></p>
<p class="m_-1402359263384047780m_-4229415893897382798gmail-p1" style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-weight:normal;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;background-color:rgb(255,255,255)"><span class="m_-1402359263384047780m_-4229415893897382798gmail-s1" style="font-variant-ligatures:no-common-ligatures"><font color="#0000ff">warning:Attempt to reference attribute (BOTTOM-TOP_GRID_DIMENSION) which is undefined</font></span></p><p class="m_-1402359263384047780m_-4229415893897382798gmail-p1" style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-weight:normal;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;background-color:rgb(255,255,255)"><span class="m_-1402359263384047780m_-4229415893897382798gmail-s1" style="font-variant-ligatures:no-common-ligatures"><font color="#0000ff"><br></font></span></p><p class="m_-1402359263384047780m_-4229415893897382798gmail-p1" style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-weight:normal;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;background-color:rgb(255,255,255)"><span class="m_-1402359263384047780m_-4229415893897382798gmail-s1" style="font-variant-ligatures:no-common-ligatures"><font color="#0000ff">The code continues like this</font></span></p>
</div><div><br></div><div> ; MAKE PLOTS</div><div> if (isvar("zoomCoords")) then</div><div> pltres@LatLonOverlay = True</div><div> end if</div><div> plot = wrf_map_overlays(a, wks, (/contour_slp/), pltres, mpres)</div><div><br></div><div> end do ; TIME LOOP </div><div>;end do ; DAY LOOP </div><div>end</div></div><div><span style="color:rgb(0,0,255);font-family:Menlo;font-size:11px;font-variant-ligatures:no-common-ligatures;background-color:rgb(255,255,255);text-decoration-style:initial;text-decoration-color:initial;float:none;display:inline">And then the error says: </span><span style="color:rgb(0,0,255);font-family:Menlo;font-size:11px;font-variant-ligatures:no-common-ligatures;background-color:rgb(255,255,255);text-decoration-style:initial;text-decoration-color:initial;float:none;display:inline"></span><span style="color:rgb(0,0,255);font-family:Menlo;font-size:11px;font-variant-ligatures:no-common-ligatures;background-color:rgb(255,255,255);text-decoration-style:initial;text-decoration-color:initial"></span><span style="color:rgb(0,0,255);font-family:Menlo;font-size:11px;font-variant-ligatures:no-common-ligatures;background-color:rgb(255,255,255);text-decoration-style:initial;text-decoration-color:initial;float:none;display:inline"></span><p class="m_-1402359263384047780m_-4229415893897382798gmail-p1" style="color:rgb(0,0,255);font-variant-numeric:normal;font-variant-east-asian:normal;font-weight:normal;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;text-decoration-style:initial;text-decoration-color:initial;margin:0px;background-color:rgb(255,255,255)"><span class="m_-1402359263384047780m_-4229415893897382798gmail-s1" style="font-variant-ligatures:no-common-ligatures">(0)<span class="m_-1402359263384047780m_-4229415893897382798gmail-Apple-tab-span" style="white-space:pre-wrap"> </span>wrf_map: Error: no MAP_PROJ attribute in input file</span></p><p class="m_-1402359263384047780m_-4229415893897382798gmail-p1" style="color:rgb(0,0,255);font-variant-numeric:normal;font-variant-east-asian:normal;font-weight:normal;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;text-decoration-style:initial;text-decoration-color:initial;margin:0px;background-color:rgb(255,255,255)"><span class="m_-1402359263384047780m_-4229415893897382798gmail-s1" style="font-variant-ligatures:no-common-ligatures">fatal:overlay: missing value as input, can't continue</span></p></div><div><br></div><div>
<span></span>
<p class="m_-1402359263384047780m_-4229415893897382798gmail-p1" style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-weight:normal;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(255,255,255)"><span class="m_-1402359263384047780m_-4229415893897382798gmail-s1" style="font-variant-ligatures:no-common-ligatures">fatal:["Execute.c":8640]:Execu<wbr>te: Error occurred at or near line 83 in file Plot_SLP_simus.ncl,</span></p><p class="m_-1402359263384047780m_-4229415893897382798gmail-p1" style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-weight:normal;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(255,255,255)"><span class="m_-1402359263384047780m_-4229415893897382798gmail-s1" style="font-variant-ligatures:no-common-ligatures"><br></span></p><p class="m_-1402359263384047780m_-4229415893897382798gmail-p1" style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-weight:normal;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(255,255,255)"><span class="m_-1402359263384047780m_-4229415893897382798gmail-s1" style="font-variant-ligatures:no-common-ligatures">line 83 corresponds to:</span></p><p class="m_-1402359263384047780m_-4229415893897382798gmail-p1" style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-weight:normal;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(255,255,255)"><span class="m_-1402359263384047780m_-4229415893897382798gmail-s1" style="font-variant-ligatures:no-common-ligatures"><br></span></p><p class="m_-1402359263384047780m_-4229415893897382798gmail-p1" style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-weight:normal;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(255,255,255)"><span class="m_-1402359263384047780m_-4229415893897382798gmail-s1" style="font-variant-ligatures:no-common-ligatures"><span style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:small;font-variant-ligatures:normal;background-color:rgb(255,255,255);text-decoration-style:initial;text-decoration-color:initial;float:none;display:inline"><span> </span>plot = wrf_map_overlays(a, wks, (/contour_slp/), pltres, mpres)</span> </span></p><p class="m_-1402359263384047780m_-4229415893897382798gmail-p1" style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-weight:normal;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(255,255,255)"><span class="m_-1402359263384047780m_-4229415893897382798gmail-s1" style="font-variant-ligatures:no-common-ligatures"><br></span></p><p class="m_-1402359263384047780m_-4229415893897382798gmail-p1" style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-weight:normal;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(255,255,255)"><span class="m_-1402359263384047780m_-4229415893897382798gmail-s1" style="font-variant-ligatures:no-common-ligatures">What I could understand is that, now, I don't have LAT and LONG coordinates, but grid points. But I don't know how to call them.</span></p><p class="m_-1402359263384047780m_-4229415893897382798gmail-p1" style="margin:0px;font-variant-numeric:normal;font-variant-east-asian:normal;font-weight:normal;font-stretch:normal;font-size:11px;line-height:normal;font-family:Menlo;color:rgb(0,0,0);background-color:rgb(255,255,255)"><span class="m_-1402359263384047780m_-4229415893897382798gmail-s1" style="font-variant-ligatures:no-common-ligatures">Thank you very much in advance.</span></p>
<br></div></div>
<br>______________________________<wbr>_________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/<wbr>mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>