| 1 | MODULE dynnxt_c1d |
|---|
| 2 | !!====================================================================== |
|---|
| 3 | !! *** MODULE dynnxt_c1d *** |
|---|
| 4 | !! Ocean dynamics: time stepping in 1D configuration |
|---|
| 5 | !!====================================================================== |
|---|
| 6 | !! History : 2.0 ! 2004-10 (C. Ethe) Original code from dynnxt.F90 |
|---|
| 7 | !! 3.0 ! 2008-04 (G.madec) Style only |
|---|
| 8 | !!---------------------------------------------------------------------- |
|---|
| 9 | #if defined key_c1d |
|---|
| 10 | !!---------------------------------------------------------------------- |
|---|
| 11 | !! 'key_c1d' 1D Configuration |
|---|
| 12 | !!---------------------------------------------------------------------- |
|---|
| 13 | !!---------------------------------------------------------------------- |
|---|
| 14 | !! dyn_nxt_c1d : update the horizontal velocity from the momentum trend |
|---|
| 15 | !!---------------------------------------------------------------------- |
|---|
| 16 | USE oce ! ocean dynamics and tracers |
|---|
| 17 | USE dom_oce ! ocean space and time domain |
|---|
| 18 | USE in_out_manager ! I/O manager |
|---|
| 19 | USE lbclnk ! lateral boundary condition (or mpp link) |
|---|
| 20 | USE prtctl ! Print control |
|---|
| 21 | |
|---|
| 22 | IMPLICIT NONE |
|---|
| 23 | PRIVATE |
|---|
| 24 | |
|---|
| 25 | PUBLIC dyn_nxt_c1d ! routine called by step.F90 |
|---|
| 26 | !!---------------------------------------------------------------------- |
|---|
| 27 | !! NEMO/C1D 3.3 , NEMO Consortium (2010) |
|---|
| 28 | !! $Id$ |
|---|
| 29 | !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) |
|---|
| 30 | !!---------------------------------------------------------------------- |
|---|
| 31 | |
|---|
| 32 | CONTAINS |
|---|
| 33 | |
|---|
| 34 | SUBROUTINE dyn_nxt_c1d ( kt ) |
|---|
| 35 | !!---------------------------------------------------------------------- |
|---|
| 36 | !! *** ROUTINE dyn_nxt_c1d *** |
|---|
| 37 | !! |
|---|
| 38 | !! ** Purpose : Compute the after horizontal velocity from the |
|---|
| 39 | !! momentum trend. |
|---|
| 40 | !! |
|---|
| 41 | !! ** Method : Apply lateral boundary conditions on the trends (ua,va) |
|---|
| 42 | !! through calls to routine lbc_lnk. |
|---|
| 43 | !! After velocity is compute using a leap-frog scheme environment: |
|---|
| 44 | !! (ua,va) = (ub,vb) + 2 rdt (ua,va) |
|---|
| 45 | !! Time filter applied on now horizontal velocity to avoid the |
|---|
| 46 | !! divergence of two consecutive time-steps and swap of dynamics |
|---|
| 47 | !! arrays to start the next time step: |
|---|
| 48 | !! (ub,vb) = (un,vn) + atfp [ (ub,vb) + (ua,va) - 2 (un,vn) ] |
|---|
| 49 | !! (un,vn) = (ua,va) |
|---|
| 50 | !! |
|---|
| 51 | !! ** Action : - Update ub,vb arrays, the before horizontal velocity |
|---|
| 52 | !! - Update un,vn arrays, the now horizontal velocity |
|---|
| 53 | !!---------------------------------------------------------------------- |
|---|
| 54 | INTEGER, INTENT( in ) :: kt ! ocean time-step index |
|---|
| 55 | !! |
|---|
| 56 | INTEGER :: jk ! dummy loop indices |
|---|
| 57 | REAL(wp) :: z2dt ! temporary scalar |
|---|
| 58 | !!---------------------------------------------------------------------- |
|---|
| 59 | |
|---|
| 60 | IF( kt == nit000 ) THEN |
|---|
| 61 | IF(lwp) WRITE(numout,*) |
|---|
| 62 | IF(lwp) WRITE(numout,*) 'dyn_nxt_c1d : time stepping on 1D configuation' |
|---|
| 63 | IF(lwp) WRITE(numout,*) '~~~~~~~' |
|---|
| 64 | ENDIF |
|---|
| 65 | |
|---|
| 66 | ! Local constant initialization |
|---|
| 67 | z2dt = 2. * rdt |
|---|
| 68 | IF( neuler == 0 .AND. kt == nit000 ) z2dt = rdt |
|---|
| 69 | |
|---|
| 70 | ! Lateral boundary conditions on ( ua, va ) |
|---|
| 71 | CALL lbc_lnk( ua, 'U', -1. ) |
|---|
| 72 | CALL lbc_lnk( va, 'V', -1. ) |
|---|
| 73 | |
|---|
| 74 | DO jk = 1, jpkm1 ! Next Velocity |
|---|
| 75 | ua(:,:,jk) = ( ub(:,:,jk) + z2dt * ua(:,:,jk) ) * umask(:,:,jk) |
|---|
| 76 | va(:,:,jk) = ( vb(:,:,jk) + z2dt * va(:,:,jk) ) * vmask(:,:,jk) |
|---|
| 77 | END DO |
|---|
| 78 | |
|---|
| 79 | DO jk = 1, jpkm1 ! Time filter and swap of dynamics arrays |
|---|
| 80 | IF( neuler == 0 .AND. kt == nit000 ) THEN ! Euler (forward) time stepping |
|---|
| 81 | ub(:,:,jk) = un(:,:,jk) |
|---|
| 82 | vb(:,:,jk) = vn(:,:,jk) |
|---|
| 83 | un(:,:,jk) = ua(:,:,jk) |
|---|
| 84 | vn(:,:,jk) = va(:,:,jk) |
|---|
| 85 | ELSE ! Leap-frog time stepping |
|---|
| 86 | ub(:,:,jk) = atfp * ( ub(:,:,jk) + ua(:,:,jk) ) + atfp1 * un(:,:,jk) |
|---|
| 87 | vb(:,:,jk) = atfp * ( vb(:,:,jk) + va(:,:,jk) ) + atfp1 * vn(:,:,jk) |
|---|
| 88 | un(:,:,jk) = ua(:,:,jk) |
|---|
| 89 | vn(:,:,jk) = va(:,:,jk) |
|---|
| 90 | ENDIF |
|---|
| 91 | END DO |
|---|
| 92 | |
|---|
| 93 | IF(ln_ctl) CALL prt_ctl( tab3d_1=un, clinfo1=' nxt_c1d - Un: ', mask1=umask, & |
|---|
| 94 | & tab3d_2=vn, clinfo2=' Vn: ' , mask2=vmask ) |
|---|
| 95 | ! |
|---|
| 96 | END SUBROUTINE dyn_nxt_c1d |
|---|
| 97 | |
|---|
| 98 | #else |
|---|
| 99 | !!---------------------------------------------------------------------- |
|---|
| 100 | !! Default key NO 1D Config |
|---|
| 101 | !!---------------------------------------------------------------------- |
|---|
| 102 | CONTAINS |
|---|
| 103 | SUBROUTINE dyn_nxt_c1d ( kt ) |
|---|
| 104 | WRITE(*,*) 'dyn_nxt_c1d: You should not have seen this print! error?', kt |
|---|
| 105 | END SUBROUTINE dyn_nxt_c1d |
|---|
| 106 | #endif |
|---|
| 107 | |
|---|
| 108 | !!====================================================================== |
|---|
| 109 | END MODULE dynnxt_c1d |
|---|
