[ncl-talk] Error for missing value

[azizul haque]

Hi NCL Users,
I am computing daily annual cycle for CMIP5 model with rcp85 and getting following error meaages

**warning:tofloat: A bad value was passed to (string) tofloat, input strings must contain numeric digits, replacing with missing value
warning:tofloat: A bad value was passed to (string) tofloat, input strings must contain numeric digits, replacing with missing value
fatal:fspan: Missing value detected in input, can't continue
fatal:["Execute.c":8567]:Execute: Error occurred at or near line 5371 in file $NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl

fatal:["Execute.c":8567]:Execute: Error occurred at or near line 5482 in file $NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl

fatal:["Execute.c":8567]:Execute: Error occurred at or near line 105 in file cmip5_dac_10m.ncl***

Please find the script herewith this mail below

Regards
Azizul
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

load
"$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl" 
load
"$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl" 
load
"$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl" 
;**************************************************************
; User Input
;***************************************************************
   begin
   lat1     = 20.0    ; Region
   lat2     = 35.0
   lon1     = 70.0    
   lon2     = 90.0    
;***************************************************************
; End User Input
;***************************************************************
   a          =addfile("/Volumes/BOULDER/Projects/CMIP5/Daily/RCP85/CCSM4/r2i1p1/pr_day_CCSM4_rcp85_r2i1p1_20800101-21001231_sa.nc","r")
 time0     = a->time
 pr0       = a->pr(:,{lat1:lat2},{lon1:lon2}) ;
(time,lat,lon)  (0,1,2)
;Convert kg to mm
 pr0       = pr0*86400.
 pr0 at units =
"mm/day"
 printVarSummary(pr0)
 printMinMax(pr0,1)
 clat0      = latRegWgt(pr0&lat, typeof(pr0), 0)
; create lat weights
 pr0Area    = wgt_areaave_Wrap(pr0, clat0, 1, 0)  ; weighted area avg
 ymdStrt0   = 20800101 ; start yyyymmdd
 ymdLast0   = 21001231       ;last
 yrStrt     = ymdStrt0/10000
 yrLast     = ymdLast0/10000
 nhar           = 4  ; no of fourier comp
 var                     = "pr"
  TIME0     = cd_calendar(time0,0)   ; type of float
  year0    = floattointeger( TIME0(:,0) )
  month0    = floattointeger( TIME0(:,1) )
  day0      = floattointeger( TIME0(:,2) )
  ddd0        = day_of_year(year0, month0, day0)
  yyyyddd0  = year0*1000 + ddd0
  nday0     = 365
  pday0    = new ( nday0, typeof(pr0),
getFillValue(pr0))
 
   do nd0  =0,nday0-2
     iday0      = ind(ddd0.eq.(nd0+1))
     pday0(nd0)         = avg(pr0Area(iday0))
     delete(iday0)                       ; may change next irruption
   end do
  pday0(nday0-1)        = 0.5*(pday0(0)+pday0(nday0-2))
  pr0AvgRegion  = avg(pr0Area)                    ; regional avg
 
;************************************************
; Create smooth annual cycle
;************************************************                        
   cf0       = ezfftf(pday0)
   cf0(:,3 ) =
0.75*cf0(:,3)                        ;
arbitrary wgts
   cf0(:,4 ) =
0.50*cf0(:,4)
   cf0(:,5 ) =
0.25*cf0(:,5)
   cf0(:,6:) =
0.0                                 ; no
contribution
   pday0_smth= ezfftb
(cf0, cf0 at xbar)                ;
reconstruct
 
;************************************************
; Create plot 
;************************************************
   ntim0                       = dimsizes(time0)
 
   PDAY0       = new( (/2,nday0/), typeof(pday0),
"No_FillValue")
   PDAY0(0,:)  = pday0
   PDAY0(1,:)  = pday0_smth
 
wks         =
gsn_open_wks("png", "Annual_Cycle_21C")
 
 
res                   = True        ; plot mods desired
   res at gsnDraw           = False
   res at gsnFrame          = False
  
   res at trXMinF           = 0           ; max value on x-axis
   res at trXMaxF           = nday0+1      ; max value on x-axis
 
   res at vpHeightF                                    = 0.4                     ; resize           
   res at vpWidthF                                      = 0.8                 
   res at tmXBMode                                     = "Explicit" 
   res at trXMinF           =  1
   res at trXMaxF           = 360
      
   res at xyLineThicknesses
= (/2.5, 2.5/)              ; make 2nd
lines thicker
   res at xyLineColors      =
(/"NavyBlue","yellow"/)          ; change line color
   res at xyMonoDashPattern = True                      ; all solid
 
   res at xyLineThicknesses = (/2.5 , 2.5/)
   res at xyLineColors      =
(/"NavyBlue","yellow"/)  ; change line color
   res at xyMonoDashPattern = True
   res at tiXAxisString     = "day
of year"     ;
 
 
; Plotting procedure
  plot=
new(2,graphic)
 
  res at gsnLeftString                          = "CCSM4" 
  res at gsnYRefLine                            =
"pr0AvgRegion"           ;
create a reference line       
  ;res at tiYAxisString             = "pr0 (mm/day)"    ; y-axis label    
  res at tiXAxisString                                       = "day
of year"     ; x      
  res at gsnRightString                        = "Areal
Mean="+sprintf("%4.2f", pr0AvgRegion)+" mm/day"      
  res at gsnCenterString                     = "regName"
  res at gsnLeftStringFontHeightF                 =
0.010
  res at gsnRightStringFontHeightF   =
0.010
  res at gsnCenterStringFontHeightF            =
0.010
  plot(0)                                              = gsn_csm_y(wks, PDAY0,
res)
 
; Paneling the resources
  resP                   = True
  resP at txString          = "Daily Annual Cycle for
Precipitation of 10 CMIP5 Models for 21 Century with RCP85"
  resP at gsnFrame          = False
  resP at txFontHeightF      = 0.0110
  resP at gsnMaximize        =
True                ; maximize 
  resP at gsnPanelFigureStrings=(/"a)","b)","c)","d)","e)","f)","g)","h)","i)","j)"/)
  resP at gsnPanelFigureStringsFontHeightF=0.006
  resP at txFontHeightF=0.01
  gsn_panel(wks,
plot,(/5,2/),resP) ; not draw as one plot
end
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20140921/7adbef8c/attachment.html