[ncl-talk] fatal:["Execute.c":8640]:Execute: Error occurred at or near line 897 in file $NCARG_ROOT/lib/ncarg/nclscripts/csm/crop.ncl

Dennis Shea shea at ucar.edu
Wed Aug 23 10:30:08 MDT 2017


Guido's response is correct.

Still, a 'trick' when encountering a malloc error is to use a small subset
of the data

Rather than:

  tasmax = a->tasmax

  tasmin = b->tasmin

  Pr_annual_total  = c->pr

  lon = a->longitude

  lat = a->latitude
  time = a->time

Use
  NT = 100                                    ; arbitrary

  tasmax = a->tasmax(0:NT,:,:)

  tasmin  = b->tasmin(0:NT,:,:)

  Pr_annual_total  = c->pr(0:NT,:,:)

  lon = a->longitude

  lat = a->latitude


Note: rather than reading the 'time' variable directly
  time = a->time(0:NT)
I suggest using the 'time' coordinate associated with the variable
   time = tasmax&time


====


On Wed, Aug 23, 2017 at 5:15 AM, Guido Cioni <guidocioni at gmail.com> wrote:

> You have to increase the maximum memory (ws maximum size) used by ncl to
> plot. See here (https://www.ncl.ucar.edu/Document/Graphics/hlures.shtml).
>
> Additionally, depending on the machine, you may need to increase the
> maximum memory available for ncl at execution time (like ulimit in Linux).
> Cheers
>
> Il 23 ago 2017 1:09 PM, "Beáta Szabó-Takács" <szabo.b at czechglobe.cz> ha
> scritto:
>
>> Dear NCL Users,
>>
>>
>>
>> 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:
>>
>> fatal:NclMalloc Failed:[errno=12]
>>
>> fatal:New: could not create new array:[errno=12]
>>
>> fatal:["Execute.c":8640]:Execute: Error occurred at or near line 897 in
>> file $NCARG_ROOT/lib/ncarg/nclscripts/csm/crop.ncl
>>
>>
>>
>> fatal:["Execute.c":8640]:Execute: Error occurred at or near line 162 in
>> file lambert_PET.ncl
>>
>>
>>
>> ; Main code.
>>
>> ;----------------------------------------------------------------------
>>
>> begin
>>
>> ;************************************************
>>
>> ; read in netCDF file
>>
>> ;************************************************
>>
>>   a = addfile("E-OBS_tasmax.nc","r")
>>
>>   b = addfile("E-OBS_tasmin.nc","r")
>>
>>   c = addfile("E-OBS_1971-2000_ys.nc","r")
>>
>>
>>
>>
>>
>>   tasmax = a->tasmax
>>
>>   tasmin = b->tasmin
>>
>>   Pr_annual_total  = c->pr
>>
>>   lon = a->longitude
>>
>>   lat = a->latitude
>>
>>   time = a->time
>>
>>   minlat             =  30                ; min lat to mask
>>
>>   maxlat             =  75                ; max lat to mask
>>
>>   minlon             = -20                ; min lon to mask
>>
>>   maxlon             =  40                ; max lon to mask
>>
>>
>>
>>   nlat = dimsizes(lat)
>>
>>   mlon = dimsizes(lon)
>>
>>
>>
>>   lat_new=conform_dims( (/nlat,mlon/), lat, 0)
>>
>>
>>
>>   time2 = calendar_decode2(time,-5)
>>
>>   doy = day_of_year(time2(:,0),time2(:,1),time2(:,2))
>>
>>    print(doy)
>>
>>     radext = radext_fao56(doy, lat_new, 0)
>>
>>     radext = where(ismissing(radext), 0, radext)
>>
>>
>>
>>     evtH_0 = refevt_hargreaves_fao56( tmin, tmax, radext, (/0,0,0/) )
>>
>>
>>
>>    print(evtH_0)
>>
>>    opt = True
>>
>>    opt at nval_crit = 12
>>
>>
>>
>>    PET_m = calculate_monthly_values (evtH_0, "sum", 0,opt)
>>
>>    print(PET_m)
>>
>>
>>
>>     PET_annual_total = month_to_annual(PET_m, 0)
>>
>> ;    Pr_annual_total = month_to_annual(Pr_m, 0)
>>
>>
>>
>>       AI = Pr_annual_total/PET_annual_total
>>
>>      AI2 = dim_avg_n(AI,0)
>>
>>   IM =  new(dimsizes(AI2),typeof(AI2))
>>
>>
>>
>>    IM = where(AI2 .lt. 0.03,1.0,IM)
>>
>>    IM = where(AI2 .ge. 0.03 .and. AI2 .lt. 0.2,2.0,IM)
>>
>>    IM = where(AI2 .ge. 0.2 .and. AI2 .lt. 0.5,3.0,IM)
>>
>>    IM = where(AI2 .ge. 0.5 .and. AI2 .lt. 0.65,4.0,IM)
>>
>>    IM = where(AI2 .ge. 0.65,5.0,IM)
>>
>>
>>
>> The error refer to the radext_fao56 in my script which is attached as
>> well. Could someone suggest a solution?
>>
>>
>>
>> Thank you for your help in advance!
>>
>>
>>
>> Best regards,
>>
>> Beata
>>
>> _______________________________________________
>> ncl-talk mailing list
>> ncl-talk at ucar.edu
>> List instructions, subscriber options, unsubscribe:
>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>
>>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at ucar.edu
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170823/8fbed4a0/attachment.html 


More information about the ncl-talk mailing list