Re: bus error?

Hi Erik,

A "bus error" could be any number of things. Do you know if this code worked under 5.1.1?

Can you help me see which line the bus error occurs on by runnining your script with the "-x" option:

ncl -x /Users/enoble/myScripts/ncl/scripts_in_use/wrf_trmm_hovmoller.ncl

Send me the last few lines of the echoed output that you should see.

Thanks,

--Mary

On May 4, 2010, at 3:51 PM, Erik Noble wrote:

> 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

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Tue May 4 16:52:48 2010