<div dir="ltr">Hi Sara,<div><br></div><div>It's throwing an error because your value for pressure at Fortran index (1,1,1) is 0 hPa. This is an unrealistic value for surface pressure. There appears to be a problem with your data. </div><div><br></div><div>In NCL, do a wrf_user_getvar(your_file, "pressure", 0) and check that the value returned for index (0,0,0) is something close to ~1000 hPa (or less if it's in the mountains). If not, something is wrong with your data, because WRF generally does not produce bogus or missing values unless something went wrong during the run (also check error logs for your WRF run). Also make sure you do this check for each of the files and times in your data set.</div><div><br></div><div>Hope this helps,</div><div><br></div><div>Bill</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Sun, Apr 15, 2018 at 10:17 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">Hi! I run my NCL script, and <pre style="color:rgb(0,0,0);font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;word-spacing:0px;text-decoration-style:initial;text-decoration-color:initial"><i>I got an error:
</i>><i>
</i>><i> Troubles finding level 100 above ground.
</i>><i> Problems first occur at ( 1, 1)
</i>><i> Surface pressure = 0.0 hPa.
</i>><i> STOP Error_in_finding_100_hPa_up</i></pre><pre style="color:rgb(0,0,0);font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;word-spacing:0px;text-decoration-style:initial;text-decoration-color:initial"><i><br></i></pre><pre style="color:rgb(0,0,0);font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;word-spacing:0px;text-decoration-style:initial;text-decoration-color:initial"><i>And this is the start of the script:</i></pre><pre style="text-align:start;text-indent:0px;text-decoration-style:initial;text-decoration-color:initial"><i><font color="#000000">load "$NCARG_ROOT/lib/ncarg/<wbr>nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/<wbr>nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/<wbr>nclscripts/wrf/WRFUserARW.ncl"
load "/Volumes/DiscoSara/ncl-6.4.0/<wbr>lib/ncarg/nclscripts/csm/gsn_<wbr>code.ncl"
load "/Volumes/DiscoSara/ncl-6.4.0/<wbr>lib/ncarg/nclscripts/csm/gsn_<wbr>csm.ncl"
load "/Volumes/DiscoSara/ncl-6.4.0/<wbr>lib/ncarg/nclscripts/csm/<wbr>contributed.ncl"
load "/Volumes/DiscoSara/ncl-6.4.0/<wbr>lib/ncarg/nclscripts/csm/shea_<wbr>util.ncl"
load "/Volumes/DiscoSara/ncl-6.4.0/<wbr>lib/ncarg/nclscripts/wrf/<wbr>WRFUserARW.ncl"<span style="font-style:normal"><br></span></font></i></pre><pre style="text-align:start;text-indent:0px;text-decoration-style:initial;text-decoration-color:initial"><i><font color="#000000">begin<br></font></i></pre><pre style="text-align:start;text-indent:0px;text-decoration-style:initial;text-decoration-color:initial"><i><font color="#000000">diri= "/Volumes/DiscoSara/"
files = systemfunc("ls " + diri + "wrfout_d02_2010*") ;para llamar a todas las salidas de golpe
a = addfiles(files,"r")
times = wrf_user_getvar(a,"times",-1) ; get all times in the files
ntimes = dimsizes(times) ;number of times in the files
do it = 0,ntimes-1,6 ; TIME LOOP
;;;;; Get the variables
;;;Presión, Latitud y Longitud
lat = wrf_user_getvar(a,"XLAT",it)
lon = wrf_user_getvar(a,"XLONG",it) ;it va de 0 a ntimes-1 cada 6 horas
slp = wrf_user_getvar(a,"slp",it)
min_slp = min(slp) ;mín de presión
SLP = ndtooned(slp) ;ndtooned: Convierte una matriz multi-D en una 1D
LAT = ndtooned(lat) ;
LON = ndtooned(lon)
ij = ind(min_slp .eq. SLP) ;ij: calcula la posición dentro del vector SLP en la que la presión es mínima (min_slp)
;ind se usa para vectores en una dimensión.
min_LAT = LAT(ij) ;
min_LON = LON(ij) ;
</font></i></pre><div><i><font color="#000000"><br></font></i></div><div><i><font color="#000000"><br></font></i></div><div><font color="#000000">I think the problem is on "slp", but I don't know how to solve it</font></div><div><font color="#000000"><br></font></div><div><font color="#000000">Thanks</font></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>