wel, de BRN shear zou het verschil moeten zijn tussen de gem.gewogen wind in de0-6km laag en de 0-500m laag. Daarvan het kwadraat en de helft
BRNshear= 0.5*(DLshear-LLshear)²
Maar m'n poging om dit te simuleren geven andere resultaten dan op vergelijkbare kaarten...Dat gewogen heb ik opgelost door gewoon de densiteitverdeling van de standard atmosfeer te nemen en dit te normaliseren (dus op 10m geeft die een factor *1, op 950hpa *0.952 enz...maar ben al niet zeker of dit ok is...
hieronder grads-poging...
'define avgLu = (u10+(u950*0.952))/2'
'define avgLv = (v10+(v950*0.952))/2'
'define avgDu = (u10+(u950*0.952)+(u900*0.912)+(u850*0.871)+(u800*0.8286)+(u750*0.787)+(u700*0.7451)+(u650*0.7016)+(u600*0.657)+(u550*0.6123)+(u500*0.5664)+(u450*0.521))/12'
'define avgDv = (v10+(v950*0.952)+(v900*0.912)+(v850*0.871)+(v800*0.8286)+(v750*0.787)+(v700*0.7451)+(v650*0.7016)+(v600*0.657)+(v550*0.6123)+(v500*0.5664)+(v450*0.521))/12'
'define magnL = mag(avgLu,avgLv)'
'define magnD = mag(avgDu,avgDv)'
'define shear = magnD-magnL'
'define brnshear = 0.5*(shear*shear)'
Quote selectie