[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

19990713: lifted parcel properties

>From: Raymond Tanabe <address@hidden>
>Organization: .
>Keywords: 199907140223.UAA08240

>Gempak Support,
>
>I am wondering is there is a way to print out certain properties ( such as
>Temperature) of a lifted parcel in using SNLIST.  I know that this is
>possible using SFLIST using LTMP but unfortunately LTMP cannot be used
>with SNLIST.  Right not, I am using SNLIST to print out various properties of
>a sounding from 1000 to 100mb in fifty mb steps.  Is there a way that I
>can also list the temp of a lifted parcel at these steps as well?  I saw
>that the formula for CAPE uses the potential temperature of a lifted
>parcel and the potential temperature of the environment in its
>calculations.  Since the potential temperature of the environment can be
>listed using THTA, is there something similar for potential temperature of
>a lifted parcel?
>
>Aloha...Ray
>
>##########################################
>Ray Tanabe - Graduate Research Assistant
>Pacific ENSO Applications Center
>c/o Department of Meteorology, HIG Room 350
>University of Hawaii - Manoa Campus
>2525 Correa Road - Honolulu, HI 96822
>Tel: 808-956-2333  Fax: 808-956-2877
>

Ray,

You should be able to use the LTMP function in SNLIST as SNPARM.
LTMP will give you the temperature Kelvin of the parcel lifted from each
level to the specified level (default is 500 mb). Since LTMP will need to
calculate THTE, levels above the Dewpoint cutoff will not be calculated. 
In the example below, I specified LTMP!500 to get the temperature of the parcel
at 500mb. If I omitted the !500, then I would get the default level which
happens to be 500mb, but you caould choose 700mb or some other level. Eg:

 SNFILE   = $UPA/19990714_upa.gem
 AREA     = @ito
 DATTIM   = 1200
 SNPARM   = tmpc;tmpk;ltmp!500
 STNDEX   = 
 LEVELS   = all
 VCOORD   = pres
 OUTPUT   = t
 MRGDAT   = y
 GEMPAK-SNLIST>r

 SNPARM = PRES;TMPC;TMPK;LTMP                                                   

 STID = ITO           STNM =    91285   TIME = 990714/1200         
 SLAT =  19.71     SLON =  -155.07   SELV =    11.0
 STIM =  1200

      PRES     TMPC     TMPK     LTMP
   1016.00    20.60   293.75   260.71
   1000.00    21.40   294.55   265.43
    982.23    20.39   293.54   265.03
    948.18    18.40   291.55   264.27
    925.00    17.00   290.15   263.76
    915.23    16.42   289.57   263.48
    882.90    14.47   287.62   262.58
    850.00    12.40   285.55   261.70
    840.00    11.60   284.75   260.37
    821.29    11.12   284.27   261.57
    809.00    10.80   283.95   262.35
    794.00     9.00   282.15   258.37
    791.81     8.12   281.27   251.01
    790.00     7.40   280.55   247.44
    780.00     9.60   282.75   250.93
    763.23     9.10   282.25   251.47
    760.00     9.00   282.15   251.57
    735.71    11.91   285.06   255.28
    735.00    12.00   285.15   255.43
    700.00    10.00   283.15   257.20
    669.00     9.40   282.55   259.99
    658.76     8.43   281.58   260.25
    611.12     3.70   276.85   261.43
    566.86    -1.03   272.12   262.54
    552.00    -2.70   270.45   262.91
    545.75    -3.30   269.85   263.19
    525.04    -5.33   267.82   264.11
    500.00    -7.90   265.25   265.25
    485.39    -9.26   263.89   266.14
    400.00   -18.10   255.05   271.84
    397.28   -18.46   254.69   271.98
    338.00   -27.10   246.05   275.18
    322.21   -30.23   242.92   275.42
    300.00   -34.90   238.25   275.69
    282.53   -38.39   234.76   276.35
    258.29   -43.60   229.55   277.22
    250.00   -45.50   227.65   277.51
    246.83   -46.19   226.96   277.68
    224.75   -51.23   221.92   278.89
    200.00   -57.50   215.65   280.19
    195.01   -58.70   214.45   280.65
    177.00   -63.30   209.85   282.34
    176.79   -63.34   209.81   282.38
    159.83   -66.51   206.64   286.25
    150.00   -68.50   204.65   288.67
    138.00   -65.70   207.45   299.68
    130.65   -65.85   207.30   304.18
    124.24   -65.99   207.16   308.38
    118.15   -66.13   207.02   312.63
    111.00   -66.30   206.85   317.99
    101.53   -69.38   203.77   321.34
    100.00   -69.90   203.25   321.91
     98.00   -70.30   202.85   323.14
     96.45   -70.42   202.73   324.43
     90.20   -70.90   202.25   329.91
     77.40   -66.10   207.05   352.83
     71.02   -66.79   206.36   360.42
     70.00   -66.90   206.25   361.71
     55.36   -62.85   210.30   394.38
     50.00   -61.10   212.05   409.40
     45.46   -60.17   212.98   422.53
     30.81   -56.36   216.79   480.68
     30.00   -56.10   217.05   484.91
     24.26   -54.28   218.87 -9999.00
     23.20   -53.90   219.25 -9999.00
     20.00   -48.10   225.05 -9999.00
     17.51   -46.45   226.70 -9999.00
     15.26   -44.74   228.41 -9999.00
     14.50   -44.10   229.05 -9999.00
     11.10   -43.38   229.77 -9999.00
     10.00   -43.10   230.05 -9999.00
      7.09   -38.87   234.28 -9999.00
      7.00   -38.70   234.45 -9999.00
      5.00   -38.90   234.25 -9999.00
      4.76   -38.28   234.87 -9999.00
      4.40   -37.30   235.85 -9999.00


There is a grid function TMST(THTE,PRES) for use with gridded data
where you can input the parcel theta-e and determine the temperature
at a given lev. Using THTE(PRES,TMPC,DWPC) at the level to lift from,
you can obtain the profile of parcel temperature in gdprof.

Steve Chiswell