Changeset 11624 for NEMO/branches/2019/dev_r11613_ENHANCE-04_namelists_as_internalfiles/src/TOP/PISCES/SED/sedini.F90
- Timestamp:
- 2019-10-01T15:07:45+02:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2019/dev_r11613_ENHANCE-04_namelists_as_internalfiles/src/TOP/PISCES/SED/sedini.F90
r11536 r11624 406 406 !!---------------------------------------------------------------------- 407 407 408 INTEGER :: numnamsed_ref = -1 !! Logical units fornamelist sediment409 INTEGER :: numnamsed_cfg = -1 !! Logical units fornamelist sediment408 CHARACTER(:), ALLOCATABLE :: numnamsed_ref !! Character buffer for reference namelist sediment 409 CHARACTER(:), ALLOCATABLE :: numnamsed_cfg !! Character buffer for configuration namelist sediment 410 410 INTEGER :: ios ! Local integer output status for namelist read 411 411 CHARACTER(LEN=20) :: clname … … 452 452 IF(lwp) WRITE(numsed,*) ' sed_init_nam : read SEDIMENT namelist' 453 453 IF(lwp) WRITE(numsed,*) ' ~~~~~~~~~~~~~~' 454 CALL ctl_opn( numnamsed_ref, TRIM( clname )//'_ref', 'OLD' , 'FORMATTED', 'SEQUENTIAL', -1, numout, .FALSE. )455 CALL ctl_opn( numnamsed_cfg, TRIM( clname )//'_cfg', 'OLD' , 'FORMATTED', 'SEQUENTIAL', -1, numout, .FALSE. )454 CALL load_nml( numnamsed_ref, TRIM( clname )//'_ref', numout, .FALSE. ) 455 CALL load_nml( numnamsed_cfg, TRIM( clname )//'_cfg', numout, .FALSE. ) 456 456 457 457 nitsed000 = nittrc000 458 458 nitsedend = nitend 459 459 ! Namelist nam_run 460 REWIND( numnamsed_ref ) ! Namelist nam_run in reference namelist : Pisces variables461 460 READ ( numnamsed_ref, nam_run, IOSTAT = ios, ERR = 901) 462 461 901 IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_run in reference namelist' ) 463 462 464 REWIND( numnamsed_cfg ) ! Namelist nam_run in reference namelist : Pisces variables465 463 READ ( numnamsed_cfg, nam_run, IOSTAT = ios, ERR = 902) 466 464 902 IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_run in configuration namelist' ) … … 474 472 IF ( ln_p5z .AND. ln_sed_2way ) CALL ctl_stop( '2 ways coupling with sediment cannot be activated with PISCES-QUOTA' ) 475 473 476 REWIND( numnamsed_ref ) ! Namelist nam_geom in reference namelist : Pisces variables477 474 READ ( numnamsed_ref, nam_geom, IOSTAT = ios, ERR = 903) 478 475 903 IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_geom in reference namelist' ) 479 476 480 REWIND( numnamsed_cfg ) ! Namelist nam_geom in reference namelist : Pisces variables481 477 READ ( numnamsed_cfg, nam_geom, IOSTAT = ios, ERR = 904) 482 478 904 IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_geom in configuration namelist' ) … … 497 493 dtsed = r2dttrc 498 494 499 REWIND( numnamsed_ref ) ! Namelist nam_trased in reference namelist : Pisces variables500 495 READ ( numnamsed_ref, nam_trased, IOSTAT = ios, ERR = 905) 501 496 905 IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_trased in reference namelist' ) 502 497 503 REWIND( numnamsed_cfg ) ! Namelist nam_trased in reference namelist : Pisces variables504 498 READ ( numnamsed_cfg, nam_trased, IOSTAT = ios, ERR = 906) 505 499 906 IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_trased in configuration namelist' ) … … 530 524 ENDIF 531 525 532 REWIND( numnamsed_ref ) ! Namelist nam_diased in reference namelist : Pisces variables533 526 READ ( numnamsed_ref, nam_diased, IOSTAT = ios, ERR = 907) 534 527 907 IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_diased in reference namelist' ) 535 528 536 REWIND( numnamsed_cfg ) ! Namelist nam_diased in reference namelist : Pisces variables537 529 READ ( numnamsed_cfg, nam_diased, IOSTAT = ios, ERR = 908) 538 530 908 IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_diased in configuration namelist' ) … … 572 564 ! Inorganic chemistry parameters 573 565 !---------------------------------- 574 REWIND( numnamsed_ref ) ! Namelist nam_inorg in reference namelist : Pisces variables575 566 READ ( numnamsed_ref, nam_inorg, IOSTAT = ios, ERR = 909) 576 567 909 IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_inorg in reference namelist' ) 577 568 578 REWIND( numnamsed_cfg ) ! Namelist nam_inorg in reference namelist : Pisces variables579 569 READ ( numnamsed_cfg, nam_inorg, IOSTAT = ios, ERR = 910) 580 570 910 IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_inorg in configuration namelist' ) … … 598 588 ! Additional parameter linked to POC/O2/No3/Po4 599 589 !---------------------------------------------- 600 REWIND( numnamsed_ref ) ! Namelist nam_poc in reference namelist : Pisces variables601 590 READ ( numnamsed_ref, nam_poc, IOSTAT = ios, ERR = 911) 602 591 911 IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_poc in reference namelist' ) 603 592 604 REWIND( numnamsed_cfg ) ! Namelist nam_poc in reference namelist : Pisces variables605 593 READ ( numnamsed_cfg, nam_poc, IOSTAT = ios, ERR = 912) 606 594 912 IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_poc in configuration namelist' ) … … 650 638 ! Bioturbation parameter 651 639 !------------------------ 652 REWIND( numnamsed_ref ) ! Namelist nam_btb in reference namelist : Pisces variables653 640 READ ( numnamsed_ref, nam_btb, IOSTAT = ios, ERR = 913) 654 641 913 IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_btb in reference namelist' ) 655 642 656 REWIND( numnamsed_cfg ) ! Namelist nam_btb in reference namelist : Pisces variables657 643 READ ( numnamsed_cfg, nam_btb, IOSTAT = ios, ERR = 914) 658 644 914 IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_btb in configuration namelist' ) … … 671 657 ! Initial value (t=0) for sediment pore water and solid components 672 658 !---------------------------------------------------------------- 673 REWIND( numnamsed_ref ) ! Namelist nam_rst in reference namelist : Pisces variables674 659 READ ( numnamsed_ref, nam_rst, IOSTAT = ios, ERR = 915) 675 660 915 IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_rst in reference namelist' ) 676 661 677 REWIND( numnamsed_cfg ) ! Namelist nam_rst in reference namelist : Pisces variables678 662 READ ( numnamsed_cfg, nam_rst, IOSTAT = ios, ERR = 916) 679 663 916 IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_rst in configuration namelist' ) … … 686 670 nn_dtsed = nn_dttrc 687 671 688 CLOSE( numnamsed_cfg )689 CLOSE( numnamsed_ref )690 672 691 673 END SUBROUTINE sed_init_nam
Note: See TracChangeset
for help on using the changeset viewer.