Changeset 10874 for NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps/src/OCE/DYN/dynldf.F90
- Timestamp:
- 2019-04-15T15:57:37+02:00 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps/src/OCE/DYN/dynldf.F90
r10806 r10874 43 43 CONTAINS 44 44 45 SUBROUTINE dyn_ldf( kt , ktlev1, ktlev2, pu_rhs, pv_rhs)45 SUBROUTINE dyn_ldf( kt ) 46 46 !!---------------------------------------------------------------------- 47 47 !! *** ROUTINE dyn_ldf *** … … 49 49 !! ** Purpose : compute the lateral ocean dynamics physics. 50 50 !!---------------------------------------------------------------------- 51 INTEGER, INTENT(in) :: kt ! ocean time-step index 52 INTEGER, INTENT(in) :: ktlev1, ktlev2 ! time level index for source terms 53 REAL(wp), INTENT(inout), DIMENSION(jpi,jpj,jpk) :: pu_rhs, pv_rhs ! momentum trends 51 INTEGER, INTENT(in) :: kt ! ocean time-step index 54 52 ! 55 53 REAL(wp), ALLOCATABLE, DIMENSION(:,:,:) :: ztrdu, ztrdv … … 60 58 IF( l_trddyn ) THEN ! temporary save of momentum trends 61 59 ALLOCATE( ztrdu(jpi,jpj,jpk) , ztrdv(jpi,jpj,jpk) ) 62 ztrdu(:,:,:) = pu_rhs(:,:,:)63 ztrdv(:,:,:) = pv_rhs(:,:,:)60 ztrdu(:,:,:) = ua(:,:,:) 61 ztrdv(:,:,:) = va(:,:,:) 64 62 ENDIF 65 63 66 64 SELECT CASE ( nldf_dyn ) ! compute lateral mixing trend and add it to the general trend 67 65 ! 68 CASE ( np_lap ) ; CALL dyn_ldf_lap( kt, ktlev1, ktlev2, uu(:,:,:,ktlev1), vv(:,:,:,ktlev1), pu_rhs, pv_rhs, 1 ) ! iso-level laplacian66 CASE ( np_lap ) ; CALL dyn_ldf_lap( kt, ub, vb, ua, va, 1 ) ! iso-level laplacian 69 67 CASE ( np_lap_i ) ; CALL dyn_ldf_iso( kt ) ! rotated laplacian 70 CASE ( np_blp ) ; CALL dyn_ldf_blp( kt, ktlev1, ktlev2, uu(:,:,:,ktlev1), vv(:,:,:,ktlev1), pu_rhs, pv_rhs) ! iso-level bi-laplacian68 CASE ( np_blp ) ; CALL dyn_ldf_blp( kt, ub, vb, ua, va ) ! iso-level bi-laplacian 71 69 ! 72 70 END SELECT 73 71 74 72 IF( l_trddyn ) THEN ! save the horizontal diffusive trends for further diagnostics 75 ztrdu(:,:,:) = pu_rhs(:,:,:) - ztrdu(:,:,:)76 ztrdv(:,:,:) = pv_rhs(:,:,:) - ztrdv(:,:,:)73 ztrdu(:,:,:) = ua(:,:,:) - ztrdu(:,:,:) 74 ztrdv(:,:,:) = va(:,:,:) - ztrdv(:,:,:) 77 75 CALL trd_dyn( ztrdu, ztrdv, jpdyn_ldf, kt ) 78 76 DEALLOCATE ( ztrdu , ztrdv )
Note: See TracChangeset
for help on using the changeset viewer.