/[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 33 by guez, Fri Apr 9 10:56:14 2010 UTC revision 64 by guez, Wed Aug 29 14:47:17 2012 UTC
# Line 22  contains Line 22  contains
22      use filtreg_m, only: filtreg      use filtreg_m, only: filtreg
23    
24      ! Arguments:      ! Arguments:
25      INTEGER, INTENT (IN) :: itau      INTEGER, INTENT(IN) :: itau
26      REAL :: ucov(ip1jmp1, llm), teta(ip1jmp1, llm), masse(ip1jmp1, llm)      REAL :: ucov(ip1jmp1, llm), masse(ip1jmp1, llm)
27      REAL :: ps(ip1jmp1), phis(ip1jmp1)      real, intent(in):: teta(ip1jmp1, llm)
28        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 :: phi(ip1jmp1, llm), bern(ip1jmp1, llm)
31      REAL :: dp(ip1jmp1)      REAL :: dp(ip1jmp1)
# Line 43  contains Line 44  contains
44    
45      !-----------------------------------------------------------------------      !-----------------------------------------------------------------------
46    
47        print *, "Call sequence information: sortvarc"
48    
49      time = real(itau) / day_step + time_0      time = real(itau) / day_step + time_0
50      dtvrs1j = dtvr/daysec      dtvrs1j = dtvr/daysec
51      rjour = real(int(itau*dtvrs1j))      rjour = real(int(itau*dtvrs1j))
# Line 57  contains Line 60  contains
60      rmsdpdt = ssum(ip1jmp1, ge, 1) - ssum(jjp1, ge, iip1)      rmsdpdt = ssum(ip1jmp1, ge, 1) - ssum(jjp1, ge, iip1)
61      rmsdpdt = daysec*1.E-2*sqrt(rmsdpdt/imjmp1)      rmsdpdt = daysec*1.E-2*sqrt(rmsdpdt/imjmp1)
62      CALL scopy(ijp1llm, bern, 1, bernf, 1)      CALL scopy(ijp1llm, bern, 1, bernf, 1)
63      CALL filtreg(bernf, jjp1, llm, -2, 2, .TRUE., 1)      CALL filtreg(bernf, jjp1, llm, -2, 2, .TRUE.)
64    
65      ! Calcul du moment  angulaire      ! Calcul du moment  angulaire
66      radsg = rad/g      radsg = rad/g
Legend:
Removed from v.33  
changed lines
  Added in v.64
  ViewVC Help
Powered by ViewVC 1.1.21