/[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 17 by guez, Tue Aug 5 13:31:32 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, intent(in):: g
36          real rconst
37        real plev(klon,klev+1),temp(klon,klev)        real plev(klon,klev+1),temp(klon,klev)
38        real ustar(klon),snstable        real ustar(klon),snstable
39        REAL zlev(klon,klev+1)        REAL zlev(klon,klev+1)
# Line 226  c....................................... Line 228  c.......................................
228  c  c
229  c   Initialisation de q2  c   Initialisation de q2
230    
231        call yamada(ngrid,dt,g,rconst,plev,temp        call yamada(ngrid,g,rconst,plev,temp
232       s   ,zlev,zlay,u,v,teta,cd,q2diag,km,kn,ustar       s   ,zlev,zlay,u,v,teta,cd,q2diag,km,kn,ustar
233       s   ,l_mix)       s   ,l_mix)
234        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 656  c  CALCUL DE LA DIFFUSION VERTICALE DE Q
656             sss=sss+plev(1,ilev-1)-plev(1,ilev+1)             sss=sss+plev(1,ilev-1)-plev(1,ilev+1)
657             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)
658          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  
659          do ilev=2,klev-1          do ilev=2,klev-1
660             sss=sss+plev(1,ilev-1)-plev(1,ilev+1)             sss=sss+plev(1,ilev-1)-plev(1,ilev+1)
661             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.17
  ViewVC Help
Powered by ViewVC 1.1.21