Re: argument mismatch issue

Mary,

In fact, it is the kf script which requires the reodering to have time, lat
and lon. The program works fine when i use olr
 short olr(time, lat, lon)

when i use other variable as below

double pr(time, lat, lon)
 I am having the error.

Here are a few the lines of the code

CS

  ; open the file
  fileName = "clm_pr_9600.nc"
   inFile = addfile( fileName, "r" )
  obsPerDay = 1

  ; find our times
  time = inFile->time
 ; lon = inFile->lon
 ; lat =inFile->lat
  pr = inFile ->pr
  filtStart = ut_inv_calendar( 1996, 01, 01, 0, 0, 0, time@units, 0)
  filtEnd = max(time)
  ;startDate = filtStart + 192 * 24
  ;endDate = filtEnd - 192 * 24
  startDate = filtStart
  endDate = filtEnd
  ; read in the unfiltered data
 print( "Reading the data..." )
   dataShort = inFile->pr({filtStart:filtEnd},:,:)
   printVarSummary(pr)
  data = doubletofloat( dataShort )
 printVarSummary(data)
 es= data

   es(time|:,lat|y,lon|:) = (/ kf_filter( data(time|:,lat|y,lon|:), \\
         obsPerDay, 2, 10, -20, -6, 90, mis, "mrg" ) /)

******************************************************************************

2011/10/17 Mary Haley <haley@ucar.edu>

> CS,
>
> The error is telling you what the problem is. You are trying to reorder
> "data" by using named dimensions "time", "lat", "lon",
> but as you can see from the printVarSummary command, "data" doesn't have
> any named dimensions.
>
> The printVarSummary of "data" makes it look like it is already ordered time
> x lat x lon (I'm assuming this because
> the first dimension is rather large, so that's likely a time dimension, and
> the second dimension is about half the
> size of the third dimension, so this looks like lat x lon).
>
> If "data" is already time x lat x lon, then don't use dimension
> reordering:
>
> es(time|:,lat|y,lon|:) = (/ kf_filter( data,obsPerDay, 2, 10, -20, -6,
> 90, mis, "mrg" ) /)
>
> I don't know why you have the reordering on the left side as well. Are you
> sure you need to do this?
>
> --Mary
>
> On Oct 16, 2011, at 10:18 AM, Coly SAR wrote:
>
> I had a progress but having an error in variable data with no specified
> time , lat ,lon
> ............................
> Variable: data
> Type: float
> Total Size: 75221244 bytes
> 18805311 values
> Number of Dimensions: 3
> Dimensions and sizes: [1827] x [73] x [141]
> Coordinates:
> Number Of Attributes: 1
> _FillValue : -9.99e+08
> (0) Filtering...dim oct 16 16:13:54 GMT 2011
> fatal:(lat) is not a named dimension in variable (data).
> fatal:Execute: Error occurred at or near line 83
>
> CS
>
> 2011/10/16 Dennis Shea <shea@ucar.edu>
>
>> This type of error can be diagnosed by doing a printVarSummary(...)
>> of kf_filter and es. Or, print(typeof(kf_filter)) and print(typeof(es))
>>
>> Perhaps, kf_filter is type double and es is type float? If so,
>>
>> es(time|:,lat|y,lon|:) = tofloat( \
>>
>> (/ kf_filter( data(time|:,lat|y,lon|:), \
>> obsPerDay, 2, 10, -20, -6, 90, mis, "mrg" ) /) )
>>
>>
>> On 10/16/11 9:08 AM, Coly SAR wrote:
>>
>>> Dear all,
>>>
>>> I am having a mismatch argument error when running a script with my data
>>> as shown below:
>>>
>>> Number of Dimensions: 3
>>> Dimensions and sizes: [time | 1827] x [lat | 73] x [lon | 141]
>>> Coordinates:
>>> time: [ 0..2629440]
>>> lat: [-0.5..35.5]
>>> lon: [-22.5..47.5]
>>> Number Of Attributes: 1
>>> _FillValue : -999000000
>>> (0) 0 72 dim oct 16 14:55:40 GMT 2011
>>> fatal:Argument type mismatch on argument (0) of (kf_filter) can not
>>> coerce
>>> fatal:Execute: Error occurred at or near line 97
>>>
>>> Hoping to have a hint to solve this urgent problem.
>>>
>>> The line of the code is below:
>>>
>>> es(time|:,lat|y,lon|:) = (/ kf_filter( data(time|:,lat|y,lon|:), \\
>>> obsPerDay, 2, 10, -20, -6, 90, mis, "mrg" ) /)
>>>
>>> Thanks,
>>>
>>> CS
>>>
>>>
>>>
>>> ______________________________**_________________
>>> ncl-talk mailing list
>>> List instructions, subscriber options, unsubscribe:
>>> http://mailman.ucar.edu/**mailman/listinfo/ncl-talk<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
>
>
>

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Mon Oct 17 11:41:38 2011