Re: missing values in Skew T plot

From: Dennis Shea <shea_at_nyahnyahspammersnyahnyah>
Date: Mon, 27 Aug 2007 21:06:09 -0600

In the skewT_PlotData function

Change

                        ; calculating thermodynamic quantities
   idx = ind (.not.ismissing(P) .and. .not.ismissing(TC) .and. \
              .not.ismissing(TDC) .and. P.ge.100. )

to

   idx = ind (.not.ismissing(P) .and. P.ge.100. )

=====
In the main script

   nlvl = 17
   ncol = 6
   TestData = asciiread (diri+fili , (/nlvl,ncol/), "float")

   p = TestData (:,0)
   z = TestData (:,1)
   tc = TestData (:,2)
   tdc = TestData (:,3)
   wspd = TestData (:,4)
   wdir = TestData (:,5)
   tdc@_FillValue = 32766 ; <=====

Good luck
D

Mateus da Silva Teixeira wrote:
> -----BEGIN PGP SIGNED MESSAGE-----
> Hash: SHA1
>
> Dear NCL users,
>
> I want to plot a Skew-T graph with the following data:
>
> (0) 1000 125 22.94 22.266 4.85489 147.619
> (1) 925 797 18.44 17.7892 5.54436 140.856
> (2) 850 1519 14.34 13.5486 6.90072 90.8304
> (3) 700 3141 6.44 5.39017 5.28015 37.304
> (4) 600 4396 0.740021 -2.48837 5.80774 357.039
> (5) 500 5840 -7.06 -26.4358 5.33105 293.198
> (6) 400 7546 -17.46 -42.1564 9.97247 263.089
> (7) 300 9633 -33.56 -41.7041 18.4329 283.81
> (8) 250 10887 -43.36 32766 23.7697 292.249
> (9) 200 12350 -55.16 32766 26.6308 292.521
> (10) 150 14140 -66.36 32766 28.1888 286.699
> (11) 100 16530 -74.06 32766 15.0283 286.626
> (12) 70 18629 -69.06 32766 6.74758 39.5881
> (13) 50 20676 -61.76 32766 12.5539 67.5206
> (14) 30 23900 -54.36 32766 16.1805 100.685
> (15) 20 26530 -49.36 32766 21.4009 89.4646
> (16) 10 31207 -36.06 32766 19.4653 65.4099
>
> in which only Td has missing values above 300 hPa (_FillValue = 32766).
> The wind and Z data are plotted in the graph, but the line of T stop in
> the last level where Td has valid data.
>
> Is this the default behavior of Skew-T plot data function?
>
> Please, see the attached figure.
>
> Thanks,
>
> Mateus
>
> -----BEGIN PGP SIGNATURE-----
> Version: GnuPG v1.4.5 (GNU/Linux)
> Comment: Using GnuPG with Mozilla - http://enigmail.mozdev.org
>
> iD8DBQFG0s96tBm/Y0ahARYRAln0AJ43utyGQuN3I3YvW3WajdZA6Uz9uQCffUyo
> aXgArMYYyryfoc6KSgvzmvs=
> =bGxK
> -----END PGP SIGNATURE-----
>
>
> ------------------------------------------------------------------------
>
>
> ------------------------------------------------------------------------
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk_at_ucar.edu
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk

_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Mon Aug 27 2007 - 21:06:09 MDT

This archive was generated by hypermail 2.2.0 : Thu Aug 30 2007 - 21:26:07 MDT