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

20051004: 20051003: Gempak Question

This archive contains answers to questions sent to Unidata support through mid-2025. Note that the archive is no longer being updated. We provide the archive for reference; many of the answers presented here remain technically correct, even if somewhat outdated. For the most up-to-date information on the use of NSF Unidata software and data services, please consult the Software Documentation first.

  • Subject: 20051004: 20051003: Gempak Question
  • Date: Tue, 04 Oct 2005 14:40:01 -0600
Ryan,

No, the arinfo just provides information about the image header. It
doesn't have any function regarding dissplay and navigation of the image.

NMAP2 and GARP have methods for reading out mouse values for images,
but for command line use, you would simply need a program wrapper to
call $GEMPAK/source/gplt/satellite/gsatpx.f
(which is what the $GEMPAK/source/gemlib/im/imgpix.f routine is doing to
get a cloud height from an IR image as called by the cloud height routine
in NMAP2).

If you build the source distribution, you can build the testim program with:
cd $GEMPAK/source/gemlib/im
make tlink

Then the testim program would be a simple lat, lon readout program such as:
> testim
  1 = IM_SIMG                 2 = IM_ISUB 
  3 = IM_DROP                 4 = IM_LUTF 
  5 = IM_QSIZ                 6 = IM_GPIX 
  7 = IM_BTOT                 8 = IM_TTOB 
  9 = IM_GTMP                10 = GG_SDEV 
 11 = GG_MAP                 12 = GG_LTLN 
 13 = GG_STTL                14 =         
 15 = IM_QTIM                16 = IM_QLUT 
 17 =                        18 = IM_QCHN 
 20 = GSBRGB                 21 = GQCLRS  
 22 = GCLEAR 
 23 = GSROAM                 24 = GROAM  
 30 = DUMP COMMON 
Select a subroutine number or type EXIT:6
 Enter IMGFIL
$SAT/GOES-10/4km/IR/IR_20051004_2000
 Enter GAREA
dset
 Enter GEMPAK time
051004/2000
 Enter RLAT, RLON
20.00, -157.00
 Enter maximum search distance
1000000
 Enter stndard sounding flag
f
 Enter moist-adiabatic cloud flag
t
 Enter pixel area, 0 -> single pixel, 5 -> 11 x 11 area 
0
 Enter pixel mode
0
 [SN 5]  Stable sounding.
 IRET =            5
 [IM 5]
 IPIX =          107
 TMPK =    276.5000

If pixel values are all you need, then this can be modified as needed.
If you need to do more advanced calibration and/or pixel geometry calculations,
then McIDAS would be a better choice.

Steve Chiswell
Unidata User Support



>From: Ryan Lyman <address@hidden>
>Organization: UCAR/Unidata
>Keywords: 200510041958.j94JwKG7029444

>
>Thanks steve for the response.  though this is certainly a useful 
>command, i'm not sure if this is quite what i was looking for...i need 
>to extract satellite data from a single grid point, lat/lon, or pixel 
>centered or close to mauna kea on the Big Island using WV and IR.  Can i 
>do that through the arinfo command, or is there another way?
>
>aloha,
>ryan
>
>
>
>Unidata Support wrote:
>
>> Ryan,
>> 
>> You can use the program arinfo to list out information about
>> a an AREA file. The display of the satellite image is done
>> through the use of SATFIL with PROJ=sat in the display programs.
>> 
>> Steve Chiswell
>> Unidata User Support
>> 
>> 
>> 
>> 
>>>From: Ryan Lyman <address@hidden>
>>>Organization: UCAR/Unidata
>>>Keywords: 200510032207.j93M7WG7004462
>> 
>> 
>>>Is there a command in gempak to exract the satellite data from a 
>>>satellite area file that comes through the LDM?  basically a command 
>>>similar to sflist or gdlist, but for satellite data.
>>>
>>>aloha,
>>>ryan lyman
>>>
>>>-- 
>>>
>>>-------------------------------------------------------
>>> Ryan Lyman                   Mauna Kea Weather Center
>>> Institute for Astronomy      Forecast Meteorologist
>>> 640 North A'ohoku Place      Tel: (808) 932-2323
>>> Hilo, Hawaii 96720-2700      Fax: (808) 933-0737
>>>
>>>             http://mkwc.ifa.hawaii.edu/
>>>-------------------------------------------------------
>>>
>> 
>> --
>> ****************************************************************************
>> Unidata User Support                                    UCAR Unidata Program
>> (303)497-8643                                                  P.O. Box 3000
>> address@hidden                                   Boulder, CO 80307
>> ----------------------------------------------------------------------------
>> Unidata WWW Service              http://my.unidata.ucar.edu/content/support 
>> ----------------------------------------------------------------------------
>> NOTE: All email exchanges with Unidata User Support are recorded in the
>> Unidata inquiry tracking system and then made publicly available
>> through the web.  If you do not want to have your interactions made
>> available in this way, you must let us know in each email you send to us.
>
>-- 
>
>-------------------------------------------------------
>  Ryan Lyman                   Mauna Kea Weather Center
>  Institute for Astronomy      Forecast Meteorologist
>  640 North A'ohoku Place      Tel: (808) 932-2323
>  Hilo, Hawaii 96720-2700      Fax: (808) 933-0737
>
>              http://mkwc.ifa.hawaii.edu/
>-------------------------------------------------------
>
--
NOTE: All email exchanges with Unidata User Support are recorded in the
Unidata inquiry tracking system and then made publicly available
through the web.  If you do not want to have your interactions made
available in this way, you must let us know in each email you send to us.