bus error?

Hi.
I am getting a "bus error.'

What are these errors related to? I am using NCL 5.2.0

progo:~ enoble$ ncl
/Users/enoble/myScripts/ncl/scripts_in_use/wrf_trmm_hovmoller.ncl
 Copyright (C) 1995-2010 - All Rights Reserved
 University Corporation for Atmospheric Research
 NCAR Command Language Version 5.2.0
 The use of this software is governed by a License Agreement.
 See http://www.ncl.ucar.edu/ for more details.

Variable: WRF_rain
Type: float
Total Size: 677600 bytes
            169400 values
Number of Dimensions: 3
Dimensions and sizes: [Time | 11] x [south_north | 110] x [west_east | 140]
Coordinates:
            south_north: [-17.51322..34.56909]
            west_east: [-34.86005..34.86005]
Number Of Attributes: 10
  Time_Label : <ARRAY of 49 elements>
  Time_string : 2006-09-02_00:00:00
  Time : <ARRAY of 49 elements>
  FieldType : 104
  MemoryOrder : XY
  description : Daily Accumulated Precipitation
  units : mm
  stagger :
  coordinates : XLONG XLAT
  _FillValue : -999

Variable: rain_trmm
Type: float
Total Size: 31276800 bytes
            7819200 values
Number of Dimensions: 3
Dimensions and sizes: [time | 30] x [latGRID | 362] x [lonGRID | 720]
Coordinates:
            latGRID: [-90..90]
            lonGRID: [-180.2504..179.2504]
Number Of Attributes: 6
  hdf_name : precipitation
  lonFlip : longitude coordinate variable has been reordered via lonFlip
  units : mm
  long_name : 24-hr accumulated precipitation
  sum_op_ncl : dim_sum_n over dimension(s): time
  _FillValue : -9999

Variable: rain_TRMM
Type: float
Total Size: 31276800 bytes
            7819200 values
Number of Dimensions: 3
Dimensions and sizes: [time | 30] x [latGRID | 362] x [lonGRID | 720]
Coordinates:
            latGRID: [-90..90]
            lonGRID: [-180.2504..179.2504]
Number Of Attributes: 6
  _FillValue : -9999
  sum_op_ncl : dim_sum_n over dimension(s): time
  long_name : 24-hr accumulated precipitation
  units : mm
  lonFlip : longitude coordinate variable has been reordered via lonFlip
  hdf_name : precipitation

Variable: trmm
Type: float
Total Size: 88000 bytes
            22000 values
Number of Dimensions: 3
Dimensions and sizes: [time | 11] x [latGRID | 20] x [lonGRID | 100]
Coordinates:
            latGRID: [5.235458..14.70914]
            lonGRID: [-19.97357..29.95827]
Number Of Attributes: 6
  hdf_name : precipitation
  lonFlip : longitude coordinate variable has been reordered via lonFlip
  units : mm
  long_name : 24-hr accumulated precipitation
  sum_op_ncl : dim_sum_n over dimension(s): time
  _FillValue : -9999

Variable: WRF_RAIN
Type: float
Total Size: 88000 bytes
            22000 values
Number of Dimensions: 3
Dimensions and sizes: [Time | 11] x [south_north | 20] x [west_east | 100]
Coordinates:
            south_north: [5.272903..14.64821]
            west_east: [-19.81254..29.84421]
Number Of Attributes: 10
  _FillValue : -999
  coordinates : XLONG XLAT
  stagger :
  units : mm
  description : Daily Accumulated Precipitation
  MemoryOrder : XY
  FieldType : 104
  Time : <ARRAY of 49 elements>
  Time_string : 2006-09-02_00:00:00
  Time_Label : <ARRAY of 49 elements>
Bus error
progo:~ enoble$

;;;;;;;;;;;;;;;;;;; Code;;;;;;;;;;;;;;;;;;;
;***************************************************************
; Load Libaries
;***************************************************************
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"
load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRF_contributed.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW_add.ncl"
;***************************************************************
; Change these parts
;***************************************************************
version = "Default"
forecast_days = "12"
dates = "9-2-2006 0z to 9-14-2006 0z"
diri = "/Users/enoble/myWorkDir/Current_work/12-Day_forecast_Default_folder/"
fili = "WRF-"+version+"_"+forecast_days+"day-forecast_rain.nc"
dirA = "/Users/enoble/myWorkDir/precipitation-wrf_work/TRMM_programs/All_daily_files/"
filA = "TRMM_Daily_PrecipTotals_05x05-GRID.nc"
;***************************************************************
; Format of ouput file
;***************************************************************
type = "x11"
;type = "pdf"
;type = "eps"
;type = "ps"
;***************************************************************
; Read in Files
;***************************************************************
; WRF
model = addfile(diri+fili,"r")
;***** wrf precipitation ********
WRF_rain = model->rain
printVarSummary(WRF_rain)
; TRMM
obs = addfile(dirA+filA,"r")
;***** TRMM precipitation ********
rain_trmm = obs->daily_accum
rain_TRMM = rain_trmm ; Check Times that you need (ALWAYS 1 DAY AHEAD)
;rain_TRMM = where(rain_TRMM.lt.0.0 , 0.0, rain_TRMM)
printVarSummary(rain_trmm)
printVarSummary(rain_TRMM)
;***************************************************************
; Adjust WRF to fit to TRMM data ; Fit each other
;***************************************************************
trmm = rain_TRMM(1:11,{5:15},{-20:30})
WRF_RAIN = WRF_rain(:,{5:15},{-20:30})
printVarSummary(trmm)
printVarSummary(WRF_RAIN)
;***************************************************************
; Hov
;***************************************************************
trmmHOV= dim_avg_n_Wrap(trmm,1)
wrfHOV = dim_avg_n_Wrap(WRF_RAIN,1)
;************************************************
; create plot
;************************************************
type = "ps" ;
create variable to hold info
type@wkColorModel = "cmyk" ; make
output cmyk
wks = gsn_open_wks(type,"TRMM_rainHov_LenBox") ; open a ps file
gsn_define_colormap(wks,"BlAqGrYeOrRe") ; choose colormap

res = True
 ; plot mods desired
;******************* Use These next set of commands if using color
res@cnFillOn = True ; turn on color fill
res@gsnSpreadColors = True ; use full range of colors
res@cnLinesOn = False ; turn off contour lines
res@cnLineLabelsOn = False ; turn off line labels
; CMYK colors
;res@gsnSpreadColorStart = 10 ; start at color 10
;res@gsnSpreadColorEnd = 65
;*******************
res@cnLevelSelectionMode = "ManualLevels" ; set manual contour levels
res@cnMinLevelValF = 0.
res@cnMaxLevelValF = 60.
res@cnLevelSpacingF = 5.

res@tiMainString = "Daily Accumulated Precipitation"
res@tiXAxisString = "Longitude"
res@gsnMajorLonSpacing = 5.
;res@tmXBTickSpacingF = 5
res@tiYAxisString = "Elapsed Time: SEP 2 to 14, 2006"
res@gsnLeftString = "Lats: 5N to 15N"
res@gsnCenterString = "TRMM 3B42"
res@gsnRightString = "Units: mm"
res@pmTickMarkDisplayMode = "Always"
plot = gsn_csm_hov(wks,trmmHOV,res)

;res@tiYAxisString = "Elapsed Time: SEP 2 to 14, 2006"
;res@gsnLeftString = "Lats: 5N to 15N"
;res@gsnCenterString = "WRF-ARW"
;res@gsnRightString = "Units: mm"

;plot=gsn_csm_hov(wks,wrfHOV,res)
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Tue May 4 15:51:39 2010