Re: zonal wind anomaly estimation

Hi Juki,

This helps tremendously, thanks! I was able to run your script and I think I see the problem. I'll write back if we come up with a work-around.

--Mary

On Mar 6, 2012, at 4:08 PM, juki juki wrote:

> Hi Mary,
>
> I have sent the data and ncl code to the mentioned ftp, hope you received it. Thank you so much.
>
> Regards,
> Juki
>
> From: Mary Haley <haley@ucar.edu>
> To: juki juki <juky_emc2@yahoo.com>
> Cc: "ncl-talk@ucar.edu" <ncl-talk@ucar.edu>
> Sent: Wednesday, March 7, 2012 7:23 AM
> Subject: Re: zonal wind anomaly estimation
>
> There's still not enough information, because the script you gave me and the line number you are reporting below don't match.
>
> Can you provide your data file so I can run the script?
>
> See the ftp information at:
>
> http://www.ncl.ucar.edu/report_bug.shtml
>
> --Mary
>
>
> On Mar 6, 2012, at 3:02 PM, juki juki wrote:
>
> > Dear Mary;
> >
> > Thanks for help.
> >
> > Besides the error of "the dimension mismatch on subscript #1", another error is:
> >
> > "near line 3826 in file contributed.ncl"
> > "near line 46" of my code (in revised code structure), that is, climo =clmDayTTL(data,yyyyddd)
> >
> > I have tried to do this:
> > printVarSummary(fin->uwnd) or printVarSummary(data), it seems O.K,
> >
> > printVarSummary(fin->uwnd
> >
> > Copyright (C) 1995-2011 - All Rights Reserved
> > University Corporation for Atmospheric Research
> > NCAR Command Language Version 6.0.0
> > The use of this software is governed by a License Agreement.
> > See http://www.ncl.ucar.edu/ for more details.
> > (0) Here we go! Wed Mar 7 06:58:13 TST 2012
> > (0) Calculate climatology... Wed Mar 7 06:58:13 TST 2012
> > (0) Reading the data...
> >
> > Variable: uwnd (file variable)
> > Type: short
> > Total Size: 130453920 bytes
> > 65226960 values
> > Number of Dimensions: 4
> > Dimensions and sizes: [time | 365] x [level | 17] x [lat | 73] x [lon | 144]
> > Coordinates:
> > time: [1770696..1779432]
> > level: [1000..10]
> > lat: [90..-90]
> > lon: [ 0..357.5]
> > Number Of Attributes: 20
> > long_name : Daily U-wind on Pressure Levels
> > valid_range : ( -32765, -1266 )
> > unpacked_valid_range : ( -140, 175 )
> > actual_range : ( -79.76001, 134.22 )
> > units : m/s
> > add_offset : 187.65
> > scale_factor : 0.01
> > missing_value : 32766
> > _FillValue : -32767
> > precision : 2
> > least_significant_digit : 1
> > GRIB_id : 33
> > GRIB_name : UGRD
> > var_desc : u-wind
> > dataset : NCEP/DOE AMIP-II Reanalysis (Reanalysis-2) Daily Averages
> > level_desc : Pressure Levels
> > statistic : Mean
> > parent_stat : Individual Obs
> > standard_name : eastward_wind
> > cell_methods : time: mean (of 4 6-hourly values in one day)
> > fatal:Dimension size mismatch on subscript #1, left-hand and right-hand side dimensions do not match
> >
> > fatal:Execute: Error occurred at or near line 3826 in file $NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl
> >
> >
> > fatal:Execute: Error occurred at or near line 46 in file annual_cycle_uwind.ncl
> >
> > printVarSummary(data), some output:
> >
> > Copyright (C) 1995-2011 - All Rights Reserved
> > University Corporation for Atmospheric Research
> > NCAR Command Language Version 6.0.0
> > The use of this software is governed by a License Agreement.
> > See http://www.ncl.ucar.edu/ for more details.
> > (0) Here we go! Wed Mar 7 06:58:33 TST 2012
> > (0) Calculate climatology... Wed Mar 7 06:58:33 TST 2012
> > (0) Reading the data...
> >
> > Variable: data
> > Type: float
> > Total Size: 433620 bytes
> > 108405 values
> > Number of Dimensions: 3
> > Dimensions and sizes: [time | 365] x [lat | 9] x [lon | 33]
> > Coordinates:
> > time: [1770696..1779432]
> > lat: [-10..10]
> > lon: [80..160]
> > Number Of Attributes: 21
> > level : 200
> > long_name : Daily U-wind on Pressure Levels
> > unpacked_valid_range : ( -140, 175 )
> > actual_range : ( -79.76001, 134.22 )
> > units : m/s
> > missing_value_original : 32766
> > _FillValue_original : -32767
> > _FillValue : -32767
> > precision : 2
> > least_significant_digit : 1
> > GRIB_id : 33
> > GRIB_name : UGRD
> > var_desc : u-wind
> > dataset : NCEP/DOE AMIP-II Reanalysis (Reanalysis-2) Daily Averages
> > level_desc : Pressure Levels
> > statistic : Mean
> > parent_stat : Individual Obs
> > standard_name : eastward_wind
> > cell_methods : time: mean (of 4 6-hourly values in one day)
> > valid_range : ( -140, 174.99 )
> > missing_value : -32767
> > fatal:Dimension size mismatch on subscript #1, left-hand and right-hand side dimensions do not match
> >
> > fatal:Execute: Error occurred at or near line 3826 in file $NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl
> >
> >
> > fatal:Execute: Error occurred at or near line 46 in file annual_cycle_uwind.ncl
> >
> > Thanks again,
> >
> > JukY
> >
> >
> >
> >
> >
> >
> >
> >
> >
> > From: Mary Haley <haley@ucar.edu>
> > To: juki juki <juky_emc2@yahoo.com>
> > Cc: "ncl-talk@ucar.edu" <ncl-talk@ucar.edu>
> > Sent: Wednesday, March 7, 2012 12:23 AM
> > Subject: Re: zonal wind anomaly estimation
> >
> > Dear Juky,
> >
> > You need to provide more information, since we can't run your script. Is there more to the error message? Does the error message provide a line number?
> >
> > My guess is that the error is with this line:
> >
> > dataShort = fin->uwnd(:,{200},{latS:latN},{lonW:lonE})
> >
> > Make sure you are subscripting this properly by first doing a printVarSummary on fin->uwnd:
> >
> > printVarSummary(fin->uwnd)
> >
> >
> > --Mary
> >
> > On Mar 5, 2012, at 4:55 PM, juki juki wrote:
> >
> > > Dear All;
> > >
> > > I want to calculate the zonal wind anomaly from daily wind data of NCAR. The code is attached. I got an error "the dimension mismatch on subscript #1". Probably you could help me to look at the error and propose possible way to improve it, I really appreciate it. Thanks for help.
> > >
> > > Regards,
> > > JuKy
> > > <annual_cycle_uwind.ncl>_______________________________________________
> > > 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
>
>
>
> _______________________________________________
> 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 Wed Mar 7 05:50:41 2012