440 |
enddo |
enddo |
441 |
|
|
442 |
do i = 1,irun |
do i = 1,irun |
443 |
robar(i) = robar(i)/(ps(i)-ple(i,Lm+1-nbase(i))) * 100.0 |
robar(i) = robar(i)/(ps(i)-ple(i,Lm+1-nbase(i)-1)) * 100.0 |
444 |
ubar(i) = ubar(i)/(ps(i)-ple(i,Lm+1-nbase(i))) |
ubar(i) = ubar(i)/(ps(i)-ple(i,Lm+1-nbase(i)-1)) |
445 |
vbar(i) = vbar(i)/(ps(i)-ple(i,Lm+1-nbase(i))) |
vbar(i) = vbar(i)/(ps(i)-ple(i,Lm+1-nbase(i)-1)) |
446 |
|
|
447 |
speed(i) = sqrt( ubar(i)*ubar(i) + vbar(i)*vbar(i) ) |
speed(i) = sqrt( ubar(i)*ubar(i) + vbar(i)*vbar(i) ) |
448 |
ang(i) = atan2(vbar(i),ubar(i)) |
ang(i) = atan2(vbar(i),ubar(i)) |