- Timestamp:
- 2020-02-27T12:27:50+01:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/UKMO/NEMO_4.0_add_pond_lids_prints/src/ICE/icedyn_rhg_evp.F90
r11081 r12475 185 185 186 186 IF( kt == nit000 .AND. lwp ) WRITE(numout,*) '-- ice_dyn_rhg_evp: EVP sea-ice rheology' 187 188 write(numout,*)'ice_dyn_rhg_evp 1: u_ice = ',u_ice(3,4) 187 189 ! 188 190 !!gm for Clem: OPTIMIZATION: I think zfmask can be computed one for all at the initialization.... … … 225 227 CALL lbc_lnk( 'icedyn_rhg_evp', zfmask, 'F', 1._wp ) 226 228 229 write(numout,*)'ice_dyn_rhg_evp 3: u_ice = ',u_ice(3,4) 230 227 231 !------------------------------------------------------------------------------! 228 232 ! 1) define some variables and initialize arrays … … 267 271 ELSE ; zkt = 0._wp 268 272 ENDIF 273 274 write(numout,*)'ice_dyn_rhg_evp 3: u_ice = ',u_ice(3,4) 269 275 ! 270 276 !------------------------------------------------------------------------------! … … 328 334 END DO 329 335 CALL lbc_lnk_multi( 'icedyn_rhg_evp', zmf, 'T', 1., zdt_m, 'T', 1. ) 336 337 write(numout,*)'ice_dyn_rhg_evp 3: u_ice = ',u_ice(3,4) 330 338 ! 331 339 ! !== Landfast ice parameterization ==! … … 368 376 IF( iom_use('tau_icebfr') ) CALL iom_put( 'tau_icebfr', tau_icebfr(:,:) ) 369 377 378 write(numout,*)'ice_dyn_rhg_evp 3: u_ice = ',u_ice(3,4) 379 370 380 !------------------------------------------------------------------------------! 371 381 ! 3) Solution of the momentum equation, iterative procedure … … 397 407 CALL lbc_lnk( 'icedyn_rhg_evp', zds, 'F', 1. ) 398 408 409 write(numout,*)'ice_dyn_rhg_evp 3: u_ice = ',u_ice(3,4) 410 399 411 DO jj = 2, jpj ! loop to jpi,jpj to avoid making a communication for zs1,zs2,zs12 400 412 DO ji = 2, jpi ! no vector loop … … 440 452 END DO 441 453 CALL lbc_lnk( 'icedyn_rhg_evp', zp_delt, 'T', 1. ) 454 455 write(numout,*)'ice_dyn_rhg_evp 3: u_ice = ',u_ice(3,4) 442 456 443 457 DO jj = 1, jpjm1 … … 488 502 END DO 489 503 END DO 504 write(numout,*)'ice_dyn_rhg_evp 3: u_ice = ',u_ice(3,4) 505 490 506 ! 491 507 ! --- Computation of ice velocity --- ! … … 558 574 & ( zmf(ji ,jj) * ( e1v(ji ,jj) * v_ice(ji ,jj) + e1v(ji ,jj-1) * v_ice(ji ,jj-1) ) & 559 575 & + zmf(ji+1,jj) * ( e1v(ji+1,jj) * v_ice(ji+1,jj) + e1v(ji+1,jj-1) * v_ice(ji+1,jj-1) ) ) 576 577 IF (ji == 3 .AND. jj == 4) THEN 578 write(numout,*)'ice_dyn_rhg_evp 6: u_ice(ji,jj), u_oce(ji,jj), v_iceU(ji,jj), v_oceU(ji,jj) = ',u_ice(ji,jj), u_oce(ji,jj), v_iceU(ji,jj), v_oceU(ji,jj) 579 write(numout,*)'ice_dyn_rhg_evp 6: zTauO, zTauU_ib(ji,jj), zvel =',zTauO, zTauU_ib(ji,jj), zvel 580 write(numout,*)'ice_dyn_rhg_evp 6: r1_e1u(ji,jj), zmf(ji,jj), e1v(ji,jj) = ',r1_e1u(ji,jj), zmf(ji,jj), e1v(ji,jj) 581 ENDIF 582 560 583 ! 561 584 ! !--- Sum of external forces (explicit solution) = F + tau_ia + Coriolis + spg + tau_io … … 580 603 & ) * zmaskU(ji,jj) 581 604 ENDIF 605 606 IF (ji == 3 .AND. jj == 4) THEN 607 write(numout,*)'ice_dyn_rhg_evp 7: u_ice(ji,jj), zmU_t(ji,jj), zTauE, zTauO = ',u_ice(ji,jj), zmU_t(ji,jj), zTauE, zTauO 608 write(numout,*)'ice_dyn_rhg_evp 7: zepsi, zTauO, zTauB, rswitch = ', zepsi, zTauO, zTauB, rswitch 609 write(numout,*)'ice_dyn_rhg_evp 7: zdtevp, rn_lfrelax, zswitchU(ji,jj), zmaskU(ji,jj) = ', zdtevp, rn_lfrelax, zswitchU(ji,jj), zmaskU(ji,jj) 610 ENDIF 611 582 612 END DO 583 613 END DO … … 608 638 & ( zmf(ji ,jj) * ( e1v(ji ,jj) * v_ice(ji ,jj) + e1v(ji ,jj-1) * v_ice(ji ,jj-1) ) & 609 639 & + zmf(ji+1,jj) * ( e1v(ji+1,jj) * v_ice(ji+1,jj) + e1v(ji+1,jj-1) * v_ice(ji+1,jj-1) ) ) 640 641 IF (ji == 3 .AND. jj == 4) THEN 642 write(numout,*)'ice_dyn_rhg_evp 8: u_ice(ji,jj), u_oce(ji,jj), v_iceU(ji,jj), v_oceU(ji,jj) = ',u_ice(ji,jj), u_oce(ji,jj), v_iceU(ji,jj), v_oceU(ji,jj) 643 write(numout,*)'ice_dyn_rhg_evp 8: zTauO, zTauU_ib(ji,jj), zvel =',zTauO, zTauU_ib(ji,jj), zvel 644 write(numout,*)'ice_dyn_rhg_evp 8: r1_e1u(ji,jj), zmf(ji,jj), e1v(ji,jj) = ',r1_e1u(ji,jj), zmf(ji,jj), e1v(ji,jj) 645 ENDIF 646 610 647 ! 611 648 ! !--- Sum of external forces (explicit solution) = F + tau_ia + Coriolis + spg + tau_io … … 630 667 & ) * zmaskU(ji,jj) 631 668 ENDIF 669 670 IF (ji == 3 .AND. jj == 4) THEN 671 write(numout,*)'ice_dyn_rhg_evp 9: u_ice(ji,jj), zmU_t(ji,jj), zTauE, zTauO = ',u_ice(ji,jj), zmU_t(ji,jj), zTauE, zTauO 672 write(numout,*)'ice_dyn_rhg_evp 9: zepsi, zTauO, zTauB, rswitch = ', zepsi, zTauO, zTauB, rswitch 673 write(numout,*)'ice_dyn_rhg_evp 9: zdtevp, rn_lfrelax, zswitchU(ji,jj), zmaskU(ji,jj) = ', zdtevp, rn_lfrelax, zswitchU(ji,jj), zmaskU(ji,jj) 674 ENDIF 632 675 END DO 633 676 END DO 634 677 CALL lbc_lnk( 'icedyn_rhg_evp', u_ice, 'U', -1. ) 678 679 write(numout,*)'ice_dyn_rhg_evp 10: u_ice = ',u_ice(3,4) 635 680 ! 636 681 #if defined key_agrif … … 706 751 CALL bdy_ice_dyn( 'V' ) 707 752 ENDIF 753 write(numout,*)'ice_dyn_rhg_evp 5: u_ice = ',u_ice(3,4) 708 754 ! 709 755 !------------------------------------------------------------------------------! … … 757 803 pstress2_i (:,:) = zs2 (:,:) 758 804 pstress12_i(:,:) = zs12(:,:) 805 write(numout,*)'ice_dyn_rhg_evp 4: u_ice = ',u_ice(3,4) 759 806 ! 760 807 … … 853 900 END DO 854 901 END DO 902 903 write(numout,*)'ice_dyn_rhg_evp 11: u_ice = ',u_ice(3,4) 855 904 856 905 CALL lbc_lnk_multi( 'icedyn_rhg_evp', zdiag_sig1 , 'T', 1., zdiag_sig2 , 'T', 1., & … … 887 936 888 937 ENDIF 938 939 write(numout,*)'ice_dyn_rhg_evp 2: u_ice = ',u_ice(3,4) 889 940 ! 890 941 END SUBROUTINE ice_dyn_rhg_evp
Note: See TracChangeset
for help on using the changeset viewer.
