conform?

From: yanling_1220 <yanling_1220_at_nyahnyahspammersnyahnyah>
Date: Mon, 31 Aug 2009 10:33:58 +0800 (CST)

Hi allúČ
 
I write a simple script(atteched) with "clmMonTLL"and "conform" function to calculate the grid difference with a long term monthly means.Though, I obatin an error following that I cannot solve even if conceptually very simple.Can any body help please?
 
Thanks in advance!
 
Regards,
sunflower
 
--------------------------------------------------------------------------------
fatal:conform: The array to be conformed must have the same number of dimensions
 as indicated by the length of the last argument
--------------------------------------------------------------------------------

load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"

begin
  in = addfile("mpi_echam5.20c3m.average_runs.monthly.1961-1-2000-12_tas_celsius_linint_1.0.nc","r")
  A = in->tas (0:359,:,:) ;(360,181, 361)
  printVarSummary(A)
  tClm = clmMonTLL( A ) ;(12,181,361)
  printVarSummary(tClm)
  
  in2= addfile("mpi_echam5_runs_average_sresa2_2001_2030_celsius_linint_1.0.nc","r")
  B = in2->tas(0:359,:,:) ;(360,181,361)
  printVarSummary(B)
  
  tconform = conform(B,tClm,(/1,2/)) ;error lies
  printVarSummary(tconform )
  
  diff = B - tconform
  printVarSummary( diff)
  
  diff_at_long_name = " grid Difference between 2001_2030 with 1961_1990_term_monthly_means " ; meta data
  diff_at_units = B_at_units
  copy_VarCoords(B, diff)
  printVarSummary(diff)

  delete(in)
  delete(A)
  delete(tClm)
  delete(in2)
  delete(B)
  delete(diff)
end
 ----------------------------------------------------------------------------
     
      

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Sun Aug 30 2009 - 20:33:58 MDT

This archive was generated by hypermail 2.2.0 : Thu Sep 03 2009 - 07:55:08 MDT