/[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 76 by guez, Fri Nov 15 18:45:49 2013 UTC revision 78 by guez, Wed Feb 5 17:51:07 2014 UTC
# Line 19  contains Line 19  contains
19      USE ener, ONLY: ang, ang0, etot, etot0, ptot, ptot0, rmsdpdt, rmsv, &      USE ener, ONLY: ang, ang0, etot, etot0, ptot, ptot0, rmsdpdt, rmsv, &
20           stot, stot0, ztot, ztot0           stot, stot0, ztot, ztot0
21      use filtreg_m, only: filtreg      use filtreg_m, only: filtreg
22        use massbarxy_m, only: massbarxy
23      USE paramet_m, ONLY: iip1, iip2, ijp1llm, ip1jm, ip1jmp1, jjp1      USE paramet_m, ONLY: iip1, iip2, ijp1llm, ip1jm, ip1jmp1, jjp1
24    
25      INTEGER, INTENT(IN):: itau      INTEGER, INTENT(IN):: itau
26      REAL, INTENT(IN):: ucov(ip1jmp1, llm)      REAL, INTENT(IN):: ucov(ip1jmp1, llm)
27      real, intent(in):: teta(ip1jmp1, llm)      REAL, INTENT(IN):: teta(ip1jmp1, llm)
28      REAL, INTENT(IN):: ps(ip1jmp1)      REAL, INTENT(IN):: ps(ip1jmp1)
29      REAL, INTENT(IN):: masse(ip1jmp1, llm)      REAL, INTENT(IN):: masse(ip1jmp1, llm)
30      REAL, INTENT (IN):: pk(ip1jmp1, llm)      REAL, INTENT(IN):: pk(ip1jmp1, llm)
31      REAL, INTENT(IN):: phis(ip1jmp1)      REAL, INTENT(IN):: phis(ip1jmp1)
32      REAL, INTENT(IN):: vorpot(ip1jm, llm)      REAL, INTENT(IN):: vorpot(ip1jm, llm)
33      REAL, intent(in):: phi(ip1jmp1, llm)      REAL, intent(in):: phi(ip1jmp1, llm)
# Line 38  contains Line 39  contains
39      REAL:: vor(ip1jm), bernf(ip1jmp1, llm), ztotl(llm)      REAL:: vor(ip1jm), bernf(ip1jmp1, llm), ztotl(llm)
40      REAL:: etotl(llm), stotl(llm), rmsvl(llm), angl(llm), ge(ip1jmp1)      REAL:: etotl(llm), stotl(llm), rmsvl(llm), angl(llm), ge(ip1jmp1)
41      REAL:: cosphi(ip1jm), omegcosp(ip1jm)      REAL:: cosphi(ip1jm), omegcosp(ip1jm)
42      REAL:: dtvrs1j, rjour, heure, radsg, radomeg      REAL dtvrs1j, rjour, heure, radsg, radomeg
43      real massebxy(ip1jm, llm)      REAL massebxy(ip1jm, llm)
44      INTEGER:: l, ij, imjmp1      INTEGER l, ij
45      REAL:: ssum      REAL ssum
46      real time      real time
47    
48      !-----------------------------------------------------------------------      !-----------------------------------------------------------------------
49    
50      print *, "Call sequence information: sortvarc"      PRINT *, "Call sequence information: sortvarc"
51    
52      time = real(itau) / day_step + time_0      time = real(itau) / day_step + time_0
53      dtvrs1j = dtvr/daysec      dtvrs1j = dtvr/daysec
54      rjour = real(int(itau*dtvrs1j))      rjour = real(int(itau*dtvrs1j))
55      heure = (itau*dtvrs1j-rjour)*24.      heure = (itau*dtvrs1j-rjour)*24.
     imjmp1 = iim*jjp1  
56      IF (abs(heure-24.)<=0.0001) heure = 0.      IF (abs(heure-24.)<=0.0001) heure = 0.
57    
58      CALL massbarxy(masse, massebxy)      CALL massbarxy(masse, massebxy)
59    
60      ! Calcul  de  rmsdpdt      ! Calcul  de  rmsdpdt
61      ge(:) = dp(:)*dp(:)      ge = dp*dp
62      rmsdpdt = ssum(ip1jmp1, ge, 1) - ssum(jjp1, ge, iip1)      rmsdpdt = ssum(ip1jmp1, ge, 1) - ssum(jjp1, ge, iip1)
63      rmsdpdt = daysec*1.E-2*sqrt(rmsdpdt/imjmp1)      rmsdpdt = daysec*1.E-2*sqrt(rmsdpdt / (iim * jjp1))
64      bernf = bern      bernf = bern
65      CALL filtreg(bernf, jjp1, llm, -2, 2, .TRUE.)      CALL filtreg(bernf, jjp1, llm, -2, 2, .TRUE.)
66    
# Line 73  contains Line 73  contains
73      END DO      END DO
74    
75      ! Calcul  de l'energie, de l'enstrophie, de l'entropie et de rmsv      ! Calcul  de l'energie, de l'enstrophie, de l'entropie et de rmsv
76    
77      DO l = 1, llm      DO l = 1, llm
78         DO ij = 1, ip1jm         DO ij = 1, ip1jm
79            vor(ij) = vorpot(ij, l)*vorpot(ij, l)*massebxy(ij, l)            vor(ij) = vorpot(ij, l)*vorpot(ij, l)*massebxy(ij, l)
Legend:
Removed from v.76  
changed lines
  Added in v.78
  ViewVC Help
Powered by ViewVC 1.1.21