Changeset 7753 for trunk/NEMOGCM/NEMO/OPA_SRC/DYN/dynzdf.F90
- Timestamp:
- 2017-03-03T12:46:59+01:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMOGCM/NEMO/OPA_SRC/DYN/dynzdf.F90
r7698 r7753 53 53 !! 54 54 INTEGER, INTENT( in ) :: kt ! ocean time-step index 55 INTEGER :: ji, jj, jk ! dummy loop indices56 55 ! 57 56 REAL(wp), POINTER, DIMENSION(:,:,:) :: ztrdu, ztrdv … … 67 66 IF( l_trddyn ) THEN ! temporary save of ta and sa trends 68 67 CALL wrk_alloc( jpi, jpj, jpk, ztrdu, ztrdv ) 69 !$OMP PARALLEL DO schedule(static) private(jk,jj,ji) 70 DO jk = 1, jpk 71 DO jj = 1, jpj 72 DO ji = 1, jpi 73 ztrdu(ji,jj,jk) = ua(ji,jj,jk) 74 ztrdv(ji,jj,jk) = va(ji,jj,jk) 75 END DO 76 END DO 77 END DO 68 ztrdu(:,:,:) = ua(:,:,:) 69 ztrdv(:,:,:) = va(:,:,:) 78 70 ENDIF 79 71 … … 86 78 87 79 IF( l_trddyn ) THEN ! save the vertical diffusive trends for further diagnostics 88 !$OMP PARALLEL DO schedule(static) private(jk,jj,ji) 89 DO jk = 1, jpk 90 DO jj = 1, jpj 91 DO ji = 1, jpi 92 ztrdu(ji,jj,jk) = ( ua(ji,jj,jk) - ub(ji,jj,jk) ) / r2dt - ztrdu(ji,jj,jk) 93 ztrdv(ji,jj,jk) = ( va(ji,jj,jk) - vb(ji,jj,jk) ) / r2dt - ztrdv(ji,jj,jk) 94 END DO 95 END DO 96 END DO 80 ztrdu(:,:,:) = ( ua(:,:,:) - ub(:,:,:) ) / r2dt - ztrdu(:,:,:) 81 ztrdv(:,:,:) = ( va(:,:,:) - vb(:,:,:) ) / r2dt - ztrdv(:,:,:) 97 82 CALL trd_dyn( ztrdu, ztrdv, jpdyn_zdf, kt ) 98 83 CALL wrk_dealloc( jpi, jpj, jpk, ztrdu, ztrdv )
Note: See TracChangeset
for help on using the changeset viewer.