Changeset 12546 for NEMO/branches/2020/dev_r12512_HPC-04_mcastril_Mixed_Precision_implementation/src/OCE/DIA/diawri.F90
- Timestamp:
- 2020-03-13T11:06:44+01:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2020/dev_r12512_HPC-04_mcastril_Mixed_Precision_implementation/src/OCE/DIA/diawri.F90
r12493 r12546 183 183 ! 184 184 END_2D 185 CALL lbc_lnk( 'diawri', z2d, 'T', 1. )185 CALL lbc_lnk( 'diawri', z2d, 'T', 1.0_wp ) 186 186 CALL iom_put( "taubot", z2d ) 187 187 ENDIF … … 237 237 & * umask(ji,jj,1) * umask(ji-1,jj,1) * vmask(ji,jj,1) * umask(ji,jj-1,1) 238 238 END_2D 239 CALL lbc_lnk( 'diawri', z2d, 'T', 1. )239 CALL lbc_lnk( 'diawri', z2d, 'T', 1.0_wp ) 240 240 CALL iom_put( "sstgrad2", z2d ) ! square of module of sst gradient 241 241 z2d(:,:) = SQRT( z2d(:,:) ) … … 269 269 & + vv(ji,jj ,jk,Kmm)**2 * e1v(ji,jj ) * e3v(ji,jj ,jk,Kmm) ) 270 270 END_3D 271 CALL lbc_lnk( 'diawri', z3d, 'T', 1. )271 CALL lbc_lnk( 'diawri', z3d, 'T', 1.0_wp ) 272 272 CALL iom_put( "eken", z3d ) ! kinetic energy 273 273 ENDIF … … 291 291 z2d(ji,jj) = z2d(ji,jj) + z3d(ji,jj,jk) * ( ts(ji,jj,jk,jp_tem,Kmm) + ts(ji+1,jj,jk,jp_tem,Kmm) ) 292 292 END_3D 293 CALL lbc_lnk( 'diawri', z2d, 'U', -1. )293 CALL lbc_lnk( 'diawri', z2d, 'U', -1.0_wp ) 294 294 CALL iom_put( "u_heattr", 0.5*rcp * z2d ) ! heat transport in i-direction 295 295 ENDIF … … 300 300 z2d(ji,jj) = z2d(ji,jj) + z3d(ji,jj,jk) * ( ts(ji,jj,jk,jp_sal,Kmm) + ts(ji+1,jj,jk,jp_sal,Kmm) ) 301 301 END_3D 302 CALL lbc_lnk( 'diawri', z2d, 'U', -1. )302 CALL lbc_lnk( 'diawri', z2d, 'U', -1.0_wp ) 303 303 CALL iom_put( "u_salttr", 0.5 * z2d ) ! heat transport in i-direction 304 304 ENDIF … … 318 318 z2d(ji,jj) = z2d(ji,jj) + z3d(ji,jj,jk) * ( ts(ji,jj,jk,jp_tem,Kmm) + ts(ji,jj+1,jk,jp_tem,Kmm) ) 319 319 END_3D 320 CALL lbc_lnk( 'diawri', z2d, 'V', -1. )320 CALL lbc_lnk( 'diawri', z2d, 'V', -1.0_wp ) 321 321 CALL iom_put( "v_heattr", 0.5*rcp * z2d ) ! heat transport in j-direction 322 322 ENDIF … … 327 327 z2d(ji,jj) = z2d(ji,jj) + z3d(ji,jj,jk) * ( ts(ji,jj,jk,jp_sal,Kmm) + ts(ji,jj+1,jk,jp_sal,Kmm) ) 328 328 END_3D 329 CALL lbc_lnk( 'diawri', z2d, 'V', -1. )329 CALL lbc_lnk( 'diawri', z2d, 'V', -1.0_wp ) 330 330 CALL iom_put( "v_salttr", 0.5 * z2d ) ! heat transport in j-direction 331 331 ENDIF … … 336 336 z2d(ji,jj) = z2d(ji,jj) + e3t(ji,jj,jk,Kmm) * ts(ji,jj,jk,jp_tem,Kmm) 337 337 END_3D 338 CALL lbc_lnk( 'diawri', z2d, 'T', -1. )338 CALL lbc_lnk( 'diawri', z2d, 'T', -1.0_wp ) 339 339 CALL iom_put( "tosmint", rho0 * z2d ) ! Vertical integral of temperature 340 340 ENDIF … … 344 344 z2d(ji,jj) = z2d(ji,jj) + e3t(ji,jj,jk,Kmm) * ts(ji,jj,jk,jp_sal,Kmm) 345 345 END_3D 346 CALL lbc_lnk( 'diawri', z2d, 'T', -1. )346 CALL lbc_lnk( 'diawri', z2d, 'T', -1.0_wp ) 347 347 CALL iom_put( "somint", rho0 * z2d ) ! Vertical integral of salinity 348 348 ENDIF
Note: See TracChangeset
for help on using the changeset viewer.