Calculating daily rainfall from 3hourly WRF output

Hi,

I have been trying to calculate daily precipitation from 3hourly WRF
history file but could not succeed. Below I am including my script
with the error message I got.

;=================================================================
; Output Frequency: 00, 03, 06, 09, 12, 15, 18, 21
; Want to generate a Rainfall time series
; In WRF rainfall variable RAINC and RAINNC are accumulated rainfall
; (rainfall at a particulat time steps is the cumulative sum of the
previous all time steps)
;**************************************************************************************************

  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/wrf/WRFUserARW.ncl"
  load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"

begin
;************************************************
; Read the file
;************************************************

     wrf_out =
addfile("/scratch/j/jcl/mkamal/DATA/Observation/DAYMET/Regrid_Rect_Curv/wrfout_d02_2002-01-01.nc","r")

     lat = wrf_out->XLAT(0,:,:)
     lon = wrf_out->XLONG(0,:,:)

     x = wrf_out->RAINC(::8,:,:) ; 8 'samples' per day

     y = wrf_out->RAINNC(::8,:,:)

     printVarSummary(x)

     printVarSummary(y)

     rainfall_tot = x+y

     printVarSummary(rainfall_tot)

;-------------------------------------------
; The following statement works fine
;-------------------------------------------

     nt = 0
     wrf_rain = rainfall_tot(nt+1,:,:)-rainfall_tot(nt,:,:)
     printVarSummary(wrf_rain)

; times = wrf_user_list_times(wrf_out)
; ntimes=dimsizes(times)

;-------------------------------------------------------------------------
; Seperating daily rainfall from accumulated rainfall (problem is here)
;-------------------------------------------------------------------------
     do nt =0,30,1
        wrf_rain = rainfall_tot(nt+1,:,:)-rainfall_tot(nt,:,:)
     end do

end

===========================================
I get the following error message:
===========================================
Variable: rainfall_tot
Type: float
Total Size: 5546520 bytes
             1386630 values
Number of Dimensions: 3
Dimensions and sizes: [31] x [210] x [213]
Coordinates:

Variable: wrf_rain
Type: float
Total Size: 178920 bytes
             44730 values
Number of Dimensions: 2
Dimensions and sizes: [210] x [213]
Coordinates:
fatal:Subscript out of range, error in subscript #0
fatal:An error occurred reading rainfall_tot
fatal:["Execute.c":7743]:Execute: Error occurred at or near line 50 in
file wrf_daily_precipitation.ncl

Where line 50 in my script is "wrf_rain =
rainfall_tot(nt+1,:,:)-rainfall_tot(nt,:,:)"
============================================================================

I look forward to hearing from you.

Thanks
Kamal

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Mon Mar 18 01:15:52 2013