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