Changeset 218 for LMDZ.3.3/trunk/libf/dyn3d/read_pstoke.F
- Timestamp:
- May 14, 2001, 5:01:36 PM (23 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ.3.3/trunk/libf/dyn3d/read_pstoke.F
r185 r218 196 196 c niveaux de pression 197 197 198 #ifdef NC_DOUBLE 199 status=NF_GET_VARA_DOUBLE(ncidp,varidpl,1,kev,pl) 200 #else 198 201 status=NF_GET_VARA_REAL(ncidp,varidpl,1,kev,pl) 202 #endif 199 203 200 204 c lecture de aire et phis … … 212 216 c 213 217 c phis 218 #ifdef NC_DOUBLE 219 status=NF_GET_VARA_DOUBLE(ncidp,varidps,start,count,phisfi2) 220 #else 214 221 status=NF_GET_VARA_REAL(ncidp,varidps,start,count,phisfi2) 222 #endif 223 224 215 225 call gr_ecrit_fi(1,kon,iim,jjm+1,phisfi2,phisfi) 216 226 217 227 c aire 228 #ifdef NC_DOUBLE 229 status=NF_GET_VARA_DOUBLE(ncidp,varidai,start,count,airefi2) 230 #else 218 231 status=NF_GET_VARA_REAL(ncidp,varidai,start,count,airefi2) 232 #endif 233 219 234 call gr_ecrit_fi(1,kon,iim,jjm+1,airefi2,airefi) 220 235 else … … 241 256 c frac_impa 242 257 258 #ifdef NC_DOUBLE 259 status=NF_GET_VARA_DOUBLE(ncidp,varidfi,start,count,frac_impa2) 260 #else 243 261 status=NF_GET_VARA_REAL(ncidp,varidfi,start,count,frac_impa2) 262 #endif 244 263 call gr_ecrit_fi(kev,kon,iim,jjm+1,frac_impa2,frac_impa) 245 264 246 265 c frac_nucl 247 266 267 #ifdef NC_DOUBLE 268 status=NF_GET_VARA_DOUBLE(ncidp,varidfn,start,count,frac_nucl2) 269 #else 248 270 status=NF_GET_VARA_REAL(ncidp,varidfn,start,count,frac_nucl2) 271 #endif 249 272 call gr_ecrit_fi(kev,kon,iim,jjm+1,frac_nucl2,frac_nucl) 250 273 251 274 c abder t 275 276 #ifdef NC_DOUBLE 277 status=NF_GET_VARA_DOUBLE(ncidp,varidt,start,count,t2) 278 #else 252 279 status=NF_GET_VARA_REAL(ncidp,varidt,start,count,t2) 280 #endif 281 253 282 call gr_ecrit_fi(kev,kon,iim,jjm+1,t2,t) 254 283 255 284 c mfu 285 286 #ifdef NC_DOUBLE 287 status=NF_GET_VARA_DOUBLE(ncidp,varidmfu,start,count,mfu2) 288 #else 256 289 status=NF_GET_VARA_REAL(ncidp,varidmfu,start,count,mfu2) 290 #endif 291 257 292 call gr_ecrit_fi(kev,kon,iim,jjm+1,mfu2,mfu) 258 293 259 294 c mfd 295 #ifdef NC_DOUBLE 296 status=NF_GET_VARA_DOUBLE(ncidp,varidmfd,start,count,mfd2) 297 #else 260 298 status=NF_GET_VARA_REAL(ncidp,varidmfd,start,count,mfd2) 299 #endif 261 300 call gr_ecrit_fi(kev,kon,iim,jjm+1,mfd2,mfd) 262 301 263 302 c en_u 303 #ifdef NC_DOUBLE 304 status=NF_GET_VARA_DOUBLE(ncidp,varidenu,start,count,en_u2) 305 #else 264 306 status=NF_GET_VARA_REAL(ncidp,varidenu,start,count,en_u2) 307 #endif 265 308 call gr_ecrit_fi(kev,kon,iim,jjm+1,en_u2,en_u) 266 309 267 310 c de_u 311 #ifdef NC_DOUBLE 312 status=NF_GET_VARA_DOUBLE(ncidp,variddeu,start,count,de_u2) 313 #else 268 314 status=NF_GET_VARA_REAL(ncidp,variddeu,start,count,de_u2) 315 #endif 269 316 call gr_ecrit_fi(kev,kon,iim,jjm+1,de_u2,de_u) 270 317 271 318 c en_d 319 #ifdef NC_DOUBLE 320 status=NF_GET_VARA_DOUBLE(ncidp,varidend,start,count,en_d2) 321 #else 272 322 status=NF_GET_VARA_REAL(ncidp,varidend,start,count,en_d2) 323 #endif 273 324 call gr_ecrit_fi(kev,kon,iim,jjm+1,en_d2,en_d) 274 325 275 326 c de_d 327 #ifdef NC_DOUBLE 328 status=NF_GET_VARA_DOUBLE(ncidp,varidded,start,count,de_d2) 329 #else 276 330 status=NF_GET_VARA_REAL(ncidp,varidded,start,count,de_d2) 331 #endif 277 332 call gr_ecrit_fi(kev,kon,iim,jjm+1,de_d2,de_d) 278 333 279 334 c coefh 280 335 print*,'LECTURE de coefh a irec =',irec 336 #ifdef NC_DOUBLE 337 status=NF_GET_VARA_DOUBLE(ncidp,varidch,start,count,coefh2) 338 #else 281 339 status=NF_GET_VARA_REAL(ncidp,varidch,start,count,coefh2) 340 #endif 341 282 342 call gr_ecrit_fi(kev,kon,iim,jjm+1,coefh2,coefh) 283 343 … … 289 349 c pyu1 290 350 print*,'LECTURE de yu1 a irec =',irec 351 #ifdef NC_DOUBLE 352 status=NF_GET_VARA_DOUBLE(ncidp,varidyu1,start,count,pyu12) 353 #else 291 354 status=NF_GET_VARA_REAL(ncidp,varidyu1,start,count,pyu12) 355 #endif 356 292 357 call gr_ecrit_fi(1,kon,iim,jjm+1,pyu12,pyu1) 293 358 294 359 c pyv1 295 360 print*,'LECTURE de yv1 a irec =',irec 361 #ifdef NC_DOUBLE 362 status=NF_GET_VARA_DOUBLE(ncidp,varidyv1,start,count,pyv12) 363 #else 296 364 status=NF_GET_VARA_REAL(ncidp,varidyv1,start,count,pyv12) 365 #endif 297 366 call gr_ecrit_fi(1,kon,iim,jjm+1,pyv12,pyv1) 298 367 299 368 c ftsol1 300 369 print*,'LECTURE de ftsol1 a irec =',irec 370 #ifdef NC_DOUBLE 371 status=NF_GET_VARA_DOUBLE(ncidp,varidfts1,start,count,ftsol12) 372 #else 301 373 status=NF_GET_VARA_REAL(ncidp,varidfts1,start,count,ftsol12) 374 #endif 302 375 call gr_ecrit_fi(1,kon,iim,jjm+1,ftsol12,ftsol1) 303 376 304 377 c ftsol2 305 378 print*,'LECTURE de ftsol2 a irec =',irec 379 #ifdef NC_DOUBLE 380 status=NF_GET_VARA_DOUBLE(ncidp,varidfts2,start,count,ftsol22) 381 #else 306 382 status=NF_GET_VARA_REAL(ncidp,varidfts2,start,count,ftsol22) 383 #endif 307 384 call gr_ecrit_fi(1,kon,iim,jjm+1,ftsol22,ftsol2) 308 385 309 386 c ftsol3 310 387 print*,'LECTURE de ftsol3 a irec =',irec 388 #ifdef NC_DOUBLE 389 status=NF_GET_VARA_DOUBLE(ncidp,varidfts3,start,count,ftsol32) 390 #else 311 391 status=NF_GET_VARA_REAL(ncidp,varidfts3,start,count,ftsol32) 392 #endif 312 393 call gr_ecrit_fi(1,kon,iim,jjm+1,ftsol32,ftsol3) 313 394 314 395 c ftsol4 396 #ifdef NC_DOUBLE 397 status=NF_GET_VARA_DOUBLE(ncidp,varidfts4,start,count,ftsol42) 398 #else 315 399 status=NF_GET_VARA_REAL(ncidp,varidfts4,start,count,ftsol42) 400 #endif 316 401 call gr_ecrit_fi(1,kon,iim,jjm+1,ftsol42,ftsol4) 317 402 318 403 c psrf1 404 #ifdef NC_DOUBLE 405 status=NF_GET_VARA_DOUBLE(ncidp,varidpsr1,start,count,psrf12) 406 #else 319 407 status=NF_GET_VARA_REAL(ncidp,varidpsr1,start,count,psrf12) 408 #endif 320 409 c call dump2d(iip1-1,jjm+1,psrf12,'PSRF1NC') 321 410 call gr_ecrit_fi(1,kon,iim,jjm+1,psrf12,psrf1) 322 411 323 412 c psrf2 413 #ifdef NC_DOUBLE 414 status=NF_GET_VARA_DOUBLE(ncidp,varidpsr2,start,count,psrf22) 415 #else 324 416 status=NF_GET_VARA_REAL(ncidp,varidpsr2,start,count,psrf22) 417 #endif 325 418 c call dump2d(iip1-1,jjm+1,psrf22,'PSRF2NC') 326 419 call gr_ecrit_fi(1,kon,iim,jjm+1,psrf22,psrf2) 327 420 328 421 c psrf3 422 #ifdef NC_DOUBLE 423 status=NF_GET_VARA_DOUBLE(ncidp,varidpsr3,start,count,psrf32) 424 #else 329 425 status=NF_GET_VARA_REAL(ncidp,varidpsr3,start,count,psrf32) 426 #endif 330 427 call gr_ecrit_fi(1,kon,iim,jjm+1,psrf32,psrf3) 331 428 332 429 c psrf4 430 #ifdef NC_DOUBLE 431 status=NF_GET_VARA_DOUBLE(ncidp,varidpsr4,start,count,psrf42) 432 #else 333 433 status=NF_GET_VARA_REAL(ncidp,varidpsr4,start,count,psrf42) 434 #endif 334 435 call gr_ecrit_fi(1,kon,iim,jjm+1,psrf42,psrf4) 335 436
Note: See TracChangeset
for help on using the changeset viewer.