Go Back  PPRuNe Forums > Flight Deck Forums > Tech Log
Reload this Page >

Accurate Temperature Correction Formula

Tech Log The very best in practical technical discussion on the web

Accurate Temperature Correction Formula

Old 11th Jan 2016, 03:38
  #1 (permalink)  
Thread Starter
 
Join Date: Jun 2001
Location: at the edge of the alps
Posts: 447
Received 0 Likes on 0 Posts
Accurate Temperature Correction Formula

I am trying to develop an Excel formula for exact temperature correction for a friend, following the "accurate" formula from PANS-OPS Section III 4.3.4.1:


Correction = trianglehPaircraft - trianglehGaircraft = (-triangleTstd/L0) * Ln[(1 + L0 * trianglehPaircraft)/(T0 + L0*hPaerodrome)]



Where:
trianglehPaircraft = Aircraft height above aerodrome (pressure) ft
trianglehGaircraft = Aircraft height above aerodrome (geopotential) ft
triangleTstd = Temperature deviation from the ISA temperature K
L0 = Standard temperature lapse rate with pressure
altitude in the first layer (sea level to tropopause)
of the ISA.(See note (b) below) K/ft
T0 = Standard temperature at sea level K
hPaerodrome = Aerodrome height (pressure) ft


I have devised the following Excel Formula:

=(-"ISADEV"/-0,0019812)*LN((1+-0,0019812*"MDH")/(288,15+-0,0019812*"APAlt"))

where "ISADEV" is the temperature deviation from ISA (usually a negative number in cold weather)
"MDH" is the minimum descent height in ft
"APAlt" is the airport altitude in ft
-0,0019812 is the lapse rate in K/ft and 288,15 is ISA/sea level in K.

This results in meaningless results.

I am challenged by not really understanding why the basic formula should work and would be grateful for any hints on where I might have erred in applying it.
Alpine Flyer is offline  
Old 11th Jan 2016, 15:04
  #2 (permalink)  
 
Join Date: May 2005
Location: France
Age: 69
Posts: 1,142
Received 3 Likes on 2 Posts
4ft per degree C dev per 1000ft above source elevation.

I don't know how to express that in terms that would fit your excel formula, sorry!
eckhard is offline  
Old 11th Jan 2016, 15:35
  #3 (permalink)  
 
Join Date: Jun 2003
Location: LSZH, oder in der Nähe
Posts: 348
Likes: 0
Received 0 Likes on 0 Posts
It all seems terribly overcomplicated


The ISA model is quite simple and straightforward...

Metric:
T = T0 - 6.5 (h/1000)
where
T0 is the ISA temperature (15C)
h is altitude in meters
T is the temperature at altitude h
Imperial:
T = T0 - 1.98 (h/1000)
where
T0 is the ISA temperature (15C)
h is altitude in feet
T is the temperature at altitude h

WeeJeem is offline  
Old 11th Jan 2016, 20:01
  #4 (permalink)  
 
Join Date: May 2005
Location: France
Age: 69
Posts: 1,142
Received 3 Likes on 2 Posts
WeeJeem,

Thanks for those formulæ but I think that Alpine Flyer is looking for a formula that expresses the error between indicated and true altitude, given a deviation from the ISA temperature.
eckhard is offline  
Old 11th Jan 2016, 20:10
  #5 (permalink)  
Thread Starter
 
Join Date: Jun 2001
Location: at the edge of the alps
Posts: 447
Received 0 Likes on 0 Posts
@eckhard: the 4% formula is the "simple" formula but limited to temperatures down to -15C. Beyond that, PANS OPS offers a table that covers lower temperatures but is limited to 5000ft and two more formulae, the more complicated of which was the reason for my post.

As probably all those formulae are but child's play for a modern Computer, IMHO it does not make sense to use the simplified versions rather than going for the most accurate calculation.

@WeeJam: I have understood lapse rate (which is 1,9812C per 1000ft to be exact ) but lapse rate is just one building block for calculating cold temperature correction.

Any other input welcome. It must be possible to solve this!
Alpine Flyer is offline  
Old 12th Jan 2016, 02:00
  #6 (permalink)  
 
Join Date: Feb 2000
Location: 500 miles from Chaikhosi, Yogistan
Posts: 4,290
Received 139 Likes on 63 Posts
There's a simpler ICAO formula that's pretty good, without the ln....

Unless your friend is seriously seriously scientific (and in which case would know how to use it anyway?) try that one.
compressor stall is offline  
Old 12th Jan 2016, 02:40
  #7 (permalink)  
 
Join Date: May 2006
Location: Australia
Posts: 494
Received 17 Likes on 7 Posts
Try this...

=(B6-B2)*((15-(B8+(0.00198*B2)))/((273+(B8+(0.00198*B2)))-(0.5*0.00198*((B6-B2)+B2))))

Where: B6 is aircraft alt (ft)
B2 is aerodrome elevation (ft)
B8 is aerodrome temperature

The formula is copied straight from my spreadsheet and should work if the correct values are assigned. The formula is an older one but should work for your purposes.

Alpha
alphacentauri is offline  
Old 16th Jan 2016, 15:25
  #8 (permalink)  
Thread Starter
 
Join Date: Jun 2001
Location: at the edge of the alps
Posts: 447
Received 0 Likes on 0 Posts
@alphacentauri:

Thanks, this seems to be the first equation from 4.3.3 that is valid for altimeter setting sources up to 10.000ft (which should be sufficient for Europe) but only 5000ft above source.

According to 8168 this equation was used to calculate the "standard table" that is also valid up to 5000ft above source.

So it can't be used for higher altitudes of rather low-lying airports in mountainous areas (such as Innsbruck or Geneva) which have minimum crossing altitudes more than 5000ft above the airport.
Alpine Flyer is offline  
Old 17th Oct 2016, 02:42
  #9 (permalink)  
 
Join Date: Oct 2016
Location: Japan
Posts: 1
Likes: 0
Received 0 Likes on 0 Posts
Cannot calculate Accurate Method

Hi everyone,

I have a simple question, cannot calculate following formula for use following example value as Ln value was minus. (excel shows "#NUM!" error)
How can I get a correct answer ?

Ex.
Airport Temp; -40C
Airport Elevation; 10ft
Aircraft Altitude; 3000ft

=(-"ISADEV"/-0,0019812)*LN((1+-0,0019812*"MDH")/(288,15+-0,0019812*"APAlt"))

where "ISADEV" is the temperature deviation from ISA (usually a negative number in cold weather)
"MDH" is the minimum descent height in ft
"APAlt" is the airport altitude in ft
-0,0019812 is the lapse rate in K/ft and 288,15 is ISA/sea level in K.
ja1234 is offline  
Old 4th Jan 2019, 11:57
  #10 (permalink)  
 
Join Date: Jan 2019
Location: Hamburg
Posts: 2
Likes: 0
Received 0 Likes on 0 Posts
Hi to all,
has anyone got a solution for this? Maybe there's a mistake within the displayed formula.
I've tried out everything (+/- values, °C iso K...) but there's never a result that makes any sense.
Has anyone access to the equation #24 of the ESDU publication Performance, Vol. 2, Item No. 77022?
Thanxx al lot so far...
Lumbanraja is offline  
Old 5th Jan 2019, 10:42
  #11 (permalink)  
 
Join Date: May 2005
Location: France
Age: 69
Posts: 1,142
Received 3 Likes on 2 Posts
Error of approx 660ft? True altitude 2340ft
eckhard is offline  
Old 7th Jan 2019, 01:21
  #12 (permalink)  
 
Join Date: Jan 2004
Location: Hiding..... in one hemisphere or another
Posts: 1,067
Received 1 Like on 1 Post
Spreadsheet

This might help. There's a link to a spreadsheet at the bottom of the article under "Further Reading"
Atlas Shrugged is offline  
Old 8th Jan 2019, 15:49
  #13 (permalink)  
 
Join Date: Jan 2011
Location: Seattle
Posts: 715
Likes: 0
Received 3 Likes on 2 Posts
Originally Posted by Alpine Flyer
I have devised the following Excel Formula:

=(-"ISADEV"/-0,0019812)*LN((1+-0,0019812*"MDH")/(288,15+-0,0019812*"APAlt"))

where "ISADEV" is the temperature deviation from ISA (usually a negative number in cold weather)
"MDH" is the minimum descent height in ft
"APAlt" is the airport altitude in ft
-0,0019812 is the lapse rate in K/ft and 288,15 is ISA/sea level in K.

This results in meaningless results.

I am challenged by not really understanding why the basic formula should work and would be grateful for any hints on where I might have erred in applying it.
A quick examination of this formula reveals that for MDH over about 500 ft, the term (1+-0,0019812*"MDH") in your spreadsheet formula turns negative. The natural log (LN) of a negative number is undefined (unless you are an electrical engineer ).
EEngr is offline  
Old 9th Jan 2019, 17:53
  #14 (permalink)  
 
Join Date: Jan 2019
Location: Hamburg
Posts: 2
Likes: 0
Received 0 Likes on 0 Posts
Good news if anyone is interested....
I found the bug, the accompaning pdf documentation from Eurocontrol folowing the link in post #12 shows the correct equation.

=(-"ISADEV"/-0,0019812)*LN(1+(-0,0019812*"MDH")/(288,15+-0,0019812*"APAlt"))

the 1+ is prior the fraction... Sombody within ICAO seems not to be aware of maths basics ;-)
Like that, it works fine...

BG
Lumbanraja is offline  
Old 25th Sep 2022, 20:45
  #15 (permalink)  
Thread Starter
 
Join Date: Jun 2001
Location: at the edge of the alps
Posts: 447
Received 0 Likes on 0 Posts
ICAO has acknowledged the error in the meantime but given the glacial pace of publication it might take a few more years for a corrected version to reach the public.
Alpine Flyer is offline  

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are Off
Refbacks are Off


Thread Tools
Search this Thread

Contact Us - Archive - Advertising - Cookie Policy - Privacy Statement - Terms of Service

Copyright © 2024 MH Sub I, LLC dba Internet Brands. All rights reserved. Use of this site indicates your consent to the Terms of Use.