Changeset 14381
- Timestamp:
- 2021-02-03T13:36:25+01:00 (3 years ago)
- Location:
- NEMO/branches/2021/dev_r14318_RK3_stage1/src/OCE
- Files:
-
- 6 edited
-
DYN/dynadv.F90 (modified) (3 diffs)
-
DYN/dynadv_cen2.F90 (modified) (4 diffs)
-
DYN/dynadv_ubs.F90 (modified) (6 diffs)
-
TRA/traadv.F90 (modified) (4 diffs)
-
oce.F90 (modified) (1 diff)
-
stpmlf.F90 (modified) (5 diffs)
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2021/dev_r14318_RK3_stage1/src/OCE/DYN/dynadv.F90
r14053 r14381 51 51 CONTAINS 52 52 53 SUBROUTINE dyn_adv( kt, Kbb, Kmm, puu, pvv, Krhs )53 SUBROUTINE dyn_adv( kt, Kbb, Kmm, puu, pvv, Krhs, pau, pav, paw ) 54 54 !!--------------------------------------------------------------------- 55 55 !! *** ROUTINE dyn_adv *** … … 64 64 !! (see dynvor module). 65 65 !!---------------------------------------------------------------------- 66 INTEGER , INTENT( in ) :: kt ! ocean time-step index 67 INTEGER , INTENT( in ) :: Kbb, Kmm, Krhs ! ocean time level indices 68 REAL(wp), DIMENSION(jpi,jpj,jpk,jpt), INTENT(inout) :: puu, pvv ! ocean velocities and RHS of momentum equation 66 INTEGER , INTENT(in ) :: kt ! ocean time-step index 67 INTEGER , INTENT(in ) :: Kbb, Kmm, Krhs ! ocean time level indices 68 REAL(wp), DIMENSION(:,:,:) , OPTIONAL, TARGET, INTENT(in ) :: pau, pav, paw ! advective velocity 69 REAL(wp), DIMENSION(jpi,jpj,jpk,jpt) , TARGET, INTENT(inout) :: puu, pvv ! ocean velocities and RHS of momentum Eq. 69 70 !!---------------------------------------------------------------------- 70 71 ! … … 73 74 SELECT CASE( n_dynadv ) !== compute advection trend and add it to general trend ==! 74 75 CASE( np_VEC_c2 ) 75 CALL dyn_keg ( kt, nn_dynkeg , Kmm, puu, pvv, Krhs )! vector form : horizontal gradient of kinetic energy76 CALL dyn_zad ( kt , Kmm, puu, pvv, Krhs )! vector form : vertical advection76 CALL dyn_keg ( kt, nn_dynkeg , Kmm, puu, pvv, Krhs ) ! vector form : horizontal gradient of kinetic energy 77 CALL dyn_zad ( kt , Kmm, puu, pvv, Krhs ) ! vector form : vertical advection 77 78 CASE( np_FLX_c2 ) 78 CALL dyn_adv_cen2( kt , Kmm, puu, pvv, Krhs )! 2nd order centered scheme79 CALL dyn_adv_cen2( kt , Kmm, puu, pvv, Krhs, pau, pav, paw ) ! 2nd order centered scheme 79 80 CASE( np_FLX_ubs ) 80 CALL dyn_adv_ubs ( kt , Kbb, Kmm, puu, pvv, Krhs )! 3rd order UBS scheme (UP3)81 CALL dyn_adv_ubs ( kt , Kbb, Kmm, puu, pvv, Krhs, pau, pav, paw ) ! 3rd order UBS scheme (UP3) 81 82 END SELECT 82 83 ! -
NEMO/branches/2021/dev_r14318_RK3_stage1/src/OCE/DYN/dynadv_cen2.F90
r13497 r14381 36 36 CONTAINS 37 37 38 SUBROUTINE dyn_adv_cen2( kt, Kmm, puu, pvv, Krhs )38 SUBROUTINE dyn_adv_cen2( kt, Kmm, puu, pvv, Krhs, pau, pav, paw ) 39 39 !!---------------------------------------------------------------------- 40 40 !! *** ROUTINE dyn_adv_cen2 *** … … 47 47 !! ** Action : (puu(:,:,:,Krhs),pvv(:,:,:,Krhs)) updated with the now vorticity term trend 48 48 !!---------------------------------------------------------------------- 49 INTEGER , INTENT( in ) :: kt ! ocean time-step index 50 INTEGER , INTENT( in ) :: Kmm, Krhs ! ocean time level indices 51 REAL(wp), DIMENSION(jpi,jpj,jpk,jpt), INTENT(inout) :: puu, pvv ! ocean velocities and RHS of momentum equation 49 INTEGER , INTENT(in ) :: kt ! ocean time-step index 50 INTEGER , INTENT(in ) :: Kmm, Krhs ! ocean time level indices 51 REAL(wp), DIMENSION(jpi,jpj,jpk,jpt) , TARGET, INTENT(inout) :: puu, pvv ! ocean velocities and RHS of momentum equation 52 REAL(wp), DIMENSION(:,:,:) , OPTIONAL, TARGET, INTENT(in ) :: pau, pav, paw ! advective velocity 52 53 ! 53 54 INTEGER :: ji, jj, jk ! dummy loop indices 54 REAL(wp), DIMENSION(jpi,jpj,jpk) :: zfu_t, zfu_f, zfu_uw, zfu 55 REAL(wp), DIMENSION(jpi,jpj,jpk) :: zfv_t, zfv_f, zfv_vw, zfv, zfw 55 REAL(wp), DIMENSION(jpi,jpj,jpk) :: zfu_t, zfu_f, zfu_uw, zfu 56 REAL(wp), DIMENSION(jpi,jpj,jpk) :: zfv_t, zfv_f, zfv_vw, zfv, zfw 57 REAL(wp), DIMENSION(:,:,:), POINTER :: zptu, zptv, zptw 56 58 !!---------------------------------------------------------------------- 57 59 ! … … 67 69 ENDIF 68 70 ! 71 IF( PRESENT( pau ) ) THEN ! RK3: advective velocity (pau,pav,paw) /= advected velocity (puu,pvv,ww) 72 zptu => pau(:,:,:) 73 zptv => pav(:,:,:) 74 zptw => paw(:,:,:) 75 ELSE ! MLF: advective velocity = (puu,pvv,ww) 76 zptu => puu(:,:,:,Kmm) 77 zptv => pvv(:,:,:,Kmm) 78 zptw => ww (:,:,: ) 79 ENDIF 80 ! 69 81 ! !== Horizontal advection ==! 70 82 ! 71 83 DO jk = 1, jpkm1 ! horizontal transport 72 zfu(:,:,jk) = 0.25_wp * e2u(:,:) * e3u(:,:,jk,Kmm) * puu(:,:,jk,Kmm)73 zfv(:,:,jk) = 0.25_wp * e1v(:,:) * e3v(:,:,jk,Kmm) * pvv(:,:,jk,Kmm)84 zfu(:,:,jk) = 0.25_wp * e2u(:,:) * e3u(:,:,jk,Kmm) * zptu(:,:,jk) 85 zfv(:,:,jk) = 0.25_wp * e1v(:,:) * e3v(:,:,jk,Kmm) * zptv(:,:,jk) 74 86 DO_2D( 1, 0, 1, 0 ) ! horizontal momentum fluxes (at T- and F-point) 75 87 zfu_t(ji+1,jj ,jk) = ( zfu(ji,jj,jk) + zfu(ji+1,jj,jk) ) * ( puu(ji,jj,jk,Kmm) + puu(ji+1,jj ,jk,Kmm) ) … … 104 116 IF( ln_linssh ) THEN ! linear free surface: advection through the surface 105 117 DO_2D( 0, 0, 0, 0 ) 106 zfu_uw(ji,jj,1) = 0.5_wp * ( e1e2t(ji,jj) * ww(ji,jj,1) + e1e2t(ji+1,jj) * ww(ji+1,jj,1) ) * puu(ji,jj,1,Kmm)107 zfv_vw(ji,jj,1) = 0.5_wp * ( e1e2t(ji,jj) * ww(ji,jj,1) + e1e2t(ji,jj+1) * ww(ji,jj+1,1) ) * pvv(ji,jj,1,Kmm)118 zfu_uw(ji,jj,1) = 0.5_wp * ( e1e2t(ji,jj) * zptw(ji,jj,1) + e1e2t(ji+1,jj) * zptw(ji+1,jj,1) ) * puu(ji,jj,1,Kmm) 119 zfv_vw(ji,jj,1) = 0.5_wp * ( e1e2t(ji,jj) * zptw(ji,jj,1) + e1e2t(ji,jj+1) * zptw(ji,jj+1,1) ) * pvv(ji,jj,1,Kmm) 108 120 END_2D 109 121 ENDIF 110 122 DO jk = 2, jpkm1 ! interior advective fluxes 111 123 DO_2D( 0, 1, 0, 1 ) ! 1/4 * Vertical transport 112 zfw(ji,jj,jk) = 0.25_wp * e1e2t(ji,jj) * ww(ji,jj,jk)124 zfw(ji,jj,jk) = 0.25_wp * e1e2t(ji,jj) * zptw(ji,jj,jk) 113 125 END_2D 114 126 DO_2D( 0, 0, 0, 0 ) -
NEMO/branches/2021/dev_r14318_RK3_stage1/src/OCE/DYN/dynadv_ubs.F90
r13497 r14381 42 42 CONTAINS 43 43 44 SUBROUTINE dyn_adv_ubs( kt, Kbb, Kmm, puu, pvv, Krhs )44 SUBROUTINE dyn_adv_ubs( kt, Kbb, Kmm, puu, pvv, Krhs, pau, pav, paw ) 45 45 !!---------------------------------------------------------------------- 46 46 !! *** ROUTINE dyn_adv_ubs *** … … 69 69 !! Reference : Shchepetkin & McWilliams, 2005, Ocean Modelling. 70 70 !!---------------------------------------------------------------------- 71 INTEGER , INTENT( in ) :: kt ! ocean time-step index 72 INTEGER , INTENT( in ) :: Kbb, Kmm, Krhs ! ocean time level indices 73 REAL(wp), DIMENSION(jpi,jpj,jpk,jpt), INTENT(inout) :: puu, pvv ! ocean velocities and RHS of momentum equation 71 INTEGER , INTENT(in ) :: kt ! ocean time-step index 72 INTEGER , INTENT(in ) :: Kbb, Kmm, Krhs ! ocean time level indices 73 REAL(wp), DIMENSION(jpi,jpj,jpk,jpt) , TARGET, INTENT(inout) :: puu, pvv ! ocean velocities and RHS of momentum equation 74 REAL(wp), DIMENSION(:,:,:) , OPTIONAL, TARGET, INTENT(in ) :: pau, pav, paw ! advective velocity 74 75 ! 75 76 INTEGER :: ji, jj, jk ! dummy loop indices … … 79 80 REAL(wp), DIMENSION(jpi,jpj,jpk,2) :: zlu_uu, zlu_uv 80 81 REAL(wp), DIMENSION(jpi,jpj,jpk,2) :: zlv_vv, zlv_vu 82 REAL(wp), DIMENSION(:,:,:), POINTER :: zptu, zptv, zptw 81 83 !!---------------------------------------------------------------------- 82 84 ! … … 101 103 zfv_vw(:,:,:) = pvv(:,:,:,Krhs) 102 104 ENDIF 105 ! 106 IF( PRESENT( pau ) ) THEN ! RK3: advective velocity (pau,pav,paw) /= advected velocity (puu,pvv,ww) 107 zptu => pau(:,:,:) 108 zptv => pav(:,:,:) 109 zptw => paw(:,:,:) 110 ELSE ! MLF: advective velocity = (puu,pvv,ww) 111 zptu => puu(:,:,:,Kmm) 112 zptv => pvv(:,:,:,Kmm) 113 zptw => ww (:,:,: ) 114 ENDIF 115 ! 103 116 ! ! =========================== ! 104 117 DO jk = 1, jpkm1 ! Laplacian of the velocity ! 105 118 ! ! =========================== ! 106 119 ! ! horizontal volume fluxes 107 zfu(:,:,jk) = e2u(:,:) * e3u(:,:,jk,Kmm) * puu(:,:,jk,Kmm)108 zfv(:,:,jk) = e1v(:,:) * e3v(:,:,jk,Kmm) * pvv(:,:,jk,Kmm)120 zfu(:,:,jk) = e2u(:,:) * e3u(:,:,jk,Kmm) * zptu(:,:,jk) 121 zfv(:,:,jk) = e1v(:,:) * e3v(:,:,jk,Kmm) * zptv(:,:,jk) 109 122 ! 110 123 DO_2D( 0, 0, 0, 0 ) ! laplacian … … 133 146 DO jk = 1, jpkm1 ! ====================== ! 134 147 ! ! horizontal volume fluxes 135 zfu(:,:,jk) = 0.25_wp * e2u(:,:) * e3u(:,:,jk,Kmm) * puu(:,:,jk,Kmm)136 zfv(:,:,jk) = 0.25_wp * e1v(:,:) * e3v(:,:,jk,Kmm) * pvv(:,:,jk,Kmm)148 zfu(:,:,jk) = 0.25_wp * e2u(:,:) * e3u(:,:,jk,Kmm) * zptu(:,:,jk) 149 zfv(:,:,jk) = 0.25_wp * e1v(:,:) * e3v(:,:,jk,Kmm) * zptv(:,:,jk) 137 150 ! 138 151 DO_2D( 1, 0, 1, 0 ) ! horizontal momentum fluxes at T- and F-point … … 195 208 IF( ln_linssh ) THEN ! constant volume : advection through the surface 196 209 DO_2D( 0, 0, 0, 0 ) 197 zfu_uw(ji,jj,1) = 0.5_wp * ( e1e2t(ji,jj) * ww(ji,jj,1) + e1e2t(ji+1,jj) * ww(ji+1,jj,1) ) * puu(ji,jj,1,Kmm)198 zfv_vw(ji,jj,1) = 0.5_wp * ( e1e2t(ji,jj) * ww(ji,jj,1) + e1e2t(ji,jj+1) * ww(ji,jj+1,1) ) * pvv(ji,jj,1,Kmm)210 zfu_uw(ji,jj,1) = 0.5_wp * ( e1e2t(ji,jj) * zptw(ji,jj,1) + e1e2t(ji+1,jj) * zptw(ji+1,jj,1) ) * puu(ji,jj,1,Kmm) 211 zfv_vw(ji,jj,1) = 0.5_wp * ( e1e2t(ji,jj) * zptw(ji,jj,1) + e1e2t(ji,jj+1) * zptw(ji,jj+1,1) ) * pvv(ji,jj,1,Kmm) 199 212 END_2D 200 213 ENDIF 201 214 DO jk = 2, jpkm1 ! interior fluxes 202 215 DO_2D( 0, 1, 0, 1 ) 203 zfw(ji,jj,jk) = 0.25_wp * e1e2t(ji,jj) * ww(ji,jj,jk)216 zfw(ji,jj,jk) = 0.25_wp * e1e2t(ji,jj) * zptw(ji,jj,jk) 204 217 END_2D 205 218 DO_2D( 0, 0, 0, 0 ) -
NEMO/branches/2021/dev_r14318_RK3_stage1/src/OCE/TRA/traadv.F90
r14189 r14381 80 80 CONTAINS 81 81 82 SUBROUTINE tra_adv( kt, Kbb, Kmm, pts, Krhs )82 SUBROUTINE tra_adv( kt, Kbb, Kmm, pts, Krhs, pau, pav, paw ) 83 83 !!---------------------------------------------------------------------- 84 84 !! *** ROUTINE tra_adv *** … … 88 88 !! ** Method : - Update (uu(:,:,:,Krhs),vv(:,:,:,Krhs)) with the advection term following nadv 89 89 !!---------------------------------------------------------------------- 90 INTEGER , INTENT(in) :: kt ! ocean time-step index 91 INTEGER , INTENT(in) :: Kbb, Kmm, Krhs ! time level indices 92 REAL(wp), DIMENSION(jpi,jpj,jpk,jpts,jpt), INTENT(inout) :: pts ! active tracers and RHS of tracer equation 90 INTEGER , INTENT(in ) :: kt ! ocean time-step index 91 INTEGER , INTENT(in ) :: Kbb, Kmm, Krhs ! time level indices 92 REAL(wp), DIMENSION(:,:,:), OPTIONAL, TARGET, INTENT(in ) :: pau, pav, paw ! advective velocity 93 REAL(wp), DIMENSION(jpi,jpj,jpk,jpts,jpt) , INTENT(inout) :: pts ! active tracers and RHS of tracer equation 93 94 ! 94 95 INTEGER :: ji, jj, jk ! dummy loop index 96 REAL(wp), DIMENSION(:,:,:), POINTER :: zptu, zptv, zptw 95 97 ! TEMP: [tiling] This change not necessary and can be A2D(nn_hls) if using XIOS (subdomain support) 96 REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, SAVE :: zuu, zvv, zww ! 3D workspace97 REAL(wp), DIMENSION(:,:,:), ALLOCATABLE ::ztrdt, ztrds98 REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, SAVE :: zuu, zvv, zww ! 3D workspace 99 REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: ztrdt, ztrds 98 100 ! TEMP: [tiling] This change not necessary after extra haloes development 99 LOGICAL :: lskip101 LOGICAL :: lskip 100 102 !!---------------------------------------------------------------------- 101 103 ! … … 119 121 ENDIF 120 122 ENDIF 123 ! 121 124 IF( .NOT. lskip ) THEN 122 ! !== effective transport ==! 125 ! !== effective advective transport ==! 126 ! 127 IF( PRESENT( pau ) ) THEN ! RK3: advective velocity (pau,pav,paw) /= advected velocity (uu,vv,ww) 128 zptu => pau(:,:,:) 129 zptv => pav(:,:,:) 130 zptw => paw(:,:,:) 131 ELSE ! MLF: advective velocity = (uu,vv,ww) 132 zptu => uu(:,:,:,Kmm) 133 zptv => vv(:,:,:,Kmm) 134 zptw => ww(:,:,: ) 135 ENDIF 136 ! 123 137 IF( ln_wave .AND. ln_sdw ) THEN 124 138 DO_3D( nn_hls, nn_hls, nn_hls, nn_hls, 1, jpkm1 ) 125 zuu(ji,jj,jk) = e2u (ji,jj) * e3u(ji,jj,jk,Kmm) * ( uu(ji,jj,jk,Kmm) + usd(ji,jj,jk) )126 zvv(ji,jj,jk) = e1v (ji,jj) * e3v(ji,jj,jk,Kmm) * ( vv(ji,jj,jk,Kmm) + vsd(ji,jj,jk) )127 zww(ji,jj,jk) = e1e2t(ji,jj) * ( ww(ji,jj,jk)+ wsd(ji,jj,jk) )139 zuu(ji,jj,jk) = e2u (ji,jj) * e3u(ji,jj,jk,Kmm) * ( zptu(ji,jj,jk) + usd(ji,jj,jk) ) 140 zvv(ji,jj,jk) = e1v (ji,jj) * e3v(ji,jj,jk,Kmm) * ( zptv(ji,jj,jk) + vsd(ji,jj,jk) ) 141 zww(ji,jj,jk) = e1e2t(ji,jj) * ( zptw(ji,jj,jk) + wsd(ji,jj,jk) ) 128 142 END_3D 129 143 ELSE 130 144 DO_3D( nn_hls, nn_hls, nn_hls, nn_hls, 1, jpkm1 ) 131 zuu(ji,jj,jk) = e2u (ji,jj) * e3u(ji,jj,jk,Kmm) * uu(ji,jj,jk,Kmm) ! eulerian transport only132 zvv(ji,jj,jk) = e1v (ji,jj) * e3v(ji,jj,jk,Kmm) * vv(ji,jj,jk,Kmm)133 zww(ji,jj,jk) = e1e2t(ji,jj) * ww(ji,jj,jk)145 zuu(ji,jj,jk) = e2u (ji,jj) * e3u(ji,jj,jk,Kmm) * zptu(ji,jj,jk) ! eulerian transport only 146 zvv(ji,jj,jk) = e1v (ji,jj) * e3v(ji,jj,jk,Kmm) * zptv(ji,jj,jk) 147 zww(ji,jj,jk) = e1e2t(ji,jj) * zptw(ji,jj,jk) 134 148 END_3D 135 149 ENDIF … … 275 289 WRITE(numout,*) ' No advection on T & S ln_traadv_OFF = ', ln_traadv_OFF 276 290 WRITE(numout,*) ' centered scheme ln_traadv_cen = ', ln_traadv_cen 277 WRITE(numout,*) ' horizontal 2nd/4th order nn_cen_h = ', nn_ fct_h278 WRITE(numout,*) ' vertical 2nd/4th order nn_cen_v = ', nn_ fct_v291 WRITE(numout,*) ' horizontal 2nd/4th order nn_cen_h = ', nn_cen_h 292 WRITE(numout,*) ' vertical 2nd/4th order nn_cen_v = ', nn_cen_v 279 293 WRITE(numout,*) ' Flux Corrected Transport scheme ln_traadv_fct = ', ln_traadv_fct 280 294 WRITE(numout,*) ' horizontal 2nd/4th order nn_fct_h = ', nn_fct_h -
NEMO/branches/2021/dev_r14318_RK3_stage1/src/OCE/oce.F90
r14064 r14381 21 21 !! dynamics and tracer fields 22 22 !! -------------------------- 23 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:,:) :: uu , vv !: horizontal velocities [m/s]24 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: ww !: vertical velocity [m/s]23 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:,:), TARGET :: uu , vv !: horizontal velocities [m/s] 24 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) , TARGET :: ww !: vertical velocity [m/s] 25 25 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: wi !: vertical vel. (adaptive-implicit) [m/s] 26 26 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: hdiv !: horizontal divergence [s-1] -
NEMO/branches/2021/dev_r14318_RK3_stage1/src/OCE/stpmlf.F90
r14239 r14381 94 94 !!---------------------------------------------------------------------- 95 95 INTEGER :: ji, jj, jk, jtile ! dummy loop indice 96 REAL(wp), ALLOCATABLE, DIMENSION(:,:,:) :: zgdept 96 REAL(wp), ALLOCATABLE, DIMENSION(:,:,:) :: zgdept 97 REAL(wp), TARGET , DIMENSION(jpi,jpj,jpk) :: zau, zav, zaw ! advective velocity 97 98 !! --------------------------------------------------------------------- 98 99 #if defined key_agrif … … 222 223 uu(:,:,:,Nrhs) = 0._wp ! set dynamics trends to zero 223 224 vv(:,:,:,Nrhs) = 0._wp 224 225 !!st 226 zau(:,:,:) = uu(:,:,:,Nnn) !!st patch for MLF will be computed in RK3 227 zav(:,:,:) = vv(:,:,:,Nnn) 228 zaw(:,:,:) = ww(:,:,: ) 229 !!st 225 230 IF( lk_asminc .AND. ln_asmiau .AND. ln_dyninc ) & 226 231 & CALL dyn_asm_inc ( kstp, Nbb, Nnn, uu, vv, Nrhs ) ! apply dynamics assimilation increment … … 230 235 & CALL Agrif_Sponge_dyn ! momentum sponge 231 236 #endif 232 CALL dyn_adv( kstp, Nbb, Nnn , uu, vv, Nrhs ) ! advection (VF or FF) ==> RHS237 CALL dyn_adv( kstp, Nbb, Nnn , uu, vv, Nrhs, zau, zav, zaw ) ! advection (VF or FF) ==> RHS 233 238 CALL dyn_vor( kstp, Nnn , uu, vv, Nrhs ) ! vorticity ==> RHS 234 239 CALL dyn_ldf( kstp, Nbb, Nnn , uu, vv, Nrhs ) ! lateral mixing … … 247 252 IF( ln_zad_Aimp ) CALL wAimp ( kstp, Nnn ) ! Adaptive-implicit vertical advection partitioning 248 253 ENDIF 249 250 254 !!st 255 zau(:,:,:) = uu(:,:,:,Nnn) !!st patch for MLF will be computed in RK3 256 zav(:,:,:) = vv(:,:,:,Nnn) 257 zaw(:,:,:) = ww(:,:,: ) 258 !!st 251 259 !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 252 260 ! cool skin … … 313 321 IF( ln_tile ) CALL dom_tile( ntsi, ntsj, ntei, ntej, ktile = jtile ) 314 322 315 CALL tra_adv ( kstp, Nbb, Nnn, ts, Nrhs ) ! hor. + vert. advection ==> RHS323 CALL tra_adv ( kstp, Nbb, Nnn, ts, Nrhs, zau, zav, zaw ) ! hor. + vert. advection ==> RHS 316 324 IF( ln_zdfmfc ) CALL tra_mfc ( kstp, Nbb, ts, Nrhs ) ! Mass Flux Convection 317 325 IF( ln_zdfosm ) THEN
Note: See TracChangeset
for help on using the changeset viewer.
