Hi all,
I have a WRF output, and I want to craet an ascii output of it, acording to a list of lat&lon.
when run my script, get this error:
fatal:Dimension (west_east) of (dusttotal) does not have an associated coordinate variable
fatal:Execute: Error occurred at or near line 35 in file pm10.ncl
I could not resolve this error, can anyone help me please?!
thanks in advance
Marzie
here is my script:
//*********************************
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"
begin
a = addfile("wrfout_d01_2010-05-18_00_00_00.nc","r")
times = wrf_user_list_times(a) ; get times in the file
ntimes = dimsizes(times)
do it = 0,ntimes-1,1 ; TIME LOOP
DUST_1 = wrf_user_getvar(a,"DUST_1",it)
DUST_2 = wrf_user_getvar(a,"DUST_2",it)
DUST_3 = wrf_user_getvar(a,"DUST_3",it)
DUST_4 = wrf_user_getvar(a,"DUST_4",it)
DUST_5 = wrf_user_getvar(a,"DUST_5",it)
dusttot=DUST_1
dusttot=((DUST_1)+(DUST_2)+(DUST_3)+(DUST_4)+(DUST_5))
dusttotal=dusttot(0,:,:)
printVarSummary(DUST_1)
printVarSummary(dusttotal)
data= asciiread("lat&lon-M.csv",(/20/),"float")
output = new((/10/),float)
output!0 = "Dust"
do i = 0,19,2
output(i/2) = dusttotal({data(i)},{data(i+1)})
end do
printVarSummary(output)
asciiwrite("dust.txt",output)
end do
end
//******************************************
and here is summary of variables:
----------------------------------------------
Variable: times
Type: string
Total Size: 4 bytes
1 values
Number of Dimensions: 1
Dimensions and sizes: [1]
Coordinates:
Number Of Attributes: 2
description : times in file
_FillValue : missing
(0) 2010-05-18_00:00:00
------------------------------------------------------
Variable: DUST_1
Type: float
Total Size: 1152000 bytes
288000 values
Number of Dimensions: 3
Dimensions and sizes: [bottom_top | 40] x [south_north | 60] x [west_east | 120]
Coordinates:
Number Of Attributes: 6
coordinates : XLONG XLAT
stagger :
units : ug/kg-dryair
description : dust size bin 1: 0.5um effective radius
MemoryOrder : XYZ
FieldType : 104
-----------------------------------------------------
Variable: dusttotal
Type: float
Total Size: 28800 bytes
7200 values
Number of Dimensions: 2
Dimensions and sizes: [south_north | 60] x [west_east | 120]
Coordinates:
Number Of Attributes: 6
coordinates : XLONG XLAT
stagger :
units : ug/kg-dryair
description : dust size bin 1: 0.5um effective radius
MemoryOrder : XYZ
FieldType : 104
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Sat Jul 7 14:04:34 2012
This archive was generated by hypermail 2.1.8 : Mon Jul 09 2012 - 10:45:32 MDT