| 1 |
module tourabs_m |
| 2 |
|
| 3 |
IMPLICIT NONE |
| 4 |
|
| 5 |
contains |
| 6 |
|
| 7 |
SUBROUTINE tourabs(ntetaSTD, vcov, ucov, vorabs) |
| 8 |
|
| 9 |
! Author: I. Musat, 28 October 2004, adaptation de tourpot |
| 10 |
! Objet : calcul de la vorticité absolue |
| 11 |
|
| 12 |
USE dimens_m, ONLY: iim, jjm |
| 13 |
USE paramet_m, ONLY: iip1, ip1jm, ip1jmp1 |
| 14 |
USE comconst, ONLY: rad |
| 15 |
USE comgeom, ONLY: cu, cv, fext, rlatv, unsairez |
| 16 |
USE filtreg_m, ONLY: filtreg |
| 17 |
USE nr_util, ONLY: pi |
| 18 |
|
| 19 |
INTEGER, intent(in):: ntetaSTD |
| 20 |
REAL, intent(in):: vcov(ip1jm, ntetaSTD), ucov(ip1jmp1, ntetaSTD) |
| 21 |
|
| 22 |
REAL, intent(out):: vorabs(ip1jm, ntetaSTD) |
| 23 |
! vorabs = Filtre(d(vcov)/dx - d(ucov)/dy) + fext |
| 24 |
|
| 25 |
! Variables locales: |
| 26 |
INTEGER l, ij, i, j |
| 27 |
REAL rot(ip1jm, ntetaSTD) |
| 28 |
|
| 29 |
!-------------------------------------------------------------------- |
| 30 |
|
| 31 |
! Calcul du rotationnel du vent puis filtrage |
| 32 |
|
| 33 |
DO l = 1, ntetaSTD |
| 34 |
DO i = 1, iip1 |
| 35 |
DO j = 1, jjm |
| 36 |
ij = i + (j - 1) * iip1 |
| 37 |
IF (ij <= ip1jm - 1) THEN |
| 38 |
IF (cv(ij) == 0. .OR. cv(ij+1) == 0. .OR. cu(ij) == 0. & |
| 39 |
.OR. cu(ij+iip1) == 0.) THEN |
| 40 |
rot(ij, l) = 0. |
| 41 |
ELSE |
| 42 |
rot(ij, l) = (vcov(ij + 1, l) / cv(ij + 1) - vcov(ij, l) & |
| 43 |
/ cv(ij)) / (2. * pi * RAD * cos(rlatv(j))) & |
| 44 |
* real(iim) + (ucov(ij + iip1, l) / cu(ij + iip1) & |
| 45 |
- ucov(ij, l) / cu(ij)) / (pi * RAD) * (real(jjm) - 1.) |
| 46 |
ENDIF |
| 47 |
ENDIF |
| 48 |
end DO |
| 49 |
end DO |
| 50 |
|
| 51 |
! correction pour rot(iip1, j, l) |
| 52 |
! rot(iip1, j, l) = rot(1, j, l) |
| 53 |
DO ij = iip1, ip1jm, iip1 |
| 54 |
rot(ij, l) = rot(ij - iim, l) |
| 55 |
end DO |
| 56 |
end DO |
| 57 |
|
| 58 |
CALL filtreg(rot, jjm, ntetaSTD, 2, 1, .FALSE.) |
| 59 |
|
| 60 |
DO l = 1, ntetaSTD |
| 61 |
DO ij = 1, ip1jm - 1 |
| 62 |
vorabs(ij, l) = (rot(ij, l) + fext(ij) * unsairez(ij)) |
| 63 |
end DO |
| 64 |
|
| 65 |
! correction pour vorabs(iip1, j, l) |
| 66 |
! vorabs(iip1, j, l)= vorabs(1, j, l) |
| 67 |
DO ij = iip1, ip1jm, iip1 |
| 68 |
vorabs(ij, l) = vorabs(ij - iim, l) |
| 69 |
end DO |
| 70 |
end DO |
| 71 |
|
| 72 |
END SUBROUTINE tourabs |
| 73 |
|
| 74 |
end module tourabs_m |
Parent Directory
| 
Revision Log