- Timestamp:
- 2020-09-14T17:40:34+02:00 (4 years ago)
- Location:
- NEMO/branches/2019/dev_r11351_fldread_with_XIOS
- Files:
-
- 2 edited
-
. (modified) (1 prop)
-
src/TOP/PISCES/SED/sedini.F90 (modified) (14 diffs)
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2019/dev_r11351_fldread_with_XIOS
- Property svn:externals
-
old new 3 3 ^/utils/build/mk@HEAD mk 4 4 ^/utils/tools@HEAD tools 5 ^/vendors/AGRIF/dev @HEADext/AGRIF5 ^/vendors/AGRIF/dev_r12970_AGRIF_CMEMS ext/AGRIF 6 6 ^/vendors/FCM@HEAD ext/FCM 7 7 ^/vendors/IOIPSL@HEAD ext/IOIPSL 8 9 # SETTE 10 ^/utils/CI/sette@13382 sette
-
- Property svn:externals
-
NEMO/branches/2019/dev_r11351_fldread_with_XIOS/src/TOP/PISCES/SED/sedini.F90
r10362 r13463 13 13 USE sedarr 14 14 USE sedadv 15 USE trc_oce, ONLY : nn_dttrc16 15 USE trcdmp_sed 17 16 USE trcdta … … 23 22 PRIVATE 24 23 24 !! * Substitutions 25 # include "do_loop_substitute.h90" 25 26 !! Module variables 26 27 REAL(wp) :: & … … 134 135 ! Determination of sediments number of points and allocate global variables 135 136 epkbot(:,:) = 0. 136 DO jj = 1, jpj 137 DO ji = 1, jpi 138 ikt = mbkt(ji,jj) 139 IF( tmask(ji,jj,ikt) == 1 ) epkbot(ji,jj) = e3t_1d(ikt) 140 gdepbot(ji,jj) = gdepw_0(ji,jj,ikt) 141 ENDDO 142 ENDDO 137 DO_2D( 1, 1, 1, 1 ) 138 ikt = mbkt(ji,jj) 139 IF( tmask(ji,jj,ikt) == 1 ) epkbot(ji,jj) = e3t_1d(ikt) 140 gdepbot(ji,jj) = gdepw_0(ji,jj,ikt) 141 END_2D 143 142 144 143 ! computation of total number of ocean points … … 248 247 ! Computation of 1D array of sediments points 249 248 indoce = 0 250 DO jj = 1, jpj 251 DO ji = 1, jpi 252 IF ( epkbot(ji,jj) > 0. ) THEN 253 indoce = indoce + 1 254 iarroce(indoce) = (jj - 1) * jpi + ji 255 ENDIF 256 END DO 257 END DO 249 DO_2D( 1, 1, 1, 1 ) 250 IF ( epkbot(ji,jj) > 0. ) THEN 251 indoce = indoce + 1 252 iarroce(indoce) = (jj - 1) * jpi + ji 253 ENDIF 254 END_2D 258 255 259 256 IF ( indoce .EQ. 0 ) THEN … … 406 403 !!---------------------------------------------------------------------- 407 404 408 INTEGER :: numnamsed_ref = -1 !! Logical units fornamelist sediment409 INTEGER :: numnamsed_cfg = -1 !! Logical units fornamelist sediment405 CHARACTER(:), ALLOCATABLE :: numnamsed_ref !! Character buffer for reference namelist sediment 406 CHARACTER(:), ALLOCATABLE :: numnamsed_cfg !! Character buffer for configuration namelist sediment 410 407 INTEGER :: ios ! Local integer output status for namelist read 411 408 CHARACTER(LEN=20) :: clname … … 452 449 IF(lwp) WRITE(numsed,*) ' sed_init_nam : read SEDIMENT namelist' 453 450 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.)451 CALL load_nml( numnamsed_ref, TRIM( clname )//'_ref', numout, lwm ) 452 CALL load_nml( numnamsed_cfg, TRIM( clname )//'_cfg', numout, lwm ) 456 453 457 454 nitsed000 = nittrc000 458 455 nitsedend = nitend 459 456 ! Namelist nam_run 460 REWIND( numnamsed_ref ) ! Namelist nam_run in reference namelist : Pisces variables461 457 READ ( numnamsed_ref, nam_run, IOSTAT = ios, ERR = 901) 462 901 IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_run in reference namelist', lwp ) 463 464 REWIND( numnamsed_cfg ) ! Namelist nam_run in reference namelist : Pisces variables 458 901 IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_run in reference namelist' ) 459 465 460 READ ( numnamsed_cfg, nam_run, IOSTAT = ios, ERR = 902) 466 902 IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_run in configuration namelist' , lwp)461 902 IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_run in configuration namelist' ) 467 462 468 463 IF (lwp) THEN … … 474 469 IF ( ln_p5z .AND. ln_sed_2way ) CALL ctl_stop( '2 ways coupling with sediment cannot be activated with PISCES-QUOTA' ) 475 470 476 REWIND( numnamsed_ref ) ! Namelist nam_geom in reference namelist : Pisces variables477 471 READ ( numnamsed_ref, nam_geom, IOSTAT = ios, ERR = 903) 478 903 IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_geom in reference namelist', lwp ) 479 480 REWIND( numnamsed_cfg ) ! Namelist nam_geom in reference namelist : Pisces variables 472 903 IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_geom in reference namelist' ) 473 481 474 READ ( numnamsed_cfg, nam_geom, IOSTAT = ios, ERR = 904) 482 904 IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_geom in configuration namelist' , lwp)475 904 IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_geom in configuration namelist' ) 483 476 484 477 IF (lwp) THEN … … 495 488 496 489 jpksedm1 = jpksed - 1 497 dtsed = r2dttrc 498 499 REWIND( numnamsed_ref ) ! Namelist nam_trased in reference namelist : Pisces variables 490 dtsed = rDt_trc 491 500 492 READ ( numnamsed_ref, nam_trased, IOSTAT = ios, ERR = 905) 501 905 IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_trased in reference namelist', lwp ) 502 503 REWIND( numnamsed_cfg ) ! Namelist nam_trased in reference namelist : Pisces variables 493 905 IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_trased in reference namelist' ) 494 504 495 READ ( numnamsed_cfg, nam_trased, IOSTAT = ios, ERR = 906) 505 906 IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_trased in configuration namelist' , lwp)496 906 IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_trased in configuration namelist' ) 506 497 507 498 DO jn = 1, jpsol … … 530 521 ENDIF 531 522 532 REWIND( numnamsed_ref ) ! Namelist nam_diased in reference namelist : Pisces variables533 523 READ ( numnamsed_ref, nam_diased, IOSTAT = ios, ERR = 907) 534 907 IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_diased in reference namelist', lwp ) 535 536 REWIND( numnamsed_cfg ) ! Namelist nam_diased in reference namelist : Pisces variables 524 907 IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_diased in reference namelist' ) 525 537 526 READ ( numnamsed_cfg, nam_diased, IOSTAT = ios, ERR = 908) 538 908 IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_diased in configuration namelist' , lwp)527 908 IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_diased in configuration namelist' ) 539 528 540 529 DO jn = 1, jpdia3dsed … … 572 561 ! Inorganic chemistry parameters 573 562 !---------------------------------- 574 REWIND( numnamsed_ref ) ! Namelist nam_inorg in reference namelist : Pisces variables575 563 READ ( numnamsed_ref, nam_inorg, IOSTAT = ios, ERR = 909) 576 909 IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_inorg in reference namelist', lwp ) 577 578 REWIND( numnamsed_cfg ) ! Namelist nam_inorg in reference namelist : Pisces variables 564 909 IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_inorg in reference namelist' ) 565 579 566 READ ( numnamsed_cfg, nam_inorg, IOSTAT = ios, ERR = 910) 580 910 IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_inorg in configuration namelist' , lwp)567 910 IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_inorg in configuration namelist' ) 581 568 582 569 IF (lwp) THEN … … 598 585 ! Additional parameter linked to POC/O2/No3/Po4 599 586 !---------------------------------------------- 600 REWIND( numnamsed_ref ) ! Namelist nam_poc in reference namelist : Pisces variables601 587 READ ( numnamsed_ref, nam_poc, IOSTAT = ios, ERR = 911) 602 911 IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_poc in reference namelist', lwp ) 603 604 REWIND( numnamsed_cfg ) ! Namelist nam_poc in reference namelist : Pisces variables 588 911 IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_poc in reference namelist' ) 589 605 590 READ ( numnamsed_cfg, nam_poc, IOSTAT = ios, ERR = 912) 606 912 IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_poc in configuration namelist' , lwp)591 912 IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_poc in configuration namelist' ) 607 592 608 593 IF (lwp) THEN … … 650 635 ! Bioturbation parameter 651 636 !------------------------ 652 REWIND( numnamsed_ref ) ! Namelist nam_btb in reference namelist : Pisces variables653 637 READ ( numnamsed_ref, nam_btb, IOSTAT = ios, ERR = 913) 654 913 IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_btb in reference namelist', lwp ) 655 656 REWIND( numnamsed_cfg ) ! Namelist nam_btb in reference namelist : Pisces variables 638 913 IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_btb in reference namelist' ) 639 657 640 READ ( numnamsed_cfg, nam_btb, IOSTAT = ios, ERR = 914) 658 914 IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_btb in configuration namelist' , lwp)641 914 IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_btb in configuration namelist' ) 659 642 660 643 IF (lwp) THEN … … 671 654 ! Initial value (t=0) for sediment pore water and solid components 672 655 !---------------------------------------------------------------- 673 REWIND( numnamsed_ref ) ! Namelist nam_rst in reference namelist : Pisces variables674 656 READ ( numnamsed_ref, nam_rst, IOSTAT = ios, ERR = 915) 675 915 IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_rst in reference namelist', lwp ) 676 677 REWIND( numnamsed_cfg ) ! Namelist nam_rst in reference namelist : Pisces variables 657 915 IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_rst in reference namelist' ) 658 678 659 READ ( numnamsed_cfg, nam_rst, IOSTAT = ios, ERR = 916) 679 916 IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_rst in configuration namelist' , lwp)660 916 IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_rst in configuration namelist' ) 680 661 681 662 IF (lwp) THEN … … 684 665 WRITE(numsed,*) ' ' 685 666 ENDIF 686 nn_dtsed = nn_dttrc 687 688 CLOSE( numnamsed_cfg ) 689 CLOSE( numnamsed_ref ) 667 nn_dtsed = 1 668 690 669 691 670 END SUBROUTINE sed_init_nam
Note: See TracChangeset
for help on using the changeset viewer.
