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

Diff of /trunk/dyn3d/Read_reanalyse/read_reanalyse.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 read_reanalyse(timestep,psi &        subroutine read_reanalyse(timestep,psi,u,v,t,q,masse,mode,nlevnc)
          ,u,v,t,q,masse,ps,mode,nlevnc)  
2    
3  !  !
4  ! $Header: /home/cvsroot/LMDZ4/libf/dyn3d/read_reanalyse.F,v 1.3 2005/04/15 12:31:21 lmdzadmin Exp $  ! $Header: /home/cvsroot/LMDZ4/libf/dyn3d/read_reanalyse.F,v 1.3 2005/04/15 12:31:21 lmdzadmin Exp $
# Line 31  Line 30 
30    
31        real, intent(in):: psi(iip1,jjp1)        real, intent(in):: psi(iip1,jjp1)
32        real u(iip1,jjp1,llm),v(iip1,jjm,llm)        real u(iip1,jjp1,llm),v(iip1,jjm,llm)
33        real t(iip1,jjp1,llm),ps(iip1,jjp1),q(iip1,jjp1,llm)        real t(iip1,jjp1,llm),q(iip1,jjp1,llm)
34        real masse(iip1,jjp1,llm),pk(iip1,jjp1,llm)        real masse(iip1,jjp1,llm),pk(iip1,jjp1,llm)
35    
36    
# Line 98  Line 97 
97              if (ncidpl.eq.-99) ncidpl=ncidQ              if (ncidpl.eq.-99) ncidpl=ncidQ
98              endif              endif
99    
 ! Pression de surface  
             if (guide_P) then  
             rcode=nf90_open('ps.nc',nf90_nowrite,ncidps)  
             rcode = nf90_inq_varid(ncidps, 'SP', varidps)  
             print*,'ncidps,varidps',ncidps,varidps  
             endif  
   
100  ! Coordonnee verticale  ! Coordonnee verticale
101              if (ncep) then              if (ncep) then
102                 print*,'Vous etes entrain de lire des donnees NCEP'                 print*,'Vous etes entrain de lire des donnees NCEP'
# Line 130  Line 122 
122        endif        endif
123    
124  ! -----------------------------------------------------------------  ! -----------------------------------------------------------------
125  !   lecture des champs u, v, T, ps  !   lecture des champs u, v, T
126  ! -----------------------------------------------------------------  ! -----------------------------------------------------------------
127    
128  !  dimensions pour les champs scalaires et le vent zonal  !  dimensions pour les champs scalaires et le vent zonal
# Line 204  Line 196 
196        count(3)=1        count(3)=1
197        count(4)=0        count(4)=0
198    
 !  Pression de surface  
 !  -------------------  
   
       if (guide_P) then  
       status=NF90_GET_VAR(ncidps,varidps,psnc,start,count)  
       call dump2d(iip1,jjp1,psnc,'PSNC COUCHE 1 ')  
       call correctbid(iim,jjp1,psnc)  
       endif  
   
   
   
 ! -----------------------------------------------------------------  
199  !  Interpollation verticale sur les niveaux modele  !  Interpollation verticale sur les niveaux modele
200  ! -----------------------------------------------------------------  ! -----------------------------------------------------------------
201        call reanalyse2nat(nlevnc,psi,unc,vnc,tnc,Qnc,psnc,pl,u,v,t,Q &        call reanalyse2nat(nlevnc,psi,unc,vnc,tnc,Qnc,psnc,pl,u,v,t,Q,masse,pk)
           ,ps,masse,pk)  
202    
203        call dump2d(iip1,jjm,v,'V COUCHE APRES ')        call dump2d(iip1,jjm,v,'V COUCHE APRES ')
204    
Legend:
Removed from v.84  
changed lines
  Added in v.85
  ViewVC Help
Powered by ViewVC 1.1.21