Re: month_to_annual function error

 Hi all,

 I have the same problem with month_to_annual function as well,
 I tried with small test from Marry, and it worked properly. But, with my
 data, it did not work and the message is :
 fatal:mask: dimension sizes of parameter 0 and parameter 1 do not match
 fatal:Execute: Error occurred at or near line 8573 in file
 $NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl

 fatal:Execute: Error occurred at or near line 40 in file
 monthtoyearlySODA.ncl

 Here is my code :load
 "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
 begin
 ;**********************************************
 ; file handling
 ;**********************************************i

 ; test

   var_raw = random_uniform(0.,2.,(/360,146,196/))
   var_raw!0 = "time"
   var_raw!1 = "lat"
   var_raw!2 = "lon"
   printVarSummary(var_raw)
   var = month_to_annual(var_raw, 1)
   printVarSummary(var)

 ; data
 a = addfile("/data03/siswanto/dataocn/TEMP_SODA_surface.nc","r") ; monthly
 data

 temp = a->TEMP(0:479,:,:,:)
 printVarSummary(temp)
 temp!0="TIME1"
 temp!1 = "DEPTH1_1"
 temp!2 = "LAT"
 temp!3 = "LON"
 temp@_FillValue=-999.

 ; annual analysis
    ;Tem = month_to_annual_weighted(196001, temp, 1)
    Tem = month_to_annual(temp, 1)
    ;Tem!0 = "year"
    ;Tem!1 = "lev"
    ;Tem!2 = "lat"
    ;Tem!3 = "lon"
    printVarSummary(Tem)
    Temp@_FillValue=-999.
    ano=dim_rmvmean_n(Tem,0) ; creating anomaly using NCL function remove
 mean flux
    annsst=dim_avg(dim_avg(ano))
   fout = addfile("TEMP_ySODA.nc","c")
   fout->TEMP = Tem

 end

 Thanks for any respons.
 Best, Siswanto

> I just tried a ver simple test.
>
>> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
>>
>> diri = "./"
>> fili = "camaqua.cam2.h1.1985-09-01-00000.nc"
>> f = addfile(diri+fili, "r")
>> x = f->LHFLX(0:23,:,:) ; read 2 years [24 months]
>> printVarSummary(x)
>>
>> xAnn = month_to_annual( x, 0)
>> printVarSummary(xAnn)
>>
>> =========
>> No problem.
>> =========
>> Variable: x
>> Type: float
>> Total Size: 786432 bytes
>> 196608 values
>> Number of Dimensions: 3
>> Dimensions and sizes: [time | 24] x [lat | 64] x [lon | 128]
>> Coordinates:
>> time: [ 0.. 23]
>> lat: [-87.86379883923262..87.86379883923262]
>> lon: [ 0..357.1875]
>> Number Of Attributes: 3
>> units : W/m2
>> long_name : Surface latent heat flux
>> cell_method : time: mean
>>
>>
>> Variable: xAnn
>> Type: float
>> Total Size: 65536 bytes
>> 16384 values
>> Number of Dimensions: 3
>> Dimensions and sizes: [year | 2] x [lat | 64] x [lon | 128]
>> Coordinates:
>> lat: [-87.86379883923262..87.86379883923262]
>> lon: [ 0..357.1875]
>> Number Of Attributes: 4
>> cell_method : time: mean
>> long_name : Surface latent heat flux
>> units : W/m2
>> NCL : month_to_annual
>>
>>
>>
>> On 6/2/10 3:23 PM, Sanjiv Kumar wrote:
>>> Hello Dennis,
>>> Below are details:
>>> printVarSummary( lht_cm20)
>>> Variable: lht_cm20
>>> Type: float
>>> Total Size: 4718592 bytes
>>> 1179648 values
>>> Number of Dimensions: 3
>>> Dimensions and sizes: [time | 144] x [lat | 64] x [lon | 128]
>>> Coordinates:
>>> time: [181071..185420]
>>> lat: [-87.8638..87.8638]
>>> lon: [ 0..357.1875]
>>> Number Of Attributes: 5
>>> long_name : runoff (mean) clm3.5
>>> units : mm/s
>>> cell_method : time: mean
>>> _FillValue : 1e+36
>>> missing_value : 1e+36
>>> and
>>> opt = 0
>>> thanks
>>> sanjiv
>>>
>>>
>>> On Wed, Jun 2, 2010 at 4:12 PM, Dennis Shea <shea@ucar.edu
>>> <mailto:shea@ucar.edu>> wrote:
>>>
>>> You must provide more information.
>>>
>>> Whenever you get an error like this, it is useful to proceed the
>>> call with
>>>
>>> printVarSummary( data )
>>>
>>> tell us what 'opt' you are using.
>>>
>>> ann = month_to_annual( data, opt)
>>>
>>>
>>> On 6/2/10 2:05 PM, Sanjiv Kumar wrote:
>>>
>>> Hello,
>>> i am trying to use month_to_annual function to calculate annual
>>> total
>>> from monthly data (var(time|:, lat|:, lon|:)), but i am getting
>>> following error message
>>> fatal:mask: dimension sizes of parameter 0 and parameter 1 do
>>> not match
>>> fatal:Execute: Error occurred at or near line 8537 in file
>>> $NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl
>>> i have no idea, why i am getting this message..any clue is
>>> highly
>>> appreciated.
>>> thanks
>>> sanjiv
>>>
>>>
>>>
>>> _______________________________________________
>>> 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
>>
>
>
> Siswanto
> Climate and Environmental Physics siswanto@climate.unibe.ch
> Physics Institute, University of Bern ph: +41 (0)31 631 48 71
> Sidlerstrasse 5, 3012 Bern, Switzerland
>
>

Siswanto
Climate and Environmental Physics siswanto@climate.unibe.ch
Physics Institute, University of Bern ph: +41 (0)31 631 48 71
Sidlerstrasse 5, 3012 Bern, Switzerland

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Fri Jun 25 07:29:41 2010