Re: WRF data fro streamlines

From: Erik Noble <enoble_at_nyahnyahspammersnyahnyah>
Date: Tue, 23 Oct 2007 18:11:31 -0400

Hi. Thank you.
Here is my code "Test-File-WRF-Streamlines.ncl" (attached as well).
Below are my errors, too (attached as well).
-Erik

*******************************************************
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/wrf/WRFUserARW.ncl"

begin

a=addfile("/Volumes/Data_and_Models/Model-Output/Athena/WRF-SOP3-
Athena_Phys-3-1-1-1.nc","r")

times = wrf_user_list_times(a)
do it = 324,324 ;;;;;Time step you want!!!!
print("Working on time: " + times(it) )

u = wrf_user_getvar(a,"ua",it) ; u averaged to mass points
v = wrf_user_getvar(a,"va",it) ; v averaged to mass points
p = wrf_user_getvar(a, "pressure",it) ; pressure
pressure = 700.
u_plane = wrf_user_intrp3d( u,p,"h",pressure,0.,False)
v_plane = wrf_user_intrp3d( v,p,"h",pressure,0.,False)
u_plane_at_lat2D = a->XLAT_U(0,:,:) ; do you need XLAT_U, or XLAT?
u_plane_at_lon2D = a->XLONG_U(0,:,:) ; do you need XLONG_U, or XLONG?
v_plane_at_lat2D = a->XLAT_V(0,:,:) ; do you need XLAT_V, or XLAT?
v_plane_at_lon2D = a->XLONG_V ; do you need XLONG_V, or XLONG?

printVarSummary(u)
printVarSummary(u_plane)
printVarSummary(a->XLAT_U)
printVarSummary(a->XLAT)

wks = gsn_open_wks("ps","test")
mpres = True ; Create map background
mpres_at_gsnMaximize = True
mpres_at_tiMainString = "Streamlines" ; some title
mpres_at_mpGeophysicalLineThicknessF = 3.0
mpres_at_gsnAddCyclic = False
plot = gsn_csm_streamline_map(wks,u_plane,v_plane,mpres)
mpres_at_mpMinLatF = -30.
mpres_at_mpMaxLatF = 20.
mpres_at_mpMinLonF = -40.
mpres_at_mpMaxLonF = 40.
plot2 = gsn_csm_streamline_map(wks,u_plane,v_plane,mpres)
end do
end

*******************************

here are my errors:
  Copyright (C) 1995-2007 - All Rights Reserved
  University Corporation for Atmospheric Research
  NCAR Command Language Version 4.3.1
  The use of this software is governed by a License Agreement.
  See http://www.ncl.ucar.edu/ for more details.

Variable: times
Type: string
Total Size: 1948 bytes
             487 values
Number of Dimensions: 1
Dimensions and sizes: [487]
Coordinates:
Number Of Attributes: 2
   description : times in file
   _FillValue : missing
(0) 2006-08-01_00:00:00
(1) 2006-08-01_03:00:00
(2) 2006-08-01_06:00:00
(3) 2006-08-01_09:00:00
(4) 2006-08-01_12:00:00
(5) 2006-08-01_15:00:00
(6) 2006-08-01_18:00:00
(7) 2006-08-01_21:00:00
(8) 2006-08-02_00:00:00
(9) 2006-08-02_03:00:00
(10) 2006-08-02_06:00:00
(11) 2006-08-02_09:00:00
(12) 2006-08-02_12:00:00
(13) 2006-08-02_15:00:00
(14) 2006-08-02_18:00:00
(15) 2006-08-02_21:00:00
(16) 2006-08-03_00:00:00
(17) 2006-08-03_03:00:00
(18) 2006-08-03_06:00:00
(19) 2006-08-03_09:00:00
(20) 2006-08-03_12:00:00
(21) 2006-08-03_15:00:00
(22) 2006-08-03_18:00:00
(23) 2006-08-03_21:00:00
(24) 2006-08-04_00:00:00
(25) 2006-08-04_03:00:00
(26) 2006-08-04_06:00:00
(27) 2006-08-04_09:00:00
(28) 2006-08-04_12:00:00
(29) 2006-08-04_15:00:00
(30) 2006-08-04_18:00:00
(31) 2006-08-04_21:00:00
(32) 2006-08-05_00:00:00
(33) 2006-08-05_03:00:00
(34) 2006-08-05_06:00:00
(35) 2006-08-05_09:00:00
(36) 2006-08-05_12:00:00
(37) 2006-08-05_15:00:00
(38) 2006-08-05_18:00:00
(39) 2006-08-05_21:00:00
(40) 2006-08-06_00:00:00
(41) 2006-08-06_03:00:00
(42) 2006-08-06_06:00:00
(43) 2006-08-06_09:00:00
(44) 2006-08-06_12:00:00
(45) 2006-08-06_15:00:00
(46) 2006-08-06_18:00:00
(47) 2006-08-06_21:00:00
(48) 2006-08-07_00:00:00
(49) 2006-08-07_03:00:00
(50) 2006-08-07_06:00:00
(51) 2006-08-07_09:00:00
(52) 2006-08-07_12:00:00
(53) 2006-08-07_15:00:00
(54) 2006-08-07_18:00:00
(55) 2006-08-07_21:00:00
(56) 2006-08-08_00:00:00
(57) 2006-08-08_03:00:00
(58) 2006-08-08_06:00:00
(59) 2006-08-08_09:00:00
(60) 2006-08-08_12:00:00
(61) 2006-08-08_15:00:00
(62) 2006-08-08_18:00:00
(63) 2006-08-08_21:00:00
(64) 2006-08-09_00:00:00
(65) 2006-08-09_03:00:00
(66) 2006-08-09_06:00:00
(67) 2006-08-09_09:00:00
(68) 2006-08-09_12:00:00
(69) 2006-08-09_15:00:00
(70) 2006-08-09_18:00:00
(71) 2006-08-09_21:00:00
(72) 2006-08-10_00:00:00
(73) 2006-08-10_03:00:00
(74) 2006-08-10_06:00:00
(75) 2006-08-10_09:00:00
(76) 2006-08-10_12:00:00
(77) 2006-08-10_15:00:00
(78) 2006-08-10_18:00:00
(79) 2006-08-10_21:00:00
(80) 2006-08-11_00:00:00
(81) 2006-08-11_03:00:00
(82) 2006-08-11_06:00:00
(83) 2006-08-11_09:00:00
(84) 2006-08-11_12:00:00
(85) 2006-08-11_15:00:00
(86) 2006-08-11_18:00:00
(87) 2006-08-11_21:00:00
(88) 2006-08-12_00:00:00
(89) 2006-08-12_03:00:00
(90) 2006-08-12_06:00:00
(91) 2006-08-12_09:00:00
(92) 2006-08-12_12:00:00
(93) 2006-08-12_15:00:00
(94) 2006-08-12_18:00:00
(95) 2006-08-12_21:00:00
(96) 2006-08-13_00:00:00
(97) 2006-08-13_03:00:00
(98) 2006-08-13_06:00:00
(99) 2006-08-13_09:00:00
(100) 2006-08-13_12:00:00
(101) 2006-08-13_15:00:00
(102) 2006-08-13_18:00:00
(103) 2006-08-13_21:00:00
(104) 2006-08-14_00:00:00
(105) 2006-08-14_03:00:00
(106) 2006-08-14_06:00:00
(107) 2006-08-14_09:00:00
(108) 2006-08-14_12:00:00
(109) 2006-08-14_15:00:00
(110) 2006-08-14_18:00:00
(111) 2006-08-14_21:00:00
(112) 2006-08-15_00:00:00
(113) 2006-08-15_03:00:00
(114) 2006-08-15_06:00:00
(115) 2006-08-15_09:00:00
(116) 2006-08-15_12:00:00
(117) 2006-08-15_15:00:00
(118) 2006-08-15_18:00:00
(119) 2006-08-15_21:00:00
(120) 2006-08-16_00:00:00
(121) 2006-08-16_03:00:00
(122) 2006-08-16_06:00:00
(123) 2006-08-16_09:00:00
(124) 2006-08-16_12:00:00
(125) 2006-08-16_15:00:00
(126) 2006-08-16_18:00:00
(127) 2006-08-16_21:00:00
(128) 2006-08-17_00:00:00
(129) 2006-08-17_03:00:00
(130) 2006-08-17_06:00:00
(131) 2006-08-17_09:00:00
(132) 2006-08-17_12:00:00
(133) 2006-08-17_15:00:00
(134) 2006-08-17_18:00:00
(135) 2006-08-17_21:00:00
(136) 2006-08-18_00:00:00
(137) 2006-08-18_03:00:00
(138) 2006-08-18_06:00:00
(139) 2006-08-18_09:00:00
(140) 2006-08-18_12:00:00
(141) 2006-08-18_15:00:00
(142) 2006-08-18_18:00:00
(143) 2006-08-18_21:00:00
(144) 2006-08-19_00:00:00
(145) 2006-08-19_03:00:00
(146) 2006-08-19_06:00:00
(147) 2006-08-19_09:00:00
(148) 2006-08-19_12:00:00
(149) 2006-08-19_15:00:00
(150) 2006-08-19_18:00:00
(151) 2006-08-19_21:00:00
(152) 2006-08-20_00:00:00
(153) 2006-08-20_03:00:00
(154) 2006-08-20_06:00:00
(155) 2006-08-20_09:00:00
(156) 2006-08-20_12:00:00
(157) 2006-08-20_15:00:00
(158) 2006-08-20_18:00:00
(159) 2006-08-20_21:00:00
(160) 2006-08-21_00:00:00
(161) 2006-08-21_03:00:00
(162) 2006-08-21_06:00:00
(163) 2006-08-21_09:00:00
(164) 2006-08-21_12:00:00
(165) 2006-08-21_15:00:00
(166) 2006-08-21_18:00:00
(167) 2006-08-21_21:00:00
(168) 2006-08-22_00:00:00
(169) 2006-08-22_03:00:00
(170) 2006-08-22_06:00:00
(171) 2006-08-22_09:00:00
(172) 2006-08-22_12:00:00
(173) 2006-08-22_15:00:00
(174) 2006-08-22_18:00:00
(175) 2006-08-22_21:00:00
(176) 2006-08-23_00:00:00
(177) 2006-08-23_03:00:00
(178) 2006-08-23_06:00:00
(179) 2006-08-23_09:00:00
(180) 2006-08-23_12:00:00
(181) 2006-08-23_15:00:00
(182) 2006-08-23_18:00:00
(183) 2006-08-23_21:00:00
(184) 2006-08-24_00:00:00
(185) 2006-08-24_03:00:00
(186) 2006-08-24_06:00:00
(187) 2006-08-24_09:00:00
(188) 2006-08-24_12:00:00
(189) 2006-08-24_15:00:00
(190) 2006-08-24_18:00:00
(191) 2006-08-24_21:00:00
(192) 2006-08-25_00:00:00
(193) 2006-08-25_03:00:00
(194) 2006-08-25_06:00:00
(195) 2006-08-25_09:00:00
(196) 2006-08-25_12:00:00
(197) 2006-08-25_15:00:00
(198) 2006-08-25_18:00:00
(199) 2006-08-25_21:00:00
(200) 2006-08-26_00:00:00
(201) 2006-08-26_03:00:00
(202) 2006-08-26_06:00:00
(203) 2006-08-26_09:00:00
(204) 2006-08-26_12:00:00
(205) 2006-08-26_15:00:00
(206) 2006-08-26_18:00:00
(207) 2006-08-26_21:00:00
(208) 2006-08-27_00:00:00
(209) 2006-08-27_03:00:00
(210) 2006-08-27_06:00:00
(211) 2006-08-27_09:00:00
(212) 2006-08-27_12:00:00
(213) 2006-08-27_15:00:00
(214) 2006-08-27_18:00:00
(215) 2006-08-27_21:00:00
(216) 2006-08-28_00:00:00
(217) 2006-08-28_03:00:00
(218) 2006-08-28_06:00:00
(219) 2006-08-28_09:00:00
(220) 2006-08-28_12:00:00
(221) 2006-08-28_15:00:00
(222) 2006-08-28_18:00:00
(223) 2006-08-28_21:00:00
(224) 2006-08-29_00:00:00
(225) 2006-08-29_03:00:00
(226) 2006-08-29_06:00:00
(227) 2006-08-29_09:00:00
(228) 2006-08-29_12:00:00
(229) 2006-08-29_15:00:00
(230) 2006-08-29_18:00:00
(231) 2006-08-29_21:00:00
(232) 2006-08-30_00:00:00
(233) 2006-08-30_03:00:00
(234) 2006-08-30_06:00:00
(235) 2006-08-30_09:00:00
(236) 2006-08-30_12:00:00
(237) 2006-08-30_15:00:00
(238) 2006-08-30_18:00:00
(239) 2006-08-30_21:00:00
(240) 2006-08-31_00:00:00
(241) 2006-08-31_03:00:00
(242) 2006-08-31_06:00:00
(243) 2006-08-31_09:00:00
(244) 2006-08-31_12:00:00
(245) 2006-08-31_15:00:00
(246) 2006-08-31_18:00:00
(247) 2006-08-31_21:00:00
(248) 2006-09-01_00:00:00
(249) 2006-09-01_03:00:00
(250) 2006-09-01_06:00:00
(251) 2006-09-01_09:00:00
(252) 2006-09-01_12:00:00
(253) 2006-09-01_15:00:00
(254) 2006-09-01_18:00:00
(255) 2006-09-01_21:00:00
(256) 2006-09-02_00:00:00
(257) 2006-09-02_03:00:00
(258) 2006-09-02_06:00:00
(259) 2006-09-02_09:00:00
(260) 2006-09-02_12:00:00
(261) 2006-09-02_15:00:00
(262) 2006-09-02_18:00:00
(263) 2006-09-02_21:00:00
(264) 2006-09-03_00:00:00
(265) 2006-09-03_03:00:00
(266) 2006-09-03_06:00:00
(267) 2006-09-03_09:00:00
(268) 2006-09-03_12:00:00
(269) 2006-09-03_15:00:00
(270) 2006-09-03_18:00:00
(271) 2006-09-03_21:00:00
(272) 2006-09-04_00:00:00
(273) 2006-09-04_03:00:00
(274) 2006-09-04_06:00:00
(275) 2006-09-04_09:00:00
(276) 2006-09-04_12:00:00
(277) 2006-09-04_15:00:00
(278) 2006-09-04_18:00:00
(279) 2006-09-04_21:00:00
(280) 2006-09-05_00:00:00
(281) 2006-09-05_03:00:00
(282) 2006-09-05_06:00:00
(283) 2006-09-05_09:00:00
(284) 2006-09-05_12:00:00
(285) 2006-09-05_15:00:00
(286) 2006-09-05_18:00:00
(287) 2006-09-05_21:00:00
(288) 2006-09-06_00:00:00
(289) 2006-09-06_03:00:00
(290) 2006-09-06_06:00:00
(291) 2006-09-06_09:00:00
(292) 2006-09-06_12:00:00
(293) 2006-09-06_15:00:00
(294) 2006-09-06_18:00:00
(295) 2006-09-06_21:00:00
(296) 2006-09-07_00:00:00
(297) 2006-09-07_03:00:00
(298) 2006-09-07_06:00:00
(299) 2006-09-07_09:00:00
(300) 2006-09-07_12:00:00
(301) 2006-09-07_15:00:00
(302) 2006-09-07_18:00:00
(303) 2006-09-07_21:00:00
(304) 2006-09-08_00:00:00
(305) 2006-09-08_03:00:00
(306) 2006-09-08_06:00:00
(307) 2006-09-08_09:00:00
(308) 2006-09-08_12:00:00
(309) 2006-09-08_15:00:00
(310) 2006-09-08_18:00:00
(311) 2006-09-08_21:00:00
(312) 2006-09-09_00:00:00
(313) 2006-09-09_03:00:00
(314) 2006-09-09_06:00:00
(315) 2006-09-09_09:00:00
(316) 2006-09-09_12:00:00
(317) 2006-09-09_15:00:00
(318) 2006-09-09_18:00:00
(319) 2006-09-09_21:00:00
(320) 2006-09-10_00:00:00
(321) 2006-09-10_03:00:00
(322) 2006-09-10_06:00:00
(323) 2006-09-10_09:00:00
(324) 2006-09-10_12:00:00
(325) 2006-09-10_15:00:00
(326) 2006-09-10_18:00:00
(327) 2006-09-10_21:00:00
(328) 2006-09-11_00:00:00
(329) 2006-09-11_03:00:00
(330) 2006-09-11_06:00:00
(331) 2006-09-11_09:00:00
(332) 2006-09-11_12:00:00
(333) 2006-09-11_15:00:00
(334) 2006-09-11_18:00:00
(335) 2006-09-11_21:00:00
(336) 2006-09-12_00:00:00
(337) 2006-09-12_03:00:00
(338) 2006-09-12_06:00:00
(339) 2006-09-12_09:00:00
(340) 2006-09-12_12:00:00
(341) 2006-09-12_15:00:00
(342) 2006-09-12_18:00:00
(343) 2006-09-12_21:00:00
(344) 2006-09-13_00:00:00
(345) 2006-09-13_03:00:00
(346) 2006-09-13_06:00:00
(347) 2006-09-13_09:00:00
(348) 2006-09-13_12:00:00
(349) 2006-09-13_15:00:00
(350) 2006-09-13_18:00:00
(351) 2006-09-13_21:00:00
(352) 2006-09-14_00:00:00
(353) 2006-09-14_03:00:00
(354) 2006-09-14_06:00:00
(355) 2006-09-14_09:00:00
(356) 2006-09-14_12:00:00
(357) 2006-09-14_15:00:00
(358) 2006-09-14_18:00:00
(359) 2006-09-14_21:00:00
(360) 2006-09-15_00:00:00
(361) 2006-09-15_03:00:00
(362) 2006-09-15_06:00:00
(363) 2006-09-15_09:00:00
(364) 2006-09-15_12:00:00
(365) 2006-09-15_15:00:00
(366) 2006-09-15_18:00:00
(367) 2006-09-15_21:00:00
(368) 2006-09-16_00:00:00
(369) 2006-09-16_03:00:00
(370) 2006-09-16_06:00:00
(371) 2006-09-16_09:00:00
(372) 2006-09-16_12:00:00
(373) 2006-09-16_15:00:00
(374) 2006-09-16_18:00:00
(375) 2006-09-16_21:00:00
(376) 2006-09-17_00:00:00
(377) 2006-09-17_03:00:00
(378) 2006-09-17_06:00:00
(379) 2006-09-17_09:00:00
(380) 2006-09-17_12:00:00
(381) 2006-09-17_15:00:00
(382) 2006-09-17_18:00:00
(383) 2006-09-17_21:00:00
(384) 2006-09-18_00:00:00
(385) 2006-09-18_03:00:00
(386) 2006-09-18_06:00:00
(387) 2006-09-18_09:00:00
(388) 2006-09-18_12:00:00
(389) 2006-09-18_15:00:00
(390) 2006-09-18_18:00:00
(391) 2006-09-18_21:00:00
(392) 2006-09-19_00:00:00
(393) 2006-09-19_03:00:00
(394) 2006-09-19_06:00:00
(395) 2006-09-19_09:00:00
(396) 2006-09-19_12:00:00
(397) 2006-09-19_15:00:00
(398) 2006-09-19_18:00:00
(399) 2006-09-19_21:00:00
(400) 2006-09-20_00:00:00
(401) 2006-09-20_03:00:00
(402) 2006-09-20_06:00:00
(403) 2006-09-20_09:00:00
(404) 2006-09-20_12:00:00
(405) 2006-09-20_15:00:00
(406) 2006-09-20_18:00:00
(407) 2006-09-20_21:00:00
(408) 2006-09-21_00:00:00
(409) 2006-09-21_03:00:00
(410) 2006-09-21_06:00:00
(411) 2006-09-21_09:00:00
(412) 2006-09-21_12:00:00
(413) 2006-09-21_15:00:00
(414) 2006-09-21_18:00:00
(415) 2006-09-21_21:00:00
(416) 2006-09-22_00:00:00
(417) 2006-09-22_03:00:00
(418) 2006-09-22_06:00:00
(419) 2006-09-22_09:00:00
(420) 2006-09-22_12:00:00
(421) 2006-09-22_15:00:00
(422) 2006-09-22_18:00:00
(423) 2006-09-22_21:00:00
(424) 2006-09-23_00:00:00
(425) 2006-09-23_03:00:00
(426) 2006-09-23_06:00:00
(427) 2006-09-23_09:00:00
(428) 2006-09-23_12:00:00
(429) 2006-09-23_15:00:00
(430) 2006-09-23_18:00:00
(431) 2006-09-23_21:00:00
(432) 2006-09-24_00:00:00
(433) 2006-09-24_03:00:00
(434) 2006-09-24_06:00:00
(435) 2006-09-24_09:00:00
(436) 2006-09-24_12:00:00
(437) 2006-09-24_15:00:00
(438) 2006-09-24_18:00:00
(439) 2006-09-24_21:00:00
(440) 2006-09-25_00:00:00
(441) 2006-09-25_03:00:00
(442) 2006-09-25_06:00:00
(443) 2006-09-25_09:00:00
(444) 2006-09-25_12:00:00
(445) 2006-09-25_15:00:00
(446) 2006-09-25_18:00:00
(447) 2006-09-25_21:00:00
(448) 2006-09-26_00:00:00
(449) 2006-09-26_03:00:00
(450) 2006-09-26_06:00:00
(451) 2006-09-26_09:00:00
(452) 2006-09-26_12:00:00
(453) 2006-09-26_15:00:00
(454) 2006-09-26_18:00:00
(455) 2006-09-26_21:00:00
(456) 2006-09-27_00:00:00
(457) 2006-09-27_03:00:00
(458) 2006-09-27_06:00:00
(459) 2006-09-27_09:00:00
(460) 2006-09-27_12:00:00
(461) 2006-09-27_15:00:00
(462) 2006-09-27_18:00:00
(463) 2006-09-27_21:00:00
(464) 2006-09-28_00:00:00
(465) 2006-09-28_03:00:00
(466) 2006-09-28_06:00:00
(467) 2006-09-28_09:00:00
(468) 2006-09-28_12:00:00
(469) 2006-09-28_15:00:00
(470) 2006-09-28_18:00:00
(471) 2006-09-28_21:00:00
(472) 2006-09-29_00:00:00
(473) 2006-09-29_03:00:00
(474) 2006-09-29_06:00:00
(475) 2006-09-29_09:00:00
(476) 2006-09-29_12:00:00
(477) 2006-09-29_15:00:00
(478) 2006-09-29_18:00:00
(479) 2006-09-29_21:00:00
(480) 2006-09-30_00:00:00
(481) 2006-09-30_03:00:00
(482) 2006-09-30_06:00:00
(483) 2006-09-30_09:00:00
(484) 2006-09-30_12:00:00
(485) 2006-09-30_15:00:00
(486) 2006-09-30_18:00:00
(0) Working on time: 2006-09-10_12:00:00

Variable: u
Type: float
Total Size: 1636308 bytes
             409077 values
Number of Dimensions: 3
Dimensions and sizes: [27] x [109] x [139]
Coordinates:
Number Of Attributes: 2
   description : u Velocity
   units : m/s

Variable: u_plane
Type: float
Total Size: 60604 bytes
             15151 values
Number of Dimensions: 2
Dimensions and sizes: [109] x [139]
Coordinates:
Number Of Attributes: 6
   lon2D : <ARRAY of 15260 elements>
   lat2D : <ARRAY of 15260 elements>
   description : u Velocity
   units : m/s
   _FillValue : -999999
   PlotLevelID : 700 hPa

Variable: XLAT_U (file variable)
Type: float
Total Size: 29726480 bytes
             7431620 values
Number of Dimensions: 3
Dimensions and sizes: [Time | 487] x [south_north | 109] x
[west_east_stag | 140]
Coordinates:
Number Of Attributes: 6
   FieldType : 104
   MemoryOrder : XY
   description : LATITUDE, SOUTH IS NEGATIVE
   units : degree_north
   stagger : X
   coordinates : XLONG_U XLAT_U

Variable: XLAT (file variable)
Type: float
Total Size: 29514148 bytes
             7378537 values
Number of Dimensions: 3
Dimensions and sizes: [Time | 487] x [south_north | 109] x [west_east
| 139]
Coordinates:
Number Of Attributes: 5
   FieldType : 104
   MemoryOrder : XY
   description : LATITUDE, SOUTH IS NEGATIVE
   units : degree_north
   stagger :
(0) check_for_y_lat_coord: Warning: Data either does not contain a
valid latitude coordinate array or doesn't contain one at all.
(0) A valid latitude coordinate array should have a 'units' attribute
equal to one of the following values:
(0) 'degrees_north' 'degrees-north' 'degree_north' 'degrees
north' 'degrees_N' 'Degrees_north' 'degree_N' 'degreeN' 'degreesN'
'deg north'
(0) check_for_lon_coord: Warning: Data either does not contain a
valid longitude coordinate array or doesn't contain one at all.
(0) A valid longitude coordinate array should have a 'units'
attribute equal to one of the following values:
(0) 'degrees_east' 'degrees-east' 'degree_east' 'degrees east'
'degrees_E' 'Degrees_east' 'degree_E' 'degreeE' 'degreesE' 'deg east'
(0) check_for_y_lat_coord: Warning: Data either does not contain a
valid latitude coordinate array or doesn't contain one at all.
(0) A valid latitude coordinate array should have a 'units' attribute
equal to one of the following values:
(0) 'degrees_north' 'degrees-north' 'degree_north' 'degrees
north' 'degrees_N' 'Degrees_north' 'degree_N' 'degreeN' 'degreesN'
'deg north'
(0) check_for_lon_coord: Warning: Data either does not contain a
valid longitude coordinate array or doesn't contain one at all.
(0) A valid longitude coordinate array should have a 'units'
attribute equal to one of the following values:
(0) 'degrees_east' 'degrees-east' 'degree_east' 'degrees east'
'degrees_E' 'Degrees_east' 'degree_E' 'degreeE' 'degreesE' 'deg east'


On Oct 23, 2007, at 5:28 PM, Adam Phillips wrote:

> 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/wrf/WRFUserARW.ncl"
>
> begin
> a=addfile("/Volumes/Data_and_Models/Model-Output/Athena/WRF-SOP3-
> Athena_Phys-3-1-1-1.nc","r")
> u = wrf_user_getvar(a,"ua",it) ; u averaged to mass points
> v = wrf_user_getvar(a,"va",it) ; v averaged to mass points
> p = wrf_user_getvar(a, "pressure",it) ; pressure
> pressure = 700.
> u_plane = wrf_user_intrp3d( u,p,"h",pressure,0.,False)
> v_plane = wrf_user_intrp3d( v,p,"h",pressure,0.,False)
> u_plane_at_lat2D = a->XLAT_U(0,:,:) ; do you need XLAT_U, or XLAT?
> u_plane_at_lon2D = a->XLONG_U(0,:,:) ; do you need XLONG_U, or XLONG?
> v_plane_at_lat2D = a->XLAT_V(0,:,:) ; do you need XLAT_V, or XLAT?
> v_plane_at_lon2D = a->XLONG_V ; do you need XLONG_V, or XLONG?
>
> printVarSummary(u)
> printVarSummary(u_plane)
> printVarSummary(a->XLAT_U)
> printVarSummary(a->XLAT)
>
> wks = gsn_open_wks("ps","test")
> mpres = True ; Create map background
> mpres_at_gsnMaximize = True
> mpres_at_tiMainString = "Streamlines" ; some tit
> mpres_at_mpGeophysicalLineThicknessF = 3.0
> mpres_at_gsnAddCyclic = False
> plot = gsn_csm_streamline_map(wks,u_plane,v_plane,mpres)
> mpres_at_mpMinLatF = -30.
> mpres_at_mpMaxLatF = 20.
> mpres_at_mpMinLonF = -40.
> mpres_at_mpMaxLonF = 40.
> plot2 = gsn_csm_streamline_map(wks,u_plane,v_plane,mpres)
> end

_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Tue Oct 23 2007 - 16:11:31 MDT

This archive was generated by hypermail 2.2.0 : Wed Oct 24 2007 - 09:14:00 MDT