Re: EOF - gaussian grid

Dennis,

The script ran, but did not generate the figure.

The problem is: "fatal:(lat) is not a dimension name in variable (clat),
could not determine dimension number
                        fatal:["Execute.c":7743]:Execute: Error occurred at
or near line 113 in file eof_NAO.ncl"

- line 113: eof_ts = eof_ts/sumWgt [where sumWgt = mln*sum(
clat({lat|latS:latN}) )]

_____________________________________________________________________________________________
This occur because clat do not have the characteristics of lat.

printVarSummary(lat):
   Variable: lat
   Type: double
   Total Size: 768 bytes
               96 values
   Number of Dimensions: 1
   Dimensions and sizes: [lat | 96]
   Coordinates:
               lat: [-88.57217..88.57217]
   Number Of Attributes: 4
     standard_name : latitude
     long_name : latitude
     units : degrees_north
     axis : Y

printVarSummary(clat) - where clat = tofloat(f->lat)
   Variable: clat
   Type: float
   Total Size: 384 bytes
               96 values
   Number of Dimensions: 1
   Dimensions and sizes: [96]
   Coordinates:
_____________________________________________________________________________________________

How I can repair this?

Renata

2013/4/17 Renata Tedeschi <rgtedeschi@gmail.com>

> Thanks for your help Dennis, now the script ran correctly.
>
> Renata
>
>
> 2013/4/17 Dennis Shea <shea@ucar.edu>
>
>> I speculate that the clat (derived fron 'lat') is type double.
>>
>> printVarSummary(lat) ; note the 'type'
>> printVarSummary(clat)
>>
>>
>> wSLP = SLP ; type float
>>
>> wSLP = SLP*conform(SLP, clat, 1)
>>
>> In the above 'conform(SLP, clat, 1)' creates an array that is
>> the same size and shape as SLP (wSLP) but is populated with
>> entries of clat which is type double. The product of a
>> float*double is a double.
>>
>> NCL is a strongly typed language it will not let you lose information.
>> Hence, the error message
>>
>>
>> "Assignment type mismatch, right hand side can't
>> be coerced to type of left hand side"
>>
>> refers to the fact you are trying to force a double to a float.
>>
>> The easiest way to address this is to concert the lat the float
>>
>> clat = tofloat(f->lat) ; convert to float
>>
>> clat = sqrt( cos(rad*clat) )
>>
>> Good luck
>>
>>
>> On 4/17/13 6:26 AM, Renata Tedeschi wrote:
>>
>>> When I use the eof_1.ncl script the following message appears:
>>>
>>> "fatal:["NclVar.c":1382]:**Assignment type mismatch, right hand side
>>> can't
>>> be coerced to type of left hand side
>>> fatal:["Execute.c":7743]:**Execute: Error occurred at or near line
>>> 91 in
>>> file eof_NAO.ncl"
>>>
>>> Line 91: wSLP = SLP*conform(SLP, clat, 1)
>>> ______________________________**______________________________**
>>> ______________________________**___
>>> - printVarSummary(SLP) [where SLP is: SLP = month_to_season (slp,
>>> season)]:
>>> Variable: SLP
>>> Type: float
>>> Total Size: 2211840 bytes
>>> 552960 values
>>> Number of Dimensions: 3
>>> Dimensions and sizes: [time | 30] x [lat | 96] x [lon | 192]
>>> Coordinates:
>>> time: [ 0..254232]
>>> lat: [-88.57217..88.57217]
>>> lon: [-180..178.125]
>>> Number Of Attributes: 6
>>> lonFlip : longitude coordinate variable has been reordered via
>>> lonFlip
>>> grid_type : gaussian
>>> long_name : DJF: ** msl SEA LEVEL PRESSURE (HPA)
>>> lev : 1000
>>> _FillValue : 9.999e+20
>>> NMO : 0
>>>
>>> - printVarSummary(wSLP) [where wSLP=SLP]:
>>> Variable: wSLP
>>> Type: float
>>> Total Size: 2211840 bytes
>>> 552960 values
>>> Number of Dimensions: 3
>>> Dimensions and sizes: [time | 30] x [lat | 96] x [lon | 192]
>>> Coordinates:
>>> time: [ 0..254232]
>>> lat: [-88.57217..88.57217]
>>> lon: [-180..178.125]
>>> Number Of Attributes: 6
>>> NMO : 0
>>> _FillValue : 9.999e+20
>>> lev : 1000
>>> long_name : DJF: ** msl SEA LEVEL PRESSURE (HPA)
>>> grid_type : gaussian
>>> lonFlip : longitude coordinate variable has been reordered via
>>> lonFlip
>>> ______________________________**______________________________**
>>> ______________________________**_______
>>>
>>> Sorry for inconvenience and thanks for your help
>>> Renata
>>>
>>>
>>>
>>> 2013/4/16 Dennis Shea <shea@ucar.edu>
>>>
>>> You never said what your problem is.
>>>> It certainly is *not* related using 'clat' rather than 'gw'
>>>>
>>>> ==
>>>> Did you get an error message?
>>>> What does the following show?
>>>>
>>>> printVarSummary(SLP)
>>>>
>>>> -
>>>> FYI: You can use the NCL function 'gaus' to generate
>>>> the gaussian weights, if desired.
>>>>
>>>> On 4/16/13 11:05 AM, Renata Tedeschi wrote:
>>>>
>>>>> Hi everyone.
>>>>>
>>>>> I'm using the eof_1.ncl script, but I have one problem. The problem
>>>>> occur
>>>>> in these commands:
>>>>>
>>>>> ;=============================**=============================
>>>>> ; create weights: sqrt(cos(lat)) [or sqrt(gw) ]
>>>>> ;=============================**=============================
>>>>> rad = 4.*atan(1.)/180.
>>>>> clat = f->lat
>>>>> clat = sqrt( cos(rad*clat) ) ; gw for gaussian
>>>>> grid
>>>>>
>>>>> ;=============================**=============================
>>>>> ; weight all observations
>>>>> ;=============================**=============================
>>>>> wSLP = SLP ; copy meta data
>>>>> wSLP = SLP*conform(SLP, clat, 1)
>>>>>
>>>>> I can believe that the problem is because my grid is gaussian, and I do
>>>>>
>>>> not
>>>>
>>>>> use gw, but I do not know what is gw.
>>>>>
>>>>> Renata Tedeschi
>>>>>
>>>>>
>>>>>
>>>>> ______________________________**_________________
>>>>> 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<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 Wed Apr 17 07:49:57 2013