Great! This is exactly what I want. It works!
Thanks alot!
Bob
Ansley Manke wrote:
Hi Bob,
Are the grids of the two input datasets exactly the same? You
could ensure that they are by regridding one to the grid of the other:yes? save/file=myTest.cdf u_velocity[d=1], v_velocity[d=2,g=u_velocity[d=1]]
Ansley Manke
Bob Yu wrote:
hi Ferreters,how I tried to combine two variables from two files into one file, I got duplicated z-aixs information:yes? set data glb_expt_01_9_3zu.nc yes? set data glb_expt_01_9_3zv.nc yes? save/file=myTest.cdf/title="my test"/clobber u_velocity[d=1], v_velocity[d=2]This is the ncdump -h myTest.cdf:netcdf myTest { dimensions: LONGITUDE = 180 ; LATITUDE = 67 ; LATITUDEedges = 68 ; Z = 33 ; Zedges = 34 ; TIME = UNLIMITED ; // (13 currently) Z1 = 33 ; <-- duplicated with Z Z1edges = 34 ; <-- duplicated with Z variables: double LONGITUDE(LONGITUDE) ; LONGITUDE:units = "degrees_east" ; LONGITUDE:modulo = 360. ; LONGITUDE:point_spacing = "even" ; LONGITUDE:AXIS = "X" ; double LATITUDE(LATITUDE) ; LATITUDE:units = "degrees_north" ; LATITUDE:point_spacing = "uneven" ; LATITUDE:AXIS = "Y" ; LATITUDE:edges = "LATITUDEedges" ; double LATITUDEedges(LATITUDEedges) ; LATITUDEedges:edges = " " ; double Z(Z) ; Z:units = "m" ; Z:positive = "down" ; Z:point_spacing = "uneven" ; Z:AXIS = "Z" ; Z:edges = "Zedges" ; double Zedges(Zedges) ; Zedges:edges = " " ; double TIME(TIME) ; TIME:units = "day since 2001-01-16 00:00:00" ; TIME:time_origin = "16-JAN-2001 00:00:00" ; TIME:calendar = "360_DAY" ; TIME:modulo = " " ; float U_VELOCITY(TIME, Z, LATITUDE, LONGITUDE) ; U_VELOCITY:missing_value = 1.267651e+30f ; U_VELOCITY:_FillValue = 1.267651e+30f ; U_VELOCITY:long_name = "Velocity U" ; U_VELOCITY:history = "From glb_expt_01_9_3zu.nc" ; U_VELOCITY:units = "cm/s" ; double Z1(Z1) ; Z1:units = "m" ; Z1:positive = "down" ; Z1:point_spacing = "uneven" ; Z1:AXIS = "Z" ; Z1:edges = "Z1edges" ; double Z1edges(Z1edges) ; Z1edges:edges = " " ; float V_VELOCITY(TIME, Z1, LATITUDE, LONGITUDE) ; V_VELOCITY:missing_value = 1.267651e+30f ; V_VELOCITY:_FillValue = 1.267651e+30f ; V_VELOCITY:long_name = "Velocity V" ; V_VELOCITY:history = "From glb_expt_01_9_3zv.nc" ; V_VELOCITY:units = "cm/s" ; // global attributes: :history = "FERRET V5.50 22-May-03" ; :title = "myTest" ; }Anyone can tell how to eliminate the duplication?Thanks-- (Bob) Yunyue YU Code 971, Goddard Space Flight Center Greenbelt, MD 20771 (301)614-6850 yuy@nemo.gsfc.nasa.gov
-- (Bob) Yunyue YU Code 971, Goddard Space Flight Center Greenbelt, MD 20771 (301)614-6850 yuy@nemo.gsfc.nasa.gov