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

20020304: 20020301: GPFRONT and Satellite Overlay



Patrick,

Are you running the script from cron? If so, most of your
environmental variables won't be in your cron environment.
You aren't source'ing Gemenviron in your script. That would
be a big problem is running from cron, or an environment that 
doesn't have Gemenviron sourced (such as a cgi script).

Check your $GEMEXE environmental variable and $PATH.

Steve Chiswell
Unidata User SUpport


>From: "Patrick O'Reilly" <patrick.oreilly@xxxxxxx>
>Organization: UCAR/Unidata
>Keywords: 200203041755.g24HtnK29806

>This is a multi-part message in MIME format.
>
>------=_NextPart_000_0018_01C1C373.809E2DF0
>Content-Type: text/plain;
>       charset="iso-8859-1"
>Content-Transfer-Encoding: 7bit
>
>Steve,
>
>I am running version 5.6.e.1.  I am not sure why I can get the image I want
>when I run GEMPAK interactively, but get different results when running a
>script.  I can plot satellite imagery fine by itself, and I can plot fronts
>and surface data/radar fine too with scripts.  If I write the script to plot
>an IR satellite image with gpmap, then to plot fronts using gpfront, I get
>the "image file missing" error, and get a satellite image, but no fronts.
>If I write the script to just run gpfront and specify a satellite image, I
>get fronts, maybe a colorbar from the satellite image, but no satellite
>image, and no error messages.  Finally, like I said, when I run gpfront
>interactively from the command line, I get what I want, IR satellite with
>fronts on a satellite projection.  I am attaching a script and the image
>created, maybe you can see where I am going wrong.  Thanks.
>
>Patrick
>
>----- Original Message -----
>From: "Unidata Support" <support@xxxxxxxxxxxxxxxx>
>To: "Patrick O'Reilly" <patrick.oreilly@xxxxxxx>
>Cc: "Unidata Support" <support-gempak@xxxxxxxxxxxxxxxx>
>Sent: Monday, March 04, 2002 11:14 AM
>Subject: 20020301: GPFRONT and Satellite Overlay
>
>
>>
>> Patrick,
>>
>> It sounds like some memory might be getting stepped on
>> if the satellite file is reported "does not exist",
>> or the projection read from the satellite image is not
>> as expected. Either are possible if the program was linked
>> against a different version of the libraries.
>>
>> What version of GEMPAK are you running? Type "version"
>> at the gpfront prompt to make sure it agrees with your
>> installed distribution.
>>
>> Steve Chiswell
>> Unidata User Support
>>
>>
>>
>>
>> >From: "Patrick O'Reilly" <patrick.oreilly@xxxxxxx>
>> >Organization: UCAR/Unidata
>> >Keywords: 200203012258.g21MwLx25795
>>
>> >This is a multi-part message in MIME format.
>> >
>> >------=_NextPart_000_0031_01C1C142.4BAF3360
>> >Content-Type: text/plain;
>> > charset="iso-8859-1"
>> >Content-Transfer-Encoding: quoted-printable
>> >
>> >Hi again,
>> >
>> >I have been fighting with gpfront to get an overlay of IR satellite and =
>> >fronts.  I already have a script running successfully creating an image =
>> >of the 6km composite radar and fronts.  When I am physically sitting at =
>> >the computer, entering commands and using gpfront, I can get the desired
>=
>> >image.  But when I run gpfront in a script, no dice.  I have tried two =
>> >approaches  (with many, many variations).
>> >
>> >I have tried running gpmap to create the satellite image, then running =
>> >gpfront to overlay the fronts, but no dice that way.  I get the =
>> >satellite image created, but get an error ( [IM -1] Image file ... does =
>> >not exist)   from gpfront and no fronts are drawn. =20
>> >
>> >The second approach is simply running gpfront and specifying a satellite
>=
>> >file.  This method produces no  error messages, but strange plots.  One =
>> >such plot is attached.  It looks like it plotted something in regards to
>=
>> >the satellite file, as the IR colorbar is there.  The clip from the =
>> >script:
>> >
>> ># Plot the Fronts/Highs/Lows
>> >
>> >gpfront<< GPFRONT_FLAG
>> >satfil=3D$GEMDATA/images/sat/GOES-8/4km/IR/$satfil
>> >map=3D31
>> >garea=3D23.0;-120.0;47.0;-65.0
>> >proj=3Dsat//1;1;1;1
>> >radfil=3D
>> >latlon=3D0
>> >panel=3D0
>> >text=3D1.5/22/HW
>> >clear=3Dyes
>> >title=3D32/-1/Fronts at $filetime
>> >asfil=3D$GEMDATA/nwx/hpc/fronts/$frontfile
>> >device=3Dgf|$tmpgif|900;650
>> >ashr=3D
>> >asatt=3D1/N
>> >
>> >run
>> >
>> >exit
>> >
>> >GPFRONT_FLAG
>> >
>> >I can see nothing wrong with this, if you do let me know.  There aren't =
>> >errors finding satellite or front files, and all the variables work =
>> >fine.  I can't understand why it works sitting entering commands running
>=
>> >gempak interactively, but when I try the SAME thing in a script, I get a
>=
>> >hard time.  I found very little in the support email archives about =
>> >gpfront.
>> >
>> >Thanks once again for your help!  I am trying to get some requested =
>> >imagery put together for posting to the web, and have about 95% of it =
>> >done.=20
>> >
>> >Patrick
>> >
>> >~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
>> >Patrick O'Reilly                               Support Scientist
>> >The STORM Project            patrick.oreilly@xxxxxxx
>> >208 Latham Hall                             ph: 319-273-3789
>> >University of Northern Iowa
>> >Cedar Falls, IA 50614            =20
>> >~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~  =20
>> >
>> >
>> >------=_NextPart_000_0031_01C1C142.4BAF3360
>> >Content-Type: text/html;
>> > charset="iso-8859-1"
>> >Content-Transfer-Encoding: quoted-printable
>> >
>> ><!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 Transitional//EN">
>> ><HTML><HEAD>
>> ><META http-equiv=3DContent-Type content=3D"text/html; =
>> >charset=3Diso-8859-1">
>> ><META content=3D"MSHTML 5.50.4522.1800" name=3DGENERATOR>
>> ><STYLE></STYLE>
>> ></HEAD>
>> ><BODY bgColor=3D#efefef>
>> ><DIV><FONT face=3D"MS Sans Serif" size=3D2>Hi again,</FONT></DIV>
>> ><DIV><FONT face=3D"MS Sans Serif" size=3D2></FONT>&nbsp;</DIV>
>> ><DIV><FONT face=3D"MS Sans Serif" size=3D2>I have been fighting with =
>> >gpfront to get=20
>> >an overlay of IR satellite and fronts.&nbsp; I already have&nbsp;a =
>> >script=20
>> >running successfully creating an image&nbsp;of the 6km composite radar =
>> >and=20
>> >fronts.&nbsp; When I am physically sitting at the computer, entering =
>> >commands=20
>> >and using gpfront, I can get&nbsp;the desired image.&nbsp; But when I =
>> >run=20
>> >gpfront in a script, no dice.&nbsp; I have tried two approaches&nbsp; =
>> >(with=20
>> >many, many variations).</FONT></DIV>
>> ><DIV><FONT face=3D"MS Sans Serif" size=3D2></FONT>&nbsp;</DIV>
>> ><DIV><FONT face=3D"MS Sans Serif" size=3D2>I have tried running gpmap to
>=
>> >create the=20
>> >satellite image, then running gpfront to overlay the fronts, but no dice
>=
>> >that=20
>> >way.&nbsp; I get the satellite image created, but get an error ( [IM -1]
>=
>> >Image=20
>> >file ... does not exist)&nbsp;&nbsp;&nbsp;from gpfront and no fronts are
>=
>> >
> >drawn.&nbsp; </FONT></DIV>
>> ><DIV><FONT face=3D"MS Sans Serif" size=3D2></FONT>&nbsp;</DIV>
>> ><DIV><FONT face=3D"MS Sans Serif" size=3D2>The second approach is simply
>=
>> >running=20
>> >gpfront and specifying a satellite file.&nbsp; This method produces =
>> >no&nbsp;=20
>> >error messages, but strange plots.&nbsp; One such plot is =
>> >attached.&nbsp; It=20
>> >looks like it plotted something in regards to the satellite file, as the
>=
>> >IR=20
>> >colorbar is there.&nbsp; The clip from the script:</FONT></DIV>
>> ><DIV><FONT face=3D"MS Sans Serif" size=3D2><BR># Plot the=20
>> >Fronts/Highs/Lows</FONT></DIV>
>> ><DIV>&nbsp;</DIV>
>> ><DIV><FONT face=3D"MS Sans Serif" size=3D2>gpfront&lt;&lt;=20
>>
>>GPFRONT_FLAG<BR>satfil=3D$GEMDATA/images/sat/GOES-8/4km/IR/$satfil<BR>map=
>>
>>=3D31<BR>garea=3D23.0;-120.0;47.0;-65.0<BR>proj=3Dsat//1;1;1;1<BR>radfil=3D
>=
>>
>><BR>latlon=3D0<BR>panel=3D0<BR>text=3D1.5/22/HW<BR>clear=3Dyes<BR>title=3D=
>> >32/-1/Fronts=20
>> >at=20
>>
>>$filetime<BR>asfil=3D$GEMDATA/nwx/hpc/fronts/$frontfile<BR>device=3Dgf|$t=
>> >mpgif|900;650<BR>ashr=3D<BR>asatt=3D1/N</FONT></DIV>
>> ><DIV>&nbsp;</DIV>
>> ><DIV><FONT face=3D"MS Sans Serif" size=3D2>run</FONT></DIV>
>> ><DIV>&nbsp;</DIV>
>> ><DIV><FONT face=3D"MS Sans Serif" size=3D2>exit</FONT></DIV>
>> ><DIV>&nbsp;</DIV>
>> ><DIV><FONT face=3D"MS Sans Serif" size=3D2>GPFRONT_FLAG<BR></FONT></DIV>
>> ><DIV><FONT face=3D"MS Sans Serif" size=3D2>I can see nothing wrong with =
>> >this, if you=20
>> >do let me know.&nbsp; There aren't errors finding satellite or front =
>> >files, and=20
>> >all the variables work fine.&nbsp; I can't understand why it works =
>> >sitting=20
>> >entering commands running gempak interactively, but when I try the SAME =
>> >thing in=20
>> >a script, I get&nbsp;a hard time.&nbsp; I found very little in the =
>> >support email=20
>> >archives about gpfront.</FONT></DIV>
>> ><DIV><FONT face=3D"MS Sans Serif" size=3D2></FONT>&nbsp;</DIV>
>> ><DIV><FONT face=3D"MS Sans Serif" size=3D2>Thanks once again for your =
>> >help!&nbsp; I=20
>> >am trying to get some requested imagery put together for posting to the =
>> >web, and=20
>> >have about 95% of it done. </FONT></DIV>
>> ><DIV><FONT face=3D"MS Sans Serif" size=3D2></FONT>&nbsp;</DIV>
>> ><DIV><FONT face=3D"MS Sans Serif" size=3D2>Patrick</FONT></DIV>
>> ><DIV><FONT face=3D"MS Sans Serif" size=3D2>&nbsp;</DIV></FONT>
>> ><DIV><FONT face=3D"MS Sans Serif"=20
>> >size=3D2>~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~<BR>Patrick=20
>>
>>O'Reilly&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp=
>>
>>;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;=
>> >&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;=20
>> >Support Scientist<BR>The STORM=20
>>
>>Project&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;=
>> > <A=20
>>
>>href=3D"mailto:patrick.oreilly@xxxxxxx";>patrick.oreilly@xxxxxxx</A><BR>20=
>> >8 Latham=20
>>
>>Hall&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nb=
>>
>>sp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbs=
>> >p;&nbsp;&nbsp;&nbsp;&nbsp;=20
>> >ph: 319-273-3789<BR>University of Northern Iowa<BR>Cedar Falls, IA=20
>>
>>50614&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&n=
>> >bsp;=20
>> ><BR>~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~&nbsp;&nbsp;=20
>> ><BR></FONT></DIV></BODY></HTML>
>> >
>> >------=_NextPart_000_0031_01C1C142.4BAF3360--
>> >
>>
>>
>****************************************************************************
><
>> Unidata User Support                                    UCAR Unidata
>> (303)497-8644                                                  P.O. Box
>> support@xxxxxxxxxxxxxxxx                                   Boulder, CO
>> --------------------------------------------------------------------------
>> Unidata WWW Service                        http://www.unidata.ucar.edu/
><
>>
>****************************************************************************
><
>
>------=_NextPart_000_0018_01C1C373.809E2DF0
>Content-Type: image/gif;
>       name="satfronts.gif"
>Content-Transfer-Encoding: base64
>Content-Disposition: attachment;
>       filename="satfronts.gif"
>
>R0lGODdhhAOKAvYAAP///wAAAP/k3P8AAAD/AAAA////AAD///8A/4tHJv+CR/+lT/+uuf9qau4s
>LIsAAM0AAO5AAP9/AM2FAP/XAO7uAH//AADNAACLABBOix6Q/wCy7gDu7olozZEs7osAixQUFBgY
>GCAgICgoKDExMTk5OUFBQUlJSVFRUVlZWV1dXWVlZW1tbXV1dX19fYaGho6OjpaWlp6enqampqqq
>qrKysrq6usLCwsrKytLS0tvb22VlAHl5AI6OAKamALq6ANLSAOfnAPv7AAAAeQAoogB5ygC+81UA
>AKoAAP9RUf+Ojnn/ebb/tgAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA
>AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA
>AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAACwA
>AAAAhAOKAgAH/oABgoOEhYaHiImKi4yNjo+QkZKTlJWWl5iZmpucnZ6foKGio6SlpqeoqaqrrK2u
>r7CxsrO0tba3uLm6u7y9vr/AwcLDxMXGx8jJysvMzc7P0NHS09TV1tfY2drb3N3e3+Dh4uPk5ebn
>6Onq6+zt7u/w8fLz9PX29/j5+vv8/f7/AAMKHEiwoMGDCBMqXMiwocOHECNKnEixosWLGDNq3Mix
>o8ePIEOKHEmypMmTKFOqXMmypcuXMGPKnEmzps2bOHPq3Mmzp8+fQIMKHUq0qNGjSJMqXcq0qdOn
>UKNKnUq1qtWrSwVo1VpoK1dCXgV09bpoK9azaI+G/RpgrVhB/m4HxT00N63duznNtv2qty/fv2/1
>GjIrGK/hwy8FE2a7+O1esYrZymUsGbHlyyQlLwbLVXNnx48HUwaNubTpjZ4hg948WTVn0qHhVj5N
>uzY/t4XpAlYduLHsz71nR4Ztu7hxeQIAKF+unPfa32Ghk30cnfr0162Pa9/+Ljnz5c6rU9/b2jV0
>8ueFj+bOvr053LPLxs/9Wzfs4e7z6+/m/TsA4tjF5ldw9Qm4Xmj07afggs/AB2B50+Em3YDiUZgg
>gxhmWEx/3z0I4Vi5iScdiAFqaOKJwziI4oos6sMhcx62KOOM57wIHo045oiOijr26KM3Njb345BE
>WsNjkUgm/rlMkP8p6eSTxBwJ5ZRU4sJkjFVmqeUpUm5ZiQErFgDLAFBeSYx/aKap5ppstunmm3DG
>KeecdNZp550GGHDnnnz26eefgAYqaJ0FFLpmoYgOCsAAAyjq6KOQRirppHE+YuYwlGaq6aacKpqn
>np2GKuqopNqJqKFsJjooo42W6uqrsMbKnKXw0XrffPFhIuuuvPb65qe+BivssG6eWoCrrBKr7LLM
>fmdpmlieB+Fo12nS7LXYZqvtttx+Z2ysyXYr7riSNtKlfYX5Vt+FlZDr7qatvivvvKUai6qsjNKr
>775tPosmls8ViOCB1vJr8MEIJ/yqvcf2Gq7CEJNra12M/uC3LsGZRKxxmvlu7PHHzH1L7MMglyys
>v/5FKzB+FhdscsKsdvzyzPoy3CzJNOdMKsod2prdxY6xO4nO+8bMaMPLxUv00srajK3MTEdN6cTP
>1VoiyxjrKjW5MQOgqnJKby22q05ri/PYaPvJM4zWVRvb2+q+nXHaZnd8LKpQ0623pAx/ve3Zewf+
>JtUiyhf0ei3PLfi1SOcb9uKQ+9n3qeMCHvnlyq19oyPpAped0ENjTizSiz4u+unFTi5y5aajvjjh
>bivSeYUBb+L6snnfrvtyqpfN+u6Xay4k5+oRBzolwMNKOti5J49537xD/67lzosN+/G4VM9r89oH
>Lr23/r5zzX33TAvfJKbkL3xv+mlP3qa9+iZrNPtEm6+yLfTXu3z+UX+f+vrzMlrX+Pey6+WKFwQM
>ld8SODP/zQmAARQg9RiIMPudiYKZWiAGPea+sUlwfBvklwHvV4sQ8g2CJjxYB/f2wRRCzILoc+Gg
>NChDeq1wcRKsYQWJJ6EY6rBPlPuhvG4ouvkJUYScg9YFj3gnFDJRW0REnRGfOC8YCoOKdHIiFkfX
>Oy1iboBbfNcIzRU77HEijHDaHxp91UX4OW+Ca7yWpZhAxzrSEWC1o9go4pgqr6mRj2RzWhR3B0dA
>EmuMsgPM5w7oCUP6p2E0dCSleocmB76xdZI8JOfs/mjHGGHtcCS0RCYf+cdRKqqNh1pd90Boyl4h
>MhGpgRsjz9hKP5aylqaipJxUuUpW4hJWc+TkHeWDncSJopaR/GWcUJlLLyavkMp0FRl7CMvdyDKU
>7TJlMqMJPmYCypmX5KavginM8EQoLsYMRSvBmck2qs5Rt9SeL8UZKnJy0pzDcU06QaFNdj7RnQDl
>pSmhSc9NvZKH1zTFKLeZNDDyL6CWjN4gB4rJgtZzk8JkgicVCbRTtJOhzCMo5CDaxTROlKIWjdVB
>0XUdPR7TkUFs0xRdF9CQ6fJNJ6VoRVOqKXt2smKxk1spDClQNc30dO6s5DuX6caCzpOn5eKhS30B
>/shE+bN0It3bTbvZVKbG1Kk7heqkLLWEspq1rNiUBR9VddWjom6pSu2qV686yqyKdVAr/UUc7wXS
>htqVbhG1aVFx2ldcPvWuiiLrWc2a1liscXmFxeph0fTXhIUvrpGNK2Ini9hAKXaxS2gsLNCoxsK6
>9U0O1Zlc0/RVQsWTm5XtLJ/ySlVJDpY58+vd4077stWylq6yZRNngztbzhHguMg9rmhfsdDbSnYA
>XexaDpfW2lS9lrio5S12/0TbXnwUoEZb4Qdji7DMem27fZouegVlqeQmd7musC1bddu1y6qXacBd
>b3rvq98+dReBRN2mscLr3Oey6rpGHe6wENxf/kDxt8F3aq97lTvN4qUiwPG0F4FNu2F3Ki21ENbb
>g0NMJwlPGI+0K9xL9/pakW24iQF1nHZJvLUR07hSUq0dIiwWN1KstcUxJfCfHDjeG7NQgEa2k4nd
>u1FQJnSoj/3f8sI75MsaOKxJ1pmNs8ym/67McwLbI2kJu7/xmjmsVg6pgrkMMiSzWU5Lfi9Qq7PP
>TyjTiWcusnWri1sQvzlqM/4zmrz8oS/D9xG/1GKeW5hKkPpZ0EwjL5fjjFx8KnJAP/MxMg+16C3b
>0p9rhvTGQv3nCldNQp/MtJi1KdNOz5jPHMOyqF9G6klzTgm4zjWum7xIJ0MZ0uat9axhJuth/ltK
>17rmdYF6vOo/L5C3+R22xoSdZEK3DZRB5SewKfto80o7YtQ28rGRvWtTi+bQkVjopv6o3Wh/m9jF
>FrW1swdTd6/q0bZ8d6TjDWlLJeHfAP83uldRb3sLCtre1jfCwk3jed8Cw+Dys28V3mZ+C9rfAQf4
>wFXxY4Mryq2BpfioLV5qzmVc40sMI6xHBfKQi/zlJXM4/lSe8EkdFa7+IbkOQQVzBmL85Bu/8BZX
>XqqZQq9Qu5V0CIHVcwyaWsW1/WfNMzVF3z3MzXH8VJ6a7nOTnzwJQUcFFYn+KiP69sx81DrXKSjz
>Eh6xwGUf4MrNvEat83ztCfx5xpusYzPa/u6HE99VbsHJ6DDaHe9OzzHUXxOwqWq7hoEXfMzainUq
>qh3xG9R7wEOJab9rDfJkxxe+E6z0/F0e84k319fBDglMh3nFKYz89ka/pkBv8PCoD6HMW+Z5UcZe
>9rOX9OTh7iaGlwz3uc885z54P97P8u8hTLOwXK3do0FUuLbfGtOTb0KHz07VsGegy4NFfdorJ8Zn
>87TYts993S9fgpZeTSxfr04KfmvqRS//PLc6XT3v7fTt535PRxZWI1R11kjih3+ip3/vg3PfsWiC
>A4ABqHzmwnzTRCL093j5E3o3o39oNn5qpn5Rs3UT6EKWYoGGYx9CFX7pw4EHc3eY5YJ+/mV+TMN+
>JUiBc+Z4JcJ4ved79KOA7kKCv0V8sWZ8EGODN5h6FYMETNiETIhi1RR2jLCBHvcuSChRClh6G3OF
>Sdh1A9iDatWCQGiFQshVVShZelOGXYiDZeGETiiFCqU9wKcxaoiFDqaFCcOFa0hAbUcLciiDEXOF
>Y5hzRrgvdbiHXriEbviEKac7gOgxdXiGCYY2h4iIfKh42QZgu/OIkKiGDFYn1qd9lWiJ+WMpi8iI
>PkRTkpiHMHhCn7iFrUiKl2gup4gEcPhrbzWICjOKj6KL9MKLspg+fTgLt8OJJqOHvfiKgRiLwUg/
>pniKt6hpz0OERIOMp7SK82KNzVg9/sMYhpczhzUIjNcoNeK4jcljKUeQjuqYjtGiYmDog4sjfZTI
>jJCCjdlIj+bIjZh4PHWhgwgYj+BIjuUYKL4oLwOZj7eDjuuojlCIHgMTjYYAkAEpippijKyIjwi5
>Owq5kEegbBi4ghrYPpRTkCCjjVVGkmSYkexDaAlygJ2gVV9jj8d4kDBGjXRIkyoZORu5kLzmfBBJ
>CHqzQMq4fjjpWjZ5kznZPSw5Fy5JS2JzlItjkjU5kcuIkUmJOZZyAFq5lVoZf+PxkOBXf1sDlVFZ
>lHNFlct4lc4zRrWSahloZ2Mpk2NjlmQGgmrZc1nJlVvpkXHzjpIgNRZ5OXQpZfJ4/pdNt5QcxWws
>2ECBKZhWaZQOaJh4l5d6eQAo1lI61myMiZJ0M5iNtlWSOZmcU5ldaW4fGYc005in45lDWDZyGZok
>1o2OlZqc+X+PuUur9ZrDkn2weTqUqZc/KZYgg5bAI5V1WUq6OX002JuiI5ujNZzEmTysKVgNOJTd
>wpvMGTyjSZrBGZIRU5imd5sEaZ3cgp3ZCTnOyVwcFJ3OM510Uogjs5zniZ6cYwT2eZ/22Z1w+Z3g
>yT/uOSfwuZvyOZ+CYyn4iZ/6+Y8KY5f+KZ53qHN/M6AEujfpGV+WNZK1iTn/CaAQWjd4OKFLY6AH
>mp+NWDMik6Ea6qAPCjHmCaJ0/oOYheOXyMMvvJScomOcgRKgvuJWCOCi2mkuI0qiiZRHmbmY7sJn
>5Mk+GyonMiOhHXg2CNCjPvo6+zhLF6KYwvkuqolBSxonMtaizDJTUTqmU1qgnFMEaJqmaOqRV/N8
>LqOlKNo9XQondEcvUzSmeFqmFFqlfOeTqEkuW7p0KvqglTc9A4SniKqnL3qmapqmfep6ILmf4kKW
>Ozen2SWCHtoxiJqoijo2ltKojpqC49GU0NctgepCllp8hScv87Opm9qp1gMAXwgcRQqWb6mgUBSn
>BJSqqgqmYdo1rvqqsBo1nwqqReCV5/Zk0ogt7MlEvNpqvoo7wBqsnDqsITqr/nwJZpGKq8vSrE/0
>rNjnpE9aftZKrBVjrMcqqufEGx7VLN5qeYNKqAbjgR1arvoCo5i5eJLKRafqrOBaezoacQxorzlT
>rKB6mRYmdk3zrlv0rwB7MPRKsAUrqzlYq961YP1pSg5LevPqgRJbsOZCBCI7siKboC8pLBiapFkX
>ryt6MMH6sVtjKSRLsibrlLwSUzZaqSzrJx8qLK4KszHLpyWqPDGpsoC0sQlmtMRSrUBrriE7syU7
>tPoDQEqbdjsrOf3KK8LatE5bsfq6C7vSWlVrtaOSsdfys1xbPpwDtVGbimWbtaOEo4JSUuLysml7
>rdgqtZxCqfSEtKQEmtqC/rZ3Wz+cMwSGe7iGW7OlqkC6akhye5I3lbOKYreDqzOWgriIq7hv2ilw
>q0x+K1hw1bmSQrmVSzNtFyJuarN7K7lx+7nnF5l8yymkW7ozc7mYm7gTMy0Jyrisq7GfC7h2uLSC
>S7u1K7Smaau4OEmi27d+m1Of1rt2MrvEazK2e7vtSGdZk6WuCL2tJIHJCHxmK7tbO72ma7zEBKky
>CgkZFLvE9bjN1FfAKyrSS77UW7i3OwRQiL6p+3mRwrDB5b1zS4TeJL/DS7/lm7fEZCC+tqzf27jM
>C66AeH2TQq15asCEay73i7/XBmb6q5/967/b5b5edZYk5UwkpRwUXMEW/gyyGHy/XtmW2eud48m+
>DQbAkvNAoFnCJQwAKRylK3zBXpuJE9KX+3sJ35tvImfDq5tTNeVOPazCP3zAFSMEVFzFVExCbpm8
>8DS2NSzCAQyeSfW6bfTEPhzFLBzEnpdP6au+fINBzkXDm+LFkOuCA+xhA/DEZmy5nGPFVozFlbHG
>6da/G/TGywsprvu3Mki3WOg+rJLCeazHCHxFBUXIDgyn5gW4JRUzFPzIkDzFfHzFeotMRAfH/Blt
>kSnGTmM01MrJZ1wWnwzKbntno1zI+1JzJyVeqjy+rFy/kRwMFkXJ3DsuSepy35NDTLvLvFwxQbDM
>zLzMmqs44gTMXDw2/oVZNi10zMgMMhXaCr+8w2IYel2lXmSazcVrLs3czM/Mv9ykw+wjwH6zqlJK
>zgVkvr48yd5MPom8Pqsqz60sO+fMzOlsxPZ8wt/MVA8Yrfy8Q70MDN1M0Pj8ilMWsAktMZzzz84c
>yut0z0/UsxN9rxVt0QENj+s8y5V8ORzd0fQyb+5YxNA80gJGy7dz0igtL5Zi0UGAsNQCyI7Q0F5E
>ytWD0DOt0GhspRyFvAwcTdLMRIUa1P0MSzadrYb2py7d0zBdREDN1Ei00Ooaw/uK1CQdzEd21Vjt
>0eYCBGZ91matMvkaljKcaF89za7jf2MtxWWB1mit1uvK1l0tyy9d/tJ0I9dzXc5DjayMp6ya6dV9
>DdZoA9iBncx1bddpvcHGo606PYUDTdV+vdj73NiOPaT+mFCkurmczaFLPdoxxzk/kNqqndrKlteV
>bdmmTdqlHdseY239aLF7TdsAi6m6LdRlsdqrfb308dqN0NuXatx4O9gYbdziitwHYynArdohnU3O
>zW31Wt0pjdrR/QPTPaPYfdAS/d3Nss2sIN591tzmLUbaHd3dHTrmzdvpvS/kTXDvDd/xnd3mst3c
>vdyxvdn3rTDzzXH1Ldb/vS3Qvd3t/Zfi7d8F7tuyo98JHsjfzeANbjABLnQTPtsVbuGc4wMe/uEe
>HuFsjN32veHu/nLhCkviBG7i2GIpIA7iIo5oGY7eLG7g9MzQKk7jNa4tltIDPv7jPh7jO01iSQ05
>Or7jPH7jekXkbx054Y3k48Q5QA7kQl7cTJ7YcA1u1w3l2dLjUx7k/C1bRY5DW87lcqTkUQdhYy44
>T27msuLlX17lsK3mTW7kZe7my4Li7Rpiax44bY7nrwLnUy7ni0Bj7GzngJ4wei7V/XXoZL7iic4r
>gk7lYd5Zjs7mJR7pyrLoWqxfl+7nma7pwzLpP07oimDoGg3qGi7qXS7lX94Dpp4IqO7QmL7qrH7m
>Wp3mDdbnYQ3pt14qlsIDwj7swh7riDDrmK3Y7kLhvz7eaM4I/mIiJtx85cme5Swa6s2+K8FO7MOO
>TYiinnxe54/u69l+UbmuCN8O7nSO5U5u6+WuSebC7d0+Temu19rr6eJe60f+7jvDOfJe7HOGKNhG
>3KdO7bfk0x7j7vx+Ms/eGgJ/aSwt0o2e76r+5wvvKJayAxq/8Rq/Ud9OxB5s8J9p7Qqz7xc/KtaW
>7iDP6PjO7nZ+5ydv7hXD8RyPJYWybFzNrb1t8TGfWA0vCA/fwSyv2zzf8+zFOTS/8fEX9ARD8Mee
>40bf6nmLG0z/kU5/CAsO81GfKRmf9DsAIKcy3Dl/siq+9bg+814P9mJypdp61Myt9WYfVbn+8T30
>2aob93hf/mJIn/aGoPLUtK1tnfeC3y83Xu+5MPiIj2PmogOM3/iMTxo3r4mJP/lqclCGf/iUn/mD
>xjmO7/iOcfmYr/minzmc3/mPPwiRf7GjL/qIBPqhv/qZbymmf/qxDPuDz+lub/uCbyk50Pu+3/vG
>jvW6T/m4f9jDn/e8//u+H/wRefyJX/xG6vxbn/zKnwPMXwjSj/jUr/zXD5TZv/s/rwvfD/7mUv3A
>X+nSRu7jj5Xh//pcp/7rr5OcY/7Wj/7DBv/xT6XnLvnvb/L5z/6AECA4SDgogIOYiChQ2Oj4CBkp
>OQlQaXmJmam5ydnp+QkaKjpKWmp6ioo5sLqa6voKGys7/ktba3uLm6u7y9vrSzppqKjIGGx8jAz5
>u8zc7BzL2vo8TV1tfY2drb3NLXosAB4unkxebtyNnq6eGT2w/g4fLz9PX5/+PZxYbM7fP2gPMCCt
>dgILGjyIMKHCat/EOfQHkd/CiRMJUryIMaPGjQHx5cOxL6LIcxxLwrNoMqXKlSxbomroMNzImcFc
>2qSG8qbOnTx7Gvx2I6jQoCFpGi3kMymunEqbOn0KlRfMmEWPWo2K1RTTrFy7ev0KAOhQoVWtGgWL
>VlW0tGzbumUpduyNsmZnvgW79a7evXzpTY1Zt27frHkHGz6M+FncsXQDR0zstDDkyZQrv/r70PFR
>yz1z/hrgDDq0aE2LhzZuJLNQas2CRttsBzu2NNe0a7vFPC5YbkFUWQewvVK2cFbAixvPWprsJMCG
>wAVY7fh4yeHDpVu/vjM50eWZnzvn/T06dozUhY8/j54jbuiSVrsPHzh9xfKx5du/f/Cbjf38959W
>Hd57/42EH0YGfFZgggoWtB587QUI4YAiLbjQgRReiOE7+vXHHyNUOShgcxI+lqFBFpaIYooMGfOh
>cy0CuE+ImqkY0Ik03oijLhtyaMOI3sUYIWs50mPjkEYemUqDPkInE3uCIbnOgQhCSWWVnnxTQ5Za
>ZrkkfL0JaWU3RYZJZphYbqlll1U5GV+Z2ozpZpxG/irpG01yYgPnnXqmeCaaNfhYZ017UpPnoIYu
>SGegEx7qTKGMPipfn2gCqmgkkC4j5aWa3hdAAbq9WClEm/bi6KimAleAp8v5yWWoop6aS6mwzgpa
>qqq2x+qfrvpD6y2y9gosYrayCOqu5gSLbLIL2norrqxSaiwh5zFbQCfUKovtodQS+2G0x453LSfh
>ZktunNQ2G4kANKzL7rrQevsbuMxaO2+59pJ5LrdfwovMtPVuMu69Ah95LrqQqNsuu+9666+t9Do8
>cMQ5FmzwIy0uHG3DqT68sTxTSgxyUxRX7AjCCdOAsbEaVyvuv+/8GnLMKo1MMmonu8tvMivTHHCU
>/h/LDHRLPOvLXM4kYcfzyB7/HHTTJSXN4gxSTy11yrvunHTHL8PsdNcHZU10d0ZTIm/WPafDtddq
>AwQxABRHTfXUVruKddbxpL123vG0jQmzYe82tiR1J30303ofXg/fm3wTt9yBC4q0y5mcjQ7eiF+e
>jeKL/83m40iVrbUmlItpOOamd6P55rrJwHrrrM8d6socs+xz6afffk3qqi/nuuuwVyp7y7pvYznu
>xvMyPGmcO+h5I8EDLDnpth9PvS/RX8li763/rujzol9P/PTVj59L8pwk2vwj3k8O/pvikw//LOaf
>n732MnAf6Pp9t49n8fH/P4r50e9T3Uqf+kBH/rvvCZBQ7wOgA0uxwN21x373M+ABI5e80WXDfw/s
>IPsSaIpvUBB/depapjyIQlNEcIDcKeCPLBivpnEwhR1cIQsnaD8PFc1zTjshDX/oCRveMF0vahLz
>xtbDBgIRhfwDBotiAMUoQnE3neNXEpeIxf2F7jLLo+IRcybDGWaReluExTekKEUvkvAsQRPjGHHX
>xBB2UUY8bKMS30g+Db7kiWicIh0fZ8c74pGMejzFGfsYAzWmz46DTCHFaqGkPwYOaD5sJABpBkkW
>wWCTnNykDj9pQEoK0pKXIxwt6LSvOsbMjaR0mt1OqclOchJIa7zKKlnZyqC9chbog+E/QDHK/kPh
>MpcxMxsIzRhLWcKgljP6xDDL9ExiDsyYZeSibpTpSV9+zhPRtFIlpWk6albTFb3UpjO7SaVvghNx
>4hTiJ76BzWVqU1rcDOae1LlOvbXTndgjYCrNWc98ClRF9RJnJq+JTWaKpxPoHKhDa6U4Yx60hf/0
>5TkfitEEZZBnE23PCz4K0o8qtE0MtWdGTwrRcWqxYB0logvnWU+TonSmkInjSvnZiW+ENKQjfRIn
>GkrToIKlkLNr6cF2CtKemiWgvZKpUEdFVFCodI/+3CFAf+pUaGb1qZBy2DGzoVOkvkCptsQqrfDJ
>1WR59avYKKdFSzortKY1WGtl6zW+4YK8/uo1r2TdDFxhBdS54mtjNlVMF+f5S7MCdquC1VNdp7oi
>3ex1r31lo2JNJdfGnuqxON0FXifLV8S25q+YZaxm3cTZzurosKJl6qgye1pNpVa1ufgsaCtrJ9Ju
>KrCxHdJsaYsLt8LQtZqCbW8fRVhMdsO2k8WtXXR7Kd4e90YSXS5rERvT3Up3uikyKDeYS1nRxvCy
>jzIud/fUTutWVWxXJa8wt3veDO1TvctpgX3va1/nEii75TVtfH073++yCL/41e+iiGsJyQzJvP81
>V4C3IVwLAhNOCs4RfBu8rH3alRrfIPB9DUyii6qFON68MIY1quENG3a9gGtvSX9W4RPL/hgefKOm
>gHXj4fyKVxSVjPGMf4wOlZpSGxEOJY9t5GMgKxkbAnwkhFnEgihLOcohcVIVgXfkEyV5yVx+hhAL
>u4xvTHnKtKyyVek2Cil9ZstdbvMvIFsPVIrNiCB2XimkxGY367l8cPYLlMdM5R/R8cq+OcWB8rzn
>RMuvz35mMZOCBJ460/MU7VCzoi+NPBUDRMyAZoGZaSkihqGi0gzGtKkhyOhGU3SHg/5i/ka9llPL
>2hbAvfFyOu1pQrRa0onVSqxnDexYgFkevdy1qCn962Ar+xRRVXV7cP3pUEf62L4m8bKvTYpzLYTT
>gI52qAldaFhbG9vkDqK2FVJsL505/nbink253/3BZhP7z93WNfPAXcJ202fc8Da10hJS5EXqe9/u
>6PeplYuQb6xg4QxfOK8lMvB9G/zgCP/JdWEacfpMHNNDtrhuGt7wh38L2QTn98bd3HEGXdzFJ295
>KMyWcBapYOY0n7nIy+HynIsC5vmRec1pfnNy6Hzos+Oox1fNXmNsCxJL3y/Rn34JGx+9PT8Hejma
>7gishxjqT/euyh3t6khovRBj5xXXu57ery+n6ja/ut+Y/nann13nD+7IypUe96znfetzb3mK1U51
>tr+r7IMgPMT77ve/C+QbbFfB4PdOdsibHfEbT3Gq0RFwuLczt5Sf+OXR7fOqP37z/s/tfL9rHee7
>B0PDnD+7HuV94s+DXjeNH704W8/11w+7wahPPdhtT03cQ133vRes7Ge/nBQof/nKB74xhY/2jRY/
>rdOfR+YfYfhOSb4fiCf+8adbfeuziPnMB77Yt3941/MP9ucNv/h1Q/7lm1/zqYI+0b2vaQx/f9uq
>n0T2s49zfadhKLd/yNceKICACYiA84d96Ddy6pd2XVaABuhSFbV6DgiAQieAdadk7uds6aKACsiA
>ejcspQeBXteBE0iBB1MsyPB/DhiAczeAS+aBvrccIZiAI9gIGagzGxiBQKaC/KcbOLiAbleCDXiE
>B3aCG8Z+p1WDNvggL8Vy0TdV/k3YWE8IhSBIhEHXgzK4frt3hUE4EYmSdPgWbl4ofWKIUlhoD99A
>hCjwSWrUYmiGhlUIhtTHhlnIghcDai9khgu1hEUlY3eoHizyhookSexWh0xIiEJ1bknhVnQmbSrj
>g3aYh5bEUkrxDSfAiZ3IiUWRGomIZYsoiP/lZJDIIp7oiWVWDKLYPZXIiJeIR/+Gir8XaZD2Qhlj
>eolWcT2xiar4iYIGI5N4Nbu4Z0bnE3K2JrhYjMbYZimXHakIjCcwh9P2hyTljDQIjToRSQVkgaOY
>jUC4S77IIiZgjudojhdjb1wYjinIc+SoG+iIjlzYL+04YygYjbbYWvYYe/jI/o3lKI/puGP8yHs/
>mI9IV40SBhUFR5CXwoH/GI8BaQL0eAwLyZANySiKxxPXB0hP4W4YOSiWl4wAGZAUeTRN8ZEgeSea
>844HGYXfaGSRcZEqKScLlIkbSZLyaJKQg5IzSZNlIosY8Q0SOZEDKZM/WZNqCJEIeY1g5JE+iZRU
>EpQZ8Q0lYJVXaZU7STZHGZVhMpVU2X8x2ZNdaSVfCZa6gZVYqZWCI5NQSZY4YpZnyZRh1zxP+ZZQ
>Epdy2R5peZVraSlceZcAln9pUZV8WQJ+qQyAGZjUpZSaGJYCN5Ymt5gX0piOqRskgJmZiZmIeUFK
>UR+TKV+VWYtzyZmO8JSw/gGalCmao9kemqmZpWlnPYmaqakgeVkSHDlJbYlotGkbtnmbLOKamQmb
>2+SZn8mb9uGbv6mP5jReSSEbx4mcq+kU3xCcm4ldzekTzwmd6JGcJkGd1TmcrYGdnWGc23kd3emd
>j1mXupls5nkc6JmelwmeRlmcs+me7wmfykma9Omc9nmfxWGFb/ENI0CgBUqg4dlr/bmb/3kYj8gZ
>uIlE7OmWDBoaN/mgLGKgBoqgoyWbKUmhonGKF7qcGNehH9qbIWoZA5qhB3qLy8iOummitIGMKYqh
>KzoCcSggCfmKshmjrrGNkNGNrRgkTblUitmjnDGOlfENIsCkTcqkXjRt/lFKiWN5pCk1o0q6cqHI
>jHRYnCyBf1UKC9UlosvhpE5aZpNIpH5FpTPzhfmpklI3pi9ZNJnhiq+2pinxpWCKCv44GUtapk8K
>I6Akpc3YpSuRp3qKanBKo7rxp4A6jHWab3dqEoeKqDv3kEA6R0O6pYpYqGyahoNZqZWgkVjKqI2a
>G6fqInR5hp2Kp20qnccpknFagYCxLzBpp6w6qZcaqvGmq4jhp38qpFamqpGKq0/Tq7sqqhBjkH2q
>nqpUrBwxg8iqQMLTi8yqGyGArdmKrRs6njwhmRoRrdIadYxGi6S6n/topNB6rJXKT41ImCyirdrK
>rRZpqOuqp24aFd8Q/q/ZOq9P8apB5qqgGqP4mq/N2pFO8a+oE7Di6jYEW7DXuq8h0K8Im7DbQKmh
>6q61AaFGExUVqw0Xi6gOixwsAgIla7IlO7EU66mW6LEYmbHAsbFO6a8C+zULy64iyxXfcLInm7Ii
>07K5Y7P3+rOYqhs7a7I967Ot+qlCO7REe67Y1bE0WxAg+6EvKx06a7QggLRK0bTUQLUMirNoEbNW
>FLUMqxNWex1Ya7RbmxRha7a4gLZpa7C5WbZvW69uK7Ykm7Vs27ZSa7cCYaE0MrbwkhV4+7epUK6C
>yyIWw58U67eHKw/VyieLWzKNy7WGC7nZJrkoMrjU5rhdm7l8dqUq/vINjIuunwu6oVsLYnojnauL
>MxugqstkrKu4umG6UAu7cSu7zaCotfu0JOqzDrq7e8OnnEu5qGG5PZG4w5sO9nohpVu5p9u3KMq8
>CruspDu3ERq8gVu9FhuurXu8qpG8Zzu63WsNlve49uG6Uzq9y2u+03Atzksh0Iu80ssTLfm+XhZR
>veu7cpp0w3W5P5q/1iNk+Nu/6XK7g5qLhBvA5TvAmQZnSXrAe5hKcpimaqq8tPvAuwBmDmy8thu9
>ggZqF2xZGWzAGwy3x+e+kwvC9csckLqqO3G9KEwLncW92KuPqRqsxMilMjxkqTu8wBW7kRK+uuaH
>3wHDYHK/1Jus/ulLwxyzaE6cHnIWpUncTEtMrlL8xPEGLPQbqFW8qeA4VECMNmrmXxSnxW4SpLQa
>xjuaFmTsPmZcau+GucfhxevYG1YMiGhRx7Ugxyb2jHD8vEVsCMNojcN6xWwhyM8gx6bXx3acpRGi
>o7f6xmlcD2ZsjLpbJndcyOtoyAzcFouMKXPscsLbxdnLsW6RKhNKEaRcyqbcK5zMG+PrFR7ayq6c
>eLBMK+tLqGzxrQphaZnMxLtMyLNsv2lhy8B8xuS2ubPCyz3sy8l8ELh8ch4MK7L8HLT8FdJcENRc
>edZ8KtjMt3jBzfYQzI4swbGMyjJ7F+VMJN5scBqszi0svscc/s2s/M7LfG38O8+/O4Vt4c7wAM/x
>zM/ETM9GbM/3bCKArGgz7MzrTLZ8kZIBfQ0DfXoOfc3FnM0JrdCV4J+y+72nfNCdjLsS3QraCdIh
>3c/+O8nrORjlkdIqbdDLkcD/3M7UEbroiy3PzKl9AdNvko0FhdHhrNHj/BY4vUH6vGzRU9APPaI2
>fdPm0T9KrWyQlc4z3R41/VaIUR0VzdDHGEEr7NT+vNWH0dXVYNHlZoW6PNZZHcLAa9blOQ1pzcy0
>9cjGIc7arBdy3ShfHciivCc8LcZxvaC04NcSeNc5ktccbdLt2Qx0XdWaTC6C7caJUdiyANnBxtZ5
>s9glDRm//uwLmT1rw6w3lE3Jlk3RvnLYP9bMnF3Uem3SjCza/gbOamPaxDoZoI0L5zx8AuzaI23M
>np3bqR0Lq32PyUXav03Tbw3VhP0Ls31pTa3cLM2Hn+y5laHbiznUtl3ULmSruG0ZxI2U83LChxOk
>FqPDoAwa2U2Wv5XcayPObKLH2Lje4u2yyF3epZ2p3THfPiUa9k2Q7v3e8N3dX9LfRToaq2DJuziq
>phPfSCyJCtzTFQrYS93gmBOJ6S3hg03hCw6L283dwL3RnnzgZSWjFT7aLCO/XXPeh9yH7HviHt7b
>CVS85l3gwbrDJGx/ICrj91dN8ow4a2zBiLzHMd6VLCvW/vot4jjqoux8oigO1ppW2yFO1gAMoJJd
>zaCa5NONwMxd1k+e2OD3eTcc5K/N2FYK5YjNbFieLbcdwzG+2Z1nw2l+J50N13Ae593X4w4O0epd
>G1Oufo1k582NpIAefYLe59fto76tsYkJF45OvnvO50sO25UB5I3embzEocjUrZaQWBf0l5FOSoMu
>XqWOREeRhKb+GKrO6q3OfQfkgs2SCSpDtq5+hvOkjrau67vO673u678O7MEu7MN+MFpN7MeO7Mmu
>7MvO7M3u7K6e688u7dNO7dVu7deO7dZOF9ya7d3u7d8O7uEu7r4e7eNu7ueO7umu7usu7dvO7u8O
>7/Eu/u/zTu/wUu71ju/5ru/7zu/6vu0v9b8iQuT9TvAFb/AHj/DR8u+1um4iPPAJD/ERL/ETT/Gz
>KoXSdkRayu0Vz/Ed7/EfL14LXxbgVuIgb/Inj/Ip3+x8OKeqqm4br/IxL/MzT/PIIPJ0quGo4W01
>z/M97/M/nz4sj+MOn/EQ/vBAj/RJr/RLXw7u/qgaXzLfzfRTT/VVb/VOT+JtPA5Hb/Vd7/Vfb/JC
>/21t7Mlgb/Znj/Y9r/FxCMaCquNpD/dxL/f07t0tz+RcP/d5r/d7b+5WNWdGj/d8L/iDT/iFb/iH
>j/iJr/iLz/iN7/iPD/mRL/mTT/mVb/mXj/mZr/mb/s/5ne/5nw/6oS/6o0/6pW/6p4/6qa/6q8/6
>re/6rw/7sS/7aa8BtW/7tR8AuF8Ot6/7vG/7hHD7heD7GgD8xN8Iul/8vE8OyD/7ze/8qs78zH8M
>v5/70W/8g0D91F/92I/80i8I3t/93i8J4j/+1//854/+imL9u2/+61/87//97X/92h//5h//8B8Y
>5J/++8//M7H+9A8IGoIaAYWGAYOHioKHjIuEiJCFjpGSk5SGmI6MlImZkoOWlZeeiqanqKmqq6yt
>rq+wsbKztLW2t7i5uru8vb6/wMHCw8TFxsGanJuQmJ/Nl42ikdCkzqCln86ky9PU3NvM0sfj5OXm
>/ufo6err7O3u7/DxssnZ3c/Wj9Gmy6GWnJndovWrR81eOG+gCspbyLChw4cQI0qcSLHiMXoIQ63C
>hlFfwID/Jn1U6LHgt34H81lcybKly5cwY8qc+a7jN1cdRxosqWzas3sfO6X8KQ0ozaNIkypdyrSp
>0182U6LK2SwqT2aISOpUKBThvoRbn4odS7as2bNo0VkdRZQQtpGeSq2NW1QcXLfXqgmlmzWt37+A
>Awse/HRtJbsDVT6qKupmXsWQqxHUKLAk4cuYM2vezLmzMaOeQ4seTbq06cxvT6tezbq169cOQcOe
>Tbu27du4c+vezbu379/AgwsfTry48ePIkytfCc68ufPn0JkGAgA7
>
>------=_NextPart_000_0018_01C1C373.809E2DF0
>Content-Type: text/plain;
>       name="satfronts.txt"
>Content-Transfer-Encoding: quoted-printable
>Content-Disposition: attachment;
>       filename="satfronts.txt"
>
>#! /bin/csh=0A=
>#=0A=
>#=0A=
># This script will plot surface fronts over the US=0A=
># for the most recent available data or user specified date.=0A=
>#=0A=
>#=0A=
>#----------------------------------------------------------------------=0A=
>#=0A=
># Verify that we have the X virtual frame buffer on this host=0A=
>#=0A=
>if ( ! -e /usr/X11R6/bin/Xvfb ) then=0A=
>    echo "This host does not have the X virtual frame buffer running\!"=0A=
>    exit=0A=
>endif=0A=
># =0A=
># set the DISPLAY variable to use the X Virtual frame buffer (display 1)=0A=
>#=0A=
>setenv DISPLAY blizzard.storm.uni.edu:1.0=0A=
>=0A=
># Set time information. This stuff should be parsed automatically=0A=
># and should not require any mods. Look down below for changing =0A=
># this script to fit additional applications.=0A=
>=0A=
># See if an "archive" time has been entered, get the parts=0A=
>=0A=
>if ( ${1}x !=3D x ) then=0A=
>=0A=
>#   check syntax=0A=
>=0A=
>   if ( `echo -n ${1} | wc -c` !=3D 8 ) then=0A=
>       echo "Usage: sfc_map [YYMMDDHH] "=0A=
>       exit=0A=
>    endif=0A=
>=0A=
>    set yy=3D`echo $1 | cut -c 1-2`=0A=
>    if ( $yy < 70 ) then=0A=
>        set cc=3D20=0A=
>    else=0A=
>        set cc=3D19=0A=
>    endif=0A=
>    set mm=3D`echo $1 | cut -c 3-4`=0A=
>    set dd=3D`echo $1 | cut -c 5-6`=0A=
>    set hh=3D`echo $1 | cut -c 7-8`=0A=
>=0A=
>    set realtime_flag=3DF=0A=
>=0A=
>else=0A=
>=0A=
>#   check on who is running this script. Only allow gempak to run this=0A=
>#   in real-time mode=0A=
>=0A=
>    set u=3D`whoami`=0A=
>    if ( $u !=3D gempak ) then=0A=
>        echo "You must be user gempak to run this script in real-time =
>mode."=0A=
>        exit=0A=
>    endif=0A=
>=0A=
>#   Otherwise, get all the date/time information for real-time=0A=
>=0A=
>    setenv TZ UTC=0A=
>    set cc=3D`date "+%Y" | cut -c 1-2`=0A=
>    set yy=3D`date "+%y"`=0A=
>    set mm=3D`date "+%m"`=0A=
>    set dd=3D`date "+%d"`=0A=
>    set HH=3D`date "+%H"`=0A=
>=0A=
>    set realtime_flag=3DT=0A=
>=0A=
>endif=0A=
>=0A=
># Get the month name=0A=
>=0A=
>if ( $mm =3D=3D "01" ) set mon=3DJAN=0A=
>if ( $mm =3D=3D "02" ) set mon=3DFEB=0A=
>if ( $mm =3D=3D "03" ) set mon=3DMAR=0A=
>if ( $mm =3D=3D "04" ) set mon=3DAPR=0A=
>if ( $mm =3D=3D "05" ) set mon=3DMAY=0A=
>if ( $mm =3D=3D "06" ) set mon=3DJUN=0A=
>if ( $mm =3D=3D "07" ) set mon=3DJUL=0A=
>if ( $mm =3D=3D "08" ) set mon=3DAUG=0A=
>if ( $mm =3D=3D "09" ) set mon=3DSEP=0A=
>if ( $mm =3D=3D "10" ) set mon=3DOCT=0A=
>if ( $mm =3D=3D "11" ) set mon=3DNOV=0A=
>if ( $mm =3D=3D "12" ) set mon=3DDEC=0A=
>=0A=
>=0A=
>#--------------------------------------------------------------------=0A=
>#=0A=
>#  Customize shell variables here...=0A=
>#=0A=
>#--------------------------------------------------------------------=0A=
>=0A=
>#  Build all the variables that use times.=0A=
>=0A=
>set frontdir=3D$GEMDATA/nwx/hpc/fronts=0A=
>=0A=
>cd $frontdir=0A=
>=0A=
>set frontfile=3D`ls *.front | tail -1`=0A=
>=0A=
>set filetime=3D`grep VALID $frontfile | cut -c 11-13`=0A=
>=0A=
># get the last satellite image=0A=
>=0A=
>set satdir=3D$GEMDATA/images/sat/GOES-8/4km/IR=0A=
>=0A=
>cd $satdir=0A=
>=0A=
>set satfil=3D`ls IR_* | tail -1`=0A=
>=0A=
>set datelabel=3D`echo $satfil | cut -c 4-11`=0A=
>=0A=
>set timelabel=3D`echo $satfil | cut -c 13-16`=0A=
>=0A=
>cd=0A=
>=0A=
>#--------------------------------------------------------------------=0A=
>#=0A=
>#  Set the output plot file=0A=
>#=0A=
>#--------------------------------------------------------------------=0A=
>=0A=
># Set the GIF file names. Real-time and archive derived files go in=0A=
># different places to avoid trampling on each other.=0A=
>=0A=
>if ( $realtime_flag =3D=3D T ) then=0A=
>    set tmpgif=3Dsatfronts.gif.tmp=0A=
>    set finalgif=3D/space/html/images/satfronts.gif=0A=
>else=0A=
>    set tmpgif=3Dfronts.gif.tmp=0A=
>    set finalgif=3Dfronts.gif=0A=
>endif=0A=
>=0A=
>=0A=
># Do minimal error checking - verify the gempak file is there.=0A=
>=0A=
>#if ( ! -e $satfil ) then=0A=
>#   echo "GEMPAK file $satfil does not exist. "=0A=
>  =0A=
>#endif=0A=
>=0A=
>=0A=
>=0A=
>#--------------------------------------------------------------------=0A=
>#=0A=
>#  The rest is changing what scripts run..=0A=
>#=0A=
>#--------------------------------------------------------------------=0A=
>=0A=
># Copy in the black color table=0A=
>=0A=
>cp ~gempak/scripts/coltbl.xwp_black coltbl.xwp=0A=
>=0A=
># Remove the temporary gif file if it exists=0A=
>=0A=
>if ( -e $tmpgif ) rm $tmpgif=0A=
>=0A=
>=0A=
># Run the gempak program(s)=0A=
>=0A=
>#-----------------------------------------------=0A=
># Plot the Fronts/Highs/Lows=0A=
>=0A=
>gpfront<< GPFRONT_FLAG=0A=
>satfil=3D$GEMDATA/images/sat/GOES-8/4km/IR/$satfil=0A=
>map=3D31=0A=
>garea=3D23.0;-120.0;55.0;-65.0=0A=
>proj=3Dsat//1;1;1;1=0A=
>radfil=3D=0A=
>latlon=3D0=0A=
>panel=3D0=0A=
>text=3D1.5/22/HW=0A=
>clear=3Dyes=0A=
>title=3D32/-1/Fronts at $filetime=0A=
>asfil=3D$GEMDATA/nwx/hpc/fronts/$frontfile=0A=
>device=3Dgf|$tmpgif|900;650=0A=
>ashr=3D=0A=
>asatt=3D1/N=0A=
>=0A=
>run=0A=
>=0A=
>exit=0A=
>=0A=
>GPFRONT_FLAG=0A=
>=0A=
>gptext<< GPTEXT_FLAG=0A=
>clear=3Dno=0A=
>txtfil=3D/export/home/gempak/scripts/storm.txt=0A=
>txtloc=3D.01;.03=0A=
>text=3D1/22/HW=0A=
>colors=3D25=0A=
>column=3D1=0A=
>=0A=
>run=0A=
>=0A=
>exit=0A=
>=0A=
>GPTEXT_FLAG=0A=
>=0A=
>gpend=0A=
>=0A=
># Move the temporary gif to its final resting place.=0A=
>=0A=
>if ( -e $tmpgif ) then=0A=
>=0A=
>    mv $tmpgif $finalgif=0A=
>else=0A=
>       echo the gif aint here dude=0A=
>=0A=
>endif=0A=
>
>------=_NextPart_000_0018_01C1C373.809E2DF0--
>