- Timestamp:
- 2021-05-11T09:51:03+02:00 (3 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2021/dev_r14116_HPC-10_mcastril_Mixed_Precision_implementation/src/ICE/icedyn_rhg_vp.F90
r14644 r14831 1438 1438 DO jj = 2, jpj - 1 1439 1439 DO ji = 2, jpi - 1 1440 z fU(ji,jj) = 0.5_wp * ( ( zs1(ji+1,jj) - zs1(ji,jj) ) * e2u(ji,jj) &1440 zFU(ji,jj) = 0.5_wp * ( ( zs1(ji+1,jj) - zs1(ji,jj) ) * e2u(ji,jj) & 1441 1441 & + ( zs2(ji+1,jj) * e2t(ji+1,jj) * e2t(ji+1,jj) - zs2(ji,jj) * e2t(ji,jj) * e2t(ji,jj) & 1442 1442 & ) * r1_e2u(ji,jj) & … … 1444 1444 & ) * 2._wp * r1_e1u(ji,jj) & 1445 1445 & ) * r1_e1e2u(ji,jj) 1446 z fV(ji,jj) = 0.5_wp * ( ( zs1(ji,jj+1) - zs1(ji,jj) ) * e1v(ji,jj) &1446 zFV(ji,jj) = 0.5_wp * ( ( zs1(ji,jj+1) - zs1(ji,jj) ) * e1v(ji,jj) & 1447 1447 & - ( zs2(ji,jj+1) * e1t(ji,jj+1) * e1t(ji,jj+1) - zs2(ji,jj) * e1t(ji,jj) * e1t(ji,jj) & 1448 1448 & ) * r1_e1v(ji,jj) & … … 1453 1453 END DO 1454 1454 1455 CALL lbc_lnk( 'icedyn_rhg_vp', z fU, 'U', -1._wp, zfV, 'V', -1._wp )1456 1457 CALL iom_put( 'intstrx' , z fU * zmsk00 ) ! Internal force term in force balance (x)1458 CALL iom_put( 'intstry' , z fV * zmsk00 ) ! Internal force term in force balance (y)1455 CALL lbc_lnk( 'icedyn_rhg_vp', zFU, 'U', -1._wp, zFV, 'V', -1._wp ) 1456 1457 CALL iom_put( 'intstrx' , zFU * zmsk00 ) ! Internal force term in force balance (x) 1458 CALL iom_put( 'intstry' , zFV * zmsk00 ) ! Internal force term in force balance (y) 1459 1459 1460 1460 ENDIF
Note: See TracChangeset
for help on using the changeset viewer.