<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-2">
<meta name="Generator" content="Microsoft Word 15 (filtered medium)">
<style><!--
/* Font Definitions */
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0cm;
margin-bottom:.0001pt;
font-size:11.0pt;
font-family:"Calibri",sans-serif;
mso-fareast-language:EN-US;}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:#0563C1;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:#954F72;
text-decoration:underline;}
span.EmailStyle17
{mso-style-type:personal-compose;
font-family:"Calibri",sans-serif;
color:windowtext;}
.MsoChpDefault
{mso-style-type:export-only;
mso-fareast-language:EN-US;}
@page WordSection1
{size:612.0pt 792.0pt;
margin:70.85pt 70.85pt 70.85pt 70.85pt;}
div.WordSection1
{page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]-->
</head>
<body lang="HU" link="#0563C1" vlink="#954F72">
<div class="WordSection1">
<p class="MsoNormal">Dear NCL Users,<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">I would like to calculate Aridity index base on Hargreaves ETo equation. I use daily maximum and minimum temperature data from 1971 to 2000 and yearly sum precipitation (30 years). The dimension of latitude is (201x464) and the day of year
is (10957) vector for radext_fao56 function. When I run the scipt I get the following error message:<o:p></o:p></p>
<p class="MsoNormal">fatal:NclMalloc Failed:[errno=12]<o:p></o:p></p>
<p class="MsoNormal">fatal:New: could not create new array:[errno=12]<o:p></o:p></p>
<p class="MsoNormal">fatal:["Execute.c":8640]:Execute: Error occurred at or near line 897 in file $NCARG_ROOT/lib/ncarg/nclscripts/csm/crop.ncl<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">fatal:["Execute.c":8640]:Execute: Error occurred at or near line 162 in file lambert_PET.ncl<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">; Main code.<o:p></o:p></p>
<p class="MsoNormal">;----------------------------------------------------------------------<o:p></o:p></p>
<p class="MsoNormal">begin<o:p></o:p></p>
<p class="MsoNormal">;************************************************<o:p></o:p></p>
<p class="MsoNormal">; read in netCDF file<o:p></o:p></p>
<p class="MsoNormal">;************************************************<o:p></o:p></p>
<p class="MsoNormal"> a = addfile("E-OBS_tasmax.nc","r")<o:p></o:p></p>
<p class="MsoNormal"> b = addfile("E-OBS_tasmin.nc","r")<o:p></o:p></p>
<p class="MsoNormal"> c = addfile("E-OBS_1971-2000_ys.nc","r")<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"> tasmax = a->tasmax<o:p></o:p></p>
<p class="MsoNormal"> tasmin = b->tasmin<o:p></o:p></p>
<p class="MsoNormal"> Pr_annual_total = c->pr<o:p></o:p></p>
<p class="MsoNormal"> lon = a->longitude<o:p></o:p></p>
<p class="MsoNormal"> lat = a->latitude<o:p></o:p></p>
<p class="MsoNormal"> time = a->time<o:p></o:p></p>
<p class="MsoNormal"> minlat = 30 ; min lat to mask<o:p></o:p></p>
<p class="MsoNormal"> maxlat = 75 ; max lat to mask<o:p></o:p></p>
<p class="MsoNormal"> minlon = -20 ; min lon to mask<o:p></o:p></p>
<p class="MsoNormal"> maxlon = 40 ; max lon to mask<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"> nlat = dimsizes(lat)<o:p></o:p></p>
<p class="MsoNormal"> mlon = dimsizes(lon)<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"> lat_new=conform_dims( (/nlat,mlon/), lat, 0)<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"> time2 = calendar_decode2(time,-5)<o:p></o:p></p>
<p class="MsoNormal"> doy = day_of_year(time2(:,0),time2(:,1),time2(:,2))<o:p></o:p></p>
<p class="MsoNormal"> print(doy)<o:p></o:p></p>
<p class="MsoNormal"> radext = radext_fao56(doy, lat_new, 0)<o:p></o:p></p>
<p class="MsoNormal"> radext = where(ismissing(radext), 0, radext)<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"> evtH_0 = refevt_hargreaves_fao56( tmin, tmax, radext, (/0,0,0/) )<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"> print(evtH_0)<o:p></o:p></p>
<p class="MsoNormal"> opt = True<o:p></o:p></p>
<p class="MsoNormal"> opt@nval_crit = 12<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"> PET_m = calculate_monthly_values (evtH_0, "sum", 0,opt)<o:p></o:p></p>
<p class="MsoNormal"> print(PET_m)<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"> PET_annual_total = month_to_annual(PET_m, 0)<o:p></o:p></p>
<p class="MsoNormal">; Pr_annual_total = month_to_annual(Pr_m, 0)<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"> AI = Pr_annual_total/PET_annual_total<o:p></o:p></p>
<p class="MsoNormal"> AI2 = dim_avg_n(AI,0)<o:p></o:p></p>
<p class="MsoNormal"> IM = new(dimsizes(AI2),typeof(AI2))<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"> IM = where(AI2 .lt. 0.03,1.0,IM)<o:p></o:p></p>
<p class="MsoNormal"> IM = where(AI2 .ge. 0.03 .and. AI2 .lt. 0.2,2.0,IM)<o:p></o:p></p>
<p class="MsoNormal"> IM = where(AI2 .ge. 0.2 .and. AI2 .lt. 0.5,3.0,IM)<o:p></o:p></p>
<p class="MsoNormal"> IM = where(AI2 .ge. 0.5 .and. AI2 .lt. 0.65,4.0,IM)<o:p></o:p></p>
<p class="MsoNormal"> IM = where(AI2 .ge. 0.65,5.0,IM)<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">The error refer to the radext_fao56 in my script which is attached as well. Could someone suggest a solution?<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Thank you for your help in advance!<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Best regards,<o:p></o:p></p>
<p class="MsoNormal">Beata<o:p></o:p></p>
</div>
</body>
</html>