Page 1 of 1

Temperature Inversion

Posted: Sun Jan 24, 2016 12:23 pm
by nemzamp
Hello guys,
I would know if someone has a script for grads for inversion ... I would need to know the temperature difference per 100m to see if there is temperature inversion , for example :
0m to 6 ° C
100m 5 ° C -1
200m 7 ° C +2
300m 4.5 ° C -2.5
400m 3 ° C -1.5
etc .
Thanks for your help

Re: Temperature Inversion

Posted: Wed Jan 27, 2016 4:40 pm
by meteoadriatic
Hello, for a single point, did this quickly:

Code: Select all

'set lev 45'
'set lon 15'
z = 1
while (z <=40)
'set z 'z
lev=subwrd(result,4)
'define tmp1=tmpprs(z='z')'
'define tmp2=tmpprs(z='z+1')'
'define inver=tmp2-tmp1'
'd inver'
inversion=subwrd(result,4)
if inversion>0
say "inversion found in lev "lev
endif
z=z+1
endwhile
Credits here:
http://gradsusr.org/pipermail/gradsusr/ ... 32309.html

Adapt more to better suit your needs.

Re: Temperature Inversion

Posted: Wed Feb 10, 2016 7:19 pm
by nemzamp
Thank you