Directory: | ./ |
---|---|
File: | dyn3d_common/tourpot.f90 |
Date: | 2022-01-11 19:19:34 |
Exec | Total | Coverage | |
---|---|---|---|
Lines: | 11 | 11 | 100.0% |
Branches: | 12 | 12 | 100.0% |
Line | Branch | Exec | Source |
---|---|---|---|
1 | 2881 | SUBROUTINE tourpot ( vcov, ucov, massebxy, vorpot ) | |
2 | ! | ||
3 | !------------------------------------------------------------------------------- | ||
4 | ! Authors: P. Le Van. | ||
5 | !------------------------------------------------------------------------------- | ||
6 | ! Purpose: Compute potential vorticity. | ||
7 | IMPLICIT NONE | ||
8 | include "dimensions.h" | ||
9 | include "paramet.h" | ||
10 | include "comgeom.h" | ||
11 | !=============================================================================== | ||
12 | ! Arguments: | ||
13 | REAL, INTENT(IN) :: vcov (ip1jm, llm) | ||
14 | REAL, INTENT(IN) :: ucov (ip1jmp1,llm) | ||
15 | REAL, INTENT(IN) :: massebxy(ip1jm, llm) | ||
16 | REAL, INTENT(OUT) :: vorpot (ip1jm, llm) | ||
17 | !=============================================================================== | ||
18 | ! Method used: | ||
19 | ! vorpot = ( Filtre( d(vcov)/dx - d(ucov)/dy ) + fext ) /psbarxy | ||
20 | !=============================================================================== | ||
21 | ! Local variables: | ||
22 | INTEGER :: l, ij | ||
23 | REAL :: rot(ip1jm,llm) | ||
24 | !=============================================================================== | ||
25 | |||
26 | !--- Wind vorticity ; correction: rot(iip1,j,l) = rot(1,j,l) | ||
27 |
2/2✓ Branch 0 taken 112359 times.
✓ Branch 1 taken 2881 times.
|
115240 | DO l=1,llm |
28 |
2/2✓ Branch 0 taken 118538745 times.
✓ Branch 1 taken 112359 times.
|
118651104 | DO ij=1,ip1jm-1 |
29 | 118651104 | rot(ij,l)=vcov(ij+1,l)-vcov(ij,l)+ucov(ij+iip1,l)-ucov(ij,l) | |
30 | END DO | ||
31 |
2/2✓ Branch 0 taken 3483129 times.
✓ Branch 1 taken 112359 times.
|
3598369 | DO ij=iip1,ip1jm,iip1; rot(ij,l)=rot(ij-iim,l); END DO |
32 | END DO | ||
33 | |||
34 | !--- Filter | ||
35 | 2881 | CALL filtreg(rot,jjm,llm,2,1,.FALSE.,1) | |
36 | |||
37 | !--- Potential vorticity ; correction: rot(iip1,j,l) = rot(1,j,l) | ||
38 |
2/2✓ Branch 0 taken 112359 times.
✓ Branch 1 taken 2881 times.
|
115240 | DO l=1,llm |
39 |
2/2✓ Branch 0 taken 118538745 times.
✓ Branch 1 taken 112359 times.
|
118651104 | DO ij=1,ip1jm-1 |
40 | 118651104 | vorpot(ij,l)=(rot(ij,l)+fext(ij))/massebxy(ij,l) | |
41 | END DO | ||
42 |
2/2✓ Branch 0 taken 3483129 times.
✓ Branch 1 taken 112359 times.
|
3598369 | DO ij=iip1,ip1jm,iip1; vorpot(ij,l)=vorpot(ij-iim,l); END DO |
43 | END DO | ||
44 | |||
45 | 2881 | END SUBROUTINE tourpot | |
46 |