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 $ |
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 |
|
|
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' |
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 |
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 |
|
|