Changeset 13295 for NEMO/trunk/src/SWE/sbcice_cice.F90
- Timestamp:
- 2020-07-10T20:24:21+02:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/trunk/src/SWE/sbcice_cice.F90
r12983 r13295 219 219 ! T point to U point 220 220 ! T point to V point 221 DO_2D _10_10221 DO_2D( 1, 0, 1, 0 ) 222 222 fr_iu(ji,jj)=0.5*(fr_i(ji,jj)+fr_i(ji+1,jj))*umask(ji,jj,1) 223 223 fr_iv(ji,jj)=0.5*(fr_i(ji,jj)+fr_i(ji,jj+1))*vmask(ji,jj,1) … … 316 316 ! x comp of wind stress (CI_1) 317 317 ! U point to F point 318 DO_2D _10_11318 DO_2D( 1, 0, 1, 1 ) 319 319 ztmp(ji,jj) = 0.5 * ( fr_iu(ji,jj) * utau(ji,jj) & 320 320 + fr_iu(ji,jj+1) * utau(ji,jj+1) ) * fmask(ji,jj,1) … … 324 324 ! y comp of wind stress (CI_2) 325 325 ! V point to F point 326 DO_2D _11_10326 DO_2D( 1, 1, 1, 0 ) 327 327 ztmp(ji,jj) = 0.5 * ( fr_iv(ji,jj) * vtau(ji,jj) & 328 328 + fr_iv(ji+1,jj) * vtau(ji+1,jj) ) * fmask(ji,jj,1) … … 339 339 qla_ice(:,:,1)= - ( emp_ice(:,:)+sprecip(:,:) ) * rLsub 340 340 ! End of temporary code 341 DO_2D _11_11341 DO_2D( 1, 1, 1, 1 ) 342 342 IF(fr_i(ji,jj).eq.0.0) THEN 343 343 DO jl=1,ncat … … 441 441 ! x comp and y comp of surface ocean current 442 442 ! U point to F point 443 DO_2D _10_11443 DO_2D( 1, 0, 1, 1 ) 444 444 ztmp(ji,jj)=0.5*(ssu_m(ji,jj)+ssu_m(ji,jj+1))*fmask(ji,jj,1) 445 445 END_2D … … 447 447 448 448 ! V point to F point 449 DO_2D _11_10449 DO_2D( 1, 1, 1, 0 ) 450 450 ztmp(ji,jj)=0.5*(ssv_m(ji,jj)+ssv_m(ji+1,jj))*fmask(ji,jj,1) 451 451 END_2D … … 471 471 ! x comp and y comp of sea surface slope (on F points) 472 472 ! T point to F point 473 DO_2D _10_10473 DO_2D( 1, 0, 1, 0 ) 474 474 ztmp(ji,jj)=0.5 * ( (zpice(ji+1,jj )-zpice(ji,jj )) * r1_e1u(ji,jj ) & 475 475 & + (zpice(ji+1,jj+1)-zpice(ji,jj+1)) * r1_e1u(ji,jj+1) ) * fmask(ji,jj,1) … … 478 478 479 479 ! T point to F point 480 DO_2D _10_10480 DO_2D( 1, 0, 1, 0 ) 481 481 ztmp(ji,jj)=0.5 * ( (zpice(ji ,jj+1)-zpice(ji ,jj)) * r1_e2v(ji ,jj) & 482 482 & + (zpice(ji+1,jj+1)-zpice(ji+1,jj)) * r1_e2v(ji+1,jj) ) * fmask(ji,jj,1) … … 507 507 ss_iou(:,:)=0.0 508 508 ! F point to U point 509 DO_2D _00_00509 DO_2D( 0, 0, 0, 0 ) 510 510 ss_iou(ji,jj) = 0.5 * ( ztmp1(ji,jj-1) + ztmp1(ji,jj) ) * umask(ji,jj,1) 511 511 END_2D … … 517 517 ! F point to V point 518 518 519 DO_2D _10_00519 DO_2D( 1, 0, 0, 0 ) 520 520 ss_iov(ji,jj) = 0.5 * ( ztmp1(ji-1,jj) + ztmp1(ji,jj) ) * vmask(ji,jj,1) 521 521 END_2D … … 601 601 CALL lbc_lnk( 'sbcice_cice', qsr , 'T', 1. ) 602 602 603 DO_2D _11_11603 DO_2D( 1, 1, 1, 1 ) 604 604 nfrzmlt(ji,jj)=MAX(nfrzmlt(ji,jj),0.0) 605 605 END_2D … … 625 625 ! T point to U point 626 626 ! T point to V point 627 DO_2D _10_10627 DO_2D( 1, 0, 1, 0 ) 628 628 fr_iu(ji,jj)=0.5*(fr_i(ji,jj)+fr_i(ji+1,jj))*umask(ji,jj,1) 629 629 fr_iv(ji,jj)=0.5*(fr_i(ji,jj)+fr_i(ji,jj+1))*vmask(ji,jj,1) … … 985 985 986 986 pn(:,:)=0.0 987 DO_2D _10_10987 DO_2D( 1, 0, 1, 0 ) 988 988 pn(ji,jj)=pc(ji+1-ji_off,jj+1-jj_off,1) 989 989 END_2D
Note: See TracChangeset
for help on using the changeset viewer.