Thanks Mary
I ran a test to see if this was doing what I need, and this is the result.
ncl 33> curl=(/ (/ (/1.,2.,3./), (/1.,2.,3./), (/1.,2.,3./) /), (/
(/1.,2.,3./)
ncl 33> print(curl)
Variable: curl
Type: float
Total Size: 72 bytes
18 values
Number of Dimensions: 3
Dimensions and sizes: [2] x [3] x [3]
Coordinates:
(0,0,0) 1
(0,0,1) 2
(0,0,2) 3
(0,1,0) 1
(0,1,1) 2
(0,1,2) 3
(0,2,0) 1
(0,2,1) 2
(0,2,2) 3
(1,0,0) 1
(1,0,1) 2
(1,0,2) 3
ncl 34> f3d = conform(curl,f,1)
ncl 35> curl = curl/(rho * f3d)
ncl 36> print(curl)
Variable: curl
Type: float
Total Size: 72 bytes
18 values
Number of Dimensions: 3
Dimensions and sizes: [2] x [3] x [3]
Coordinates:
(0,0,0) 9.765625e-05
(0,0,1) 0.0001953125
(0,0,2) 0.0002929688
(0,1,0) 9.765625e-05
(0,1,1) 0.0001953125
(0,1,2) 0.0002929688
(0,2,0) 9.765625e-05
(0,2,1) 0.0001953125
(0,2,2) 0.0002929688
(1,0,0) 9.765625e-05
(1,0,1) 0.0001953125
(1,0,2) 0.0002929688
--More--
Thanks for your help. I had been playing with the conform_dims function,
but did not have it quite right yet.
Sam
On 06/23/2013 02:28 PM, Mary Haley wrote:
> Hi Sam,
>
> You can use the "conform" or "conform_dims" function to conform the
> "f" array to be the same size as the "curl" array. You don't need to
> conform scalars.
>
> http://www.ncl.ucar.edu/Document/Functions/Built-in/conform_dims.shtml
>
>
> curl = curl/(rho * conform(curl,1,f)) ; "f" conforms to the 1-th
> (second) dimesion of "curl"
>
> If you want to be more explicit about what you're doing:
>
>
> f3d = conform(curl,1,f) ; "f" conforms to the 1-th (second)
Mary - I think parameters 2 and 3 are reversed in the above statement,
so they are reversed in the test I presented above.
> dimesion of "curl"
> curl = curl/(rho * f3d)
> delete(f3d) ; not necessary
>
>
> --Mary
>
>
-- Sam McClatchie, Supervisory oceanographer, Fisheries oceanography Southwest Fisheries Science Center, NOAA, 8901 La Jolla Shores Dr. La Jolla, CA 92037-1509, U.S.A. email <Sam.McClatchie@noaa.gov> Office: 858 546 7083, Cellular: 858 752 8495 Research home page <www.fishocean.info> /\ ...>><xX(¡> //// \\\\ <¡)Xx><< ///// \\\\\\ ><(((¡> >><(((¡> ...>><xX(¡>O<¡)Xx><<
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Mon Jun 24 10:26:07 2013
This archive was generated by hypermail 2.1.8 : Mon Jul 01 2013 - 12:35:42 MDT