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

Diff of /trunk/dyn3d/sortvarc.f

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

revision 57 by guez, Mon Jan 30 12:54:02 2012 UTC revision 68 by guez, Wed Nov 14 16:59:30 2012 UTC
# Line 27  contains Line 27  contains
27      real, intent(in):: teta(ip1jmp1, llm)      real, intent(in):: teta(ip1jmp1, llm)
28      REAL, INTENT(IN):: ps(ip1jmp1), phis(ip1jmp1)      REAL, INTENT(IN):: ps(ip1jmp1), phis(ip1jmp1)
29      REAL :: vorpot(ip1jm, llm)      REAL :: vorpot(ip1jm, llm)
30      REAL :: phi(ip1jmp1, llm), bern(ip1jmp1, llm)      REAL, intent(in):: phi(ip1jmp1, llm)
31        real bern(ip1jmp1, llm)
32      REAL :: dp(ip1jmp1)      REAL :: dp(ip1jmp1)
33      REAL, INTENT (IN):: time_0      REAL, INTENT (IN):: time_0
34      REAL, INTENT (IN):: pk(ip1jmp1, llm)      REAL, INTENT (IN):: pk(ip1jmp1, llm)
# Line 60  contains Line 61  contains
61      rmsdpdt = ssum(ip1jmp1, ge, 1) - ssum(jjp1, ge, iip1)      rmsdpdt = ssum(ip1jmp1, ge, 1) - ssum(jjp1, ge, iip1)
62      rmsdpdt = daysec*1.E-2*sqrt(rmsdpdt/imjmp1)      rmsdpdt = daysec*1.E-2*sqrt(rmsdpdt/imjmp1)
63      CALL scopy(ijp1llm, bern, 1, bernf, 1)      CALL scopy(ijp1llm, bern, 1, bernf, 1)
64      CALL filtreg(bernf, jjp1, llm, -2, 2, .TRUE., 1)      CALL filtreg(bernf, jjp1, llm, -2, 2, .TRUE.)
65    
66      ! Calcul du moment  angulaire      ! Calcul du moment  angulaire
67      radsg = rad/g      radsg = rad/g
Legend:
Removed from v.57  
changed lines
  Added in v.68
  ViewVC Help
Powered by ViewVC 1.1.21