Changeset 14219 for NEMO/branches/2020/dev_r14116_HPC-04_mcastril_Mixed_Precision_implementation_final/src/ICE/icedyn_rhg_vp.F90
- Timestamp:
- 2020-12-18T18:52:57+01:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2020/dev_r14116_HPC-04_mcastril_Mixed_Precision_implementation_final/src/ICE/icedyn_rhg_vp.F90
r14072 r14219 460 460 END DO 461 461 462 CALL lbc_lnk( 'icedyn_rhg_vp', zds, 'F', 1. ) ! MV TEST could be un-necessary according to Gurvan462 CALL lbc_lnk( 'icedyn_rhg_vp', zds, 'F', 1._wp ) ! MV TEST could be un-necessary according to Gurvan 463 463 CALL iom_put( 'zds' , zds ) ! MV DEBUG 464 464 … … 506 506 END DO 507 507 508 CALL lbc_lnk_multi( 'icedyn_rhg_vp', zp_deltastar_t , 'T', 1. , zzt , 'T', 1., zet, 'T', 1.)508 CALL lbc_lnk_multi( '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 … … 526 526 END DO 527 527 528 CALL lbc_lnk( 'icedyn_rhg_vp', zef, 'F', 1. )528 CALL lbc_lnk( 'icedyn_rhg_vp', zef, 'F', 1._wp ) 529 529 CALL iom_put( 'zef' , zef ) ! MV DEBUG 530 530 IF( lwp ) WRITE(numout,*) ' outer loop 1c i_out : ', i_out … … 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. , zCwV, 'V', -1.)570 CALL lbc_lnk_multi( 'icedyn_rhg_vp', zCorU, 'U', -1. , zCorV, 'V', -1.)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 ) 571 571 572 572 CALL iom_put( 'zCwU' , zCwU ) ! MV DEBUG … … 628 628 END DO 629 629 630 CALL lbc_lnk( 'icedyn_rhg_vp', zs12_rhsu, 'F', 1. )631 CALL lbc_lnk( 'icedyn_rhg_vp', zs12_rhsv, 'F', 1. )630 CALL lbc_lnk( 'icedyn_rhg_vp', zs12_rhsu, 'F', 1._wp ) 631 CALL lbc_lnk( 'icedyn_rhg_vp', zs12_rhsv, 'F', 1._wp ) 632 632 633 633 CALL iom_put( 'zs12_rhsu' , zs12_rhsu ) ! MV DEBUG … … 674 674 END DO 675 675 676 CALL lbc_lnk_multi( 'icedyn_rhg_vp', zrhsu, 'U', -1. , zrhsv, 'V', -1.)677 CALL lbc_lnk_multi( 'icedyn_rhg_vp', zmU_t, 'U', -1. , zmV_t, 'V', -1.)678 CALL lbc_lnk_multi( 'icedyn_rhg_vp', ztaux_oi_rhsu, 'U', -1. , ztauy_oi_rhsv, 'V', -1.)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) 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. , zAV , 'V', 1.)782 CALL lbc_lnk_multi( 'icedyn_rhg_vp', zBU , 'U', 1. , zBV , 'V', 1.)783 CALL lbc_lnk_multi( 'icedyn_rhg_vp', zCU , 'U', 1. , zCV , 'V', 1.)784 CALL lbc_lnk_multi( 'icedyn_rhg_vp', zDU , 'U', 1. , zDV , 'V', 1.)785 CALL lbc_lnk_multi( 'icedyn_rhg_vp', zEU , 'U', 1. , zEV , 'V', 1.)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 ) 786 786 787 787 CALL iom_put( 'zAU' , zAU ) ! MV DEBUG … … 867 867 END DO 868 868 869 CALL lbc_lnk( 'icedyn_rhg_vp', zFU, 'U', 1. )869 CALL lbc_lnk( 'icedyn_rhg_vp', zFU, 'U', 1._wp ) 870 870 871 871 !--------------- … … 885 885 END DO 886 886 887 CALL lbc_lnk_multi( 'icedyn_rhg_vp', zFU_prime, 'U', 1. , zBU_prime, 'U', 1.)887 CALL lbc_lnk_multi( 'icedyn_rhg_vp', zFU_prime, 'U', 1._wp, zBU_prime, 'U', 1._wp ) 888 888 889 889 !----------------------------- … … 965 965 END DO 966 966 967 CALL lbc_lnk _multi( 'icedyn_rhg_vp', zFV, 'V', 1.)967 CALL lbc_lnk( 'icedyn_rhg_vp', zFV, 'V', 1._wp) 968 968 969 969 !--------------- … … 983 983 END DO 984 984 985 CALL lbc_lnk_multi( 'icedyn_rhg_vp', zFV_prime, 'V', 1. , zBV_prime, 'V', 1.)985 CALL lbc_lnk_multi( '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. , v_ice, 'V', -1.)1022 CALL lbc_lnk_multi( '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. , zFV , 'V', 1.)1113 CALL lbc_lnk_multi( 'icedyn_rhg_vp', zBU_prime , 'U', 1. , zBV_prime , 'V', 1.)1114 CALL lbc_lnk_multi( 'icedyn_rhg_vp', zFU_prime , 'U', 1. , zFV_prime , 'V', 1.)1115 CALL lbc_lnk_multi( 'icedyn_rhg_vp', zCU_prime , 'U', 1. , zCV_prime , 'V', 1.)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 ) 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. , v_ice, 'V', -1.)1127 CALL lbc_lnk_multi( '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. , v_ice, 'V', -1.)1163 CALL lbc_lnk_multi( '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. , pdivu_i, 'T', 1., pdelta_i, 'T', 1.)1224 CALL lbc_lnk_multi( '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. , zs2, 'T', 1., zs12, 'T', 1.)1251 CALL lbc_lnk_multi( 'icedyn_rhg_vp', zs1, 'T', 1._wp, zs2, 'T', 1._wp, zs12, 'T', 1._wp ) 1252 1252 1253 1253 ENDIF … … 1268 1268 END DO 1269 1269 1270 CALL lbc_lnk( 'icedyn_rhg_vp', zs12f, 'F', 1. )1270 CALL lbc_lnk( 'icedyn_rhg_vp', zs12f, 'F', 1._wp ) 1271 1271 1272 1272 ENDIF … … 1307 1307 1308 1308 ! 1309 CALL lbc_lnk_multi( 'icedyn_rhg_vp', ztaux_oi, 'U', -1. , ztauy_oi, 'V', -1., ztaux_ai, 'U', -1., ztauy_ai, 'V', -1.) !, &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 1310 ! & ztaux_bi, 'U', -1., ztauy_bi, 'V', -1. ) 1311 1311 ! … … 1348 1348 END DO 1349 1349 1350 CALL lbc_lnk_multi( 'icedyn_rhg_vp', zsig_I, 'T', 1. , zsig_II, 'T', 1.)1350 CALL lbc_lnk_multi( '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. , zsig2_p, 'T', 1.)1395 CALL lbc_lnk_multi( '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. , zspgV, 'V', -1., &1426 & zCorU, 'U', -1. , zCorV, 'V', -1.)1425 CALL lbc_lnk_multi( '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. , zfV, 'V', -1.)1455 CALL lbc_lnk_multi( '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. , zdiag_ymtrp_ice, 'V', -1., &1488 & zdiag_xmtrp_snw, 'U', -1. , zdiag_ymtrp_snw, 'V', -1., &1489 & zdiag_xatrp , 'U', -1. , zdiag_yatrp , 'V', -1.)1487 CALL lbc_lnk_multi( '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.