override dimension reduction

From: Jeanne Thibeault <jeanne.thibeault_at_nyahnyahspammersnyahnyah>
Date: Mon, 21 Sep 2009 15:56:20 -0400

Hello NCL,

Is there an easy way to prevent dimension reduction when calculating
the average of a 3D variable? Or is there a way to convert the
resulting 2D variable back to a 3D variable? A portion of my script
is below.

  B1cdd8099CNRM = cddCNRM_19012099B1(79:98,:,:)
  copy_VarAtts(cdd20CNRM,B1cdd8099CNRM)
  printVarSummary(B1cdd8099CNRM)

nlatCNRM = 64
nlonCNRM = 128
ntime = 1
  meanB1cdd8099CNRM = dim_avg_Wrap(B1cdd8099CNRM(lat|:, lon|:, time|:))
  printVarSummary(meanB1cdd8099CNRM)
  B1CDD8099CNRMmean = meanB1cdd8099CNRM(0:0,:,:) ******this line is
producing an error
B1CDD8099CNRMmean = new((/ntime,nlatCNRM, nlonCNRM/), typeof
(meanB1cdd8099CNRM), meanB1cdd8099CNRM@_FillValue)

fatal:Number of subscripts on rhs do not match number of dimensions
of variable,
(3) Subscripts used, (2) Subscripts expected

Thanks for your help,
Jeanne

Variable: B1cdd8099CNRM
Type: float
Total Size: 655360 bytes
             163840 values
Number of Dimensions: 3
Dimensions and sizes: [time | 20] x [lat | 64] x [lon | 128]
Coordinates:
             time: [44012..50951]
             lat: [-87.86380004882812..87.86380004882812]
             lon: [ 0..357.1875]
Number Of Attributes: 5
   original_name : pxcdd
   missing_value : 1e+20
   _FillValue : 1e+20
   units : days
   long_name : Maximum Number of Consecutive Dry Days

Variable: meanB1cdd8099CNRM
Type: float
Total Size: 32768 bytes
             8192 values
Number of Dimensions: 2
Dimensions and sizes: [lat | 64] x [lon | 128]
Coordinates:
             lat: [-87.86380004882812..87.86380004882812]
             lon: [ 0..357.1875]
Number Of Attributes: 6
   long_name : Maximum Number of Consecutive Dry Days
   units : days
   _FillValue : 1e+20
   missing_value : 1e+20
   original_name : pxcdd
   average_op_ncl : dim_avg over dimension: time

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Mon Sep 21 2009 - 13:56:20 MDT

This archive was generated by hypermail 2.2.0 : Thu Sep 24 2009 - 13:55:09 MDT