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...
Denk dat dit wel OK is, het zou nog wel een stukje beter kunnen door te werken met geopotentiele hoogtes ipv de standaardatmosfeer, maar dat zal hooguit enkele procenten verschil maken denk ik.
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)'
Ziet er ook goed uit. Dacht eerst even dat het misging omdat je hier alleen snelheidsschering meeneemt, maar dit is juist.
Ik weet niet met wat voor kaartjes je het vergelijkt, maar misschien dat zij de BRN shear verkeerd berekenen? Nemen zij misschien wel de directionele schering mee?
Groeten,
Bas
Quote selectie