Changeset 15292
- Timestamp:
- 2021-09-27T17:35:36+02:00 (21 months ago)
- File:
-
- 1 edited
-
NEMO/trunk/src/ICE/icedyn_rhg_vp.F90 (modified) (15 diffs)
Legend:
- Unmodified
- Added
- Removed
-
NEMO/trunk/src/ICE/icedyn_rhg_vp.F90
r15014 r15292 218 218 219 219 !!---------------------------------------------------------------------------------------------------------------------- 220 ! DEBUG put all forcing terms to zero221 ! air-ice drag222 utau_ice(:,:) = 0._wp223 vtau_ice(:,:) = 0._wp224 ! coriolis225 ff_t(:,:) = 0._wp226 ! ice-ocean drag227 rn_cio = 0._wp228 ! ssh229 ! done line 330 !!! dont forget to act there230 ! END DEBUG220 !!$ ! DEBUG put all forcing terms to zero 221 !!$ ! air-ice drag 222 !!$ utau_ice(:,:) = 0._wp 223 !!$ vtau_ice(:,:) = 0._wp 224 !!$ ! coriolis 225 !!$ ff_t(:,:) = 0._wp 226 !!$ ! ice-ocean drag 227 !!$ rn_cio = 0._wp 228 !!$ ! ssh 229 !!$ ! done line 330 !!! dont forget to act there 230 !!$ ! END DEBUG 231 231 232 232 IF( kt == nit000 .AND. lwp ) WRITE(numout,*) '-- ice_dyn_rhg_vp: VP sea-ice rheology (LSR solver)' … … 240 240 241 241 ! for diagnostics and convergence tests 242 DO_2D( 1, 1, 1, 1)242 DO_2D( nn_hls, nn_hls, nn_hls, nn_hls ) 243 243 zmsk00(ji,jj) = MAX( 0._wp , SIGN( 1._wp , at_i(ji,jj) - epsi06 ) ) ! 1 if ice , 0 if no ice 244 244 END_2D 245 245 IF( nn_rhg_chkcvg > 0 ) THEN 246 DO_2D( 1, 1, 1, 1)246 DO_2D( nn_hls, nn_hls, nn_hls, nn_hls ) 247 247 zmsk15(ji,jj) = MAX( 0._wp , SIGN( 1._wp , at_i(ji,jj) - 0.15_wp ) ) ! 1 if 15% ice, 0 if less 248 248 END_2D … … 329 329 ! non-embedded sea ice: use ocean surface for slope calculation 330 330 zsshdyn(:,:) = ice_var_sshdyn( ssh_m, snwice_mass, snwice_mass_b) 331 zsshdyn(:,:) = 0._wp ! DEBUG CAREFUL !!!332 333 zmt(:,:) = rhos * vt_s(:,:) + rhoi * vt_i(:,:) ! Snow and ice mass at T-point334 zmf(:,:) = zmt(:,:) * ff_t(:,:) ! Coriolis factor at T points (m*f)335 336 DO jj = 2, jpj - 1337 DO ji = 2, jpi - 1338 339 ! Ice fraction at U-V points 340 za_iU(ji,jj) = 0.5_wp * ( at_i(ji,jj) * e1e2t(ji,jj) + at_i(ji+1,jj) * e1e2t(ji+1,jj) ) * r1_e1e2u(ji,jj) * umask(ji,jj,1)341 za_iV(ji,jj) = 0.5_wp * ( at_i(ji,jj) * e1e2t(ji,jj) + at_i(ji,jj+1) * e1e2t(ji,jj+1) ) * r1_e1e2v(ji,jj) * vmask(ji,jj,1)342 343 ! Snow and ice mass at U-V points344 zm1 = zmt(ji,jj)345 zm2 = zmt(ji+1,jj)346 zm3 = zmt(ji,jj+1)347 348 349 350 351 352 353 354 355 356 357 358 359 360 361 362 363 364 365 366 367 368 369 370 371 372 373 374 375 376 377 378 379 331 !!$ zsshdyn(:,:) = 0._wp ! DEBUG CAREFUL !!! 332 333 DO_2D( nn_hls, nn_hls, nn_hls, nn_hls ) 334 zm1 = ( rhos * vt_s(ji,jj) + rhoi * vt_i(ji,jj) ) ! Ice/snow mass at U-V points 335 zmf (ji,jj) = zm1 * ff_t(ji,jj) ! Coriolis at T points (m*f) 336 END_2D 337 338 DO_2D( nn_hls-1, nn_hls-1, nn_hls-1, nn_hls-1 ) 339 340 ! Ice fraction at U-V points 341 za_iU(ji,jj) = 0.5_wp * ( at_i(ji,jj) * e1e2t(ji,jj) + at_i(ji+1,jj) * e1e2t(ji+1,jj) ) * r1_e1e2u(ji,jj) * umask(ji,jj,1) 342 za_iV(ji,jj) = 0.5_wp * ( at_i(ji,jj) * e1e2t(ji,jj) + at_i(ji,jj+1) * e1e2t(ji,jj+1) ) * r1_e1e2v(ji,jj) * vmask(ji,jj,1) 343 344 ! Ice/snow mass at U-V points 345 zm1 = ( rhos * vt_s(ji ,jj ) + rhoi * vt_i(ji ,jj ) ) 346 zm2 = ( rhos * vt_s(ji+1,jj ) + rhoi * vt_i(ji+1,jj ) ) 347 zm3 = ( rhos * vt_s(ji ,jj+1) + rhoi * vt_i(ji ,jj+1) ) 348 zmassU = 0.5_wp * ( zm1 * e1e2t(ji,jj) + zm2 * e1e2t(ji+1,jj) ) * r1_e1e2u(ji,jj) * umask(ji,jj,1) 349 zmassV = 0.5_wp * ( zm1 * e1e2t(ji,jj) + zm3 * e1e2t(ji,jj+1) ) * r1_e1e2v(ji,jj) * vmask(ji,jj,1) 350 351 ! Mass per unit area divided by time step 352 zmassU_t(ji,jj) = zmassU * r1_Dt_ice 353 zmassV_t(ji,jj) = zmassV * r1_Dt_ice 354 355 ! Acceleration term contribution to RHS (depends on velocity at previous time step) 356 zmU_t(ji,jj) = zmassU_t(ji,jj) * u_ice(ji,jj) 357 zmV_t(ji,jj) = zmassV_t(ji,jj) * v_ice(ji,jj) 358 359 ! Ocean currents at U-V points 360 v_oceU(ji,jj) = 0.25_wp * ( v_oce(ji,jj) + v_oce(ji,jj-1) + v_oce(ji+1,jj) + v_oce(ji+1,jj-1) ) * umask(ji,jj,1) 361 u_oceV(ji,jj) = 0.25_wp * ( u_oce(ji,jj) + u_oce(ji-1,jj) + u_oce(ji,jj+1) + u_oce(ji-1,jj+1) ) * vmask(ji,jj,1) 362 363 ! Wind stress 364 ztaux_ai(ji,jj) = za_iU(ji,jj) * utau_ice(ji,jj) 365 ztauy_ai(ji,jj) = za_iV(ji,jj) * vtau_ice(ji,jj) 366 367 ! Force due to sea surface tilt(- m*g*GRAD(ssh)) 368 zspgU(ji,jj) = - zmassU * grav * ( zsshdyn(ji+1,jj) - zsshdyn(ji,jj) ) * r1_e1u(ji,jj) 369 zspgV(ji,jj) = - zmassV * grav * ( zsshdyn(ji,jj+1) - zsshdyn(ji,jj) ) * r1_e2v(ji,jj) 370 371 ! Mask for ice presence (1) or absence (0) 372 zmsk00x(ji,jj) = 1._wp - MAX( 0._wp, SIGN( 1._wp, -zmassU ) ) ! 0 if no ice 373 zmsk00y(ji,jj) = 1._wp - MAX( 0._wp, SIGN( 1._wp, -zmassV ) ) ! 0 if no ice 374 375 ! Mask for lots of ice (1) or little ice (0) 376 IF ( zmassU <= zmmin .AND. za_iU(ji,jj) <= zamin ) THEN ; zmsk01x(ji,jj) = 0._wp 377 ELSE ; zmsk01x(ji,jj) = 1._wp ; ENDIF 378 IF ( zmassV <= zmmin .AND. za_iV(ji,jj) <= zamin ) THEN ; zmsk01y(ji,jj) = 0._wp 379 ELSE ; zmsk01y(ji,jj) = 1._wp ; ENDIF 380 380 381 381 ! MV TEST DEBUG 382 IF ( ( zmt(ji,jj) <= zmmin .OR. zmt(ji+1,jj) <= zmmin ) .AND. &383 & ( at_i(ji,jj) <= zamin .OR. at_i(ji+1,jj) <= zamin ) ) THEN ; zmsk01x(ji,jj) = 0._wp384 ELSE ; zmsk01x(ji,jj) = 1._wp ; ENDIF385 386 IF ( ( zmt(ji,jj) <= zmmin .OR. zmt(ji,jj+1) <= zmmin ) .AND. &387 & ( at_i(ji,jj) <= zamin .OR. at_i(ji,jj+1) <= zamin ) ) THEN ; zmsk01y(ji,jj) = 0._wp388 ELSE ; zmsk01y(ji,jj) = 1._wp ; ENDIF382 !!$ IF ( ( zmt(ji,jj) <= zmmin .OR. zmt(ji+1,jj) <= zmmin ) .AND. & 383 !!$ & ( at_i(ji,jj) <= zamin .OR. at_i(ji+1,jj) <= zamin ) ) THEN ; zmsk01x(ji,jj) = 0._wp 384 !!$ ELSE ; zmsk01x(ji,jj) = 1._wp ; ENDIF 385 !!$ 386 !!$ IF ( ( zmt(ji,jj) <= zmmin .OR. zmt(ji,jj+1) <= zmmin ) .AND. & 387 !!$ & ( at_i(ji,jj) <= zamin .OR. at_i(ji,jj+1) <= zamin ) ) THEN ; zmsk01y(ji,jj) = 0._wp 388 !!$ ELSE ; zmsk01y(ji,jj) = 1._wp ; ENDIF 389 389 ! END MV TEST DEBUG 390 390 391 END DO 392 END DO 393 394 CALL iom_put( 'zmsk00x' , zmsk00x ) ! MV DEBUG 395 CALL iom_put( 'zmsk00y' , zmsk00y ) ! MV DEBUG 396 CALL iom_put( 'zmsk01x' , zmsk01x ) ! MV DEBUG 397 CALL iom_put( 'zmsk01y' , zmsk01y ) ! MV DEBUG 398 CALL iom_put( 'ztaux_ai' , ztaux_ai ) ! MV DEBUG 399 CALL iom_put( 'ztauy_ai' , ztauy_ai ) ! MV DEBUG 400 CALL iom_put( 'zspgU' , zspgU ) ! MV DEBUG 401 CALL iom_put( 'zspgV' , zspgV ) ! MV DEBUG 402 391 END_2D 392 393 !!$ CALL iom_put( 'zmsk00x' , zmsk00x ) ! MV DEBUG 394 !!$ CALL iom_put( 'zmsk00y' , zmsk00y ) ! MV DEBUG 395 !!$ CALL iom_put( 'zmsk01x' , zmsk01x ) ! MV DEBUG 396 !!$ CALL iom_put( 'zmsk01y' , zmsk01y ) ! MV DEBUG 397 !!$ CALL iom_put( 'ztaux_ai' , ztaux_ai ) ! MV DEBUG 398 !!$ CALL iom_put( 'ztauy_ai' , ztauy_ai ) ! MV DEBUG 399 !!$ CALL iom_put( 'zspgU' , zspgU ) ! MV DEBUG 400 !!$ CALL iom_put( 'zspgV' , zspgV ) ! MV DEBUG 401 !!$ 403 402 !------------------------------------------------------------------------------! 404 403 ! … … 437 436 438 437 ! --- divergence, tension & shear (Appendix B of Hunke & Dukowicz, 2002) --- ! 439 DO jj = 1, jpj - 1 ! loops start at 1 since there is no boundary condition (lbc_lnk) at i=1 and j=1 for F points 440 DO ji = 1, jpi - 1 441 442 ! shear at F points 443 zds(ji,jj) = ( ( zu_c(ji,jj+1) * r1_e1u(ji,jj+1) - zu_c(ji,jj) * r1_e1u(ji,jj) ) * e1f(ji,jj) * e1f(ji,jj) & 444 & + ( zv_c(ji+1,jj) * r1_e2v(ji+1,jj) - zv_c(ji,jj) * r1_e2v(ji,jj) ) * e2f(ji,jj) * e2f(ji,jj) & 445 & ) * r1_e1e2f(ji,jj) * fimask(ji,jj) 446 447 END DO 448 END DO 438 DO_2D( nn_hls, nn_hls-1, nn_hls, nn_hls-1 ) 439 440 ! shear at F points 441 zds(ji,jj) = ( ( zu_c(ji,jj+1) * r1_e1u(ji,jj+1) - zu_c(ji,jj) * r1_e1u(ji,jj) ) * e1f(ji,jj) * e1f(ji,jj) & 442 & + ( zv_c(ji+1,jj) * r1_e2v(ji+1,jj) - zv_c(ji,jj) * r1_e2v(ji,jj) ) * e2f(ji,jj) * e2f(ji,jj) & 443 & ) * r1_e1e2f(ji,jj) * fimask(ji,jj) 444 445 END_2D 449 446 450 447 CALL lbc_lnk( 'icedyn_rhg_vp', zds, 'F', 1. ) ! MV TEST could be un-necessary according to Gurvan 451 CALL iom_put( 'zds' , zds ) ! MV DEBUG448 !!$ CALL iom_put( 'zds' , zds ) ! MV DEBUG 452 449 453 450 IF( lwp ) WRITE(numout,*) ' outer loop 1a i_out : ', i_out … … 496 493 CALL lbc_lnk( 'icedyn_rhg_vp', zp_deltastar_t , 'T', 1. , zzt , 'T', 1., zet, 'T', 1. ) 497 494 498 CALL iom_put( 'zzt' , zzt ) ! MV DEBUG499 CALL iom_put( 'zet' , zet ) ! MV DEBUG500 CALL iom_put( 'zp_deltastar_t', zp_deltastar_t ) ! MV DEBUG495 !!$ CALL iom_put( 'zzt' , zzt ) ! MV DEBUG 496 !!$ CALL iom_put( 'zet' , zet ) ! MV DEBUG 497 !!$ CALL iom_put( 'zp_deltastar_t', zp_deltastar_t ) ! MV DEBUG 501 498 502 499 IF( lwp ) WRITE(numout,*) ' outer loop 1b i_out : ', i_out … … 515 512 516 513 CALL lbc_lnk( 'icedyn_rhg_vp', zef, 'F', 1. ) 517 CALL iom_put( 'zef' , zef ) ! MV DEBUG514 !!$ CALL iom_put( 'zef' , zef ) ! MV DEBUG 518 515 IF( lwp ) WRITE(numout,*) ' outer loop 1c i_out : ', i_out 519 516 … … 558 555 CALL lbc_lnk( 'icedyn_rhg_vp', zCorU, 'U', -1., zCorV, 'V', -1. ) 559 556 560 CALL iom_put( 'zCwU' , zCwU ) ! MV DEBUG561 CALL iom_put( 'zCwV' , zCwV ) ! MV DEBUG562 CALL iom_put( 'zCorU' , zCorU ) ! MV DEBUG563 CALL iom_put( 'zCorV' , zCorV ) ! MV DEBUG557 !!$ CALL iom_put( 'zCwU' , zCwU ) ! MV DEBUG 558 !!$ CALL iom_put( 'zCwV' , zCwV ) ! MV DEBUG 559 !!$ CALL iom_put( 'zCorU' , zCorU ) ! MV DEBUG 560 !!$ CALL iom_put( 'zCorV' , zCorV ) ! MV DEBUG 564 561 565 562 IF( lwp ) WRITE(numout,*) ' outer loop 1f i_out : ', i_out … … 591 588 END DO 592 589 593 CALL iom_put( 'zs1_rhsu' , zs1_rhsu ) ! MV DEBUG594 CALL iom_put( 'zs2_rhsu' , zs2_rhsu ) ! MV DEBUG595 CALL iom_put( 'zs1_rhsv' , zs1_rhsv ) ! MV DEBUG596 CALL iom_put( 'zs2_rhsv' , zs2_rhsv ) ! MV DEBUG590 !!$ CALL iom_put( 'zs1_rhsu' , zs1_rhsu ) ! MV DEBUG 591 !!$ CALL iom_put( 'zs2_rhsu' , zs2_rhsu ) ! MV DEBUG 592 !!$ CALL iom_put( 'zs1_rhsv' , zs1_rhsv ) ! MV DEBUG 593 !!$ CALL iom_put( 'zs2_rhsv' , zs2_rhsv ) ! MV DEBUG 597 594 598 595 ! a priori, no lbclnk, because rhsu is only used in the inner domain … … 619 616 CALL lbc_lnk( 'icedyn_rhg_vp', zs12_rhsv, 'F', 1. ) 620 617 621 CALL iom_put( 'zs12_rhsu' , zs12_rhsu ) ! MV DEBUG622 CALL iom_put( 'zs12_rhsv' , zs12_rhsv ) ! MV DEBUG618 !!$ CALL iom_put( 'zs12_rhsu' , zs12_rhsu ) ! MV DEBUG 619 !!$ CALL iom_put( 'zs12_rhsv' , zs12_rhsv ) ! MV DEBUG 623 620 624 621 ! a priori, no lbclnk, because rhsu are only used in the inner domain … … 644 641 END DO 645 642 646 CALL iom_put( 'zf_rhsu' , zf_rhsu ) ! MV DEBUG647 CALL iom_put( 'zf_rhsv' , zf_rhsv ) ! MV DEBUG643 !!$ CALL iom_put( 'zf_rhsu' , zf_rhsu ) ! MV DEBUG 644 !!$ CALL iom_put( 'zf_rhsv' , zf_rhsv ) ! MV DEBUG 648 645 649 646 !--------------------------- … … 666 663 CALL lbc_lnk( 'icedyn_rhg_vp', ztaux_oi_rhsu, 'U', -1., ztauy_oi_rhsv, 'V', -1.) 667 664 668 CALL iom_put( 'zmU_t' , zmU_t ) ! MV DEBUG669 CALL iom_put( 'zmV_t' , zmV_t ) ! MV DEBUG670 CALL iom_put( 'ztaux_oi_rhsu' , ztaux_oi_rhsu ) ! MV DEBUG671 CALL iom_put( 'ztauy_oi_rhsv' , ztauy_oi_rhsv ) ! MV DEBUG672 CALL iom_put( 'zrhsu' , zrhsu ) ! MV DEBUG673 CALL iom_put( 'zrhsv' , zrhsv ) ! MV DEBUG674 665 !!$ CALL iom_put( 'zmU_t' , zmU_t ) ! MV DEBUG 666 !!$ CALL iom_put( 'zmV_t' , zmV_t ) ! MV DEBUG 667 !!$ CALL iom_put( 'ztaux_oi_rhsu' , ztaux_oi_rhsu ) ! MV DEBUG 668 !!$ CALL iom_put( 'ztauy_oi_rhsv' , ztauy_oi_rhsv ) ! MV DEBUG 669 !!$ CALL iom_put( 'zrhsu' , zrhsu ) ! MV DEBUG 670 !!$ CALL iom_put( 'zrhsv' , zrhsv ) ! MV DEBUG 671 !!$ 675 672 ! inner domain calculations -> no lbclnk 676 673 … … 773 770 CALL lbc_lnk( 'icedyn_rhg_vp', zEU , 'U', 1., zEV , 'V', 1. ) 774 771 775 CALL iom_put( 'zAU' , zAU ) ! MV DEBUG776 CALL iom_put( 'zBU' , zBU ) ! MV DEBUG777 CALL iom_put( 'zCU' , zCU ) ! MV DEBUG778 CALL iom_put( 'zDU' , zDU ) ! MV DEBUG779 CALL iom_put( 'zEU' , zEU ) ! MV DEBUG780 CALL iom_put( 'zAV' , zAV ) ! MV DEBUG781 CALL iom_put( 'zBV' , zBV ) ! MV DEBUG782 CALL iom_put( 'zCV' , zCV ) ! MV DEBUG783 CALL iom_put( 'zDV' , zDV ) ! MV DEBUG784 CALL iom_put( 'zEV' , zEV ) ! MV DEBUG772 !!$ CALL iom_put( 'zAU' , zAU ) ! MV DEBUG 773 !!$ CALL iom_put( 'zBU' , zBU ) ! MV DEBUG 774 !!$ CALL iom_put( 'zCU' , zCU ) ! MV DEBUG 775 !!$ CALL iom_put( 'zDU' , zDU ) ! MV DEBUG 776 !!$ CALL iom_put( 'zEU' , zEU ) ! MV DEBUG 777 !!$ CALL iom_put( 'zAV' , zAV ) ! MV DEBUG 778 !!$ CALL iom_put( 'zBV' , zBV ) ! MV DEBUG 779 !!$ CALL iom_put( 'zCV' , zCV ) ! MV DEBUG 780 !!$ CALL iom_put( 'zDV' , zDV ) ! MV DEBUG 781 !!$ CALL iom_put( 'zEV' , zEV ) ! MV DEBUG 785 782 786 783 !------------------------------------------------------------------------------! … … 1103 1100 CALL lbc_lnk( 'icedyn_rhg_vp', zCU_prime , 'U', 1., zCV_prime , 'V', 1. ) 1104 1101 1105 CALL iom_put( 'zFU' , zFU ) ! MV DEBUG1106 CALL iom_put( 'zBU_prime' , zBU_prime ) ! MV DEBUG1107 CALL iom_put( 'zCU_prime' , zCU_prime ) ! MV DEBUG1108 CALL iom_put( 'zFU_prime' , zFU_prime ) ! MV DEBUG1109 1110 CALL iom_put( 'zFV' , zFV ) ! MV DEBUG1111 CALL iom_put( 'zBV_prime' , zBV_prime ) ! MV DEBUG1112 CALL iom_put( 'zCV_prime' , zCV_prime ) ! MV DEBUG1113 CALL iom_put( 'zFV_prime' , zFV_prime ) ! MV DEBUG1102 !!$ CALL iom_put( 'zFU' , zFU ) ! MV DEBUG 1103 !!$ CALL iom_put( 'zBU_prime' , zBU_prime ) ! MV DEBUG 1104 !!$ CALL iom_put( 'zCU_prime' , zCU_prime ) ! MV DEBUG 1105 !!$ CALL iom_put( 'zFU_prime' , zFU_prime ) ! MV DEBUG 1106 !!$ 1107 !!$ CALL iom_put( 'zFV' , zFV ) ! MV DEBUG 1108 !!$ CALL iom_put( 'zBV_prime' , zBV_prime ) ! MV DEBUG 1109 !!$ CALL iom_put( 'zCV_prime' , zCV_prime ) ! MV DEBUG 1110 !!$ CALL iom_put( 'zFV_prime' , zFV_prime ) ! MV DEBUG 1114 1111 1115 1112 CALL lbc_lnk( 'icedyn_rhg_vp', u_ice, 'U', -1., v_ice, 'V', -1. ) 1116 1113 1117 IF ( lwp ) WRITE(numout,*) ' We are about to output uice_dbg '1118 IF( iom_use('uice_dbg' ) ) CALL iom_put( 'uice_dbg' , u_ice ) ! ice velocity u after solver1119 IF( iom_use('vice_dbg' ) ) CALL iom_put( 'vice_dbg' , v_ice ) ! ice velocity v after solver1114 !!$ IF ( lwp ) WRITE(numout,*) ' We are about to output uice_dbg ' 1115 !!$ IF( iom_use('uice_dbg' ) ) CALL iom_put( 'uice_dbg' , u_ice ) ! ice velocity u after solver 1116 !!$ IF( iom_use('vice_dbg' ) ) CALL iom_put( 'vice_dbg' , v_ice ) ! ice velocity v after solver 1120 1117 1121 1118 !------------------------------------------------------------------------------! … … 1134 1131 ENDIF 1135 1132 1136 ! MV DEBUG test - replace ice velocity by ocean current to give the model the means to go ahead1137 DO jj = 2, jpj - 11138 DO ji = 2, jpi - 11139 1140 u_ice(ji,jj) = zmsk00x(ji,jj) &1141 & * ( zmsk01x(ji,jj) * u_oce(ji,jj) * 0.01_wp &1142 + ( 1._wp - zmsk01x(ji,jj) ) * u_oce(ji,jj) * 0.01_wp )1143 1144 v_ice(ji,jj) = zmsk00y(ji,jj) &1145 & * ( zmsk01y(ji,jj) * v_oce(ji,jj) * 0.01_wp &1146 + ( 1._wp - zmsk01y(ji,jj) ) * v_oce(ji,jj) * 0.01_wp )1147 1148 END DO1149 END DO1150 1151 CALL lbc_lnk( 'icedyn_rhg_vp', u_ice, 'U', -1., v_ice, 'V', -1. )1152 1153 IF ( lwp ) WRITE(numout,*) ' Velocity replaced '1133 !!$ ! MV DEBUG test - replace ice velocity by ocean current to give the model the means to go ahead 1134 !!$ DO jj = 2, jpj - 1 1135 !!$ DO ji = 2, jpi - 1 1136 !!$ 1137 !!$ u_ice(ji,jj) = zmsk00x(ji,jj) & 1138 !!$ & * ( zmsk01x(ji,jj) * u_oce(ji,jj) * 0.01_wp & 1139 !!$ + ( 1._wp - zmsk01x(ji,jj) ) * u_oce(ji,jj) * 0.01_wp ) 1140 !!$ 1141 !!$ v_ice(ji,jj) = zmsk00y(ji,jj) & 1142 !!$ & * ( zmsk01y(ji,jj) * v_oce(ji,jj) * 0.01_wp & 1143 !!$ + ( 1._wp - zmsk01y(ji,jj) ) * v_oce(ji,jj) * 0.01_wp ) 1144 !!$ 1145 !!$ END DO 1146 !!$ END DO 1147 !!$ 1148 !!$ CALL lbc_lnk( 'icedyn_rhg_vp', u_ice, 'U', -1., v_ice, 'V', -1. ) 1149 !!$ 1150 !!$ IF ( lwp ) WRITE(numout,*) ' Velocity replaced ' 1154 1151 1155 1152 ! END DEBUG … … 1365 1362 zfac = zp_deltastar_t(ji,jj) 1366 1363 zsig1 = zfac * ( pdivu_i(ji,jj) - zdeltastar_t(ji,jj) ) 1367 zsig1 = 0._wp !!! FUCKING DEBUG TEST !!!1364 !!$ zsig1 = 0._wp !!! FUCKING DEBUG TEST !!! 1368 1365 zsig2 = zfac * z1_ecc2 * zten_i(ji,jj) 1369 1366 zsig12 = zfac * z1_ecc2 * pshear_i(ji,jj)
Note: See TracChangeset
for help on using the changeset viewer.
