- Timestamp:
- 2017-09-07T20:08:11+02:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3/iceitd.F90
r8505 r8512 38 38 PUBLIC ice_itd_reb ! called in iceerr 39 39 40 ! ** ice-thickness distribution namelist (namice itd) **40 ! ** ice-thickness distribution namelist (namice_itd) ** 41 41 REAL(wp) :: rn_himean ! mean thickness of the domain (used to compute the distribution) 42 42 … … 79 79 !!------------------------------------------------------------------ 80 80 81 IF( kt == nit000 .AND. lwp) THEN 82 WRITE(numout,*) 83 WRITE(numout,*) 'ice_itd_rem : Remapping the ice thickness distribution' 84 WRITE(numout,*) '~~~~~~~~~~~~~~~' 85 ENDIF 81 IF( kt == nit000 .AND. lwp ) WRITE(numout,*) '-- ice_itd_rem: remapping ice thickness distribution' 86 82 87 83 IF( ln_limdiachk ) CALL ice_cons_hsm(0, 'iceitd_rem', zvi_b, zsmv_b, zei_b, zfw_b, zfs_b, zft_b) … … 545 541 546 542 547 SUBROUTINE ice_itd_reb 543 SUBROUTINE ice_itd_reb( kt ) 548 544 !!------------------------------------------------------------------ 549 545 !! *** ROUTINE ice_itd_reb *** … … 555 551 !! to the neighboring category 556 552 !!------------------------------------------------------------------ 553 INTEGER , INTENT (in) :: kt ! Ocean time step 557 554 INTEGER :: ji, jj, jl ! dummy loop indices 558 555 ! … … 561 558 !!------------------------------------------------------------------ 562 559 ! 560 IF( kt == nit000 .AND. lwp ) WRITE(numout,*) '-- ice_itd_reb: rebining ice thickness distribution' 561 563 562 jdonor(:,:) = 0 564 563 zdaice(:,:) = 0._wp … … 645 644 !! ** Purpose : Initializes the ice thickness distribution 646 645 !! ** Method : ... 647 !! ** input : Namelist namice itd646 !! ** input : Namelist namice_itd 648 647 !!------------------------------------------------------------------- 649 648 INTEGER :: jl ! dummy loop index … … 651 650 REAL(wp) :: zhmax, znum, zden, zalpha ! - - 652 651 !! 653 NAMELIST/namice itd/ rn_himean654 !!------------------------------------------------------------------ 655 ! 656 REWIND( numnam_ice_ref ) ! Namelist namice itd in reference namelist : Parameters for ice657 READ ( numnam_ice_ref, namice itd, IOSTAT = ios, ERR = 901)658 901 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namice itd in reference namelist', lwp )659 660 REWIND( numnam_ice_cfg ) ! Namelist namice itd in configuration namelist : Parameters for ice661 READ ( numnam_ice_cfg, namice itd, IOSTAT = ios, ERR = 902 )662 902 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namice itd in configuration namelist', lwp )663 IF(lwm) WRITE ( numoni, namice itd )652 NAMELIST/namice_itd/ rn_himean 653 !!------------------------------------------------------------------ 654 ! 655 REWIND( numnam_ice_ref ) ! Namelist namice_itd in reference namelist : Parameters for ice 656 READ ( numnam_ice_ref, namice_itd, IOSTAT = ios, ERR = 901) 657 901 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namice_itd in reference namelist', lwp ) 658 659 REWIND( numnam_ice_cfg ) ! Namelist namice_itd in configuration namelist : Parameters for ice 660 READ ( numnam_ice_cfg, namice_itd, IOSTAT = ios, ERR = 902 ) 661 902 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namice_itd in configuration namelist', lwp ) 662 IF(lwm) WRITE ( numoni, namice_itd ) 664 663 ! 665 664 IF(lwp) THEN ! control print … … 667 666 WRITE(numout,*) 'ice_itd_init : Initialization of ice cat distribution ' 668 667 WRITE(numout,*) '~~~~~~~~~~~~' 669 WRITE(numout,*) ' Namelist namice run: '668 WRITE(numout,*) ' Namelist namice_itd : ' 670 669 WRITE(numout,*) ' mean ice thickness in the domain rn_himean = ', rn_himean 671 670 ENDIF
Note: See TracChangeset
for help on using the changeset viewer.