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

Diff of /trunk/dyn3d/tourpot.f

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

revision 135 by guez, Thu Apr 30 14:22:32 2015 UTC revision 137 by guez, Wed May 6 15:51:03 2015 UTC
# Line 13  contains Line 13  contains
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.135  
changed lines
  Added in v.137

  ViewVC Help
Powered by ViewVC 1.1.21