FW: passing values

Dear all
Sorry I post it again, nobody answers me yet. I try to fix my error, it did not show the severe error anymore only warning message that
Somebody please give me idea.
Thanks
Jam
.
(66) 290.64
(67) 289.83
warning:A bad value was passed to stringtointeger, input strings must contain numeric digits, replacing with missing value
warning:A bad value was passed to stringtointeger, input strings must contain numeric digits, replacing with missing value
warning:A bad value was passed to stringtointeger, input strings must contain numeric digits, replacing with missing value
warning:A bad value was passed to stringtointeger, input strings must contain numeric digits, replacing with missing value
warning:A bad value was passed to stringtointeger, input strings must contain numeric digits, replacing with missing value
warning:A bad value was passed to stringtointeger, input strings must contain numeric digits, replacing with missing value
warning:A bad value was passed to stringtointeger, input strings must contain numeric digits, replacing with missing value
warning:A bad value was passed to stringtointeger, input strings must contain numeric digits, replacing with missing value
warning:A bad value was passed to stringtointeger, input strings must contain numeric digits, replacing with missing value
warning:A bad value was passed to stringtointeger, input strings must contain numeric digits, replacing with missing value
warning:A bad value was passed to stringtointeger, input strings must contain numeric digits, replacing with missing value
warning:A bad value was passed to stringtointeger, input strings must contain numeric digits, replacing with missing value
warning:A bad value was passed to stringtointeger, input strings must contain numeric digits, replacing with missing value
warning:A bad value was passed to stringtointeger, input strings must contain numeric digits, replacing with missing value
warning:A bad value was passed to stringtointeger, input strings must contain numeric digits, replacing with missing value
warning:A bad value was passed to stringtointeger, input strings must contain numeric digits, replacing with missing value
warning:A bad value was passed to stringtointeger, input strings must contain numeric digits, replacing with missing value
warning:A bad value was passed to stringtointeger, input strings must contain numeric digits, replacing with missing value

Variable: ntim
Type: integer
Total Size: 4 bytes
            1 values
Number of Dimensions: 1
Dimensions and sizes: [1]
Coordinates:
(0) 68
warning:lbLabelFont is not a valid resource in ihop_site2_xy.PlotManager at this time
warning:mpShapeMode is not a valid resource in ihop_site2_xy at this time

This is my code
;************************************************
; open file and read in data
;************************************************

diri = "/ptmp/ucharuso/jam/WRF3.1_RC/WRFV3/run/" ; input directory
fili = systemfunc("cd "+diri +"; ls *wrfout_d02*")
fili = fili +".nc" ; explicitly add file extension
f = addfiles(diri+fili, "r") ; refers to *all* files
print (f)
 ListSetType (f, "cat") ; concatenate or "merge" (default)

  T = f[:]->T2(:,109,82)
  T@_FillValue = 9.96921e+36
  U = f[:]->U10(:,109,82)
  V = f[:]->V10(:,109,82)
  WIND = sqrt(U*U + V*V)
  TH = f[:]->TH2(:,109,82)
  PS = f[:]->PSFC(:,109,82)
 PS = PS*.01
  TS = f[:]->TSK(:,109,82)
  RANC = f[:]->RAINC(:,109,82)
  RAIN = f[:]->RAINNC(:,109,82)
  TOTRAIN = RANC + RAIN
  QF = f[:]->QFX(:,109,82)
  HF = f[:]->HFX(:,109,82)
  L = f[:]->LH(:,109,82)
 TM = f[:]->TMN(:,109,82)
 PBL = f[:]->PBLH(:,109,82)

lat = f[:]->XLAT(0,109,82)
lon = f[:]->XLONG(0,109,82)
 SM = f[:]->SMOIS(:,0,109,82)
 ST = f[:]->TSLB(:,0,109,82)
print(T)
; lat = 37.3781 and Lon=-98.1636 for ihop site 5.

;print(lat)
;print(lon)
ivg_usgs = f[:]->IVGTYP(:,109,82)
;print(ivg_usgs)
diri2 = "/ptmp/ucharuso/jam/LLJ/WRF_LDAS_USGS/" ; input directory
fili2 = systemfunc("cd "+diri2 +"; ls *wrfout_d02*")
fili2 = fili2 +".nc" ; explicitly add file extension
f2 = addfiles(diri2+fili2, "r") ; refers to *all* files
print (f2)
 ListSetType (f2, "cat") ; concatenate or "merge" (default)

  T1 = f2[:]->T2(:,109,82)
  U1 = f2[:]->U10(:,109,82)
  V1 = f2[:]->V10(:,109,82)
  WIND1 = sqrt(U1*U1 + V1*V1)
  TH1 = f2[:]->TH2(:,109,82)
  PS1 = f2[:]->PSFC(:,109,82)
 PS1 = PS1*.01
  TS1 = f2[:]->TSK(:,109,82)
  RANC1 = f2[:]->RAINC(:,109,82)
  RAIN1 = f2[:]->RAINNC(:,109,82)
  TOTRAIN1 = RANC1 + RAIN1
  QF1 = f2[:]->QFX(:,109,82)
  HF1 = f2[:]->HFX(:,109,82)
  L1 = f2[:]->LH(:,109,82)
 TM1 = f2[:]->TMN(:,109,82)
 PBL1 = f2[:]->PBLH(:,109,82)
 T1A = T1(0:67)
 print(T1A)
lat1 = f2[:]->XLAT(0,109,82)
lon1 = f2[:]->XLONG(0,109,82)
 SM1 = f2[:]->SMOIS(:,0,109,82)
 ST1 = f2[:]->TSLB(:,0,109,82)
ivg_modis = f2[:]->IVGTYP(:,109,82)
;print(ivg_modis)
 g1 = asciiread("IHOP_site21.txt",(/146,18/), "float")

  T11 = g1(:,5)+273.15
  T12 = T11(0:145:2)
  T13 = T12(0:67)
  T11@_FillValue = 9999
  T12@_FillValue = 9999
        HF11 = g1(:,16)
  HF11@_FillValue = 9999
  HF12 = HF11(0:145:2)
 print(T13)
  L11 = g1(:,15)
  L11@_FillValue = 9999
  L12 = L11(0:145:2)
  WIND11 = g1(:,7)
  WIND11@_FillValue = 9999
  WIND12 = WIND11(0:145:2)
  PS11 = g1(:,6)
  PS12 = PS11(0:145:2)
  PS11@_FillValue = 9999
 ; MR11 = g1(:,3)
 ;SFCTEMP11 = g1(:,5)+273.15
; print(HF12)
  ntim= dimsizes(T11)
 HF12@_FillValue = 9999
 L12@_FillValue = 9999

;print("usgs="+L+"modis="+L1+"obs="+L11)

;print("usgs="+ivg_usgs+"modis="+ivg_modis)

;************************************************
; Read all the times. Convert to units for plotting
;************************************************

 Times = chartostring(f[:]->Times) ; Times is of type character
 Time = WRF_Times2Udunits_c(f[:]->Times, 0) ; convert to "hours since"
 ntim = dimsizes(Time) ; # of times on file
print(ntim)
 wks = gsn_open_wks("ps" ,"ihop_site2") ; ps,pdf,x11,ncgm,eps
   plot = new(1,graphic) ; create a plot array
  res1 = True ; plot mods desired
  res1@gsnDraw = False ; don't draw
  res1@gsnFrame = False ; don't advance frame

   res1@lbLabelFont = "times-roman"
   res1@tmXBLabelFont = "Times-Roman"
   res1@tiYAxisFont = "Times-Roman"
   res1@tiXAxisFont = "Times-Roman"
   res1@tiMainFont = "Times-Roman"
   res1@tiYAxisFontHeightF = 0.020
   res1@tiXAxisFontHeightF = 0.020
  res1@mpShapeMode = "FreeAspect"
  res1@vpWidthF = 0.7
  res1@vpHeightF = 0.4

; res1@tiMainString ="Walnut River(KS)" ; "+ lat+"N "+fabs(lon)+"E"
  res1@tiXAxisString = "UTC (2-5 June 2002)" ;Time@units ; label bottom axis with units
  res1@tiYAxisString = "Temp at 2m (K)" ; axis string
  res1@xyDashPattern = 1 ; Make curves all solid
  res1@xyLineThicknesses = 1
  res1@xyMarkLineMode = "MarkLines" ; Markers *and* lines
  res1@xyMarkers = 16 ; 3 different markers
  res1@xyMarkerSizeF = 0.015
  res1@trYMaxF = 310 ; axis max
  res1@trYMinF = 285
  res1@tmYROn = False ; Don't draw labels for
  res1@tmYRLabelsOn = False ; right Y axes, b/c left
  res1@tmXTOn = False
  res1@tmXTLabelsOn = False
res1@tmXBLabelFontHeightF = 0.020
res1@tmYLLabelFontHeightF = 0.020
  res1@tmYLMode = "Manual"
  res1@tmYLTickStartF = 285
    res1@tmYLTickEndF = 310
  res1@tmYLTickSpacingF = 4
  res1@tmXBMode = "Explicit" ; Define your own tick mark labels.
  res1@tmXBMinorOn = False ; No minor tick marks.
  res1@tmXBValues = ispan(0,67,1); Values from 0 to 12.
 ; res1@tmXBLabels = (/"00UTC","12UTC","00UTC","12UTC","00UTC","12UTC"/)

;res1@gsnMaximize = True ; uncomment to maximize size

  res2 = True
  res2@gsnDraw = False ; don't draw
    res2@gsnFrame = False ; don't advance frame
; res2@gsnMaximize = True ; uncomment to maximize size
  res2@xyDashPattern = 1 ; Make curves all solid
  res2@xyLineThicknesses = 1
  res2@xyMarkLineMode = "MarkLines" ; Markers *and* lines
  res2@xyMarkers = 4 ; 3 different markers
   res2@xyMarkerSizeF = 0.015
  res2@trYMaxF = 310 ; axis max
  res2@trYMinF = 285
  res2@tmYROn = False ; Don't draw labels for
  res2@tmYRLabelsOn = False ; right Y axes, b/c left
  res2@tmXTOn = False
  res2@tmXTLabelsOn = False

   res3 = True
   res3@gsnDraw = False ; don't draw
     res3@gsnFrame = False ; don't advance frame
; res3@gsnMaximize = True ; uncomment to maximize size
   res3@xyDashPattern = 0 ; Make curves all solid
   res3@xyLineThicknesses = 1
; res3@xyMarkLineMode = "MarkLines" ; Markers *and* lines
; res3@xyMarkers = 4 ; 3 different markers
   res3@xyMarkerSizeF = 0.015
   res3@trYMaxF = 310 ; axis max
   res3@trYMinF = 285
   res3@tmYROn = False ; Don't draw labels for
   res3@tmYRLabelsOn = False ; right Y axes, b/c left
   res3@tmXTOn = False
   res3@tmXTLabelsOn = False

plot(0) = gsn_csm_xy3(wks,Time,T,T1A,T13,res1,res2,res3)

end

                                               
_________________________________________________________________
Windows Live Hotmail: Your friends can get your Facebook updates, right from Hotmail®.
http://www.microsoft.com/middleeast/windows/windowslive/see-it-in-action/social-network-basics.aspx?ocid=PID23461::T:WLMTAGL:ON:WL:en-xm:SI_SB_4:092009

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk