Changes between Version 7 and Version 8 of 2021WP/KNL-01_Sibylle_RK3_stage1

Timestamp:

2021-03-01T18:51:01+01:00 (3 years ago)

Author:

techene

Comment:

--

2021WP/KNL-01_Sibylle_RK3_stage1

@@ -86 +86 @@
 }}}
 
-3D velocity added as an input argument of advective routines 
+3D velocity added as an input argument of advective routines passed through dyn_adv
 {{{
 OCE
  |--DYN
      |-- dynadv.F90
-         SUBROUTINE dyn_adv( kt, Kbb, Kmm, puu, pvv, Krhs, pau, pav, paw ) 
+         SUBROUTINE dyn_adv( kt, Kbb, Kmm, puu, pvv, Krhs, pau, pav, paw )
+         ...
+         CALL dyn_adv_cen2( kt                , Kmm, puu, pvv, Krhs, pau, pav, paw )   ! 2nd order centered scheme 
+         CALL dyn_adv_ubs ( kt           , Kbb, Kmm, puu, pvv, Krhs, pau, pav, paw )   ! 3rd order UBS      scheme (UP3) 
+     |-- dynadv_cen2.F90
+         SUBROUTINE dyn_adv_cen2( kt, Kmm, puu, pvv, Krhs, pau, pav, paw )  
+         ...       
+         IF( PRESENT( pau ) ) THEN     ! RK3: advective velocity (pau,pav,paw) /= advected velocity (puu,pvv,ww) 
+            zptu => pau(:,:,:) 
+            ...
+         zfu(:,:,jk) = 0.25_wp * e2u(:,:) * e3u(:,:,jk,Kmm) * zptu(:,:,jk) 
+     |-- dynadv_ubs.F90
+         SUBROUTINE dyn_adv_ubs( kt, Kbb, Kmm, puu, pvv, Krhs, pau, pav, paw ) 
+         ...
+         IF( PRESENT( pau ) ) THEN     ! RK3: advective velocity (pau,pav,paw) /= advected velocity (puu,pvv,ww) 
+            zptu => pau(:,:,:) 
+            ...
+         zfu(:,:,jk) = e2u(:,:) * e3u(:,:,jk,Kmm) * zptu(:,:,jk) 
+ |--TRA  
+     |-- traadv.F90
+         SUBROUTINE tra_adv( kt, Kbb, Kmm, pts, Krhs, pau, pav, paw ) 
+         ...
+         IF( PRESENT( pau ) ) THEN     ! RK3: advective velocity (pau,pav,paw) /= advected velocity (puu,pvv,ww) 
+            zptu => pau(:,:,:) 
+            ...
+         zuu(ji,jj,jk) = e2u  (ji,jj) * e3u(ji,jj,jk,Kmm) * ( zptu(ji,jj,jk) + usd(ji,jj,jk) ) 
+
 }}}
+
 
 === Documentation updates