1 |
subroutine reanalyse2nat(nlevnc,psi & |
subroutine reanalyse2nat(nlevnc,psi,unc,vnc,tnc,qnc,psnc,pl,u,v,t,q ,masse,pk) |
|
,unc,vnc,tnc,qnc,psnc,pl,u,v,t,q & |
|
|
,ps,masse,pk) |
|
2 |
|
|
3 |
! Inversion Nord/sud de la grille + interpollation sur les niveaux |
! Inversion Nord/sud de la grille + interpollation sur les niveaux |
4 |
! verticaux du modele. |
! verticaux du modele. |
17 |
|
|
18 |
|
|
19 |
integer nlevnc |
integer nlevnc |
20 |
real psi(iip1,jjp1) |
real, intent(in):: psi(iip1,jjp1) |
21 |
real u(iip1,jjp1,llm),v(iip1,jjm,llm) |
real u(iip1,jjp1,llm),v(iip1,jjm,llm) |
22 |
real t(iip1,jjp1,llm),ps(iip1,jjp1),q(iip1,jjp1,llm) |
real t(iip1,jjp1,llm), q(iip1,jjp1,llm) |
23 |
|
|
24 |
real pl(nlevnc) |
real pl(nlevnc) |
25 |
real unc(iip1,jjp1,nlevnc),vnc(iip1,jjm,nlevnc) |
real unc(iip1,jjp1,nlevnc),vnc(iip1,jjm,nlevnc) |
94 |
|
|
95 |
! ----------------------------------------------------------------- |
! ----------------------------------------------------------------- |
96 |
|
|
|
if (guide_P) then |
|
|
do j=1,jjp1 |
|
|
do i=1,iim |
|
|
ps(i,j)=psnc(i,jjp1+1-j) |
|
|
enddo |
|
|
ps(iip1,j)=ps(1,j) |
|
|
enddo |
|
|
endif |
|
|
|
|
|
|
|
|
! ----------------------------------------------------------------- |
|
97 |
call pres2lev(unc,zu,nlevnc,llm,pl,plunc,iip1,jjp1) |
call pres2lev(unc,zu,nlevnc,llm,pl,plunc,iip1,jjp1) |
98 |
call pres2lev(vnc,zv,nlevnc,llm,pl,plvnc,iip1,jjm ) |
call pres2lev(vnc,zv,nlevnc,llm,pl,plvnc,iip1,jjm ) |
99 |
call pres2lev(tnc,zt,nlevnc,llm,pl,plsnc,iip1,jjp1) |
call pres2lev(tnc,zt,nlevnc,llm,pl,plsnc,iip1,jjp1) |
100 |
call pres2lev(qnc,zq,nlevnc,llm,pl,plsnc,iip1,jjp1) |
call pres2lev(qnc,zq,nlevnc,llm,pl,plsnc,iip1,jjp1) |
101 |
|
|
|
! call dump2d(iip1,jjp1,ps,'PS ') |
|
|
! call dump2d(iip1,jjp1,psu,'PS ') |
|
|
! call dump2d(iip1,jjm,psv,'PS ') |
|
102 |
! Inversion Nord/Sud |
! Inversion Nord/Sud |
103 |
do l=1,llm |
do l=1,llm |
104 |
do j=1,jjp1 |
do j=1,jjp1 |