Changeset 14072 for NEMO/trunk/src/OCE/DYN/dynzad.F90
- Timestamp:
- 2020-12-04T08:48:38+01:00 (3 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/trunk/src/OCE/DYN/dynzad.F90
r14007 r14072 7 7 !! NEMO 0.5 ! 2002-07 (G. Madec) Free form, F90 8 8 !!---------------------------------------------------------------------- 9 9 10 10 !!---------------------------------------------------------------------- 11 11 !! dyn_zad : vertical advection momentum trend … … 25 25 IMPLICIT NONE 26 26 PRIVATE 27 27 28 28 PUBLIC dyn_zad ! routine called by dynadv.F90 29 29 … … 41 41 !!---------------------------------------------------------------------- 42 42 !! *** ROUTINE dynzad *** 43 !! 44 !! ** Purpose : Compute the now vertical momentum advection trend and 43 !! 44 !! ** Purpose : Compute the now vertical momentum advection trend and 45 45 !! add it to the general trend of momentum equation. 46 46 !! … … 73 73 74 74 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) 75 ALLOCATE( ztrdu(jpi,jpj,jpk) , ztrdv(jpi,jpj,jpk) ) 76 ztrdu(:,:,:) = puu(:,:,:,Krhs) 77 ztrdv(:,:,:) = pvv(:,:,:,Krhs) 78 78 ENDIF 79 79 80 80 DO jk = 2, jpkm1 ! Vertical momentum advection at level w and u- and v- vertical 81 81 DO_2D( 0, 1, 0, 1 ) ! vertical fluxes … … 111 111 ztrdv(:,:,:) = pvv(:,:,:,Krhs) - ztrdv(:,:,:) 112 112 CALL trd_dyn( ztrdu, ztrdv, jpdyn_zad, kt, Kmm ) 113 DEALLOCATE( ztrdu, ztrdv ) 113 DEALLOCATE( ztrdu, ztrdv ) 114 114 ENDIF 115 115 ! ! Control print
Note: See TracChangeset
for help on using the changeset viewer.