/[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 3 by guez, Wed Feb 27 13:16:39 2008 UTC revision 20 by guez, Wed Oct 15 16:19:57 2008 UTC
# Line 17  c -------------------------------------- Line 17  c --------------------------------------
17        use comvert        use comvert
18        use comgeom        use comgeom
19        use guide_m        use guide_m
20          use netcdf
21    
22        IMPLICIT NONE        IMPLICIT NONE
23    
24  c common  c common
# Line 69  c -------------------------------------- Line 71  c --------------------------------------
71    
72  c Vent zonal  c Vent zonal
73              if (guide_u) then              if (guide_u) then
74              ncidu=NCOPN('u.nc',NCNOWRIT,rcode)              rcode=nf90_open('u.nc',nf90_nowrite,ncidu)
75              varidu=NCVID(ncidu,'UWND',rcode)              rcode = nf90_inq_varid(ncidu, 'UWND', varidu)
76              print*,'ncidu,varidu',ncidu,varidu              print*,'ncidu,varidu',ncidu,varidu
77              if (ncidpl.eq.-99) ncidpl=ncidu              if (ncidpl.eq.-99) ncidpl=ncidu
78              endif              endif
79    
80  c Vent meridien  c Vent meridien
81              if (guide_v) then              if (guide_v) then
82              ncidv=NCOPN('v.nc',NCNOWRIT,rcode)              rcode=nf90_open('v.nc',nf90_nowrite,ncidv)
83              varidv=NCVID(ncidv,'VWND',rcode)              rcode = nf90_inq_varid(ncidv, 'VWND', varidv)
84              print*,'ncidv,varidv',ncidv,varidv              print*,'ncidv,varidv',ncidv,varidv
85              if (ncidpl.eq.-99) ncidpl=ncidv              if (ncidpl.eq.-99) ncidpl=ncidv
86              endif              endif
87    
88  c Temperature  c Temperature
89              if (guide_T) then              if (guide_T) then
90              ncidt=NCOPN('T.nc',NCNOWRIT,rcode)              rcode=nf90_open('T.nc',nf90_nowrite,ncidt)
91              varidt=NCVID(ncidt,'AIR',rcode)              rcode = nf90_inq_varid(ncidt, 'AIR', varidt)
92              print*,'ncidt,varidt',ncidt,varidt              print*,'ncidt,varidt',ncidt,varidt
93              if (ncidpl.eq.-99) ncidpl=ncidt              if (ncidpl.eq.-99) ncidpl=ncidt
94              endif              endif
95    
96  c Humidite  c Humidite
97              if (guide_Q) then              if (guide_Q) then
98              ncidQ=NCOPN('hur.nc',NCNOWRIT,rcode)              rcode=nf90_open('hur.nc',nf90_nowrite,ncidQ)
99              varidQ=NCVID(ncidQ,'RH',rcode)              rcode = nf90_inq_varid(ncidQ, 'RH', varidQ)
100              print*,'ncidQ,varidQ',ncidQ,varidQ              print*,'ncidQ,varidQ',ncidQ,varidQ
101              if (ncidpl.eq.-99) ncidpl=ncidQ              if (ncidpl.eq.-99) ncidpl=ncidQ
102              endif              endif
103    
104  c Pression de surface  c Pression de surface
105              if (guide_P) then              if (guide_P) then
106              ncidps=NCOPN('ps.nc',NCNOWRIT,rcode)              rcode=nf90_open('ps.nc',nf90_nowrite,ncidps)
107              varidps=NCVID(ncidps,'SP',rcode)              rcode = nf90_inq_varid(ncidps, 'SP', varidps)
108              print*,'ncidps,varidps',ncidps,varidps              print*,'ncidps,varidps',ncidps,varidps
109              endif              endif
110    
111  c Coordonnee verticale  c Coordonnee verticale
112              if (ncep) then              if (ncep) then
113                 print*,'Vous etes entrain de lire des donnees NCEP'                 print*,'Vous etes entrain de lire des donnees NCEP'
114                 varidpl=NCVID(ncidpl,'LEVEL',rcode)                 rcode = nf90_inq_varid(ncidpl, 'LEVEL', varidpl)
115              else              else
116                 print*,'Vous etes entrain de lire des donnees ECMWF'                 print*,'Vous etes entrain de lire des donnees ECMWF'
117                 varidpl=NCVID(ncidpl,'PRESSURE',rcode)                 rcode = nf90_inq_varid(ncidpl, 'PRESSURE', varidpl)
118              endif              endif
119              print*,'ncidu,varidpl',ncidu,varidpl              print*,'ncidu,varidpl',ncidu,varidpl
120        endif        endif
Legend:
Removed from v.3  
changed lines
  Added in v.20
  ViewVC Help
Powered by ViewVC 1.1.21