- Timestamp:
- 2019-08-22T11:49:08+02:00 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/UKMO/dev_r10037_GPU/src/OCE/DYN/dynspg.F90
r10843 r11467 71 71 !! period is used to prevent the divergence of odd and even time step. 72 72 !!---------------------------------------------------------------------- 73 USE scoce, ONLY : zpice => scr2D1, ztrdu => scr1, ztrdv => scr2 73 74 INTEGER, INTENT(in) :: kt ! ocean time-step index 74 75 ! 75 76 INTEGER :: ji, jj, jk ! dummy loop indices 76 77 REAL(wp) :: z2dt, zg_2, zintp, zgrau0r, zld ! local scalars 77 REAL(wp), ALLOCATABLE, DIMENSION(:,:) :: zpice78 REAL(wp), ALLOCATABLE, DIMENSION(:,:,:) :: ztrdu, ztrdv79 78 !!---------------------------------------------------------------------- 80 79 ! … … 82 81 ! 83 82 IF( l_trddyn ) THEN ! temporary save of ta and sa trends 84 ALLOCATE( ztrdu(jpi,jpj,jpk) , ztrdv(jpi,jpj,jpk) )85 83 ztrdu(:,:,:) = ua(:,:,:) 86 84 ztrdv(:,:,:) = va(:,:,:) … … 134 132 ! 135 133 IF( ln_ice_embd ) THEN !== embedded sea ice: Pressure gradient due to snow-ice mass ==! 136 ALLOCATE( zpice(jpi,jpj) )137 134 zintp = REAL( MOD( kt-1, nn_fsbc ) ) / REAL( nn_fsbc ) 138 135 zgrau0r = - grav * r1_rau0 … … 144 141 END DO 145 142 END DO 146 DEALLOCATE( zpice )147 143 ENDIF 148 144 ! … … 169 165 ztrdv(:,:,:) = va(:,:,:) - ztrdv(:,:,:) 170 166 CALL trd_dyn( ztrdu, ztrdv, jpdyn_spg, kt ) 171 DEALLOCATE( ztrdu , ztrdv )172 167 ENDIF 173 168 ! ! print mean trends (used for debugging)
Note: See TracChangeset
for help on using the changeset viewer.