/[lmdze]/trunk/phylmd/vdif_kcay.f
ViewVC logotype

Diff of /trunk/phylmd/vdif_kcay.f

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

revision 3 by guez, Wed Feb 27 13:16:39 2008 UTC revision 12 by guez, Mon Jul 21 16:05:07 2008 UTC
# Line 31  c kn : diffusivite turbulente des scalai Line 31  c kn : diffusivite turbulente des scalai
31  c      (en sortie : la valeur a la fin du pas de temps)  c      (en sortie : la valeur a la fin du pas de temps)
32  c  c
33  c.......................................................................  c.......................................................................
34        REAL dt,g,rconst        REAL, intent(in):: dt
35          real g,rconst
36        real plev(klon,klev+1),temp(klon,klev)        real plev(klon,klev+1),temp(klon,klev)
37        real ustar(klon),snstable        real ustar(klon),snstable
38        REAL zlev(klon,klev+1)        REAL zlev(klon,klev+1)
# Line 226  c....................................... Line 227  c.......................................
227  c  c
228  c   Initialisation de q2  c   Initialisation de q2
229    
230        call yamada(ngrid,dt,g,rconst,plev,temp        call yamada(ngrid,g,rconst,plev,temp
231       s   ,zlev,zlay,u,v,teta,cd,q2diag,km,kn,ustar       s   ,zlev,zlay,u,v,teta,cd,q2diag,km,kn,ustar
232       s   ,l_mix)       s   ,l_mix)
233        if (first.and.1.eq.1) then        if (first.and.1.eq.1) then
# Line 654  c  CALCUL DE LA DIFFUSION VERTICALE DE Q Line 655  c  CALCUL DE LA DIFFUSION VERTICALE DE Q
655             sss=sss+plev(1,ilev-1)-plev(1,ilev+1)             sss=sss+plev(1,ilev-1)-plev(1,ilev+1)
656             sssq=sssq+(plev(1,ilev-1)-plev(1,ilev+1))*q2(1,ilev)             sssq=sssq+(plev(1,ilev-1)-plev(1,ilev+1))*q2(1,ilev)
657          enddo          enddo
 c        print*,'Q2moy avant',sssq/sss  
 c       print*,'Q2q20 ',(q2(1,ilev),ilev=1,10)  
 c       print*,'Q2km0 ',(km(1,ilev),ilev=1,10)  
 c   ! C'est quoi ca qu'etait dans l'original???  
 c       do igrid=1,ngrid  
 c          q2(igrid,1)=10.  
 c       enddo  
 c        q2s=q2  
 c       do iii=1,10  
 c       call vdif_q2(dt,g,rconst,plev,temp,km,q2)  
 c       do ilev=1,klev+1  
 c          write(iii+49,*) q2(1,ilev),zlev(1,ilev)  
 c       enddo  
 c       enddo  
 c       stop  
 c       do ilev=1,klev  
 c          print*,zlev(1,ilev),q2s(1,ilev),q2(1,ilev)  
 c       enddo  
 c        q2s=q2-q2s  
 c       do ilev=1,klev  
 c          print*,q2s(1,ilev),zlev(1,ilev)  
 c       enddo  
658          do ilev=2,klev-1          do ilev=2,klev-1
659             sss=sss+plev(1,ilev-1)-plev(1,ilev+1)             sss=sss+plev(1,ilev-1)-plev(1,ilev+1)
660             sssq=sssq+(plev(1,ilev-1)-plev(1,ilev+1))*q2(1,ilev)             sssq=sssq+(plev(1,ilev-1)-plev(1,ilev+1))*q2(1,ilev)
Legend:
Removed from v.3  
changed lines
  Added in v.12
  ViewVC Help
Powered by ViewVC 1.1.21