Changeset 14644 for NEMO/branches/2020/dev_r14116_HPC-04_mcastril_Mixed_Precision_implementation_final/src/ICE/icedyn_rhg_vp.F90
- Timestamp:
- 2021-03-26T15:33:49+01:00 (3 years ago)
- Location:
- NEMO/branches/2020/dev_r14116_HPC-04_mcastril_Mixed_Precision_implementation_final
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2020/dev_r14116_HPC-04_mcastril_Mixed_Precision_implementation_final
- Property svn:externals
-
old new 9 9 10 10 # SETTE 11 ^/utils/CI/sette _wave@13990sette11 ^/utils/CI/sette@14244 sette
-
- Property svn:externals
-
NEMO/branches/2020/dev_r14116_HPC-04_mcastril_Mixed_Precision_implementation_final/src/ICE/icedyn_rhg_vp.F90
r14219 r14644 461 461 462 462 CALL lbc_lnk( 'icedyn_rhg_vp', zds, 'F', 1._wp ) ! MV TEST could be un-necessary according to Gurvan 463 CALL iom_put( 'zds' , zds ) ! MV DEBUG463 CALL iom_put( 'zds' , zds ) ! MV DEBUG 464 464 465 465 IF( lwp ) WRITE(numout,*) ' outer loop 1a i_out : ', i_out … … 506 506 END DO 507 507 508 CALL lbc_lnk _multi( 'icedyn_rhg_vp', zp_deltastar_t , 'T', 1._wp , zzt , 'T', 1._wp, zet, 'T', 1._wp )508 CALL lbc_lnk( 'icedyn_rhg_vp', zp_deltastar_t , 'T', 1._wp , zzt , 'T', 1._wp, zet, 'T', 1._wp ) 509 509 510 510 CALL iom_put( 'zzt' , zzt ) ! MV DEBUG … … 527 527 528 528 CALL lbc_lnk( 'icedyn_rhg_vp', zef, 'F', 1._wp ) 529 CALL iom_put( 'zef' 529 CALL iom_put( 'zef' , zef ) ! MV DEBUG 530 530 IF( lwp ) WRITE(numout,*) ' outer loop 1c i_out : ', i_out 531 531 … … 567 567 IF( lwp ) WRITE(numout,*) ' outer loop 1d i_out : ', i_out 568 568 569 CALL lbc_lnk _multi( 'icedyn_rhg_vp', zCwU , 'U', -1._wp, zCwV, 'V', -1._wp )570 CALL lbc_lnk _multi( 'icedyn_rhg_vp', zCorU, 'U', -1._wp, zCorV, 'V', -1._wp )569 CALL lbc_lnk( 'icedyn_rhg_vp', zCwU , 'U', -1._wp, zCwV, 'V', -1._wp ) 570 CALL lbc_lnk( 'icedyn_rhg_vp', zCorU, 'U', -1._wp, zCorV, 'V', -1._wp ) 571 571 572 572 CALL iom_put( 'zCwU' , zCwU ) ! MV DEBUG … … 674 674 END DO 675 675 676 CALL lbc_lnk _multi( 'icedyn_rhg_vp', zrhsu, 'U', -1._wp, zrhsv, 'V', -1._wp)677 CALL lbc_lnk _multi( 'icedyn_rhg_vp', zmU_t, 'U', -1._wp, zmV_t, 'V', -1._wp)678 CALL lbc_lnk _multi( 'icedyn_rhg_vp', ztaux_oi_rhsu, 'U', -1._wp, ztauy_oi_rhsv, 'V', -1._wp)676 CALL lbc_lnk( 'icedyn_rhg_vp', zrhsu, 'U', -1._wp, zrhsv, 'V', -1._wp) 677 CALL lbc_lnk( 'icedyn_rhg_vp', zmU_t, 'U', -1._wp, zmV_t, 'V', -1._wp) 678 CALL lbc_lnk( 'icedyn_rhg_vp', ztaux_oi_rhsu, 'U', -1., ztauy_oi_rhsv, 'V', -1.) 679 679 680 680 CALL iom_put( 'zmU_t' , zmU_t ) ! MV DEBUG … … 779 779 END DO 780 780 781 CALL lbc_lnk _multi( 'icedyn_rhg_vp', zAU , 'U', 1._wp, zAV , 'V', 1._wp )782 CALL lbc_lnk _multi( 'icedyn_rhg_vp', zBU , 'U', 1._wp, zBV , 'V', 1._wp )783 CALL lbc_lnk _multi( 'icedyn_rhg_vp', zCU , 'U', 1._wp, zCV , 'V', 1._wp )784 CALL lbc_lnk _multi( 'icedyn_rhg_vp', zDU , 'U', 1._wp, zDV , 'V', 1._wp )785 CALL lbc_lnk _multi( 'icedyn_rhg_vp', zEU , 'U', 1._wp, zEV , 'V', 1._wp )781 CALL lbc_lnk( 'icedyn_rhg_vp', zAU , 'U', 1._wp, zAV , 'V', 1._wp ) 782 CALL lbc_lnk( 'icedyn_rhg_vp', zBU , 'U', 1._wp, zBV , 'V', 1._wp ) 783 CALL lbc_lnk( 'icedyn_rhg_vp', zCU , 'U', 1._wp, zCV , 'V', 1._wp ) 784 CALL lbc_lnk( 'icedyn_rhg_vp', zDU , 'U', 1._wp, zDV , 'V', 1._wp ) 785 CALL lbc_lnk( 'icedyn_rhg_vp', zEU , 'U', 1._wp, zEV , 'V', 1._wp ) 786 786 787 787 CALL iom_put( 'zAU' , zAU ) ! MV DEBUG … … 885 885 END DO 886 886 887 CALL lbc_lnk _multi( 'icedyn_rhg_vp', zFU_prime, 'U', 1._wp, zBU_prime, 'U', 1._wp )887 CALL lbc_lnk( 'icedyn_rhg_vp', zFU_prime, 'U', 1._wp, zBU_prime, 'U', 1._wp ) 888 888 889 889 !----------------------------- … … 983 983 END DO 984 984 985 CALL lbc_lnk _multi( 'icedyn_rhg_vp', zFV_prime, 'V', 1._wp, zBV_prime, 'V', 1._wp )985 CALL lbc_lnk( 'icedyn_rhg_vp', zFV_prime, 'V', 1._wp, zBV_prime, 'V', 1._wp ) 986 986 987 987 !----------------------------- … … 1020 1020 ENDIF ! ll_v_iterate 1021 1021 1022 CALL lbc_lnk _multi( 'icedyn_rhg_vp', u_ice, 'U', -1._wp, v_ice, 'V', -1._wp )1022 CALL lbc_lnk( 'icedyn_rhg_vp', u_ice, 'U', -1._wp, v_ice, 'V', -1._wp ) 1023 1023 1024 1024 !-------------------------------------------------------------------------------------- … … 1110 1110 IF ( lwp ) WRITE(numout,*) ' We are out of outer loop ' 1111 1111 1112 CALL lbc_lnk _multi( 'icedyn_rhg_vp', zFU , 'U', 1._wp, zFV , 'V', 1._wp)1113 CALL lbc_lnk _multi( 'icedyn_rhg_vp', zBU_prime , 'U', 1._wp, zBV_prime , 'V', 1._wp )1114 CALL lbc_lnk _multi( 'icedyn_rhg_vp', zFU_prime , 'U', 1._wp, zFV_prime , 'V', 1._wp )1115 CALL lbc_lnk _multi( 'icedyn_rhg_vp', zCU_prime , 'U', 1._wp, zCV_prime , 'V', 1._wp )1112 CALL lbc_lnk( 'icedyn_rhg_vp', zFU , 'U', 1., zFV , 'V', 1. ) 1113 CALL lbc_lnk( 'icedyn_rhg_vp', zBU_prime , 'U', 1._wp, zBV_prime , 'V', 1._wp ) 1114 CALL lbc_lnk( 'icedyn_rhg_vp', zFU_prime , 'U', 1._wp, zFV_prime , 'V', 1._wp ) 1115 CALL lbc_lnk( 'icedyn_rhg_vp', zCU_prime , 'U', 1._wp, zCV_prime , 'V', 1._wp ) 1116 1116 1117 1117 CALL iom_put( 'zFU' , zFU ) ! MV DEBUG … … 1125 1125 CALL iom_put( 'zFV_prime' , zFV_prime ) ! MV DEBUG 1126 1126 1127 CALL lbc_lnk _multi( 'icedyn_rhg_vp', u_ice, 'U', -1._wp, v_ice, 'V', -1._wp )1127 CALL lbc_lnk( 'icedyn_rhg_vp', u_ice, 'U', -1._wp, v_ice, 'V', -1._wp ) 1128 1128 1129 1129 IF ( lwp ) WRITE(numout,*) ' We are about to output uice_dbg ' … … 1161 1161 END DO 1162 1162 1163 CALL lbc_lnk _multi( 'icedyn_rhg_vp', u_ice, 'U', -1._wp, v_ice, 'V', -1._wp )1163 CALL lbc_lnk( 'icedyn_rhg_vp', u_ice, 'U', -1._wp, v_ice, 'V', -1._wp ) 1164 1164 1165 1165 IF ( lwp ) WRITE(numout,*) ' Velocity replaced ' … … 1222 1222 IF ( lwp ) WRITE(numout,*) ' Deformation recalculated ' 1223 1223 1224 CALL lbc_lnk _multi( 'icedyn_rhg_vp', pshear_i, 'T', 1._wp, pdivu_i, 'T', 1._wp, pdelta_i, 'T', 1._wp )1224 CALL lbc_lnk( 'icedyn_rhg_vp', pshear_i, 'T', 1._wp, pdivu_i, 'T', 1._wp, pdelta_i, 'T', 1._wp ) 1225 1225 1226 1226 !------------------------------------------------------------------------------! … … 1249 1249 END DO 1250 1250 1251 CALL lbc_lnk _multi( 'icedyn_rhg_vp', zs1, 'T', 1._wp, zs2, 'T', 1._wp, zs12, 'T', 1._wp )1251 CALL lbc_lnk( 'icedyn_rhg_vp', zs1, 'T', 1._wp, zs2, 'T', 1._wp, zs12, 'T', 1._wp ) 1252 1252 1253 1253 ENDIF … … 1307 1307 1308 1308 ! 1309 CALL lbc_lnk _multi( 'icedyn_rhg_vp', ztaux_oi, 'U', -1._wp, ztauy_oi, 'V', -1._wp, ztaux_ai, 'U', -1._wp, ztauy_ai, 'V', -1._wp) !, &1310 ! & ztaux_bi, 'U', -1., ztauy_bi, 'V', -1.)1309 CALL lbc_lnk( 'icedyn_rhg_vp', ztaux_oi, 'U', -1._wp, ztauy_oi, 'V', -1._wp, ztaux_ai, 'U', -1._wp, ztauy_ai, 'V', -1. ) !, & 1310 ! & ztaux_bi, 'U', -1.)_wp, ztauy_bi, 'V', -1._wp ) 1311 1311 ! 1312 1312 CALL iom_put( 'utau_oi' , ztaux_oi * zmsk00 ) … … 1348 1348 END DO 1349 1349 1350 CALL lbc_lnk _multi( 'icedyn_rhg_vp', zsig_I, 'T', 1._wp, zsig_II, 'T', 1._wp)1350 CALL lbc_lnk( 'icedyn_rhg_vp', zsig_I, 'T', 1._wp, zsig_II, 'T', 1._wp) 1351 1351 1352 1352 IF( iom_use('normstr') ) CALL iom_put( 'normstr' , zsig_I(:,:) * zmsk00(:,:) ) ! Normal stress … … 1393 1393 IF ( lwp ) WRITE(numout,*) 'Some shitty stress work done' 1394 1394 ! 1395 CALL lbc_lnk _multi( 'icedyn_rhg_vp', zsig1_p, 'T', 1._wp, zsig2_p, 'T', 1._wp)1395 CALL lbc_lnk( 'icedyn_rhg_vp', zsig1_p, 'T', 1._wp, zsig2_p, 'T', 1._wp) 1396 1396 ! 1397 1397 IF ( lwp ) WRITE(numout,*) ' Beauaaaarflblbllll ' … … 1423 1423 END DO 1424 1424 ! 1425 CALL lbc_lnk _multi( 'icedyn_rhg_vp', zspgU, 'U', -1._wp, zspgV, 'V', -1._wp, &1426 & 1425 CALL lbc_lnk( 'icedyn_rhg_vp', zspgU, 'U', -1._wp, zspgV, 'V', -1._wp, & 1426 & zCorU, 'U', -1._wp, zCorV, 'V', -1._wp ) 1427 1427 ! 1428 1428 CALL iom_put( 'dssh_dx' , zspgU * zmsk00 ) ! Sea-surface tilt term in force balance (x) … … 1453 1453 END DO 1454 1454 1455 CALL lbc_lnk _multi( 'icedyn_rhg_vp', zfU, 'U', -1._wp, zfV, 'V', -1._wp )1455 CALL lbc_lnk( 'icedyn_rhg_vp', zfU, 'U', -1._wp, zfV, 'V', -1._wp ) 1456 1456 1457 1457 CALL iom_put( 'intstrx' , zfU * zmsk00 ) ! Internal force term in force balance (x) … … 1485 1485 END DO 1486 1486 1487 CALL lbc_lnk _multi( 'icedyn_rhg_vp', zdiag_xmtrp_ice, 'U', -1._wp, zdiag_ymtrp_ice, 'V', -1._wp, &1488 & 1489 & 1487 CALL lbc_lnk( 'icedyn_rhg_vp', zdiag_xmtrp_ice, 'U', -1._wp, zdiag_ymtrp_ice, 'V', -1._wp, & 1488 & zdiag_xmtrp_snw, 'U', -1._wp, zdiag_ymtrp_snw, 'V', -1._wp, & 1489 & zdiag_xatrp , 'U', -1._wp, zdiag_yatrp , 'V', -1._wp ) 1490 1490 1491 1491 CALL iom_put( 'xmtrpice' , zdiag_xmtrp_ice ) ! X-component of sea-ice mass transport (kg/s)
Note: See TracChangeset
for help on using the changeset viewer.