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