/[lmdze]/trunk/libf/dyn3d/Read_reanalyse/read_reanalyse.f90
ViewVC logotype

Diff of /trunk/libf/dyn3d/Read_reanalyse/read_reanalyse.f90

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

revision 60 by guez, Wed Apr 13 12:29:18 2011 UTC revision 61 by guez, Fri Apr 20 14:58:43 2012 UTC
# Line 24  Line 24 
24  ! common  ! common
25  ! ------  ! ------
26    
       include "netcdf.inc"  
   
   
27  ! arguments  ! arguments
28  ! ---------  ! ---------
29        integer nlevnc        integer nlevnc
# Line 123  Line 120 
120  ! Niveaux de pression  ! Niveaux de pression
121        print*,'WARNING!!! Il n y a pas de test de coherence'        print*,'WARNING!!! Il n y a pas de test de coherence'
122        print*,'sur le nombre de niveaux verticaux dans le fichier nc'        print*,'sur le nombre de niveaux verticaux dans le fichier nc'
123        status=NF_GET_VARA_REAL(ncidpl,varidpl,1,nlevnc,pl)        status=NF90_GET_VAR(ncidpl,varidpl,pl)
124  !  passage en pascal  !  passage en pascal
125        pl(:)=100.*pl(:)        pl(:)=100.*pl(:)
126        if (first) then        if (first) then
# Line 161  Line 158 
158    
159        if (guide_u) then        if (guide_u) then
160        print*,'avant la lecture de UNCEP nd de niv:',nlevnc        print*,'avant la lecture de UNCEP nd de niv:',nlevnc
161        status=NF_GET_VARA_REAL(ncidu,varidu,start,count,unc)        status=NF90_GET_VAR(ncidu,varidu,unc,start,count)
162  !     call dump2d(iip1,jjp1,unc,'VENT NCEP   ')  !     call dump2d(iip1,jjp1,unc,'VENT NCEP   ')
163  !     call dump2d(iip1,40,unc(1,1,nlevnc),'VENT NCEP   ')  !     call dump2d(iip1,40,unc(1,1,nlevnc),'VENT NCEP   ')
164        print*,'WARNING!!! Correction bidon pour palier a un '        print*,'WARNING!!! Correction bidon pour palier a un '
# Line 176  Line 173 
173        print*,'ncidt=',ncidt,'varidt=',varidt,'start=',start        print*,'ncidt=',ncidt,'varidt=',varidt,'start=',start
174        print*,'count=',count        print*,'count=',count
175        if (guide_T) then        if (guide_T) then
176        status=NF_GET_VARA_REAL(ncidt,varidt,start,count,tnc)        status=NF90_GET_VAR(ncidt,varidt,tnc,start,count)
177        call dump2d(iip1,jjp1,tnc,'TNC COUCHE 1 AAA ')        call dump2d(iip1,jjp1,tnc,'TNC COUCHE 1 AAA ')
178        call correctbid(iim,jjp1*nlevnc,tnc)        call correctbid(iim,jjp1*nlevnc,tnc)
179        call dump2d(iip1,jjp1,tnc,'TNC COUCHE 1 BBB ')        call dump2d(iip1,jjp1,tnc,'TNC COUCHE 1 BBB ')
# Line 186  Line 183 
183  !  --------  !  --------
184    
185        if (guide_Q) then        if (guide_Q) then
186        status=NF_GET_VARA_REAL(ncidQ,varidQ,start,count,Qnc)        status=NF90_GET_VAR(ncidQ,varidQ,Qnc,start,count)
187        call correctbid(iim,jjp1*nlevnc,Qnc)        call correctbid(iim,jjp1*nlevnc,Qnc)
188        call dump2d(iip1,jjp1,Qnc,'QNC COUCHE 1 ')        call dump2d(iip1,jjp1,Qnc,'QNC COUCHE 1 ')
189        endif        endif
# Line 196  Line 193 
193  !  -------------  !  -------------
194    
195        if (guide_v) then        if (guide_v) then
196        status=NF_GET_VARA_REAL(ncidv,varidv,start,count,vnc)        status=NF90_GET_VAR(ncidv,varidv,vnc,start,count)
197        call correctbid(iim,jjm*nlevnc,vnc)        call correctbid(iim,jjm*nlevnc,vnc)
198        call dump2d(iip1,jjm,vnc,'VNC COUCHE 1 ')        call dump2d(iip1,jjm,vnc,'VNC COUCHE 1 ')
199        endif        endif
# Line 211  Line 208 
208  !  -------------------  !  -------------------
209    
210        if (guide_P) then        if (guide_P) then
211        status=NF_GET_VARA_REAL(ncidps,varidps,start,count,psnc)        status=NF90_GET_VAR(ncidps,varidps,psnc,start,count)
212        call dump2d(iip1,jjp1,psnc,'PSNC COUCHE 1 ')        call dump2d(iip1,jjp1,psnc,'PSNC COUCHE 1 ')
213        call correctbid(iim,jjp1,psnc)        call correctbid(iim,jjp1,psnc)
214        endif        endif
Legend:
Removed from v.60  
changed lines
  Added in v.61
  ViewVC Help
Powered by ViewVC 1.1.21