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

RE: hurricane plots



Don and Tom,

> -----Original Message-----
> From: Don Murray [mailto:address@hidden 
> Sent: Friday, October 24, 2003 8:50 PM
> 
> How do you decode the positions?
> 
> Don

I use a Mcbasi script to parse the text bulletin - see below.

Tom Y. might be interested too, so I'll cc him.

REM ? Name:    HURPLT.MCB (version 1.0)
REM ?
REM ? Purpose: plot forecast hurricane positions on a map
REM ?
REM ? Invocation:
REM ?          RUN lat lon blowup FILE=HURPLT.MCB
REM ?               |   |    |         |__ McBASI script name
REM ?               |   |    |_____________image blowup factor
REM ?               |   |__________________longitude to center image
REM ?               |_____________________ latitude to center image
REM ?                                      
REM ?
REM   Written by Bill Fingerhut
REM   History: 030914 - Created for Unidata McIDAS-X Version 2003
REM ? ----------

REM   draw the map
    KEYIN "EG B"
    KEYIN "GU REST GRAPHIC"
    INPUT A$
    IF A$ <> "" GOTO 50
    A$ = "31"
    B$ = "79"
    C$ = "1"
    GOTO 75
50  INPUT B$
    INPUT C$
75 KEYIN "IMGDISP RTIMAGES/GE-IR LATLON=";A$;" ";B$;" MAG=";C$;" EU=SAB"

REM   retrieve the text bulletin
    KEYIN "WXTLIST WMO=WTNT4? WSTN=KNHC DTIME=6 NUM=6 DEV=F HURPLT.DAT"
    KEYIN "MAP H 6 LALO 4 INT=5 5"
    KEYIN "BAR SU=IRTEMP"

REM   open the file
    OPEN "HURPLT.DAT",10,80
100 INPUT @10,L$
    IF MID$(L$,1,4) <> EOF$ GOTO 150

REM   end of file read
    KEYIN "ECHO {end of file read"
    GOTO 900

REM   locate data
150 T$=MID$(L$,4,5)
    IF T$<>"TIAL " GOTO 200

REM   found initial position
    T$="HR VT"
    GOTO 500

REM   test for data
200 IF T$<>"HR VT" GOTO 100

REM   process data
500 D$=MID$(L$,14,2)
    H$=MID$(L$,17,2)
    Y$=MID$(L$,22,5)
    X$=MID$(L$,29,5)
    S$=MID$(L$,38,3)
    KEYIN "ECHO {D";D$;"DH";H$;"HY";Y$;"YX";X$;"XS";S$;"S"

REM   calculate plotting positions
    PTABLE STR$(VAL(Y$)+.5),"LP"
    PTABLE STR$(VAL(Y$)-.5),"LM"

REM   plot data
    KEYIN "PC T 1 1"
    KEYIN "PC E ";Y$;" ";X$"
    KEYIN "ZA 5 14 C TCY=5 N"
    KEYIN "PC E #LP ";X$;""
    KEYIN "ZA 5 12 L WID=2 {";S$;" kt"
    KEYIN "PC E #LM ";X$;""
    KEYIN "ZA 5 12 L WID=2 { ";D$;"/";H$;" GMT"
    GOTO 100

900 CLOSE 10
    KEYIN "FRMSAVE X REMOTE"


Bill

< Bill Fingerhut, Professor           PHONE: 802-626-6257 >
< Meteorology Dept                    FAX:   802-626-9770 >
< Lyndon State College                                    >
< Lyndonville, Vt 05851                                   >
<                                                         >
< EMAIL:     address@hidden            >
<            address@hidden                  >
< WWW:       http://apollo.lsc.vsc.edu/                   >
<                                                         >
< disclaimer: I know nothing - I only work here.          >


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.