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

[McIDAS #BFX-595521]: Question about the stats of the EF5 Hackleburg Tornadic Supercell at 20:10

Hi Rod,

re:
> Question 1 if I only just changed the time and kept the same area file
> number. Is that possible?

Yes as long as you make sure that you change the TIME= values consistently
throughout the BATCH file.

NOTE: The output that you sent was very confusing given the extra lines, etc.
Also, some parts of the listing looked like more than one thing (e.g., BATCH
job) were running at the same time.

re:
> BATCH EFFRAD7.BAT
> SF 1
> ERASE F 1 7
> Erased image frame(s) 1-7
> Erased graphic frame(s) 1-7
> ERASE: Done
> 
> IMGCOPY G13/IR27APR11 LOC/AREAS.3003* TIME=20:10* LAT=30:00 86:30 SIZE=320 320
> Beginning Image Data transfer, bytes= 208288
> Transferring AREA data outbound, bytes= 208448
> IMGCOPY: G13/IR27APR11.3 copied to LOC/AREAS.3003
> IMGCOPY: done
> 
> IMGDISP LOC/AREAS.3003 1 LAT=30:00 86:30 EU=IMAGE SF=YES
> Beginning Image Data transfer, bytes= 143328
> IMGDISP: loaded frame  1
> IMGDISP: done
> MAP STATE=MS AL GA COUNTY=ALL MCOL=8
> EU: Restoring IMAGE.ET to frame(s)= 1
> EU: Done
> MAP: Completed frame 1
> MAP FILE=OUTVHPOL MCOL=1
> MAP: Completed frame 1
> 
> IMGCOPY LOC/AREAS.3003 LOC/AREAS.3103 SIZE=SAME
> Beginning Image Data transfer, bytes= 208448
> Transferring AREA data outbound, bytes= 208528
> IMGCOPY: LOC/AREAS.3003 copied to LOC/AREAS.3103
> IMGCOPY: done
> 
> IMGREMAP G13/ERAD27APR11 LOC/AREAS.3103 *TIME=20:10*
> Beginning Image Data transfer, bytes= 108532
> IMGREMAP: transformations complete ... begin data move
> Transferring AREA data outbound, bytes= 106208
> IMGREMAP: Done...
> 
> IMGLIST LOC/AREAS.3003 FORM=ALL
> Image file directory listing for:LOC/AREAS
> Pos Satellite/         Date       Time      Center      Res (km) Image_Size
>     sensor                                 Lat  Lon    Lat   Lon
> --- -------------  ------------  --------  ---- ----  ----- ----- ------------
> 3003  G-13 IMG      27 APR 11117  20:10:00    30   86
> Band: 4    10.7 um IR Surface/Cloud-top Temp         5.09  2.43   320 x 320
> proj:    0 created: 2013040 220443  memo: RT GVAR
> type:GVAR     cal type:RAW
> offsets:  data=    3328 navigation=  256 calibration= 2816 auxiliary= 0
> doc length:   0   cal length:   0   lev length:   0 PREFIX=   0
> valcod:          0 zcor:  0 avg-smp: N
> start yyddd: 2011117  start time:201020  start scan:  382
> lcor: 4341  ecor: 12849  bytes per pixel: 2  ss:180
> Resolution Factors (base=1):   Line=    4.0   Element=    4.0
> IMGLIST: done
> 
> IMGCHA LOC/AREAS.3103 SS=180 STYPE=VISR
> Directory Block Changes
> Changed Sensor Source from: 180  to: 180
> Changed Source Type from: PRD  to: VISR
> Transferring AREA data outbound, bytes= 816
> Image file directory listing for:LOC/AREAS
> 
> Pos Satellite/         Date       Time      Center      Res (km) Image_Size
>     sensor                                 Lat  Lon    Lat   Lon
> --- -------------  ------------  --------  ---- ----  ----- ----- ------------
> 3103  G-13 IMG      27 APR 11117  20:10:00    30   86
> Band: 1    No Information Available                  5.09  2.43   320 x 320
> proj:    0 created: 2013040 220443  memo: RT GVAR
> type:VISR     cal type:BRIT
> offsets:  data=    3328 navigation=  256 calibration= 2816 auxiliary= 0
> doc length:   0   cal length:   0   lev length:   0 PREFIX=   0
> valcod:          0 zcor:  0 avg-smp: N
> start yyddd: 2011117  start time:201020  start scan:  382
> lcor: 4341  ecor: 12849  bytes per pixel: 1  ss:180
> Resolution Factors (base=1):   Line=    4.0   Element=    4.0
> IMGLIST: done
> 
> IMGCHA: Done ...
> 
> IMGDISP LOC/AREAS.3103 2 LAT=30:00 86:30 EU=CLW3 SF=YES
> Beginning Image Data transfer, bytes= 143056
> IMGDISP: loaded frame  2
> IMGDISP: done
> MAP LALO -2 INT=5 5 STATE=MS AL GA COUNTY=ALL MCOL=8
> EU: Restoring CLW3.ET to frame(s)= 2
> EU: Done
> MAP: Completed frame 2
> MAP FILE=OUTVHPOL MCOL=1
> MAP: Completed frame 2
> 
> IMGFILT LOC/AREAS.3103 LOC/AREAS.3003 LOC/AREAS.3203 FILTER=DIS 0 255 185 255 > 0
> Beginning Image Data transfer, bytes= 105776
> Beginning Image Data transfer, bytes= 106048
> Minimum value of BRIT = 0
> Maximum value of BRIT = 255
> Apply DISCRIMINATE Filter
> Apply CONTRAST Stretch
> Transferring AREA data outbound, bytes= 311248
> IMGFILT: done
> 
> IMGCHA LOC/AREAS.3203 OTYPE=PRD
> Directory Block Changes
> Changed Original Source Type from:       to: PRD
> Transferring AREA data outbound, bytes= 1056
> Image file directory listing for:LOC/AREAS
> 
> Pos Satellite/         Date       Time      Center      Res (km) Image_Size
>     sensor                                 Lat  Lon    Lat   Lon
> --- -------------  ------------  --------  ---- ----  ----- ----- ------------
> 3203  G-13 IMG      27 APR 11117  20:10:00    26   82
> Band: 1    No Information Available                  4.80  2.37   480 x 640
> proj:    0 created: 2013040 221220  memo: RT GVAR
> type:PRD      cal type:BRIT
> offsets:  data=    3328 navigation=  256 calibration= 2816 auxiliary= 0
> doc length:   0   cal length:   0   lev length:   0 PREFIX=   0
> valcod:          0 zcor:  0 avg-smp: N
> start yyddd: 2011117  start time:201020  start scan:  382
> lcor: 4341  ecor: 12849  bytes per pixel: 1  ss:180
> Resolution Factors (base=1):   Line=    4.0   Element=    4.0
> IMGLIST: done
> 
> Current Calibration Block Segments In File:
> Output         Input
> Product  Product  Brit  Brit
> Num Type low      high     low   high  Units Scale Missing Val
> --- ---- -------- -------- ----- ----- ----- ----- -----------
> 1: RICE    -1095      515     0   255  ERAD    10 -2139062144
> 
> IMGDISP LOC/AREAS.3203 3 LAT=30:00 86:30 EU=CLW3 SF=YES
> Beginning Image Data transfer, bytes= 185088
> IMGDISP: loaded frame  3
> IMGDISP: done
> MAP LALO -2 INT=5 5 STATE=MS AL GA COUNTY=ALL MCOL=8
> EU: Restoring CLW3.ET to frame(s)= 3
> EU: Done
> MAP: Completed frame 3
> MAP FILE=OUTVHPOL MCOL=1
> MAP: Completed frame 3
> 
> IMGPROBE LIST POINT MODE=N LATLON=*34.494762 88.234978 (Upper left of box)*
> Image Name           Day      Nominal Time   Scan Time    Band
> ----------------      -------    ------------   ---------    ----
> LOC/AREAS.3203      27 Apr 11117   20:10:00      MISSING        1
> File     Nominal  Image     RAW         BRIT       RICE
> Lat/Lon      Line/Element  Line/Element
> ERAD
> 34:28:44 / 88:16:04    65/  119    4601/13325            226         226 33.2
> IMGPROBE: Done
> 
> IMGPROBE LIST POINT MODE=N LATLON=*34.509640 87.485329 (Upper Right of box) *
> 
> Image Name           Day      Nominal Time   Scan Time    Band
> ----------------      -------    ------------   ---------    ----
> LOC/AREAS.3203      27 Apr 11117   20:10:00      MISSING        1
> 
> File     Nominal  Image     RAW         BRIT    RICE
> Lat/Lon      Line/Element  Line/Element
> ERAD
> 34:30:36 / 87:29:44    64/  148    4597/13441            229         229 35.1
> IMGPROBE: Done
> 
> IMGPROBE LIST POINT MODE=N LATLON=*34.129467 88.284866 (Lower left of box)*
> Image Name           Day      Nominal Time   Scan Time    Band
> ----------------      -------    ------------   ---------    ----
> LOC/AREAS.3203      27 Apr 11117   20:10:00      MISSING        1
> File     Nominal  Image     RAW         BRIT    RICE
> 
> Lat/Lon      Line/Element  Line/Element
> ERAD
> 34:08:18 / 88:18:40    72/  115    4629/13309            229         229 35.1
> IMGPROBE: Done
> 
> IMGPROBE LIST POINT MODE=N LATLON=*34.121307 87.521309 (Lower right of box) *
> Image Name           Day      Nominal Time   Scan Time    Band
> ----------------      -------    ------------   ---------    ----
> LOC/AREAS.3203      27 Apr 11117   20:10:00      MISSING        1
> File     Nominal  Image     RAW         BRIT    RICE
> Lat/Lon      Line/Element  Line/Element
> ERAD
> 34:07:13 / 87:32:02    72/  144    4629/13425            233         233
> 37.6
> IMGPROBE: Done
> 
> IMGPROBE STAT BOX RICE OUTLINE=YES 6 MODE=N* LINELE=65 119 F* PLACE=ULEFT 
> SIZE=8 33 EXCLUDE=-109.57 -6.77


re:
> *My commentary: The position of the box using the file number

File number?  The 'F' in LINELE=65 119 F stands for file coordinates.  Are
you sure that this is what you want?

re:
> place initially at the lower left. I took the difference of the range of
> lines and elements which came out to the following*
> 
> *Line: 72-64 = 8 *
> *Elements: 148-115=33  *
> 
> *Hence the SIZE=8 33.  *

Whatever works...

re:
> Name:LOC/AREAS.3203
> 
> 
> Calb  Pixels  Minimum Maximum  Mean    Stn. Dev.
> 
> 
> ---- -------- ------- ------- ------- -----------
> 
> 
> RICE      264  -109.5    42.7    31.5      16.082


> *Now I dont understand why the minimum value is -109.5. When I am
> just interested in the range of effective radius of ice ranging from 6 .77
> micrometers to 56 micrometers.

Given that you are working in File coordinates, it may be that you are
not sampling the portion of the image that you think you are sampling.

re:
> Does this affect the "mean" statistic?*

Does what affect the "mean" statistics?

re:
> *
> *I have another question. Now I was working on finding the maximum and
> minimum effective radius of ice location at 20:24 for the
> nominal satellite time 20:25 but Bob suggested to use the parallax command
> where I can use the cloud top height  from a air sounding and subtract 2 KM
> from that height to find the secondary point while using cloud to ground
> command to find the point. I use the lat and long max and min effective
> radius of ice locations in the lat lon command of the PLAX command. However
> I was able to get an accurate minimum effective radius of ice. However for
> the maximum effective radius of ice it wasn't located at the maximum
> reflectivity on radar. However with the lag time of the nominal time of
> the satellite you suggested...where I went to the next radar scan
> 20:28...the maximum effective radius of ice location strangely was at the
> debris ball which does have the highest reflectivity out of the entire
> supercell. So I am confused why it is showing this. If you want me to post
> pictures I can do that...at your request.*

I don't know what you are asking (and I'm not too sure of what you described
above).

re:
> *
> *One last question...is there any way I can find the actual scan time of
> where the region of the supercell storm of interest is scanned?

The GVARINFO command will list the time of a GVAR image as long as the
image is in its original form (meaning that it still contains GVAR
line prefixes).  The effective radius images do not have this information
as they are products,  but the IR and VIS images you have do.  Unfortunately,
GVARINFO does not understand ADDE datasets; it is an old command that
only works with AREA files directly:

HELP GVARINFO
GVARINFO -- Displays information from a GVAR area line
   GVARINFO area_number line_number
Parameters:
   area_number | area number containing GVAR data
   line_number | line number in area coordinates; use T for table of
                 calibration changes
Keyword:
   BUG=201     | Run program in debug mode                       (def=101)
----------

This can be overcome in a couple of different ways:

- link the original image to an AREA file name

  Example:

cd $MCDATA
rm -f AREA2000 AREA2001
ln -s archive/27apr2011/2011117201000i13.c01 AREA2000            <- GVAR VIS 
image
ln -s archive/27apr2011/goes13.2011.117.201020.BAND_04 AREA2001  <- GVAR IR 
image

re:
> Is it in the IMGLIST FORM=ALL command?*

No.

re:
> *Thanks for your help! Looking forward to hear you response.*

No worries.

Cheers,

Tom
--
****************************************************************************
Unidata User Support                                    UCAR Unidata Program
(303) 497-8642                                                 P.O. Box 3000
address@hidden                                   Boulder, CO 80307
----------------------------------------------------------------------------
Unidata HomePage                       http://www.unidata.ucar.edu
****************************************************************************


Ticket Details
===================
Ticket ID: BFX-595521
Department: Support McIDAS
Priority: Normal
Status: Closed