/[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 68 by guez, Wed Nov 14 16:59:30 2012 UTC revision 69 by guez, Mon Feb 18 16:33:12 2013 UTC
# Line 11  contains Line 11  contains
11      ! Author: P. Le Van      ! Author: P. Le Van
12      ! Objet : sortie des variables de contrôle      ! Objet : sortie des variables de contrôle
13    
14        USE comconst, ONLY: daysec, dtvr, g, omeg, rad
15        USE comgeom, ONLY: aire, cu, rlatu
16      USE conf_gcm_m, ONLY: day_step      USE conf_gcm_m, ONLY: day_step
17      USE dimens_m, ONLY : iim, jjm, llm      USE dimens_m, ONLY: iim, jjm, llm
18      USE paramet_m, ONLY : iip1, iip2, ijp1llm, ip1jm, ip1jmp1, jjp1      USE dynetat0_m, ONLY: day_ini
19      USE comconst, ONLY : daysec, dtvr, g, omeg, rad      USE ener, ONLY: ang, ang0, etot, etot0, ptot, ptot0, rmsdpdt, rmsv, &
     USE comgeom, ONLY : aire, cu, rlatu  
     USE dynetat0_m, ONLY : day_ini  
     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 paramet_m, ONLY: iip1, iip2, ijp1llm, ip1jm, ip1jmp1, jjp1
23    
24      ! Arguments:      INTEGER, INTENT(IN):: itau
25      INTEGER, INTENT(IN) :: itau      REAL, INTENT(IN):: ucov(ip1jmp1, llm)
     REAL :: ucov(ip1jmp1, llm), masse(ip1jmp1, llm)  
26      real, intent(in):: teta(ip1jmp1, llm)      real, intent(in):: teta(ip1jmp1, llm)
27      REAL, INTENT(IN):: ps(ip1jmp1), phis(ip1jmp1)      REAL, INTENT(IN):: ps(ip1jmp1)
28      REAL :: vorpot(ip1jm, llm)      REAL, INTENT(IN):: masse(ip1jmp1, llm)
29        REAL, INTENT (IN):: pk(ip1jmp1, llm)
30        REAL, INTENT(IN):: phis(ip1jmp1)
31        REAL, INTENT(IN):: vorpot(ip1jm, llm)
32      REAL, intent(in):: phi(ip1jmp1, llm)      REAL, intent(in):: phi(ip1jmp1, llm)
33      real bern(ip1jmp1, llm)      real, intent(in):: bern(ip1jmp1, llm)
34      REAL :: dp(ip1jmp1)      REAL, intent(in):: dp(ip1jmp1)
35      REAL, INTENT (IN):: time_0      REAL, INTENT (IN):: time_0
     REAL, INTENT (IN):: pk(ip1jmp1, llm)  
36    
37      ! Local:      ! Local:
38      REAL :: vor(ip1jm), bernf(ip1jmp1, llm), ztotl(llm)      REAL:: vor(ip1jm), bernf(ip1jmp1, llm), ztotl(llm)
39      REAL :: etotl(llm), stotl(llm), rmsvl(llm), angl(llm), ge(ip1jmp1)      REAL:: etotl(llm), stotl(llm), rmsvl(llm), angl(llm), ge(ip1jmp1)
40      REAL :: cosphi(ip1jm), omegcosp(ip1jm)      REAL:: cosphi(ip1jm), omegcosp(ip1jm)
41      REAL :: dtvrs1j, rjour, heure, radsg, radomeg      REAL:: dtvrs1j, rjour, heure, radsg, radomeg
42      real massebxy(ip1jm, llm)      real massebxy(ip1jm, llm)
43      INTEGER :: l, ij, imjmp1      INTEGER:: l, ij, imjmp1
44      REAL :: ssum      REAL:: ssum
45      real time      real time
46    
47      !-----------------------------------------------------------------------      !-----------------------------------------------------------------------
# Line 60  contains Line 61  contains
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/imjmp1)
64      CALL scopy(ijp1llm, bern, 1, bernf, 1)      bernf = bern
65      CALL filtreg(bernf, jjp1, llm, -2, 2, .TRUE.)      CALL filtreg(bernf, jjp1, llm, -2, 2, .TRUE.)
66    
67      ! Calcul du moment  angulaire      ! Calcul du moment  angulaire
Legend:
Removed from v.68  
changed lines
  Added in v.69
  ViewVC Help
Powered by ViewVC 1.1.21