Changeset 1146 for LMDZ4/trunk/libf/dyn3dpar/read_reanalyse.F
- Timestamp:
- Apr 9, 2009, 12:11:35 PM (15 years ago)
- Location:
- LMDZ4/trunk
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ4/trunk
-
Property
svn:mergeinfo
set to
/LMDZ4/branches/LMDZ4-dev merged eligible
-
Property
svn:mergeinfo
set to
-
LMDZ4/trunk/libf/dyn3dpar/read_reanalyse.F
r1122 r1146 11 11 12 12 USE parallel 13 use netcdf 13 14 c ----------------------------------------------------------------- 14 15 c Declarations … … 72 73 print *,'Vous êtes entrain de lire des données sur 73 74 . niveaux modèle' 74 ncidpl=NCOPN('apbp.nc',NCNOWRIT,rcode)75 varidap=NCVID(ncidpl,'AP',rcode)76 varidbp=NCVID(ncidpl,'BP',rcode)75 rcode=nf90_open('apbp.nc',nf90_nowrite,ncidpl) 76 rcode = nf90_inq_varid(ncidpl, 'AP', varidap) 77 rcode = nf90_inq_varid(ncidpl, 'BP', varidbp) 77 78 print*,'ncidpl,varidap',ncidpl,varidap 78 79 endif … … 80 81 c Vent zonal 81 82 if (guide_u) then 82 ncidu=NCOPN('u.nc',NCNOWRIT,rcode)83 varidu=NCVID(ncidu,'UWND',rcode)83 rcode=nf90_open('u.nc',nf90_nowrite,ncidu) 84 rcode = nf90_inq_varid(ncidu, 'UWND', varidu) 84 85 print*,'ncidu,varidu',ncidu,varidu 85 86 if (ncidpl.eq.-99) ncidpl=ncidu … … 88 89 c Vent meridien 89 90 if (guide_v) then 90 ncidv=NCOPN('v.nc',NCNOWRIT,rcode)91 varidv=NCVID(ncidv,'VWND',rcode)91 rcode=nf90_open('v.nc',nf90_nowrite,ncidv) 92 rcode = nf90_inq_varid(ncidv, 'VWND', varidv) 92 93 print*,'ncidv,varidv',ncidv,varidv 93 94 if (ncidpl.eq.-99) ncidpl=ncidv … … 96 97 c Temperature 97 98 if (guide_T) then 98 ncidt=NCOPN('T.nc',NCNOWRIT,rcode)99 varidt=NCVID(ncidt,'AIR',rcode)99 rcode=nf90_open('T.nc',nf90_nowrite,ncidt) 100 rcode = nf90_inq_varid(ncidt, 'AIR', varidt) 100 101 print*,'ncidt,varidt',ncidt,varidt 101 102 if (ncidpl.eq.-99) ncidpl=ncidt … … 104 105 c Humidite 105 106 if (guide_Q) then 106 ncidQ=NCOPN('hur.nc',NCNOWRIT,rcode)107 varidQ=NCVID(ncidQ,'RH',rcode)107 rcode=nf90_open('hur.nc',nf90_nowrite,ncidQ) 108 rcode = nf90_inq_varid(ncidQ, 'RH', varidQ) 108 109 print*,'ncidQ,varidQ',ncidQ,varidQ 109 110 if (ncidpl.eq.-99) ncidpl=ncidQ … … 112 113 c Pression de surface 113 114 if ((guide_P).OR.(guide_modele)) then 114 ncidps=NCOPN('ps.nc',NCNOWRIT,rcode)115 varidps=NCVID(ncidps,'SP',rcode)115 rcode=nf90_open('ps.nc',nf90_nowrite,ncidps) 116 rcode = nf90_inq_varid(ncidps, 'SP', varidps) 116 117 print*,'ncidps,varidps',ncidps,varidps 117 118 endif … … 119 120 c Coordonnee verticale 120 121 if (.not.guide_modele) then 121 if (ncep) then122 print*,'Vous etes entrain de lire des donnees NCEP'123 varidpl=NCVID(ncidpl,'LEVEL',rcode)124 else125 print*,'Vous etes entrain de lire des donnees ECMWF'126 varidpl=NCVID(ncidpl,'PRESSURE',rcode)127 endif128 print*,'ncidpl,varidpl',ncidpl,varidpl122 if (ncep) then 123 print*,'Vous etes entrain de lire des donnees NCEP' 124 rcode = nf90_inq_varid(ncidpl, 'LEVEL', varidpl) 125 else 126 print*,'Vous etes entrain de lire des donnees ECMWF' 127 rcode = nf90_inq_varid(ncidpl, 'PRESSURE', varidpl) 128 endif 129 print*,'ncidpl,varidpl',ncidpl,varidpl 129 130 endif 130 131 print*,'ncidu,varidpl',ncidu,varidpl
Note: See TracChangeset
for help on using the changeset viewer.