/[lmdze]/trunk/Sources/dyn3d/tourpot.f
ViewVC logotype

Diff of /trunk/Sources/dyn3d/tourpot.f

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 82 by guez, Wed Mar 5 14:57:53 2014 UTC revision 107 by guez, Thu Sep 11 15:09:15 2014 UTC
# Line 11  contains Line 11  contains
11      ! Author: P. Le Van      ! Author: P. Le Van
12      ! Objet : calcul du tourbillon potentiel      ! Objet : calcul du tourbillon potentiel
13    
     USE dimens_m, ONLY: iim, jjm, llm  
14      USE comgeom, ONLY: fext_2d      USE comgeom, ONLY: fext_2d
15        USE dimens_m, ONLY: iim, jjm, llm
16      use filtreg_m, only: filtreg      use filtreg_m, only: filtreg
17    
18      REAL, intent(in):: vcov(:, :, :) ! (iim + 1, jjm, llm)      REAL, intent(in):: vcov(:, :, :) ! (iim + 1, jjm, llm)
# Line 23  contains Line 23  contains
23      ! = (Filtre(d(vcov)/dx - d(ucov)/dy) + fext) / massebxy      ! = (Filtre(d(vcov)/dx - d(ucov)/dy) + fext) / massebxy
24    
25      ! Local:      ! Local:
26    
27      REAL rot(iim + 1, jjm, llm)      REAL rot(iim + 1, jjm, llm)
28      ! relative vorticity multiplied by cell area, in m2 s-1      ! relative vorticity multiplied by cell area, in m2 s-1
29    
# Line 36  contains Line 37  contains
37           = vcov(i + 1, j, :) - vcov(i, j, :) + ucov(i, j + 1, :) - ucov(i, j, :)           = vcov(i + 1, j, :) - vcov(i, j, :) + ucov(i, j + 1, :) - ucov(i, j, :)
38      rot(iim + 1, :, :) = rot(1, :, :)      rot(iim + 1, :, :) = rot(1, :, :)
39    
40      CALL filtreg(rot, jjm, llm, 2, 1, .FALSE.)      CALL filtreg(rot, direct = .true., intensive = .true.)
41    
42      forall (l = 1: llm) vorpot(:iim, :, l) &      forall (l = 1: llm) vorpot(:iim, :, l) &
43           = (rot(:iim, :, l) + fext_2d(:iim, :)) / massebxy(:iim, :, l)           = (rot(:iim, :, l) + fext_2d(:iim, :)) / massebxy(:iim, :, l)
Legend:
Removed from v.82  
changed lines
  Added in v.107
  ViewVC Help
Powered by ViewVC 1.1.21