- Timestamp:
- 2016-05-26T11:08:07+02:00 (8 years ago)
- Location:
- branches/UKMO/dev_r5518_v3.4_asm_nemovar_community/NEMOGCM/TOOLS/REBUILD_NEMO
- Files:
-
- 1 deleted
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/dev_r5518_v3.4_asm_nemovar_community/NEMOGCM/TOOLS/REBUILD_NEMO/src/rebuild_nemo.f90
r6617 r6625 200 200 WRITE(numerr,*) 'Attribute DOMAIN_number_total is : ', ndomain_file 201 201 WRITE(numerr,*) 'Number of files specified in namelist is: ', ndomain 202 STOP 9202 STOP 203 203 ENDIF 204 204 … … 268 268 WRITE(numerr,*) 'Attribute DOMAIN_local_sizes is : ', local_sizes 269 269 WRITE(numerr,*) 'Dimensions to be rebuilt are of size : ', outdimlens(rebuild_dims(1)), outdimlens(rebuild_dims(2)) 270 STOP 9270 STOP 271 271 ENDIF 272 272 … … 384 384 SELECT CASE( xtype ) 385 385 CASE( NF90_BYTE ) 386 globaldata_0d_i1 = 0387 386 CALL check_nf90( nf90_get_var( ncid, jv, globaldata_0d_i1 ) ) 388 387 CASE( NF90_SHORT ) 389 globaldata_0d_i2 = 0390 388 CALL check_nf90( nf90_get_var( ncid, jv, globaldata_0d_i2 ) ) 391 389 CASE( NF90_INT ) 392 globaldata_0d_i4 = 0393 390 CALL check_nf90( nf90_get_var( ncid, jv, globaldata_0d_i4 ) ) 394 391 CASE( NF90_FLOAT ) 395 globaldata_0d_sp = 0.396 392 CALL check_nf90( nf90_get_var( ncid, jv, globaldata_0d_sp ) ) 397 393 CASE( NF90_DOUBLE ) 398 globaldata_0d_dp = 0.399 394 CALL check_nf90( nf90_get_var( ncid, jv, globaldata_0d_dp ) ) 400 395 CASE DEFAULT 401 396 WRITE(numerr,*) 'Unknown nf90 type: ', xtype 402 STOP 9397 STOP 403 398 END SELECT 404 399 … … 408 403 CASE( NF90_BYTE ) 409 404 ALLOCATE(globaldata_1d_i1(indimlens(dimids(1)))) 410 globaldata_1d_i1(:) = 0411 405 CALL check_nf90( nf90_get_var( ncid, jv, globaldata_1d_i1 ) ) 412 406 CASE( NF90_SHORT ) 413 407 ALLOCATE(globaldata_1d_i2(indimlens(dimids(1)))) 414 globaldata_1d_i2(:) = 0415 408 CALL check_nf90( nf90_get_var( ncid, jv, globaldata_1d_i2 ) ) 416 409 CASE( NF90_INT ) 417 410 ALLOCATE(globaldata_1d_i4(indimlens(dimids(1)))) 418 globaldata_1d_i4(:) = 0419 411 CALL check_nf90( nf90_get_var( ncid, jv, globaldata_1d_i4 ) ) 420 412 CASE( NF90_FLOAT ) 421 413 ALLOCATE(globaldata_1d_sp(indimlens(dimids(1)))) 422 globaldata_1d_sp(:) = 0.423 414 CALL check_nf90( nf90_get_var( ncid, jv, globaldata_1d_sp ) ) 424 415 CASE( NF90_DOUBLE ) 425 416 ALLOCATE(globaldata_1d_dp(indimlens(dimids(1)))) 426 globaldata_1d_dp(:) = 0.427 417 CALL check_nf90( nf90_get_var( ncid, jv, globaldata_1d_dp ) ) 428 418 CASE DEFAULT 429 419 WRITE(numerr,*) 'Unknown nf90 type: ', xtype 430 STOP 9420 STOP 431 421 END SELECT 432 422 … … 436 426 CASE( NF90_BYTE ) 437 427 ALLOCATE(globaldata_2d_i1(indimlens(dimids(1)),indimlens(dimids(2)))) 438 globaldata_2d_i1(:,:) = 0439 428 CALL check_nf90( nf90_get_var( ncid, jv, globaldata_2d_i1 ) ) 440 429 CASE( NF90_SHORT ) 441 430 ALLOCATE(globaldata_2d_i2(indimlens(dimids(1)),indimlens(dimids(2)))) 442 globaldata_2d_i2(:,:) = 0443 431 CALL check_nf90( nf90_get_var( ncid, jv, globaldata_2d_i2 ) ) 444 432 CASE( NF90_INT ) 445 433 ALLOCATE(globaldata_2d_i4(indimlens(dimids(1)),indimlens(dimids(2)))) 446 globaldata_2d_i4(:,:) = 0447 434 CALL check_nf90( nf90_get_var( ncid, jv, globaldata_2d_i4 ) ) 448 435 CASE( NF90_FLOAT ) 449 436 ALLOCATE(globaldata_2d_sp(indimlens(dimids(1)),indimlens(dimids(2)))) 450 globaldata_2d_sp(:,:) = 0.451 437 CALL check_nf90( nf90_get_var( ncid, jv, globaldata_2d_sp ) ) 452 438 CASE( NF90_DOUBLE ) 453 439 ALLOCATE(globaldata_2d_dp(indimlens(dimids(1)),indimlens(dimids(2)))) 454 globaldata_2d_dp(:,:) = 0.455 440 CALL check_nf90( nf90_get_var( ncid, jv, globaldata_2d_dp ) ) 456 441 CASE DEFAULT 457 442 WRITE(numerr,*) 'Unknown nf90 type: ', xtype 458 STOP 9443 STOP 459 444 END SELECT 460 445 … … 465 450 ALLOCATE(globaldata_3d_i1(indimlens(dimids(1)),indimlens(dimids(2)), & 466 451 & indimlens(dimids(3)))) 467 globaldata_3d_i1(:,:,:) = 0468 452 CALL check_nf90( nf90_get_var( ncid, jv, globaldata_3d_i1 ) ) 469 453 CASE( NF90_SHORT ) 470 454 ALLOCATE(globaldata_3d_i2(indimlens(dimids(1)),indimlens(dimids(2)), & 471 455 & indimlens(dimids(3)))) 472 globaldata_3d_i2(:,:,:) = 0473 456 CALL check_nf90( nf90_get_var( ncid, jv, globaldata_3d_i2 ) ) 474 457 CASE( NF90_INT ) 475 458 ALLOCATE(globaldata_3d_i4(indimlens(dimids(1)),indimlens(dimids(2)), & 476 459 & indimlens(dimids(3)))) 477 globaldata_3d_i4(:,:,:) = 0478 460 CALL check_nf90( nf90_get_var( ncid, jv, globaldata_3d_i4 ) ) 479 461 CASE( NF90_FLOAT ) 480 462 ALLOCATE(globaldata_3d_sp(indimlens(dimids(1)),indimlens(dimids(2)), & 481 463 & indimlens(dimids(3)))) 482 globaldata_3d_sp(:,:,:) = 0.483 464 CALL check_nf90( nf90_get_var( ncid, jv, globaldata_3d_sp ) ) 484 465 CASE( NF90_DOUBLE ) 485 466 ALLOCATE(globaldata_3d_dp(indimlens(dimids(1)),indimlens(dimids(2)), & 486 467 & indimlens(dimids(3)))) 487 globaldata_3d_dp(:,:,:) = 0.488 468 CALL check_nf90( nf90_get_var( ncid, jv, globaldata_3d_dp ) ) 489 469 CASE DEFAULT 490 470 WRITE(numerr,*) 'Unknown nf90 type: ', xtype 491 STOP 9471 STOP 492 472 END SELECT 493 473 … … 498 478 ALLOCATE(globaldata_4d_i1(indimlens(dimids(1)),indimlens(dimids(2)), & 499 479 & indimlens(dimids(3)),ntchunk)) 500 globaldata_4d_i1(:,:,:,:) = 0501 480 CALL check_nf90( nf90_get_var( ncid, jv, globaldata_4d_i1, start=(/1,1,1,nt/) ) ) 502 481 CASE( NF90_SHORT ) 503 482 ALLOCATE(globaldata_4d_i2(indimlens(dimids(1)),indimlens(dimids(2)), & 504 483 & indimlens(dimids(3)),ntchunk)) 505 globaldata_4d_i2(:,:,:,:) = 0506 484 CALL check_nf90( nf90_get_var( ncid, jv, globaldata_4d_i2, start=(/1,1,1,nt/) ) ) 507 485 CASE( NF90_INT ) 508 486 ALLOCATE(globaldata_4d_i4(indimlens(dimids(1)),indimlens(dimids(2)), & 509 487 & indimlens(dimids(3)),ntchunk)) 510 globaldata_4d_i4(:,:,:,:) = 0511 488 CALL check_nf90( nf90_get_var( ncid, jv, globaldata_4d_i4, start=(/1,1,1,nt/) ) ) 512 489 CASE( NF90_FLOAT ) 513 490 ALLOCATE(globaldata_4d_sp(indimlens(dimids(1)),indimlens(dimids(2)), & 514 491 & indimlens(dimids(3)),ntchunk)) 515 globaldata_4d_sp(:,:,:,:) = 0.516 492 CALL check_nf90( nf90_get_var( ncid, jv, globaldata_4d_sp, start=(/1,1,1,nt/) ) ) 517 493 CASE( NF90_DOUBLE ) 518 494 ALLOCATE(globaldata_4d_dp(indimlens(dimids(1)),indimlens(dimids(2)), & 519 495 & indimlens(dimids(3)),ntchunk)) 520 globaldata_4d_dp(:,:,:,:) = 0.521 496 CALL check_nf90( nf90_get_var( ncid, jv, globaldata_4d_dp, start=(/1,1,1,nt/) ) ) 522 497 CASE DEFAULT 523 498 WRITE(numerr,*) 'Unknown nf90 type: ', xtype 524 STOP 9499 STOP 525 500 END SELECT 526 501 … … 542 517 CASE( NF90_BYTE ) 543 518 ALLOCATE(globaldata_1d_i1(outdimlens(dimids(1)))) 544 globaldata_1d_i1(:) = 0545 519 CASE( NF90_SHORT ) 546 520 ALLOCATE(globaldata_1d_i2(outdimlens(dimids(1)))) 547 globaldata_1d_i2(:) = 0548 521 CASE( NF90_INT ) 549 522 ALLOCATE(globaldata_1d_i4(outdimlens(dimids(1)))) 550 globaldata_1d_i4(:) = 0551 523 CASE( NF90_FLOAT ) 552 524 ALLOCATE(globaldata_1d_sp(outdimlens(dimids(1)))) 553 globaldata_1d_sp(:) = 0.554 525 CASE( NF90_DOUBLE ) 555 526 ALLOCATE(globaldata_1d_dp(outdimlens(dimids(1)))) 556 globaldata_1d_dp(:) = 0.557 527 CASE DEFAULT 558 528 WRITE(numerr,*) 'Unknown nf90 type: ', xtype 559 STOP 9529 STOP 560 530 END SELECT 561 531 … … 565 535 CASE( NF90_BYTE ) 566 536 ALLOCATE(globaldata_2d_i1(outdimlens(dimids(1)),outdimlens(dimids(2)))) 567 globaldata_2d_i1(:,:) = 0568 537 CASE( NF90_SHORT ) 569 538 ALLOCATE(globaldata_2d_i2(outdimlens(dimids(1)),outdimlens(dimids(2)))) 570 globaldata_2d_i2(:,:) = 0571 539 CASE( NF90_INT ) 572 540 ALLOCATE(globaldata_2d_i4(outdimlens(dimids(1)),outdimlens(dimids(2)))) 573 globaldata_2d_i4(:,:) = 0574 541 CASE( NF90_FLOAT ) 575 542 ALLOCATE(globaldata_2d_sp(outdimlens(dimids(1)),outdimlens(dimids(2)))) 576 globaldata_2d_sp(:,:) = 0.577 543 CASE( NF90_DOUBLE ) 578 544 ALLOCATE(globaldata_2d_dp(outdimlens(dimids(1)),outdimlens(dimids(2)))) 579 globaldata_2d_dp(:,:) = 0.580 545 CASE DEFAULT 581 546 WRITE(numerr,*) 'Unknown nf90 type: ', xtype 582 STOP 9547 STOP 583 548 END SELECT 584 549 … … 589 554 ALLOCATE(globaldata_3d_i1(outdimlens(dimids(1)),outdimlens(dimids(2)), & 590 555 & outdimlens(dimids(3)))) 591 globaldata_3d_i1(:,:,:) = 0592 556 CASE( NF90_SHORT ) 593 557 ALLOCATE(globaldata_3d_i2(outdimlens(dimids(1)),outdimlens(dimids(2)), & 594 558 & outdimlens(dimids(3)))) 595 globaldata_3d_i2(:,:,:) = 0596 559 CASE( NF90_INT ) 597 560 ALLOCATE(globaldata_3d_i4(outdimlens(dimids(1)),outdimlens(dimids(2)), & 598 561 & outdimlens(dimids(3)))) 599 globaldata_3d_i4(:,:,:) = 0600 562 CASE( NF90_FLOAT ) 601 563 ALLOCATE(globaldata_3d_sp(outdimlens(dimids(1)),outdimlens(dimids(2)), & 602 564 & outdimlens(dimids(3)))) 603 globaldata_3d_sp(:,:,:) = 0.604 565 CASE( NF90_DOUBLE ) 605 566 ALLOCATE(globaldata_3d_dp(outdimlens(dimids(1)),outdimlens(dimids(2)), & 606 567 & outdimlens(dimids(3)))) 607 globaldata_3d_dp(:,:,:) = 0.608 568 CASE DEFAULT 609 569 WRITE(numerr,*) 'Unknown nf90 type: ', xtype 610 STOP 9570 STOP 611 571 END SELECT 612 572 … … 617 577 ALLOCATE(globaldata_4d_i1(outdimlens(dimids(1)),outdimlens(dimids(2)), & 618 578 & outdimlens(dimids(3)),ntchunk)) 619 globaldata_4d_i1(:,:,:,:) = 0620 579 CASE( NF90_SHORT ) 621 580 ALLOCATE(globaldata_4d_i2(outdimlens(dimids(1)),outdimlens(dimids(2)), & 622 581 & outdimlens(dimids(3)),ntchunk)) 623 globaldata_4d_i2(:,:,:,:) = 0624 582 CASE( NF90_INT ) 625 583 ALLOCATE(globaldata_4d_i4(outdimlens(dimids(1)),outdimlens(dimids(2)), & 626 584 & outdimlens(dimids(3)),ntchunk)) 627 globaldata_4d_i4(:,:,:,:) = 0628 585 CASE( NF90_FLOAT ) 629 586 ALLOCATE(globaldata_4d_sp(outdimlens(dimids(1)),outdimlens(dimids(2)), & 630 587 & outdimlens(dimids(3)),ntchunk)) 631 globaldata_4d_sp(:,:,:,:) = 0.632 588 CASE( NF90_DOUBLE ) 633 589 ALLOCATE(globaldata_4d_dp(outdimlens(dimids(1)),outdimlens(dimids(2)), & 634 590 & outdimlens(dimids(3)),ntchunk)) 635 globaldata_4d_dp(:,:,:,:) = 0.636 591 CASE DEFAULT 637 592 WRITE(numerr,*) 'Unknown nf90 type: ', xtype 638 STOP 9593 STOP 639 594 END SELECT 640 595 ELSE 641 596 WRITE(numerr,*) 'ERROR! : A netcdf variable has more than 4 dimensions which is not taken into account' 642 STOP 9597 STOP 643 598 ENDIF 644 599 … … 1012 967 IF( nthreads == 1 .AND. istop /= nf90_noerr ) THEN 1013 968 WRITE(numerr,*) '*** NEMO rebuild failed! ***' 1014 STOP 9969 STOP 1015 970 ENDIF 1016 971 … … 1021 976 IF( istop /= nf90_noerr ) THEN 1022 977 WRITE(numerr,*) '*** NEMO rebuild failed! ***' 1023 STOP 9978 STOP 1024 979 ENDIF 1025 980 … … 1095 1050 CASE DEFAULT 1096 1051 WRITE(numerr,*) 'Unknown nf90 type: ', xtype 1097 STOP 91052 STOP 1098 1053 END SELECT 1099 1054 … … 1118 1073 CASE DEFAULT 1119 1074 WRITE(numerr,*) 'Unknown nf90 type: ', xtype 1120 STOP 91075 STOP 1121 1076 END SELECT 1122 1077 … … 1141 1096 CASE DEFAULT 1142 1097 WRITE(numerr,*) 'Unknown nf90 type: ', xtype 1143 STOP 91098 STOP 1144 1099 END SELECT 1145 1100 … … 1191 1146 WRITE(numerr,*) "*** NEMO rebuild failed ***" 1192 1147 WRITE(numerr,*) 1193 STOP 91148 STOP 1194 1149 ENDIF 1195 1150 ENDIF
Note: See TracChangeset
for help on using the changeset viewer.