Changeset 15628
- Timestamp:
- 2022-01-04T22:41:46+01:00 (3 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2021/dev_r14318_RK3_stage1/src/OCE/DOM/istate.F90
r15626 r15628 142 142 #endif 143 143 ! 144 ! Initialize "now" barotropic velocities:145 ! Do it whatever the free surface method, these arrays being used eventually146 !147 !!gm the use of umask & vmask is not necessary below as uu(:,:,:,Kmm), vv(:,:,:,Kmm), uu(:,:,:,Kbb), vv(:,:,:,Kbb) are always masked148 #if ! defined key_RK3149 uu_b(:,:,Kmm) = 0._wp ; vv_b(:,:,Kmm) = 0._wp150 DO_3D( nn_hls, nn_hls, nn_hls, nn_hls, 1, jpkm1 )151 uu_b(ji,jj,Kmm) = uu_b(ji,jj,Kmm) + e3u(ji,jj,jk,Kmm) * uu(ji,jj,jk,Kmm) * umask(ji,jj,jk)152 vv_b(ji,jj,Kmm) = vv_b(ji,jj,Kmm) + e3v(ji,jj,jk,Kmm) * vv(ji,jj,jk,Kmm) * vmask(ji,jj,jk)153 END_3D154 uu_b(:,:,Kmm) = uu_b(:,:,Kmm) * r1_hu(:,:,Kmm)155 vv_b(:,:,Kmm) = vv_b(:,:,Kmm) * r1_hv(:,:,Kmm)156 #endif157 144 ! 158 145 #if defined key_RK3 … … 169 156 uu_b(:,:,Kbb) = uu_b(:,:,Kbb) * r1_hu(:,:,Kbb) 170 157 vv_b(:,:,Kbb) = vv_b(:,:,Kbb) * r1_hv(:,:,Kbb) 171 ! 158 ! 159 #if defined key_RK3 160 ENDIF 161 #endif 162 ! Initialize "now" barotropic velocities: 163 ! Do it whatever the free surface method, these arrays being used eventually 164 ! 165 !!gm the use of umask & vmask is not necessary below as uu(:,:,:,Kmm), vv(:,:,:,Kmm), uu(:,:,:,Kbb), vv(:,:,:,Kbb) are always masked 166 #if ! defined key_RK3 167 uu_b(:,:,Kmm) = 0._wp ; vv_b(:,:,Kmm) = 0._wp 168 DO_3D( nn_hls, nn_hls, nn_hls, nn_hls, 1, jpkm1 ) 169 uu_b(ji,jj,Kmm) = uu_b(ji,jj,Kmm) + e3u(ji,jj,jk,Kmm) * uu(ji,jj,jk,Kmm) * umask(ji,jj,jk) 170 vv_b(ji,jj,Kmm) = vv_b(ji,jj,Kmm) + e3v(ji,jj,jk,Kmm) * vv(ji,jj,jk,Kmm) * vmask(ji,jj,jk) 171 END_3D 172 uu_b(:,:,Kmm) = uu_b(:,:,Kmm) * r1_hu(:,:,Kmm) 173 vv_b(:,:,Kmm) = vv_b(:,:,Kmm) * r1_hv(:,:,Kmm) 174 #else 175 ! 176 IF( .NOT. ln_rstart ) THEN 172 177 uu_b(:,:,Kmm) = uu_b(:,:,Kbb) ! Kmm value set to Kbb for initialisation in Agrif_Regrid 173 178 vv_b(:,:,Kmm) = vv_b(:,:,Kbb) 174 !175 #if defined key_RK3176 179 ENDIF 177 180 #endif
Note: See TracChangeset
for help on using the changeset viewer.