[ncl-talk] Could not Coerce data into double

[Jacob Alberto Garcia]

I found that it took too much memory hence the error.
But what can I do so that ncl gets the values only in the lowest model
layer, and then continues to calculate the average through time?

Thank you,
Jacob

On Wed, Aug 24, 2016 at 5:10 PM, Jacob Alberto Garcia <
jacob_garcia at dlsu.edu.ph> wrote:

> Good day!
>
> After running this script, I get this error wherein it cannot continue:
>
> fatal:NclMalloc Failed:[errno=12]
> fatal:dim_avg_n: Could not coerce input data to double, can't continue
> fatal:["Execute.c":8578]:Execute: Error occurred at or near line 2335 in
> file /usr/local/lib/ncarg/nclscripts/csm/contributed.ncl
>
>
> This is my script
>
> ===============================================
>
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
> ;load "./WRFUserARW.ncl"
>
>
>   diri = "./"
>   jake = (/ "d03_DEC_YSU_F.nc", "d03_JAN_YSU_F.nc", "d03_FEB_YSU_F.nc" /)
>    j = addfile("d03_DEC_YSU_F.nc", "r")
>   pthi = diri+jake
> fa  = addfiles(pthi, "r")
> pm25 = fa[:]->PM2_5_DRY
> pm25_avg = dim_avg_n_Wrap(pm25, 0)
> printVarSummary(pm25_avg)   ; south_north | ??] x [west_east | ??]
>
>
>
> ; We generate plots, but what kind do we prefer?
>   type = "pdf"
> ; type = "pdf"
> ; type = "ps"
> ; type = "ncgm"
>
> wks = gsn_open_wks(type,"PM25")
>
>
>
>   gsn_define_colormap(wks, "WhBlGrYeRe")
>
>
>
>
>   pltres = True
>   mpres = True
>   res = True
>   mpres at mpDataBaseVersion = "HighRes"
>   mpres at mpGridAndLimbOn = False
>   res at gsnMaximize = True
>   res at cnFillOn = True
>   res at gsnSpreadColors = True
>   res at cnLevelSelectionMode = "ManualLevels"
>  ; res at cnMinLevelValF  = -1.8            ; set the minimum contour level
>   ;res at cnMaxLevelValF  = 1.8            ; set the maximum contour
> level
>   ;res at cnLevelSpacingF = 0.1
>   res at cnLinesOn = False
>   res at cnLineLabelsOn = False
>   mpres at mpGeophysicalLineColor = "Black"
>
>
>
>
> plot    = gsn_csm_contour_map(wks,pm25_avg,res)
>
>      contour_DIFF = wrf_contour(j,wks,pm25_avg,res)
>
>      plot = wrf_map_overlays(j,wks,(/contour_DIFF/),pltres,mpres)
>
> ==================================================
>
> I am not sure anymore as to what part im going wrong. Hoping for your help
>
> Jacob
>
>
>
>
>
>

-- 
 <http://www.dlsu.edu.ph> <https://www.facebook.com/DLSU.Manila.100> 
<http://instagram.com/dlsu> <https://twitter.com/dlsumanila> 
<http://dlsumanila.tumblr.com/> <http://iblog.lasalle.ph/>

DISCLAIMER AND CONFIDENTIALITY NOTICE 
The information contained in this e-mail, including those in its 
attachments, is confidential and intended only for the person(s) or 
entity(ies) to which it is addressed. If you are not an intended recipient, 
you must not read, copy, store, disclose, distribute this message, or act 
in reliance upon the information contained in it. If you received this 
e-mail in error, please contact the sender and delete the material from any 
computer or system. Any views expressed in this message are those of the 
individual sender and may not necessarily reflect the views of De La Salle 
University. 
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20160824/07db074e/attachment.html