Re: zonal wind anomaly estimation

From: juki juki <juky_emc2_at_nyahnyahspammersnyahnyah>
Date: Tue Mar 06 2012 - 16:08:07 MST

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_at_ucar.edu> To: juki juki <juky_emc2_at_yahoo.com> Cc: "ncl-talk_at_ucar.edu" <ncl-talk_at_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_at_ucar.edu> > To: juki juki <juky_emc2_at_yahoo.com> > Cc: "ncl-talk_at_ucar.edu" <ncl-talk_at_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
Received on Tue Mar 6 13:08:18 2012

This archive was generated by hypermail 2.1.8 : Tue Mar 13 2012 - 14:00:14 MDT