Changeset 13295 for NEMO/trunk/src/ICE
- Timestamp:
- 2020-07-10T20:24:21+02:00 (4 years ago)
- Location:
- NEMO/trunk/src/ICE
- Files:
-
- 16 edited
-
icealb.F90 (modified) (1 diff)
-
icecor.F90 (modified) (2 diffs)
-
icectl.F90 (modified) (9 diffs)
-
icedyn.F90 (modified) (2 diffs)
-
icedyn_adv_pra.F90 (modified) (16 diffs)
-
icedyn_adv_umx.F90 (modified) (54 diffs)
-
icedyn_rdgrft.F90 (modified) (4 diffs)
-
icedyn_rhg_evp.F90 (modified) (18 diffs)
-
iceistate.F90 (modified) (4 diffs)
-
iceitd.F90 (modified) (3 diffs)
-
icesbc.F90 (modified) (1 diff)
-
icethd.F90 (modified) (4 diffs)
-
icethd_do.F90 (modified) (2 diffs)
-
iceupdate.F90 (modified) (3 diffs)
-
icevar.F90 (modified) (9 diffs)
-
icewri.F90 (modified) (3 diffs)
Legend:
- Unmodified
- Added
- Removed
-
NEMO/trunk/src/ICE/icealb.F90
r12377 r13295 119 119 ! 120 120 DO jl = 1, jpl 121 DO_2D _11_11121 DO_2D( 1, 1, 1, 1 ) 122 122 ! !--- Specific snow, ice and pond fractions (for now, we prevent melt ponds and snow at the same time) 123 123 IF( ph_snw(ji,jj,jl) == 0._wp ) THEN -
NEMO/trunk/src/ICE/icecor.F90
r13226 r13295 88 88 zzc = rhoi * r1_Dt_ice 89 89 DO jl = 1, jpl 90 DO_2D _11_1190 DO_2D( 1, 1, 1, 1 ) 91 91 zsal = sv_i(ji,jj,jl) 92 92 sv_i(ji,jj,jl) = MIN( MAX( rn_simin*v_i(ji,jj,jl) , sv_i(ji,jj,jl) ) , rn_simax*v_i(ji,jj,jl) ) … … 106 106 ! !----------------------------------------------------- 107 107 IF( kn == 2 ) THEN ! Ice drift case: Corrections to avoid wrong values ! 108 DO_2D _00_00108 DO_2D( 0, 0, 0, 0 ) 109 109 IF ( at_i(ji,jj) == 0._wp ) THEN ! what to do if there is no ice 110 110 IF ( at_i(ji+1,jj) == 0._wp ) u_ice(ji ,jj) = 0._wp ! right side -
NEMO/trunk/src/ICE/icectl.F90
r13286 r13295 368 368 cl_alname(ialert_id) = ' Incompat vol and con ' ! name of the alert 369 369 DO jl = 1, jpl 370 DO_2D _11_11370 DO_2D( 1, 1, 1, 1 ) 371 371 IF( v_i(ji,jj,jl) /= 0._wp .AND. a_i(ji,jj,jl) == 0._wp ) THEN 372 372 WRITE(numout,*) ' ALERTE 2 : Incompatible volume and concentration ' … … 380 380 cl_alname(ialert_id) = ' Very thick ice ' ! name of the alert 381 381 jl = jpl 382 DO_2D _11_11382 DO_2D( 1, 1, 1, 1 ) 383 383 IF( h_i(ji,jj,jl) > 50._wp ) THEN 384 384 WRITE(numout,*) ' ALERTE 3 : Very thick ice' … … 391 391 ialert_id = 4 ! reference number of this alert 392 392 cl_alname(ialert_id) = ' Very fast ice ' ! name of the alert 393 DO_2D _11_11393 DO_2D( 1, 1, 1, 1 ) 394 394 IF( MAX( ABS( u_ice(ji,jj) ), ABS( v_ice(ji,jj) ) ) > 2. .AND. & 395 395 & at_i(ji,jj) > 0._wp ) THEN … … 403 403 ialert_id = 5 ! reference number of this alert 404 404 cl_alname(ialert_id) = ' High salt flux ' ! name of the alert 405 DO_2D _11_11405 DO_2D( 1, 1, 1, 1 ) 406 406 IF( ABS( sfx (ji,jj) ) > 1.0e-2 ) THEN ! = 1 psu/day for 1m ocean depth 407 407 WRITE(numout,*) ' ALERTE 5 : High salt flux' … … 414 414 ialert_id = 6 ! reference number of this alert 415 415 cl_alname(ialert_id) = ' Ice on continents ' ! name of the alert 416 DO_2D _11_11416 DO_2D( 1, 1, 1, 1 ) 417 417 IF( tmask(ji,jj,1) <= 0._wp .AND. at_i(ji,jj) > 0._wp ) THEN 418 418 WRITE(numout,*) ' ALERTE 6 : Ice on continents' … … 427 427 cl_alname(ialert_id) = ' Very fresh ice ' ! name of the alert 428 428 DO jl = 1, jpl 429 DO_2D _11_11429 DO_2D( 1, 1, 1, 1 ) 430 430 IF( s_i(ji,jj,jl) < 0.1 .AND. a_i(ji,jj,jl) > 0._wp ) THEN 431 431 WRITE(numout,*) ' ALERTE 7 : Very fresh ice' … … 439 439 ialert_id = 8 ! reference number of this alert 440 440 cl_alname(ialert_id) = ' fnsolar very big ' ! name of the alert 441 DO_2D _11_11441 DO_2D( 1, 1, 1, 1 ) 442 442 IF( ABS( qns(ji,jj) ) > 1500._wp .AND. at_i(ji,jj) > 0._wp ) THEN 443 443 ! … … 454 454 cl_alname(ialert_id) = ' Very old ice ' ! name of the alert 455 455 DO jl = 1, jpl 456 DO_2D _11_11456 DO_2D( 1, 1, 1, 1 ) 457 457 IF ( ( ( ABS( o_i(ji,jj,jl) ) > rDt_ice ) .OR. & 458 458 ( ABS( o_i(ji,jj,jl) ) < 0._wp) ) .AND. & … … 470 470 inb_alp(ialert_id) = 0 471 471 DO jl = 1, jpl 472 DO_3D _11_11(1, nlay_i )472 DO_3D( 1, 1, 1, 1, 1, nlay_i ) 473 473 ztmelts = -rTmlt * sz_i(ji,jj,jk,jl) + rt0 474 474 IF( t_i(ji,jj,jk,jl) > ztmelts .AND. v_i(ji,jj,jl) > 1.e-10 & -
NEMO/trunk/src/ICE/icedyn.F90
r13226 r13295 126 126 ! CFL = 0.5 at a distance from the bound of 1/6 of the basin length 127 127 ! Then for dx = 2m and dt = 1s => rn_uice = u (1/6th) = 1m/s 128 DO_2D _11_11128 DO_2D( 1, 1, 1, 1 ) 129 129 zcoefu = ( REAL(jpiglo+1)*0.5 - REAL(ji+nimpp-1) ) / ( REAL(jpiglo+1)*0.5 - 1. ) 130 130 zcoefv = ( REAL(jpjglo+1)*0.5 - REAL(jj+njmpp-1) ) / ( REAL(jpjglo+1)*0.5 - 1. ) … … 155 155 156 156 ALLOCATE( zdivu_i(jpi,jpj) ) 157 DO_2D _00_00157 DO_2D( 0, 0, 0, 0 ) 158 158 zdivu_i(ji,jj) = ( e2u(ji,jj) * u_ice(ji,jj) - e2u(ji-1,jj) * u_ice(ji-1,jj) & 159 159 & + e1v(ji,jj) * v_ice(ji,jj) - e1v(ji,jj-1) * v_ice(ji,jj-1) ) * r1_e1e2t(ji,jj) -
NEMO/trunk/src/ICE/icedyn_adv_pra.F90
r13286 r13295 102 102 ! --- Record max of the surrounding 9-pts ice thick. (for call Hbig) --- ! 103 103 DO jl = 1, jpl 104 DO_2D _00_00104 DO_2D( 0, 0, 0, 0 ) 105 105 zhip_max(ji,jj,jl) = MAX( epsi20, ph_ip(ji,jj,jl), ph_ip(ji+1,jj ,jl), ph_ip(ji ,jj+1,jl), & 106 106 & ph_ip(ji-1,jj ,jl), ph_ip(ji ,jj-1,jl), & … … 250 250 ! derive open water from ice concentration 251 251 zati2(:,:) = SUM( pa_i(:,:,:), dim=3 ) 252 DO_2D _00_00252 DO_2D( 0, 0, 0, 0 ) 253 253 pato_i(ji,jj) = pato_i(ji,jj) - ( zati2(ji,jj) - zati1(ji,jj) ) & !--- open water 254 254 & - ( zudy(ji,jj) - zudy(ji-1,jj) + zvdx(ji,jj) - zvdx(ji,jj-1) ) * r1_e1e2t(ji,jj) * zdt … … 305 305 ! 306 306 ! Limitation of moments. 307 DO_2D _00_11307 DO_2D( 0, 0, 1, 1 ) 308 308 ! Initialize volumes of boxes (=area if adv_x first called, =psm otherwise) 309 309 psm (ji,jj,jl) = MAX( pcrh * e1e2t(ji,jj) + ( 1.0 - pcrh ) * psm(ji,jj,jl) , epsi20 ) … … 325 325 326 326 ! Calculate fluxes and moments between boxes i<-->i+1 327 DO_2D _00_11327 DO_2D( 0, 0, 1, 1 ) 328 328 zbet(ji,jj) = MAX( 0._wp, SIGN( 1._wp, put(ji,jj) ) ) 329 329 zalf = MAX( 0._wp, put(ji,jj) ) * pdt / psm(ji,jj,jl) … … 350 350 END_2D 351 351 352 DO_2D _00_10352 DO_2D( 0, 0, 1, 0 ) 353 353 zalf = MAX( 0._wp, -put(ji,jj) ) * pdt / psm(ji+1,jj,jl) 354 354 zalg (ji,jj) = zalf … … 369 369 END_2D 370 370 371 DO_2D _00_00371 DO_2D( 0, 0, 0, 0 ) 372 372 zbt = zbet(ji-1,jj) 373 373 zbt1 = 1.0 - zbet(ji-1,jj) … … 383 383 384 384 ! Put the temporary moments into appropriate neighboring boxes. 385 DO_2D _00_00385 DO_2D( 0, 0, 0, 0 ) 386 386 zbt = zbet(ji-1,jj) 387 387 zbt1 = 1.0 - zbet(ji-1,jj) … … 403 403 END_2D 404 404 405 DO_2D _00_00405 DO_2D( 0, 0, 0, 0 ) 406 406 zbt = zbet(ji,jj) 407 407 zbt1 = 1.0 - zbet(ji,jj) … … 462 462 ! 463 463 ! Limitation of moments. 464 DO_2D _11_00464 DO_2D( 1, 1, 0, 0 ) 465 465 ! Initialize volumes of boxes (=area if adv_x first called, =psm otherwise) 466 466 psm(ji,jj,jl) = MAX( pcrh * e1e2t(ji,jj) + ( 1.0 - pcrh ) * psm(ji,jj,jl) , epsi20 ) … … 482 482 483 483 ! Calculate fluxes and moments between boxes j<-->j+1 484 DO_2D _11_00484 DO_2D( 1, 1, 0, 0 ) 485 485 zbet(ji,jj) = MAX( 0._wp, SIGN( 1._wp, pvt(ji,jj) ) ) 486 486 zalf = MAX( 0._wp, pvt(ji,jj) ) * pdt / psm(ji,jj,jl) … … 507 507 END_2D 508 508 ! 509 DO_2D _10_00509 DO_2D( 1, 0, 0, 0 ) 510 510 zalf = MAX( 0._wp, -pvt(ji,jj) ) * pdt / psm(ji,jj+1,jl) 511 511 zalg (ji,jj) = zalf … … 527 527 528 528 ! Readjust moments remaining in the box. 529 DO_2D _00_00529 DO_2D( 0, 0, 0, 0 ) 530 530 zbt = zbet(ji,jj-1) 531 531 zbt1 = ( 1.0 - zbet(ji,jj-1) ) … … 541 541 542 542 ! Put the temporary moments into appropriate neighboring boxes. 543 DO_2D _00_00543 DO_2D( 0, 0, 0, 0 ) 544 544 zbt = zbet(ji,jj-1) 545 545 zbt1 = 1.0 - zbet(ji,jj-1) … … 562 562 END_2D 563 563 564 DO_2D _00_00564 DO_2D( 0, 0, 0, 0 ) 565 565 zbt = zbet(ji,jj) 566 566 zbt1 = 1.0 - zbet(ji,jj) … … 618 618 DO jl = 1, jpl 619 619 620 DO_2D _11_11620 DO_2D( 1, 1, 1, 1 ) 621 621 IF ( pv_i(ji,jj,jl) > 0._wp ) THEN 622 622 ! … … 684 684 ! -- check snow load -- ! 685 685 DO jl = 1, jpl 686 DO_2D _11_11686 DO_2D( 1, 1, 1, 1 ) 687 687 IF ( pv_i(ji,jj,jl) > 0._wp ) THEN 688 688 ! -
NEMO/trunk/src/ICE/icedyn_adv_umx.F90
r13226 r13295 107 107 ! --- Record max of the surrounding 9-pts ice thick. (for call Hbig) --- ! 108 108 DO jl = 1, jpl 109 DO_2D _00_00109 DO_2D( 0, 0, 0, 0 ) 110 110 zhip_max(ji,jj,jl) = MAX( epsi20, ph_ip(ji,jj,jl), ph_ip(ji+1,jj ,jl), ph_ip(ji ,jj+1,jl), & 111 111 & ph_ip(ji-1,jj ,jl), ph_ip(ji ,jj-1,jl), & … … 150 150 ! 151 151 ! --- define velocity for advection: u*grad(H) --- ! 152 DO_2D _00_00152 DO_2D( 0, 0, 0, 0 ) 153 153 IF ( pu_ice(ji,jj) * pu_ice(ji-1,jj) <= 0._wp ) THEN ; zcu_box(ji,jj) = 0._wp 154 154 ELSEIF( pu_ice(ji,jj) > 0._wp ) THEN ; zcu_box(ji,jj) = pu_ice(ji-1,jj) … … 183 183 IF( .NOT. ALLOCATED(jmsk_small) ) ALLOCATE( jmsk_small(jpi,jpj,jpl) ) 184 184 DO jl = 1, jpl 185 DO_2D _10_10185 DO_2D( 1, 0, 1, 0 ) 186 186 zvi_cen = 0.5_wp * ( pv_i(ji+1,jj,jl) + pv_i(ji,jj,jl) ) 187 187 IF( zvi_cen < epsi06) THEN ; imsk_small(ji,jj,jl) = 0 … … 332 332 !== Open water area ==! 333 333 zati2(:,:) = SUM( pa_i(:,:,:), dim=3 ) 334 DO_2D _00_00334 DO_2D( 0, 0, 0, 0 ) 335 335 pato_i(ji,jj) = pato_i(ji,jj) - ( zati2(ji,jj) - zati1(ji,jj) ) & 336 336 & - ( zudy(ji,jj) - zudy(ji-1,jj) + zvdx(ji,jj) - zvdx(ji,jj-1) ) * r1_e1e2t(ji,jj) * zdt … … 441 441 IF( pamsk == 0._wp ) THEN 442 442 DO jl = 1, jpl 443 DO_2D _10_10443 DO_2D( 1, 0, 1, 0 ) 444 444 IF( ABS( pu(ji,jj) ) > epsi10 ) THEN 445 445 zfu_ho (ji,jj,jl) = zfu_ho (ji,jj,jl) * puc (ji,jj,jl) / pu(ji,jj) … … 463 463 ! thus we calculate the upstream solution and apply a limiter again 464 464 DO jl = 1, jpl 465 DO_2D _00_00465 DO_2D( 0, 0, 0, 0 ) 466 466 ztra = - ( zfu_ups(ji,jj,jl) - zfu_ups(ji-1,jj,jl) + zfv_ups(ji,jj,jl) - zfv_ups(ji,jj-1,jl) ) 467 467 ! … … 484 484 IF( PRESENT( pua_ho ) ) THEN 485 485 DO jl = 1, jpl 486 DO_2D _10_10486 DO_2D( 1, 0, 1, 0 ) 487 487 pua_ho (ji,jj,jl) = zfu_ho (ji,jj,jl) ; pva_ho (ji,jj,jl) = zfv_ho (ji,jj,jl) 488 488 pua_ups(ji,jj,jl) = zfu_ups(ji,jj,jl) ; pva_ups(ji,jj,jl) = zfv_ups(ji,jj,jl) … … 494 494 ! --------------------------------- 495 495 DO jl = 1, jpl 496 DO_2D _00_00496 DO_2D( 0, 0, 0, 0 ) 497 497 ztra = - ( zfu_ho(ji,jj,jl) - zfu_ho(ji-1,jj,jl) + zfv_ho(ji,jj,jl) - zfv_ho(ji,jj-1,jl) ) 498 498 ! … … 528 528 ! 529 529 DO jl = 1, jpl 530 DO_2D _10_10530 DO_2D( 1, 0, 1, 0 ) 531 531 pfu_ups(ji,jj,jl) = MAX( pu(ji,jj), 0._wp ) * pt(ji,jj,jl) + MIN( pu(ji,jj), 0._wp ) * pt(ji+1,jj,jl) 532 532 pfv_ups(ji,jj,jl) = MAX( pv(ji,jj), 0._wp ) * pt(ji,jj,jl) + MIN( pv(ji,jj), 0._wp ) * pt(ji,jj+1,jl) … … 539 539 ! 540 540 DO jl = 1, jpl !-- flux in x-direction 541 DO_2D _10_10541 DO_2D( 1, 0, 1, 0 ) 542 542 pfu_ups(ji,jj,jl) = MAX( pu(ji,jj), 0._wp ) * pt(ji,jj,jl) + MIN( pu(ji,jj), 0._wp ) * pt(ji+1,jj,jl) 543 543 END_2D … … 545 545 ! 546 546 DO jl = 1, jpl !-- first guess of tracer from u-flux 547 DO_2D _00_00547 DO_2D( 0, 0, 0, 0 ) 548 548 ztra = - ( pfu_ups(ji,jj,jl) - pfu_ups(ji-1,jj,jl) ) & 549 549 & + ( pu (ji,jj ) - pu (ji-1,jj ) ) * pt(ji,jj,jl) * (1.-pamsk) … … 555 555 ! 556 556 DO jl = 1, jpl !-- flux in y-direction 557 DO_2D _10_10557 DO_2D( 1, 0, 1, 0 ) 558 558 pfv_ups(ji,jj,jl) = MAX( pv(ji,jj), 0._wp ) * zpt(ji,jj,jl) + MIN( pv(ji,jj), 0._wp ) * zpt(ji,jj+1,jl) 559 559 END_2D … … 563 563 ! 564 564 DO jl = 1, jpl !-- flux in y-direction 565 DO_2D _10_10565 DO_2D( 1, 0, 1, 0 ) 566 566 pfv_ups(ji,jj,jl) = MAX( pv(ji,jj), 0._wp ) * pt(ji,jj,jl) + MIN( pv(ji,jj), 0._wp ) * pt(ji,jj+1,jl) 567 567 END_2D … … 569 569 ! 570 570 DO jl = 1, jpl !-- first guess of tracer from v-flux 571 DO_2D _00_00571 DO_2D( 0, 0, 0, 0 ) 572 572 ztra = - ( pfv_ups(ji,jj,jl) - pfv_ups(ji,jj-1,jl) ) & 573 573 & + ( pv (ji,jj ) - pv (ji,jj-1 ) ) * pt(ji,jj,jl) * (1.-pamsk) … … 579 579 ! 580 580 DO jl = 1, jpl !-- flux in x-direction 581 DO_2D _10_10581 DO_2D( 1, 0, 1, 0 ) 582 582 pfu_ups(ji,jj,jl) = MAX( pu(ji,jj), 0._wp ) * zpt(ji,jj,jl) + MIN( pu(ji,jj), 0._wp ) * zpt(ji+1,jj,jl) 583 583 END_2D … … 589 589 ! 590 590 DO jl = 1, jpl !-- after tracer with upstream scheme 591 DO_2D _00_00591 DO_2D( 0, 0, 0, 0 ) 592 592 ztra = - ( pfu_ups(ji,jj,jl) - pfu_ups(ji-1,jj ,jl) & 593 593 & + pfv_ups(ji,jj,jl) - pfv_ups(ji ,jj-1,jl) ) & … … 628 628 ! 629 629 DO jl = 1, jpl 630 DO_2D _10_10630 DO_2D( 1, 0, 1, 0 ) 631 631 pfu_ho(ji,jj,jl) = 0.5_wp * pu(ji,jj) * ( pt(ji,jj,jl) + pt(ji+1,jj ,jl) ) 632 632 pfv_ho(ji,jj,jl) = 0.5_wp * pv(ji,jj) * ( pt(ji,jj,jl) + pt(ji ,jj+1,jl) ) … … 646 646 ! 647 647 DO jl = 1, jpl !-- flux in x-direction 648 DO_2D _10_10648 DO_2D( 1, 0, 1, 0 ) 649 649 pfu_ho(ji,jj,jl) = 0.5_wp * pu(ji,jj) * ( pt(ji,jj,jl) + pt(ji+1,jj,jl) ) 650 650 END_2D … … 653 653 654 654 DO jl = 1, jpl !-- first guess of tracer from u-flux 655 DO_2D _00_00655 DO_2D( 0, 0, 0, 0 ) 656 656 ztra = - ( pfu_ho(ji,jj,jl) - pfu_ho(ji-1,jj,jl) ) & 657 657 & + ( pu (ji,jj ) - pu (ji-1,jj ) ) * pt(ji,jj,jl) * (1.-pamsk) … … 663 663 664 664 DO jl = 1, jpl !-- flux in y-direction 665 DO_2D _10_10665 DO_2D( 1, 0, 1, 0 ) 666 666 pfv_ho(ji,jj,jl) = 0.5_wp * pv(ji,jj) * ( zpt(ji,jj,jl) + zpt(ji,jj+1,jl) ) 667 667 END_2D … … 672 672 ! 673 673 DO jl = 1, jpl !-- flux in y-direction 674 DO_2D _10_10674 DO_2D( 1, 0, 1, 0 ) 675 675 pfv_ho(ji,jj,jl) = 0.5_wp * pv(ji,jj) * ( pt(ji,jj,jl) + pt(ji,jj+1,jl) ) 676 676 END_2D … … 679 679 ! 680 680 DO jl = 1, jpl !-- first guess of tracer from v-flux 681 DO_2D _00_00681 DO_2D( 0, 0, 0, 0 ) 682 682 ztra = - ( pfv_ho(ji,jj,jl) - pfv_ho(ji,jj-1,jl) ) & 683 683 & + ( pv (ji,jj ) - pv (ji,jj-1 ) ) * pt(ji,jj,jl) * (1.-pamsk) … … 689 689 ! 690 690 DO jl = 1, jpl !-- flux in x-direction 691 DO_2D _10_10691 DO_2D( 1, 0, 1, 0 ) 692 692 pfu_ho(ji,jj,jl) = 0.5_wp * pu(ji,jj) * ( zpt(ji,jj,jl) + zpt(ji+1,jj,jl) ) 693 693 END_2D … … 737 737 ! !-- advective form update in zpt --! 738 738 DO jl = 1, jpl 739 DO_2D _00_00739 DO_2D( 0, 0, 0, 0 ) 740 740 zpt(ji,jj,jl) = ( pt(ji,jj,jl) - ( pubox(ji,jj ) * ( zt_u(ji,jj,jl) - zt_u(ji-1,jj,jl) ) * r1_e1t (ji,jj) & 741 741 & + pt (ji,jj,jl) * ( pu (ji,jj ) - pu (ji-1,jj ) ) * r1_e1e2t(ji,jj) & … … 764 764 ! !-- advective form update in zpt --! 765 765 DO jl = 1, jpl 766 DO_2D _00_00766 DO_2D( 0, 0, 0, 0 ) 767 767 zpt(ji,jj,jl) = ( pt(ji,jj,jl) - ( pvbox(ji,jj ) * ( zt_v(ji,jj,jl) - zt_v(ji,jj-1,jl) ) * r1_e2t (ji,jj) & 768 768 & + pt (ji,jj,jl) * ( pv (ji,jj ) - pv (ji,jj-1 ) ) * r1_e1e2t(ji,jj) & … … 846 846 ! 847 847 DO jl = 1, jpl 848 DO_2D _10_10848 DO_2D( 1, 0, 1, 0 ) 849 849 pt_u(ji,jj,jl) = 0.5_wp * umask(ji,jj,1) * ( pt(ji+1,jj,jl) + pt(ji,jj,jl) & 850 850 & - SIGN( 1._wp, pu(ji,jj) ) * ( pt(ji+1,jj,jl) - pt(ji,jj,jl) ) ) … … 855 855 ! 856 856 DO jl = 1, jpl 857 DO_2D _10_10857 DO_2D( 1, 0, 1, 0 ) 858 858 zcu = pu(ji,jj) * r1_e2u(ji,jj) * pdt * r1_e1u(ji,jj) 859 859 pt_u(ji,jj,jl) = 0.5_wp * umask(ji,jj,1) * ( pt(ji+1,jj,jl) + pt(ji,jj,jl) & … … 865 865 ! 866 866 DO jl = 1, jpl 867 DO_2D _10_10867 DO_2D( 1, 0, 1, 0 ) 868 868 zcu = pu(ji,jj) * r1_e2u(ji,jj) * pdt * r1_e1u(ji,jj) 869 869 zdx2 = e1u(ji,jj) * e1u(ji,jj) … … 879 879 ! 880 880 DO jl = 1, jpl 881 DO_2D _10_10881 DO_2D( 1, 0, 1, 0 ) 882 882 zcu = pu(ji,jj) * r1_e2u(ji,jj) * pdt * r1_e1u(ji,jj) 883 883 zdx2 = e1u(ji,jj) * e1u(ji,jj) … … 893 893 ! 894 894 DO jl = 1, jpl 895 DO_2D _10_10895 DO_2D( 1, 0, 1, 0 ) 896 896 zcu = pu(ji,jj) * r1_e2u(ji,jj) * pdt * r1_e1u(ji,jj) 897 897 zdx2 = e1u(ji,jj) * e1u(ji,jj) … … 914 914 IF( ll_neg ) THEN 915 915 DO jl = 1, jpl 916 DO_2D _10_10916 DO_2D( 1, 0, 1, 0 ) 917 917 IF( pt_u(ji,jj,jl) < 0._wp .OR. ( imsk_small(ji,jj,jl) == 0 .AND. pamsk == 0. ) ) THEN 918 918 pt_u(ji,jj,jl) = 0.5_wp * umask(ji,jj,1) * ( pt(ji+1,jj,jl) + pt(ji,jj,jl) & … … 924 924 ! !-- High order flux in i-direction --! 925 925 DO jl = 1, jpl 926 DO_2D _10_10926 DO_2D( 1, 0, 1, 0 ) 927 927 pfu_ho(ji,jj,jl) = pu(ji,jj) * pt_u(ji,jj,jl) 928 928 END_2D … … 957 957 ! !-- Laplacian in j-direction --! 958 958 DO jl = 1, jpl 959 DO_2D _10_00959 DO_2D( 1, 0, 0, 0 ) 960 960 ztv1(ji,jj,jl) = ( pt(ji,jj+1,jl) - pt(ji,jj,jl) ) * r1_e2v(ji,jj) * vmask(ji,jj,1) 961 961 END_2D 962 DO_2D _00_00962 DO_2D( 0, 0, 0, 0 ) 963 963 ztv2(ji,jj,jl) = ( ztv1(ji,jj,jl) - ztv1(ji,jj-1,jl) ) * r1_e2t(ji,jj) 964 964 END_2D … … 968 968 ! !-- BiLaplacian in j-direction --! 969 969 DO jl = 1, jpl 970 DO_2D _10_00970 DO_2D( 1, 0, 0, 0 ) 971 971 ztv3(ji,jj,jl) = ( ztv2(ji,jj+1,jl) - ztv2(ji,jj,jl) ) * r1_e2v(ji,jj) * vmask(ji,jj,1) 972 972 END_2D 973 DO_2D _00_00973 DO_2D( 0, 0, 0, 0 ) 974 974 ztv4(ji,jj,jl) = ( ztv3(ji,jj,jl) - ztv3(ji,jj-1,jl) ) * r1_e2t(ji,jj) 975 975 END_2D … … 982 982 CASE( 1 ) !== 1st order central TIM ==! (Eq. 21) 983 983 DO jl = 1, jpl 984 DO_2D _10_10984 DO_2D( 1, 0, 1, 0 ) 985 985 pt_v(ji,jj,jl) = 0.5_wp * vmask(ji,jj,1) * ( pt(ji,jj+1,jl) + pt(ji,jj,jl) & 986 986 & - SIGN( 1._wp, pv(ji,jj) ) * ( pt(ji,jj+1,jl) - pt(ji,jj,jl) ) ) … … 990 990 CASE( 2 ) !== 2nd order central TIM ==! (Eq. 23) 991 991 DO jl = 1, jpl 992 DO_2D _10_10992 DO_2D( 1, 0, 1, 0 ) 993 993 zcv = pv(ji,jj) * r1_e1v(ji,jj) * pdt * r1_e2v(ji,jj) 994 994 pt_v(ji,jj,jl) = 0.5_wp * vmask(ji,jj,1) * ( pt(ji,jj+1,jl) + pt(ji,jj,jl) & … … 999 999 CASE( 3 ) !== 3rd order central TIM ==! (Eq. 24) 1000 1000 DO jl = 1, jpl 1001 DO_2D _10_101001 DO_2D( 1, 0, 1, 0 ) 1002 1002 zcv = pv(ji,jj) * r1_e1v(ji,jj) * pdt * r1_e2v(ji,jj) 1003 1003 zdy2 = e2v(ji,jj) * e2v(ji,jj) … … 1012 1012 CASE( 4 ) !== 4th order central TIM ==! (Eq. 27) 1013 1013 DO jl = 1, jpl 1014 DO_2D _10_101014 DO_2D( 1, 0, 1, 0 ) 1015 1015 zcv = pv(ji,jj) * r1_e1v(ji,jj) * pdt * r1_e2v(ji,jj) 1016 1016 zdy2 = e2v(ji,jj) * e2v(ji,jj) … … 1025 1025 CASE( 5 ) !== 5th order central TIM ==! (Eq. 29) 1026 1026 DO jl = 1, jpl 1027 DO_2D _10_101027 DO_2D( 1, 0, 1, 0 ) 1028 1028 zcv = pv(ji,jj) * r1_e1v(ji,jj) * pdt * r1_e2v(ji,jj) 1029 1029 zdy2 = e2v(ji,jj) * e2v(ji,jj) … … 1046 1046 IF( ll_neg ) THEN 1047 1047 DO jl = 1, jpl 1048 DO_2D _10_101048 DO_2D( 1, 0, 1, 0 ) 1049 1049 IF( pt_v(ji,jj,jl) < 0._wp .OR. ( jmsk_small(ji,jj,jl) == 0 .AND. pamsk == 0. ) ) THEN 1050 1050 pt_v(ji,jj,jl) = 0.5_wp * vmask(ji,jj,1) * ( ( pt(ji,jj+1,jl) + pt(ji,jj,jl) ) & … … 1056 1056 ! !-- High order flux in j-direction --! 1057 1057 DO jl = 1, jpl 1058 DO_2D _10_101058 DO_2D( 1, 0, 1, 0 ) 1059 1059 pfv_ho(ji,jj,jl) = pv(ji,jj) * pt_v(ji,jj,jl) 1060 1060 END_2D … … 1092 1092 ! -------------------------------------------------- 1093 1093 DO jl = 1, jpl 1094 DO_2D _10_101094 DO_2D( 1, 0, 1, 0 ) 1095 1095 pfu_ho(ji,jj,jl) = pfu_ho(ji,jj,jl) - pfu_ups(ji,jj,jl) 1096 1096 pfv_ho(ji,jj,jl) = pfv_ho(ji,jj,jl) - pfv_ups(ji,jj,jl) … … 1109 1109 1110 1110 DO jl = 1, jpl 1111 DO_2D _00_001111 DO_2D( 0, 0, 0, 0 ) 1112 1112 zti_ups(ji,jj,jl)= pt_ups(ji+1,jj ,jl) 1113 1113 ztj_ups(ji,jj,jl)= pt_ups(ji ,jj+1,jl) … … 1117 1117 1118 1118 DO jl = 1, jpl 1119 DO_2D _00_001119 DO_2D( 0, 0, 0, 0 ) 1120 1120 IF ( pfu_ho(ji,jj,jl) * ( pt_ups(ji+1,jj ,jl) - pt_ups(ji,jj,jl) ) <= 0._wp .AND. & 1121 1121 & pfv_ho(ji,jj,jl) * ( pt_ups(ji ,jj+1,jl) - pt_ups(ji,jj,jl) ) <= 0._wp ) THEN … … 1146 1146 DO jl = 1, jpl 1147 1147 1148 DO_2D _11_111148 DO_2D( 1, 1, 1, 1 ) 1149 1149 IF ( pt(ji,jj,jl) <= 0._wp .AND. pt_ups(ji,jj,jl) <= 0._wp ) THEN 1150 1150 zbup(ji,jj) = -zbig … … 1162 1162 END_2D 1163 1163 1164 DO_2D _00_001164 DO_2D( 0, 0, 0, 0 ) 1165 1165 ! 1166 1166 zup = MAX( zbup(ji,jj), zbup(ji-1,jj), zbup(ji+1,jj), zbup(ji,jj-1), zbup(ji,jj+1) ) ! search max/min in neighbourhood … … 1199 1199 ! --------------------------------- 1200 1200 DO jl = 1, jpl 1201 DO_2D _10_101201 DO_2D( 1, 0, 1, 0 ) 1202 1202 zau = MIN( 1._wp , zbetdo(ji,jj,jl) , zbetup(ji+1,jj,jl) ) 1203 1203 zbu = MIN( 1._wp , zbetup(ji,jj,jl) , zbetdo(ji+1,jj,jl) ) … … 1210 1210 END_2D 1211 1211 1212 DO_2D _10_101212 DO_2D( 1, 0, 1, 0 ) 1213 1213 zav = MIN( 1._wp , zbetdo(ji,jj,jl) , zbetup(ji,jj+1,jl) ) 1214 1214 zbv = MIN( 1._wp , zbetup(ji,jj,jl) , zbetdo(ji,jj+1,jl) ) … … 1244 1244 ! 1245 1245 DO jl = 1, jpl 1246 DO_2D _00_001246 DO_2D( 0, 0, 0, 0 ) 1247 1247 zslpx(ji,jj,jl) = ( pt(ji+1,jj,jl) - pt(ji,jj,jl) ) * umask(ji,jj,1) 1248 1248 END_2D … … 1251 1251 1252 1252 DO jl = 1, jpl 1253 DO_2D _00_001253 DO_2D( 0, 0, 0, 0 ) 1254 1254 uCFL = pdt * ABS( pu(ji,jj) ) * r1_e1e2t(ji,jj) 1255 1255 … … 1335 1335 ! 1336 1336 DO jl = 1, jpl 1337 DO_2D _00_001337 DO_2D( 0, 0, 0, 0 ) 1338 1338 zslpy(ji,jj,jl) = ( pt(ji,jj+1,jl) - pt(ji,jj,jl) ) * vmask(ji,jj,1) 1339 1339 END_2D … … 1342 1342 1343 1343 DO jl = 1, jpl 1344 DO_2D _00_001344 DO_2D( 0, 0, 0, 0 ) 1345 1345 vCFL = pdt * ABS( pv(ji,jj) ) * r1_e1e2t(ji,jj) 1346 1346 … … 1436 1436 DO jl = 1, jpl 1437 1437 1438 DO_2D _11_111438 DO_2D( 1, 1, 1, 1 ) 1439 1439 IF ( pv_i(ji,jj,jl) > 0._wp ) THEN 1440 1440 ! … … 1502 1502 ! -- check snow load -- ! 1503 1503 DO jl = 1, jpl 1504 DO_2D _11_111504 DO_2D( 1, 1, 1, 1 ) 1505 1505 IF ( pv_i(ji,jj,jl) > 0._wp ) THEN 1506 1506 ! -
NEMO/trunk/src/ICE/icedyn_rdgrft.F90
r13226 r13295 161 161 npti = 0 ; nptidx(:) = 0 162 162 ipti = 0 ; iptidx(:) = 0 163 DO_2D _11_11163 DO_2D( 1, 1, 1, 1 ) 164 164 IF ( at_i(ji,jj) > epsi10 ) THEN 165 165 npti = npti + 1 … … 774 774 ! !--------------------------------------------------! 775 775 CASE( 1 ) !--- Spatial smoothing 776 DO_2D _00_00776 DO_2D( 0, 0, 0, 0 ) 777 777 IF ( SUM( a_i(ji,jj,:) ) > 0._wp ) THEN 778 778 zworka(ji,jj) = ( 4.0 * strength(ji,jj) & … … 785 785 END_2D 786 786 787 DO_2D _00_00787 DO_2D( 0, 0, 0, 0 ) 788 788 strength(ji,jj) = zworka(ji,jj) 789 789 END_2D … … 796 796 ENDIF 797 797 ! 798 DO_2D _00_00798 DO_2D( 0, 0, 0, 0 ) 799 799 IF ( SUM( a_i(ji,jj,:) ) > 0._wp ) THEN 800 800 itframe = 1 ! number of time steps for the running mean -
NEMO/trunk/src/ICE/icedyn_rhg_evp.F90
r13286 r13295 181 181 !------------------------------------------------------------------------------! 182 182 ! ocean/land mask 183 DO_2D _10_10183 DO_2D( 1, 0, 1, 0 ) 184 184 zfmask(ji,jj) = tmask(ji,jj,1) * tmask(ji+1,jj,1) * tmask(ji,jj+1,1) * tmask(ji+1,jj+1,1) 185 185 END_2D … … 188 188 ! Lateral boundary conditions on velocity (modify zfmask) 189 189 zwf(:,:) = zfmask(:,:) 190 DO_2D _00_00190 DO_2D( 0, 0, 0, 0 ) 191 191 IF( zfmask(ji,jj) == 0._wp ) THEN 192 192 zfmask(ji,jj) = rn_ishlat * MIN( 1._wp , MAX( zwf(ji+1,jj), zwf(ji,jj+1), zwf(ji-1,jj), zwf(ji,jj-1) ) ) … … 254 254 zsshdyn(:,:) = ice_var_sshdyn( ssh_m, snwice_mass, snwice_mass_b) 255 255 256 DO_2D _00_00256 DO_2D( 0, 0, 0, 0 ) 257 257 258 258 ! ice fraction at U-V points … … 305 305 ! 306 306 IF( ln_landfast_L16 ) THEN !-- Lemieux 2016 307 DO_2D _00_00307 DO_2D( 0, 0, 0, 0 ) 308 308 ! ice thickness at U-V points 309 309 zvU = 0.5_wp * ( vt_i(ji,jj) * e1e2t(ji,jj) + vt_i(ji+1,jj) * e1e2t(ji+1,jj) ) * r1_e1e2u(ji,jj) * umask(ji,jj,1) … … 322 322 ! 323 323 ELSE !-- no landfast 324 DO_2D _00_00324 DO_2D( 0, 0, 0, 0 ) 325 325 ztaux_base(ji,jj) = 0._wp 326 326 ztauy_base(ji,jj) = 0._wp … … 345 345 346 346 ! --- divergence, tension & shear (Appendix B of Hunke & Dukowicz, 2002) --- ! 347 DO_2D _10_10347 DO_2D( 1, 0, 1, 0 ) 348 348 349 349 ! shear at F points … … 355 355 CALL lbc_lnk( 'icedyn_rhg_evp', zds, 'F', 1.0_wp ) 356 356 357 DO_2D _01_01357 DO_2D( 0, 1, 0, 1 ) 358 358 359 359 ! shear**2 at T points (doc eq. A16) … … 397 397 CALL lbc_lnk( 'icedyn_rhg_evp', zp_delt, 'T', 1.0_wp ) 398 398 399 DO_2D _10_10399 DO_2D( 1, 0, 1, 0 ) 400 400 401 401 ! alpha & beta for aEVP … … 415 415 416 416 ! --- Ice internal stresses (Appendix C of Hunke and Dukowicz, 2002) --- ! 417 DO_2D _00_00417 DO_2D( 0, 0, 0, 0 ) 418 418 ! !--- U points 419 419 zfU(ji,jj) = 0.5_wp * ( ( zs1(ji+1,jj) - zs1(ji,jj) ) * e2u(ji,jj) & … … 443 443 IF( MOD(jter,2) == 0 ) THEN ! even iterations 444 444 ! 445 DO_2D _00_00445 DO_2D( 0, 0, 0, 0 ) 446 446 ! !--- tau_io/(v_oce - v_ice) 447 447 zTauO = zaV(ji,jj) * zrhoco * SQRT( ( v_ice (ji,jj) - v_oce (ji,jj) ) * ( v_ice (ji,jj) - v_oce (ji,jj) ) & … … 492 492 IF( ln_bdy ) CALL bdy_ice_dyn( 'V' ) 493 493 ! 494 DO_2D _00_00494 DO_2D( 0, 0, 0, 0 ) 495 495 ! !--- tau_io/(u_oce - u_ice) 496 496 zTauO = zaU(ji,jj) * zrhoco * SQRT( ( u_ice (ji,jj) - u_oce (ji,jj) ) * ( u_ice (ji,jj) - u_oce (ji,jj) ) & … … 543 543 ELSE ! odd iterations 544 544 ! 545 DO_2D _00_00545 DO_2D( 0, 0, 0, 0 ) 546 546 ! !--- tau_io/(u_oce - u_ice) 547 547 zTauO = zaU(ji,jj) * zrhoco * SQRT( ( u_ice (ji,jj) - u_oce (ji,jj) ) * ( u_ice (ji,jj) - u_oce (ji,jj) ) & … … 592 592 IF( ln_bdy ) CALL bdy_ice_dyn( 'U' ) 593 593 ! 594 DO_2D _00_00594 DO_2D( 0, 0, 0, 0 ) 595 595 ! !--- tau_io/(v_oce - v_ice) 596 596 zTauO = zaV(ji,jj) * zrhoco * SQRT( ( v_ice (ji,jj) - v_oce (ji,jj) ) * ( v_ice (ji,jj) - v_oce (ji,jj) ) & … … 658 658 ! 4) Recompute delta, shear and div (inputs for mechanical redistribution) 659 659 !------------------------------------------------------------------------------! 660 DO_2D _10_10660 DO_2D( 1, 0, 1, 0 ) 661 661 662 662 ! shear at F points … … 667 667 END_2D 668 668 669 DO_2D _00_00669 DO_2D( 0, 0, 0, 0 ) 670 670 671 671 ! tension**2 at T points … … 706 706 ! 5) diagnostics 707 707 !------------------------------------------------------------------------------! 708 DO_2D _11_11708 DO_2D( 1, 1, 1, 1 ) 709 709 zmsk00(ji,jj) = MAX( 0._wp , SIGN( 1._wp , at_i(ji,jj) - epsi06 ) ) ! 1 if ice, 0 if no ice 710 710 END_2D … … 735 735 ALLOCATE( zsig1(jpi,jpj) , zsig2(jpi,jpj) , zsig3(jpi,jpj) ) 736 736 ! 737 DO_2D _00_00737 DO_2D( 0, 0, 0, 0 ) 738 738 zdum1 = ( zmsk00(ji-1,jj) * pstress12_i(ji-1,jj) + zmsk00(ji ,jj-1) * pstress12_i(ji ,jj-1) + & ! stress12_i at T-point 739 739 & zmsk00(ji ,jj) * pstress12_i(ji ,jj) + zmsk00(ji-1,jj-1) * pstress12_i(ji-1,jj-1) ) & … … 786 786 & zdiag_xmtrp_snw(jpi,jpj) , zdiag_ymtrp_snw(jpi,jpj) , zdiag_xatrp(jpi,jpj) , zdiag_yatrp(jpi,jpj) ) 787 787 ! 788 DO_2D _00_00788 DO_2D( 0, 0, 0, 0 ) 789 789 ! 2D ice mass, snow mass, area transport arrays (X, Y) 790 790 zfac_x = 0.5 * u_ice(ji,jj) * e2u(ji,jj) * zmsk00(ji,jj) -
NEMO/trunk/src/ICE/iceistate.F90
r13237 r13295 288 288 ! select ice covered grid points 289 289 npti = 0 ; nptidx(:) = 0 290 DO_2D _11_11290 DO_2D( 1, 1, 1, 1 ) 291 291 IF ( zht_i_ini(ji,jj) > 0._wp ) THEN 292 292 npti = npti + 1 … … 338 338 CALL ice_var_salprof ! for sz_i 339 339 DO jl = 1, jpl 340 DO_2D _11_11340 DO_2D( 1, 1, 1, 1 ) 341 341 v_i (ji,jj,jl) = h_i(ji,jj,jl) * a_i(ji,jj,jl) 342 342 v_s (ji,jj,jl) = h_s(ji,jj,jl) * a_i(ji,jj,jl) … … 346 346 ! 347 347 DO jl = 1, jpl 348 DO_3D _11_11(1, nlay_s )348 DO_3D( 1, 1, 1, 1, 1, nlay_s ) 349 349 t_s(ji,jj,jk,jl) = zts_3d(ji,jj,jl) 350 350 e_s(ji,jj,jk,jl) = zswitch(ji,jj) * v_s(ji,jj,jl) * r1_nlay_s * & … … 354 354 ! 355 355 DO jl = 1, jpl 356 DO_3D _11_11(1, nlay_i )356 DO_3D( 1, 1, 1, 1, 1, nlay_i ) 357 357 t_i (ji,jj,jk,jl) = zti_3d(ji,jj,jl) 358 358 ztmelts = - rTmlt * sz_i(ji,jj,jk,jl) + rt0 ! melting temperature in K -
NEMO/trunk/src/ICE/iceitd.F90
r13226 r13295 98 98 ! 99 99 npti = 0 ; nptidx(:) = 0 100 DO_2D _11_11100 DO_2D( 1, 1, 1, 1 ) 101 101 IF ( at_i(ji,jj) > epsi10 ) THEN 102 102 npti = npti + 1 … … 611 611 ! !--------------------------------------- 612 612 npti = 0 ; nptidx(:) = 0 613 DO_2D _11_11613 DO_2D( 1, 1, 1, 1 ) 614 614 IF( a_i(ji,jj,jl) > 0._wp .AND. v_i(ji,jj,jl) > (a_i(ji,jj,jl) * hi_max(jl)) ) THEN 615 615 npti = npti + 1 … … 650 650 ! !----------------------------------------- 651 651 npti = 0 ; nptidx(:) = 0 652 DO_2D _11_11652 DO_2D( 1, 1, 1, 1 ) 653 653 IF( a_i(ji,jj,jl+1) > 0._wp .AND. v_i(ji,jj,jl+1) <= (a_i(ji,jj,jl+1) * hi_max(jl)) ) THEN 654 654 npti = npti + 1 -
NEMO/trunk/src/ICE/icesbc.F90
r13226 r13295 82 82 IF( ln_mixcpl) THEN ! Case of a mixed Bulk/Coupled formulation 83 83 CALL sbc_cpl_ice_tau( zutau_ice , zvtau_ice ) 84 DO_2D _00_0084 DO_2D( 0, 0, 0, 0 ) 85 85 utau_ice(ji,jj) = utau_ice(ji,jj) * xcplmask(ji,jj,0) + zutau_ice(ji,jj) * ( 1. - xcplmask(ji,jj,0) ) 86 86 vtau_ice(ji,jj) = vtau_ice(ji,jj) * xcplmask(ji,jj,0) + zvtau_ice(ji,jj) * ( 1. - xcplmask(ji,jj,0) ) -
NEMO/trunk/src/ICE/icethd.F90
r13226 r13295 109 109 zu_io(:,:) = u_ice(:,:) - ssu_m(:,:) 110 110 zv_io(:,:) = v_ice(:,:) - ssv_m(:,:) 111 DO_2D _00_00111 DO_2D( 0, 0, 0, 0 ) 112 112 zfric(ji,jj) = rn_cio * ( 0.5_wp * & 113 113 & ( zu_io(ji,jj) * zu_io(ji,jj) + zu_io(ji-1,jj) * zu_io(ji-1,jj) & … … 115 115 END_2D 116 116 ELSE ! if no ice dynamics => transmit directly the atmospheric stress to the ocean 117 DO_2D _00_00117 DO_2D( 0, 0, 0, 0 ) 118 118 zfric(ji,jj) = r1_rho0 * SQRT( 0.5_wp * & 119 119 & ( utau(ji,jj) * utau(ji,jj) + utau(ji-1,jj) * utau(ji-1,jj) & … … 126 126 ! Partial computation of forcing for the thermodynamic sea ice model 127 127 !--------------------------------------------------------------------! 128 DO_2D _11_11128 DO_2D( 1, 1, 1, 1 ) 129 129 rswitch = tmask(ji,jj,1) * MAX( 0._wp , SIGN( 1._wp , at_i(ji,jj) - epsi10 ) ) ! 0 if no ice 130 130 ! … … 196 196 ! select ice covered grid points 197 197 npti = 0 ; nptidx(:) = 0 198 DO_2D _11_11198 DO_2D( 1, 1, 1, 1 ) 199 199 IF ( a_i(ji,jj,jl) > epsi10 ) THEN 200 200 npti = npti + 1 -
NEMO/trunk/src/ICE/icethd_do.F90
r13226 r13295 145 145 zgamafr = 0.03 146 146 ! 147 DO_2D _00_00147 DO_2D( 0, 0, 0, 0 ) 148 148 IF ( qlead(ji,jj) < 0._wp .AND. tau_icebfr(ji,jj) == 0._wp ) THEN ! activated if cooling and no landfast 149 149 ! -- Wind stress -- ! … … 202 202 ! Identify grid points where new ice forms 203 203 npti = 0 ; nptidx(:) = 0 204 DO_2D _11_11204 DO_2D( 1, 1, 1, 1 ) 205 205 IF ( qlead(ji,jj) < 0._wp .AND. tau_icebfr(ji,jj) == 0._wp ) THEN 206 206 npti = npti + 1 -
NEMO/trunk/src/ICE/iceupdate.F90
r13286 r13295 113 113 ENDIF 114 114 115 DO_2D _11_11115 DO_2D( 1, 1, 1, 1 ) 116 116 117 117 ! Solar heat flux reaching the ocean = zqsr (W.m-2) … … 332 332 ! 333 333 IF( MOD( kt-1, nn_fsbc ) == 0 ) THEN !== Ice time-step only ==! (i.e. surface module time-step) 334 DO_2D _00_00334 DO_2D( 0, 0, 0, 0 ) 335 335 ! ! 2*(U_ice-U_oce) at T-point 336 336 zu_t = u_ice(ji,jj) + u_ice(ji-1,jj) - u_oce(ji,jj) - u_oce(ji-1,jj) … … 351 351 ! !== every ocean time-step ==! 352 352 ! 353 DO_2D _00_00353 DO_2D( 0, 0, 0, 0 ) 354 354 ! ice area at u and v-points 355 355 zat_u = ( at_i(ji,jj) * tmask(ji,jj,1) + at_i (ji+1,jj ) * tmask(ji+1,jj ,1) ) & -
NEMO/trunk/src/ICE/icevar.F90
r13226 r13295 243 243 zlay_i = REAL( nlay_i , wp ) ! number of layers 244 244 DO jl = 1, jpl 245 DO_3D _11_11(1, nlay_i )245 DO_3D( 1, 1, 1, 1, 1, nlay_i ) 246 246 IF ( v_i(ji,jj,jl) > epsi20 ) THEN !--- icy area 247 247 ! … … 347 347 z1_dS = 1._wp / ( zsi1 - zsi0 ) 348 348 DO jl = 1, jpl 349 DO_2D _11_11349 DO_2D( 1, 1, 1, 1 ) 350 350 zalpha(ji,jj,jl) = MAX( 0._wp , MIN( ( zsi1 - s_i(ji,jj,jl) ) * z1_dS , 1._wp ) ) 351 351 ! ! force a constant profile when SSS too low (Baltic Sea) … … 356 356 ! Computation of the profile 357 357 DO jl = 1, jpl 358 DO_3D _11_11(1, nlay_i )358 DO_3D( 1, 1, 1, 1, 1, nlay_i ) 359 359 ! ! linear profile with 0 surface value 360 360 zs0 = z_slope_s(ji,jj,jl) * ( REAL(jk,wp) - 0.5_wp ) * h_i(ji,jj,jl) * r1_nlay_i … … 486 486 ! Zap ice energy and use ocean heat to melt ice 487 487 !----------------------------------------------------------------- 488 DO_3D _11_11(1, nlay_i )488 DO_3D( 1, 1, 1, 1, 1, nlay_i ) 489 489 ! update exchanges with ocean 490 490 hfx_res(ji,jj) = hfx_res(ji,jj) - (1._wp - zswitch(ji,jj) ) * e_i(ji,jj,jk,jl) * r1_Dt_ice ! W.m-2 <0 … … 493 493 END_3D 494 494 ! 495 DO_3D _11_11(1, nlay_s )495 DO_3D( 1, 1, 1, 1, 1, nlay_s ) 496 496 ! update exchanges with ocean 497 497 hfx_res(ji,jj) = hfx_res(ji,jj) - (1._wp - zswitch(ji,jj) ) * e_s(ji,jj,jk,jl) * r1_Dt_ice ! W.m-2 <0 … … 503 503 ! zap ice and snow volume, add water and salt to ocean 504 504 !----------------------------------------------------------------- 505 DO_2D _11_11505 DO_2D( 1, 1, 1, 1 ) 506 506 ! update exchanges with ocean 507 507 sfx_res(ji,jj) = sfx_res(ji,jj) + (1._wp - zswitch(ji,jj) ) * sv_i(ji,jj,jl) * rhoi * r1_Dt_ice … … 574 574 ! zap ice energy and send it to the ocean 575 575 !---------------------------------------- 576 DO_3D _11_11(1, nlay_i )576 DO_3D( 1, 1, 1, 1, 1, nlay_i ) 577 577 IF( pe_i(ji,jj,jk,jl) < 0._wp .OR. pa_i(ji,jj,jl) <= 0._wp ) THEN 578 578 hfx_res(ji,jj) = hfx_res(ji,jj) - pe_i(ji,jj,jk,jl) * z1_dt ! W.m-2 >0 … … 581 581 END_3D 582 582 ! 583 DO_3D _11_11(1, nlay_s )583 DO_3D( 1, 1, 1, 1, 1, nlay_s ) 584 584 IF( pe_s(ji,jj,jk,jl) < 0._wp .OR. pa_i(ji,jj,jl) <= 0._wp ) THEN 585 585 hfx_res(ji,jj) = hfx_res(ji,jj) - pe_s(ji,jj,jk,jl) * z1_dt ! W.m-2 <0 … … 591 591 ! zap ice and snow volume, add water and salt to ocean 592 592 !----------------------------------------------------- 593 DO_2D _11_11593 DO_2D( 1, 1, 1, 1 ) 594 594 IF( pv_i(ji,jj,jl) < 0._wp .OR. pa_i(ji,jj,jl) <= 0._wp ) THEN 595 595 wfx_res(ji,jj) = wfx_res(ji,jj) + pv_i (ji,jj,jl) * rhoi * z1_dt -
NEMO/trunk/src/ICE/icewri.F90
r13226 r13295 71 71 72 72 ! tresholds for outputs 73 DO_2D _11_1173 DO_2D( 1, 1, 1, 1 ) 74 74 zmsk00(ji,jj) = MAX( 0._wp , SIGN( 1._wp , at_i(ji,jj) - epsi06 ) ) ! 1 if ice , 0 if no ice 75 75 zmsk05(ji,jj) = MAX( 0._wp , SIGN( 1._wp , at_i(ji,jj) - 0.05_wp ) ) ! 1 if 5% ice , 0 if less … … 78 78 END_2D 79 79 DO jl = 1, jpl 80 DO_2D _11_1180 DO_2D( 1, 1, 1, 1 ) 81 81 zmsk00l(ji,jj,jl) = MAX( 0._wp , SIGN( 1._wp , a_i(ji,jj,jl) - epsi06 ) ) 82 82 zmsksnl(ji,jj,jl) = MAX( 0._wp , SIGN( 1._wp , v_s(ji,jj,jl) - epsi06 ) ) … … 130 130 ! 131 131 IF( iom_use('icevel') .OR. iom_use('fasticepres') ) THEN ! module of ice velocity 132 DO_2D _00_00132 DO_2D( 0, 0, 0, 0 ) 133 133 z2da = u_ice(ji,jj) + u_ice(ji-1,jj) 134 134 z2db = v_ice(ji,jj) + v_ice(ji,jj-1)
Note: See TracChangeset
for help on using the changeset viewer.
