Changeset 753 for LMDZ4/branches/V3_test/libf/dyn3dpar/read_reanalyse.F
- Timestamp:
- Mar 1, 2007, 11:28:53 AM (18 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ4/branches/V3_test/libf/dyn3dpar/read_reanalyse.F
r709 r753 10 10 c mode=1 variabels GCM 11 11 12 USE parallel 12 13 c ----------------------------------------------------------------- 13 14 c Declarations … … 42 43 integer ncidpl 43 44 integer varidpl,ncidQ,varidQ 44 save ncid u,varidu,ncidv,varidv,ncidt,varidt,ncidps,varidps45 save ncidpl,ncidu,varidu,ncidv,varidv,ncidt,varidt,ncidps,varidps 45 46 save varidpl,ncidQ,varidQ 46 47 … … 55 56 logical first 56 57 save first 58 INTEGER ierr 57 59 58 60 data first/.true./ … … 69 71 c Vent zonal 70 72 if (guide_u) then 71 ncidu=NCOPN('u.nc',NCNOWRIT,rcode)72 varidu=NCVID(ncidu,'UWND',rcode)73 print*,'ncidu,varidu',ncidu,varidu74 if (ncidpl.eq.-99) ncidpl=ncidu73 ncidu=NCOPN('u.nc',NCNOWRIT,rcode) 74 varidu=NCVID(ncidu,'UWND',rcode) 75 print*,'ncidu,varidu',ncidu,varidu 76 if (ncidpl.eq.-99) ncidpl=ncidu 75 77 endif 76 78 77 79 c Vent meridien 78 80 if (guide_v) then 79 ncidv=NCOPN('v.nc',NCNOWRIT,rcode) 80 varidv=NCVID(ncidv,'VWND',rcode) 81 print*,'ncidv,varidv',ncidv,varidv 82 if (ncidpl.eq.-99) ncidpl=ncidu 83 endif 81 ncidv=NCOPN('v.nc',NCNOWRIT,rcode) 82 varidv=NCVID(ncidv,'VWND',rcode) 83 print*,'ncidv,varidv',ncidv,varidv 84 if (ncidpl.eq.-99) ncidpl=ncidu 85 endif 86 84 87 85 88 c Temperature 86 89 if (guide_T) then 87 ncidt=NCOPN('T.nc',NCNOWRIT,rcode)88 varidt=NCVID(ncidt,'AIR',rcode)89 print*,'ncidt,varidt',ncidt,varidt90 if (ncidpl.eq.-99) ncidpl=ncidu91 endif 92 90 ncidt=NCOPN('T.nc',NCNOWRIT,rcode) 91 varidt=NCVID(ncidt,'AIR',rcode) 92 print*,'ncidt,varidt',ncidt,varidt 93 if (ncidpl.eq.-99) ncidpl=ncidu 94 endif 95 93 96 c Humidite 94 97 if (guide_Q) then 95 ncidQ=NCOPN('hur.nc',NCNOWRIT,rcode)96 varidQ=NCVID(ncidQ,'RH',rcode)97 print*,'ncidQ,varidQ',ncidQ,varidQ98 if (ncidpl.eq.-99) ncidpl=ncidu99 endif 100 98 ncidQ=NCOPN('hur.nc',NCNOWRIT,rcode) 99 varidQ=NCVID(ncidQ,'RH',rcode) 100 print*,'ncidQ,varidQ',ncidQ,varidQ 101 if (ncidpl.eq.-99) ncidpl=ncidu 102 endif 103 101 104 c Pression de surface 102 105 if (guide_P) then 103 ncidps=NCOPN('ps.nc',NCNOWRIT,rcode)104 varidps=NCVID(ncidps,'SP',rcode)105 print*,'ncidps,varidps',ncidps,varidps106 endif 107 106 ncidps=NCOPN('ps.nc',NCNOWRIT,rcode) 107 varidps=NCVID(ncidps,'SP',rcode) 108 print*,'ncidps,varidps',ncidps,varidps 109 endif 110 108 111 c Coordonnee vertcale 109 112 if (ncep) then … … 134 137 status=NF_GET_VARA_REAL(ncidpl,varidpl,1,nlevnc,pl) 135 138 #endif 139 136 140 c passage en pascal 137 141 pl(:)=100.*pl(:) … … 160 164 count(4)=1 161 165 166 167 c mise a zero des tableaux 168 c ------------------------ 169 unc(:,:,:)=0. 170 vnc(:,:,:)=0. 171 tnc(:,:,:)=0. 172 Qnc(:,:,:)=0. 173 162 174 c Vent zonal 163 175 c ---------- 164 176 165 177 if (guide_u) then 166 print*,'avant la lecture de UNCEP nd de niv:',nlevnc 178 print*,'avant la lecture de UNCEP nd de niv:',nlevnc 179 167 180 #ifdef NC_DOUBLE 168 status=NF_GET_VARA_DOUBLE(ncidu,varidu,start,count,unc)181 status=NF_GET_VARA_DOUBLE(ncidu,varidu,start,count,unc) 169 182 #else 170 status=NF_GET_VARA_REAL(ncidu,varidu,start,count,unc)183 status=NF_GET_VARA_REAL(ncidu,varidu,start,count,unc) 171 184 #endif 172 185 c call dump2d(iip1,jjp1,unc,'VENT NCEP ') 173 186 c call dump2d(iip1,40,unc(1,1,nlevnc),'VENT NCEP ') 174 print*,'WARNING!!! Correction bidon pour palier a un '175 print*,'probleme dans la creation des fichiers nc'176 call correctbid(iim,jjp1*nlevnc,unc)177 call dump2d(iip1,jjp1,unc,'UNC COUCHE 1 ')187 print*,'WARNING!!! Correction bidon pour palier a un ' 188 print*,'probleme dans la creation des fichiers nc' 189 call correctbid(iim,jjp1*nlevnc,unc) 190 call dump2d(iip1,jjp1,unc,'UNC COUCHE 1 ') 178 191 endif 179 192 … … 185 198 if (guide_T) then 186 199 #ifdef NC_DOUBLE 187 status=NF_GET_VARA_DOUBLE(ncidt,varidt,start,count,tnc)200 status=NF_GET_VARA_DOUBLE(ncidt,varidt,start,count,tnc) 188 201 #else 189 status=NF_GET_VARA_REAL(ncidt,varidt,start,count,tnc)202 status=NF_GET_VARA_REAL(ncidt,varidt,start,count,tnc) 190 203 #endif 191 call dump2d(iip1,jjp1,tnc,'TNC COUCHE 1 AAA ')192 call correctbid(iim,jjp1*nlevnc,tnc)193 call dump2d(iip1,jjp1,tnc,'TNC COUCHE 1 BBB ')204 call dump2d(iip1,jjp1,tnc,'TNC COUCHE 1 AAA ') 205 call correctbid(iim,jjp1*nlevnc,tnc) 206 call dump2d(iip1,jjp1,tnc,'TNC COUCHE 1 BBB ') 194 207 endif 195 208 … … 199 212 if (guide_Q) then 200 213 #ifdef NC_DOUBLE 201 status=NF_GET_VARA_DOUBLE(ncidQ,varidQ,start,count,Qnc)214 status=NF_GET_VARA_DOUBLE(ncidQ,varidQ,start,count,Qnc) 202 215 #else 203 status=NF_GET_VARA_REAL(ncidQ,varidQ,start,count,Qnc)216 status=NF_GET_VARA_REAL(ncidQ,varidQ,start,count,Qnc) 204 217 #endif 205 call correctbid(iim,jjp1*nlevnc,Qnc)206 call dump2d(iip1,jjp1,Qnc,'QNC COUCHE 1 ')218 call correctbid(iim,jjp1*nlevnc,Qnc) 219 call dump2d(iip1,jjp1,Qnc,'QNC COUCHE 1 ') 207 220 endif 208 221 … … 213 226 if (guide_v) then 214 227 #ifdef NC_DOUBLE 215 status=NF_GET_VARA_DOUBLE(ncidv,varidv,start,count,vnc)228 status=NF_GET_VARA_DOUBLE(ncidv,varidv,start,count,vnc) 216 229 #else 217 status=NF_GET_VARA_REAL(ncidv,varidv,start,count,vnc)230 status=NF_GET_VARA_REAL(ncidv,varidv,start,count,vnc) 218 231 #endif 219 call correctbid(iim,jjm*nlevnc,vnc)220 call dump2d(iip1,jjm,vnc,'VNC COUCHE 1 ')232 call correctbid(iim,jjm*nlevnc,vnc) 233 call dump2d(iip1,jjm,vnc,'VNC COUCHE 1 ') 221 234 endif 222 235 … … 232 245 if (guide_P) then 233 246 #ifdef NC_DOUBLE 234 status=NF_GET_VARA_DOUBLE(ncidps,varidps,start,count,psnc)247 status=NF_GET_VARA_DOUBLE(ncidps,varidps,start,count,psnc) 235 248 #else 236 status=NF_GET_VARA_REAL(ncidps,varidps,start,count,psnc)249 status=NF_GET_VARA_REAL(ncidps,varidps,start,count,psnc) 237 250 #endif 238 call dump2d(iip1,jjp1,psnc,'PSNC COUCHE 1 ')239 call correctbid(iim,jjp1,psnc)251 call dump2d(iip1,jjp1,psnc,'PSNC COUCHE 1 ') 252 call correctbid(iim,jjp1,psnc) 240 253 endif 241 254
Note: See TracChangeset
for help on using the changeset viewer.
