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) |
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 |
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) |