/[lmdze]/trunk/dyn3d/tourpot.f90
ViewVC logotype

Diff of /trunk/dyn3d/tourpot.f90

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

trunk/Sources/dyn3d/tourpot.f revision 134 by guez, Wed Apr 29 15:47:56 2015 UTC trunk/dyn3d/tourpot.f revision 254 by guez, Mon Feb 5 10:39:38 2018 UTC
# Line 6  contains Line 6  contains
6    
7    SUBROUTINE tourpot(vcov, ucov, massebxy, vorpot)    SUBROUTINE tourpot(vcov, ucov, massebxy, vorpot)
8    
9      ! From LMDZ4/libf/dyn3d/tourpot.F, version 1.1.1.1 2004/05/19 12:53:06      ! From LMDZ4/libf/dyn3d/tourpot.F, version 1.1.1.1, 2004/05/19 12:53:06
10    
11      ! Author: P. Le Van      ! Author: P. Le Van
12      ! Objet : calcul du tourbillon potentiel      ! Objet : calcul du tourbillon potentiel
13    
14      USE comgeom, ONLY: fext_2d      USE comgeom, ONLY: fext_2d
15      USE dimens_m, ONLY: iim, jjm, llm      USE dimens_m, ONLY: iim, jjm, llm
16      use filtreg_m, only: filtreg      use filtreg_v_m, only: filtreg_v
17    
18      REAL, intent(in):: vcov(:, :, :) ! (iim + 1, jjm, llm)      REAL, intent(in):: vcov(:, :, :) ! (iim + 1, jjm, llm)
19      REAL, intent(in):: ucov(:, :, :) ! (iim + 1, jjm + 1, llm)      REAL, intent(in):: ucov(:, :, :) ! (iim + 1, jjm + 1, llm)
# Line 37  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, direct = .true., intensive = .true.)      CALL filtreg_v(rot, 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.134  
changed lines
  Added in v.254
  ViewVC Help
Powered by ViewVC 1.1.21