PROGRAM WeatherCalc

!---------------------------------------------------------------------------
! This program calculates the heat index, hurricane scale, potential temp
! temperature conversion between farhenheit and celsius, wind chill,
! windspeed conversion between knots and mph.  Variables are:
!
! Author: Meteorologist Andy Kren
! Date: May 3rd, 2008
!---------------------------------------------------------------------------

	IMPLICIT NONE

    REAL :: TempH, RH, Heat_Index, ptemp, Celsius, p0, p1, Maxtemp, Kelvin, Celsius2, TempF, TempC, Fahrenheit
    REAL :: WindChill, WindSpeed, Knots, Mph, WindSpeed_mph, WindSpeed_knots, Temperature, dpf, Tc, Tdc, E, Es
    INTEGER :: Knots_or_Mph
    CHARACTER :: ForCresponse, Response, FinalChoice


   PRINT *
   PRINT *, "Weather Calculator"
   PRINT *, "  By: Andy Kren   "
   PRINT *, "------------------" 

  DO
   ! Ask user what he wants to calculate
   PRINT *, "Enter 1-6 corresponding to what you want to calculate:"
   PRINT *
   PRINT *, "1. Potential Temperature"
   PRINT *, "2. Temperature Conversion"
   PRINT *, "3. Wind Speed Conversion"
   PRINT *, "4. Wind Chill"
   PRINT *, "5. Heat Index"
   PRINT *, "6. Hurricane Scale"
   PRINT *
   READ *, Response

   DO
     IF (Response == "1") THEN
        PRINT *
        PRINT *, "Enter the level (mb) at which you wish to determine Theta:"
        READ *, p1

        PRINT *, "Enter the surface pressure (i.e. around 1030mb or less):"
        READ *, p0

        PRINT *, "Enter the temperature in Celsius at level to mix down:"
        READ *, Celsius

        !convert temperature to kelvin
        Kelvin = Celsius + 273

        ! calculate potential temperature
        ptemp = (Kelvin)*((p0 / p1)**(.286))

        ! convert potential temp to fahrenheit
        Celsius2 = ptemp - 273
        Maxtemp = (1.8 * Celsius2) + 32

        PRINT *, "The Max Temperature (or Potential Temp) = ", MaxTemp
        EXIT

     ELSE IF (Response == "2") THEN
        PRINT *
        PRINT *, "Enter (1) to convert from F to C or (2) for vice versa:"
        READ *, ForCresponse
        IF (ForCresponse == "1") THEN
           PRINT *, "Enter the Fahrenheit Temperature"
           READ *, TempF
           ! Calculate the Celsius temp.
           Celsius = (0.555555556) * (TempF - 32) 
           ! Print out celsius temp
           100 FORMAT (1X, A28, F10.2)
           PRINT 100, "The Celsius Temperature is:", Celsius
        ELSE IF (ForCresponse == "2") THEN
           PRINT *, "Enter the Celsius Temperature"
           READ *, TempC
           ! Calculate the Fahrenheit
           Fahrenheit = (1.8 * TempC) + 32
           ! Print out Fahrenheit temp
           110 FORMAT (1X, A31, F10.2)
           PRINT 110, "The Fahrenheit Temperature is:", Fahrenheit
        ELSE
           PRINT *, "THAT IS NOT A VALID RESPONSE. RE-RUN PROGRAM."
           STOP
        END IF
        EXIT

     ELSE IF (Response == "3") THEN
        PRINT *
        PRINT *, "Enter a (1) to convert from knots to mph or (2) for vice versa:"
        READ *, Knots_or_Mph
        IF (Knots_or_Mph == 1) THEN
           ! calculate the knots to mph conversion and ask if wants to reconvert
           PRINT *, "Enter the wind speed in knots:"
           READ *, Knots
           WindSpeed_mph = Knots * 1.15
           130 FORMAT(1X, A24, F12.1, T39, A3)
           PRINT 130, "The Wind Speed in MPH = ", WindSpeed_mph, "mph"
           PRINT *
        ELSE IF (Knots_or_Mph == 2) THEN
           ! calculate the mph to knots conversion and ask if wants to reconvert
           PRINT *, "Enter the wind speed in mph:"
           READ *, Mph
           WindSpeed_knots = Mph / (1.15)
           140 FORMAT(1X, A26, F12.1, T41, A5)
           PRINT 140, "The Wind Speed in Knots = ", WindSpeed_knots, "knots"
        ELSE 
           PRINT *, "***Enter a valid number (1) or (2)***"
           STOP
        END IF
        EXIT

     ELSE IF (Response == "4") THEN
        PRINT *
        PRINT *, "Enter the Surface Temperature (degrees F)"
        READ *, Temperature
        PRINT *, "Enter the Wind Speed (mph)"
        READ *, WindSpeed
        IF (WindSpeed < 0) STOP "*Re-run program. Wind speed not negative*"
        WindChill = 35.74 + (0.6215 * Temperature) - (35.75 * (WindSpeed**&
        &0.16)) + (0.4275 * Temperature * (WindSpeed ** 0.16))

        200 FORMAT(1X, A, 1X, F10.1)
        IF (WindChill > 40) THEN
           PRINT 200, "The Wind Chill is:", WindChill
           PRINT *, "    This weather is no problem!"
        ELSE IF (WindChill > 20) THEN
           PRINT 200, "The Wind Chill is:", WindChill
           PRINT *
           PRINT *, "    That's Chilly!"
        ELSE IF (WindChill > 10) THEN
           PRINT 200, "The Wind Chill is:", WindChill
           PRINT *
           PRINT *, "    Pretty Nippy Outside!"
        ELSE IF (WindChill > 0) THEN
           PRINT 200, "The Wind Chill is:", WindChill
           PRINT *
           PRINT *, "    That's Pretty Cold my Friend!!"
        ELSE IF (WindChill < 0) THEN
           PRINT 200, "The Wind Chill is:", WindChill
           PRINT *
           PRINT *, "    Bitterly Freezing Cold!!!"
        END IF
        PRINT *
        EXIT

     ELSE IF (Response == "5") THEN
        DO
          PRINT *
          PRINT *, "Enter the air Temperature (Fahrenheit >= 80): "
          READ *, TempH
          IF (TempH >= 80) EXIT
        END DO
        
        PRINT *, "Enter the Dewpoint Temperature in Fahrenheit: "
        READ *, dpf
        
        ! convert temp and dewpoint to celsius
        Tc = (0.555555556) * (TempH - 32) 
        Tdc = (0.555555556) * (dpf - 32) 

        ! calculate vapor pressure and saturation vapor pressure (E, Es)
        E = 6.11 * (10.0**((7.5*Tdc)/(237.7+Tdc)))
        Es = 6.11 * (10.0**((7.5*Tc)/(237.7+Tc)))

        ! calulate RH to be used in finding the heat index:
        RH = (E / Es) * 100

         ! calculate heat index
         Heat_Index = (-42.379) + (2.04901523 * TempH) + (10.14333127 * RH) &
         &- (0.22475541 * TempH * RH) - (6.83783e-3 * (TempH ** 2)) &
         &- (5.481717e-2 * (RH ** 2)) + (1.22874e-3 * (TempH ** 2) * RH) &
         &+ (8.5282e-4 * TempH * (RH ** 2)) - (1.99e-6 * (TempH ** 2) * (RH ** 2))

         ! print out results
         PRINT *
         PRINT 190, "----------------------------------------"
         PRINT 170, "For your given Temp: ", TempH
         PRINT 240, "For your given Dewpoint: ", dpf
         PRINT *
         PRINT 180, "Heat Index = ", Heat_Index
         PRINT 230, "Relative Humidity = ", RH, "%"
         PRINT 190, "----------------------------------------"
         PRINT *
         170 FORMAT(1X, A21, F6.2)
         240 FORMAT(1X, A25, F6.2)
         180 FORMAT(1X, A13, F6.1)
         190 FORMAT(1X, A40)   
         230 FORMAT(1X, A20, F6.2, A1)

        EXIT

     ELSE IF (Response == "6") THEN
        PRINT *
        PRINT *, "Saffir Simpson Hurricane Scale"
        PRINT *, "------------------------------"
        
          PRINT *
          PRINT *, "Enter the Wind Speed of the Storm (mph):"
          READ *, WindSpeed

          ! Determine the category of the storm
          IF (WindSpeed <= 0) THEN
             PRINT *, "Storm cannot have a 0 or negative wind-RE-RUN PROGRAM"
             STOP
          ELSE IF (WindSpeed < 38) THEN
             PRINT *
             PRINT *, "This disturbance is too weak! No organized storm."
          ELSE IF (WindSpeed == 38) THEN
             PRINT *
             PRINT *, "A Tropical Depression"
          ELSE IF (WindSpeed >= 39.AND.WindSpeed <= 73) THEN
             PRINT *
             PRINT *, "A Tropical Storm"
          ELSE IF (WindSpeed >= 74.AND.WindSpeed <= 95) THEN
             PRINT *
             PRINT *, "A Category 1 Hurricane"
          ELSE IF (WindSpeed >= 96.AND.WindSpeed <= 110) THEN
             PRINT *
             PRINT *, "A Category 2 Hurricane"
          ELSE IF (Windspeed >= 111.AND.WindSpeed <= 130) THEN
             PRINT *
             PRINT *, "A Category 3 Hurricane"
          ELSE IF (WindSpeed >= 131.AND.WindSpeed <= 155) THEN
             PRINT *
             PRINT *, "A Category 4 Hurricane"
          ELSE 
             PRINT *
             PRINT *, "A Category 5 Hurricane"
          END IF
          EXIT  
     ELSE 
        PRINT *, "Enter a Valid Integer between 1 and 6, and re-run program"
        EXIT
     END IF
    END DO
     ! ask user if he wants to re-calcluate another quantity
     PRINT *
     PRINT *, "Do you want to Calculate another thing? (Y or N)?"
     READ *, FinalChoice
     IF (FinalChoice == "N".OR.FinalChoice == "n") THEN
       EXIT
     END IF
  END DO

END PROGRAM WeatherCalc<!--'"</title></head>-->

<script type="text/javascript">
//OwnerIQ
var __oiq_pct = 50;
if( __oiq_pct>=100 || Math.floor(Math.random()*100/(100-__oiq_pct)) > 0 ) {
var _oiqq = _oiqq || [];
_oiqq.push(['oiq_addPageBrand','Lycos']);
_oiqq.push(['oiq_addPageCat','Internet > Websites']);
_oiqq.push(['oiq_addPageLifecycle','Intend']);
_oiqq.push(['oiq_doTag']);
(function() {
var oiq = document.createElement('script'); oiq.type = 'text/javascript'; oiq.async = true;
oiq.src = document.location.protocol + '//px.owneriq.net/stas/s/lycosn.js';
var s = document.getElementsByTagName('script')[0]; s.parentNode.insertBefore(oiq, s);
})();
}
//Google Analytics
var _gaq = _gaq || [];
_gaq.push(['_setAccount','UA-21402695-19']);
_gaq.push(['_setDomainName','tripod.com']);
_gaq.push(['_setCustomVar',1,'member_name','krenfamily',3]);
_gaq.push(['_trackPageview']);
(function() {
var ga = document.createElement('script'); ga.type = 'text/javascript'; ga.async = true;
ga.src = ('https:' == document.location.protocol ? 'https://ssl' : 'http://www') + '.google-analytics.com/ga.js';
var s = document.getElementsByTagName('script')[0]; s.parentNode.insertBefore(ga, s);
})();
//Lycos Init
function getReferrer() {
var all= this.document.cookie;
if (all== '') return false;
var cookie_name = 'REFERRER=';
var start = all.lastIndexOf(cookie_name);
if (start == -1) return false;
start += cookie_name.length;
var end = all.indexOf(';', start);
if (end == -1) end = all.length;
return all.substring(start, end);
}
function getQuery() {
var rfr = getReferrer();
if (rfr == '') return false;
var q = extractQuery(rfr, 'yahoo.com', 'p=');
if (q) return q;
q = extractQuery(rfr, '', 'q=');
return q ? q : "";
}
function extractQuery(full, site, q_param) {
var start = full.lastIndexOf(site);
if (start == -1) return false;
start = full.lastIndexOf(q_param);
if (start == -1) return false;
start += q_param.length;
var end = full.indexOf('&', start);
if (end == -1) end = full.length;
return unescape(full.substring(start, end)).split(" ").join("+");
}
function generateHref(atag, template){
atag.href=template.replace('_MYURL_', window.location.href.replace('http://', '')).replace('_MYTITLE_','Check%20out%20this%20Tripod%20Member%20site!'); 
}
var lycos_ad = Array();
var lycos_onload_timer;
var cm_role = "live";
var cm_host = "tripod.lycos.com";
var cm_taxid = "/memberembedded";
var tripod_member_name = "krenfamily";
var tripod_member_page = "krenfamily/sitebuildercontent/sitebuilderfiles/weathercalc.f90";
var tripod_ratings_hash = "1713518595:e84f2d5df2bbda7ca888d09b10333056";

var lycos_ad_category = null;

var lycos_ad_remote_addr = "209.202.244.9";
var lycos_ad_www_server = "www.tripod.lycos.com";
var lycos_ad_track_small = "";
var lycos_ad_track_served = "";
var lycos_search_query = getQuery();
</script>

<script type="text/javascript" src="https://scripts.lycos.com/catman/init.js"></script>

<script type='text/javascript'>
 var googletag = googletag || {};
 googletag.cmd = googletag.cmd || [];
 (function() {
   var gads = document.createElement('script');
   gads.async = true;
   gads.type = 'text/javascript';
   var useSSL = 'https:' == document.location.protocol;
   gads.src = (useSSL ? 'https:' : 'http:') +
     '//www.googletagservices.com/tag/js/gpt.js';
   var node = document.getElementsByTagName('script')[0];
   node.parentNode.insertBefore(gads, node);
 })();
</script>

<script type='text/javascript'>
 googletag.cmd.push(function() {
   googletag.defineSlot('/95963596/TRI_300X250_dfp', [300, 250], 'div-gpt-ad-1450204159126-0').addService(googletag.pubads());
   googletag.defineSlot('/95963596/TRI_above_728x90_dfp', [728, 90], 'div-gpt-ad-1450204159126-1').addService(googletag.pubads());
   googletag.defineSlot('/95963596/TRI_below_728x90_dfp', [728, 90], 'div-gpt-ad-1450204159126-2').addService(googletag.pubads());
   googletag.pubads().enableSingleRequest();
   googletag.enableServices();
 });
</script>


<script type="text/javascript"> 
(function(isV)
{
    if( !isV )
    {
        return;
    }
    var adMgr = new AdManager();
    var lycos_prod_set = adMgr.chooseProductSet();
    var slots = ["leaderboard", "leaderboard2", "toolbar_image", "toolbar_text", "smallbox", "top_promo", "footer2", "slider"];
    var adCat = this.lycos_ad_category;
    adMgr.setForcedParam('page', (adCat && adCat.dmoz) ? adCat.dmoz : 'member');
    if (this.lycos_search_query)
    {
        adMgr.setForcedParam("keyword", this.lycos_search_query);
    } 
    else if(adCat && adCat.find_what)
    {
        adMgr.setForcedParam('keyword', adCat.find_what);
    }
    
    for (var s in slots)
    {
        var slot = slots[s];
        if (adMgr.isSlotAvailable(slot))
        {
            this.lycos_ad[slot] = adMgr.getSlot(slot);
        }
    }

    adMgr.renderHeader();
    adMgr.renderFooter();
}((function() {

var w = 0, h = 0, minimumThreshold = 300;

if (top == self)
{
    return true;
}
if (typeof(window.innerWidth) == 'number' )
{
    w = window.innerWidth;
    h = window.innerHeight;
}
else if (document.documentElement && (document.documentElement.clientWidth || document.documentElement.clientHeight))
{
    w = document.documentElement.clientWidth;
    h = document.documentElement.clientHeight;
}
else if (document.body && (document.body.clientWidth || document.body.clientHeight))
{
    w = document.body.clientWidth;
    h = document.body.clientHeight;
}
return ((w > minimumThreshold) && (h > minimumThreshold));
}())));




window.onload = function()
{
    var f = document.getElementById("FooterAd");
    var b = document.getElementsByTagName("body")[0];
    b.appendChild(f);
    f.style.display = "block";
    document.getElementById('lycosFooterAdiFrame').src = '/adm/ad/footerAd.iframe.html';
    


    
    // DOM Inj Ad
    (function(isTrellix)
    {
        var e = document.createElement('iframe');
        e.style.border = '0';
        e.style.margin = 0;
        e.style.display = 'block';
        e.style.cssFloat = 'right';
        e.style.height = '254px';
        e.style.overflow = 'hidden';
        e.style.padding = 0;
        e.style.width = '300px';


        var isBlokedByDomain = function( href )
        {
            var blockedDomains = [
                "ananyaporn13000.tripod.com",
                "xxxpornxxx.tripod.com"
            ];
            var flag = false;
            
            for( var i=0; i<blockedDomains.length; i++ )
            {
                if( href.search( blockedDomains[ i ] ) >= 0 )
                {
                    flag = true;
                }
            }
            return flag;
        }

        var getMetaContent = function( metaName )
        {
            var metas = document.getElementsByTagName('meta');
            for (i=0; i<metas.length; i++)
            { 
                if( metas[i].getAttribute("name") == metaName )
                { 
                    return metas[i].getAttribute("content"); 
                } 
            }
            return false;
        }
        
        var getCommentNodes = function(regexPattern)
        {
            var nodes = {};
            var nodesA = [];
            var preferredNodesList = ['a', 'c', 'b'];
        
            (function getNodesThatHaveComments(n, pattern)
            {
                if (n.hasChildNodes())
                {
                    if (n.tagName === 'IFRAME')
                    {
                        return false;
                    }
                    for (var i = 0; i < n.childNodes.length; i++)
                    {
                        if ((n.childNodes[i].nodeType === 8) && (pattern.test(n.childNodes[i].nodeValue)))
                        {
                            var areaName = pattern.exec(n.childNodes[i].nodeValue)[1];
                            nodes[areaName] = n;
                        }
                        else if (n.childNodes[i].nodeType === 1)
                        {
                            getNodesThatHaveComments(n.childNodes[i], pattern);
                        }
                    }
                }
            }(document.body, regexPattern));

            for (var i in preferredNodesList)
            {
                if (nodes[preferredNodesList[i]])
                {
                    if( isTrellix && nodes[preferredNodesList[i]].parentNode.parentNode.parentNode.parentNode )
                    {
                        nodesA.push(nodes[preferredNodesList[i]].parentNode.parentNode.parentNode.parentNode);
                    }
                    else
                    {
                        nodesA.push( nodes[preferredNodesList[i]] );
                    }
                }
            }
            return nodesA;
        }
        
        
        var properNode = null;
        var areaNodes = getCommentNodes( new RegExp( '^area Type="area_(\\w+)"' ) );

        for (var i = 0; i < areaNodes.length; i++)
        {
            var a = parseInt(getComputedStyle(areaNodes[i]).width);
            if ((a >= 300) && (a <= 400))
            {
                properNode = areaNodes[i];
                break;
            }
        }


        var propertyName = getMetaContent("property") || false;
        if( isTrellix && (properNode) )
        {
            e.src = '/adm/ad/injectAd.iframe.html';
            properNode.insertBefore(e, properNode.firstChild);
        }
        else if( isTrellix && !( properNode ) ) // Slap the ad eventhought there is no alocated slot
        {
            e.src = '/adm/ad/injectAd.iframe.html';
            e.style.cssFloat = 'none';
            var cdiv = document.createElement('div');
            cdiv.style = "width:300px;margin:10px auto;";
            cdiv.appendChild( e );
            b.insertBefore(cdiv, b.lastChild);
        }
        else if( !isBlokedByDomain( location.href ) )
        {
            var injF = document.createElement('iframe');
            injF.style.border = '0';
            injF.style.margin = 0;
            injF.style.display = 'block';
            injF.style.cssFloat = 'none';
            injF.style.height = '254px';
            injF.style.overflow = 'hidden';
            injF.style.padding = 0;
            injF.style.width = '300px';
            injF.src = '/adm/ad/injectAd.iframe.html';

            if( b && ( !isTrellix || ( typeof isTrellix == "undefined" ) ) ) // All other tripod props
            {
                var cdiv = document.createElement('div');
                cdiv.style = "width:300px;margin:10px auto;";
                cdiv.appendChild( injF );
                b.insertBefore(cdiv, b.lastChild);
            } 
        }
  }( document.isTrellix ));
}

</script>

<div id="tb_container" style="background:#DFDCCF; border-bottom:1px solid #393939; position:relative; z-index:999999999!important">
<!--form name="search" onSubmit="return searchit()" id='header_search' >
<input type="text" placeholder="Search" size=30 name="search2" value="">
<input type="button" value="Go!" onClick="searchit()">
</form>
<style>
form#header_search {
    width: 916px;
    margin: 0 auto 8px;
    position: relative;
}


form#header_search input {
    height: 40px;
    font-size: 14px;
    line-height: 40px;
    padding: 0 8px;
    box-sizing: border-box;
    background: #F4F2E9;
    border: 1px solid #BBB8B8;
    transition: background-color 300ms ease-out,
                color 300ms ease;
}

form#header_search input[type="text"] {
    width: 100%;
}
form#header_search input[type="text"]:focus {
    border-color: #A2D054;
    background-color: #fff;
    box-shadow: 0 0px 12px -4px #A2D054;
}



form#header_search input[type="button"] {
    position: absolute;
    top: 1px;
    right: 1px;
    opacity: 1;
    background: #DFDCCF;
    color: #463734;
    width: 125px;
    cursor: pointer;
    height: 38px;
    border: none;
}
form#header_search input[type="text"]:focus ~ input[type='button']:hover,
form#header_search input[type='button']:hover {
    background-color: #A5CE56;
    color: #fff;
}
form#header_search input[type="text"]:focus ~ input[type='button'] {
    background-color: #52AEDF;
    color: #fff;
}

</style>

<script>
function searchit(){
    
    // determine environment 
    var search_env 
    if (lycos_ad_www_server.indexOf(".pd.") > -1) {
       search_env = 'http://search51.pd.lycos.com/a/';
    } else if (lycos_ad_www_server.indexOf(".qa.") > -1) {
       search_env = 'http://search51.qa.lycos.com/a/';
    } else {
       search_env = 'http://search51.lycos.com/a/';
    }

var search_term = encodeURIComponent(document.search.search2.value)
var search_url = search_env+search_term;
window.open(search_url);

return false
}
</script-->
<style>
    .adCenterClass{margin:0 auto}
</style>
<div id="tb_ad" class="adCenterClass" style="display:block!important; overflow:hidden; width:916px;">

<div id="ad_container" style="display:block!important; float:left; width:728px ">
<script type="text/javascript">
if (typeof lycos_ad !== "undefined" && "leaderboard" in lycos_ad) {
  document.write(lycos_ad['leaderboard']);
}
</script>
</div>
</div>
</div>
<script type="text/javascript">
if (typeof lycos_ad !== "undefined" && "slider" in lycos_ad) {
  document.write(lycos_ad['slider']);
}
</script> <!-- added 7/22 -->
<div id="FooterAd" style="background:#DFDCCF; border-top:1px solid #393939; clear:both; display:none; width:100%!important; position:relative; z-index:999999!important; height:90px!important"> 
<div class="adCenterClass" style="display:block!important; overflow:hidden; width:916px;">
<div id="footerAd_container" style="display:block!important; float:left; width:728px">
<iframe id="lycosFooterAdiFrame" style="border:0; display:block; float:left; height:96px; overflow:hidden; padding:0; width:750px"></iframe>
</div>
</div>
</div>