Changeset 14219 for NEMO/branches/2020/dev_r14116_HPC-04_mcastril_Mixed_Precision_implementation_final/src/OCE/DYN/dynzad.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/OCE/DYN/dynzad.F90
r14072 r14219 31 31 # include "do_loop_substitute.h90" 32 32 # include "domzgr_substitute.h90" 33 # include "single_precision_substitute.h90" 33 34 !!---------------------------------------------------------------------- 34 35 !! NEMO/OCE 4.0 , NEMO Consortium (2018) … … 56 57 INTEGER , INTENT( in ) :: kt ! ocean time-step inedx 57 58 INTEGER , INTENT( in ) :: Kmm, Krhs ! ocean time level indices 58 REAL( wp), DIMENSION(jpi,jpj,jpk,jpt), INTENT(inout) :: puu, pvv ! ocean velocities and RHS of momentum equation59 REAL(dp), DIMENSION(jpi,jpj,jpk,jpt), INTENT(inout) :: puu, pvv ! ocean velocities and RHS of momentum equation 59 60 ! 60 61 INTEGER :: ji, jj, jk ! dummy loop indices … … 62 63 REAL(wp), DIMENSION(jpi,jpj) :: zww 63 64 REAL(wp), DIMENSION(jpi,jpj,jpk) :: zwuw, zwvw 64 REAL( wp), DIMENSION(:,:,:), ALLOCATABLE :: ztrdu, ztrdv65 REAL(dp), DIMENSION(:,:,:), ALLOCATABLE :: ztrdu, ztrdv 65 66 !!---------------------------------------------------------------------- 66 67 ! … … 73 74 74 75 IF( l_trddyn ) THEN ! Save puu(:,:,:,Krhs) and pvv(:,:,:,Krhs) trends 75 ALLOCATE( ztrdu(jpi,jpj,jpk) , ztrdv(jpi,jpj,jpk) ) 76 ztrdu(:,:,:) = puu(:,:,:,Krhs) 77 ztrdv(:,:,:) = pvv(:,:,:,Krhs) 76 ALLOCATE( ztrdu(jpi,jpj,jpk) , ztrdv(jpi,jpj,jpk) ) 77 ztrdu(:,:,:) = puu(:,:,:,Krhs) 78 ztrdv(:,:,:) = pvv(:,:,:,Krhs) 78 79 ENDIF 79 80 … … 114 115 ENDIF 115 116 ! ! Control print 116 IF(sn_cfctl%l_prtctl) CALL prt_ctl( tab3d_1= puu(:,:,:,Krhs), clinfo1=' zad - Ua: ', mask1=umask, &117 & tab3d_2= pvv(:,:,:,Krhs), clinfo2= ' Va: ', mask2=vmask, clinfo3='dyn' )117 IF(sn_cfctl%l_prtctl) CALL prt_ctl( tab3d_1=CASTWP(puu(:,:,:,Krhs)), clinfo1=' zad - Ua: ', mask1=umask, & 118 & tab3d_2=CASTWP(pvv(:,:,:,Krhs)), clinfo2= ' Va: ', mask2=vmask, clinfo3='dyn' ) 118 119 ! 119 120 IF( ln_timing ) CALL timing_stop('dyn_zad')
Note: See TracChangeset
for help on using the changeset viewer.