/[lmdze]/trunk/dyn3d/Read_reanalyse/reanalyse2nat.f
ViewVC logotype

Diff of /trunk/dyn3d/Read_reanalyse/reanalyse2nat.f

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

revision 84 by guez, Thu Mar 6 15:12:00 2014 UTC revision 85 by guez, Thu Mar 6 17:35:22 2014 UTC
# Line 1  Line 1 
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.
# Line 21  subroutine reanalyse2nat(nlevnc,psi & Line 19  subroutine reanalyse2nat(nlevnc,psi &
19    integer nlevnc    integer nlevnc
20    real, intent(in):: 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)
# Line 96  subroutine reanalyse2nat(nlevnc,psi & Line 94  subroutine reanalyse2nat(nlevnc,psi &
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
Legend:
Removed from v.84  
changed lines
  Added in v.85
  ViewVC Help
Powered by ViewVC 1.1.21