Changes between Version 28 and Version 29 of ticket/0863

Timestamp:

2013-01-08T11:43:19+01:00 (10 years ago)

Author:

mlelod

Comment:

--

ticket/0863

@@ -153 +153 @@
 
 }}}
+
+== To be done ==
+
+=== Compatibility with een vorticity scheme ===
+The vertical scale factor at f points should remain constant in time if the f-point belongs to an "earth cell".
+{{{{
+      IF( kt == nit000 .OR. lk_vvl ) THEN      ! reciprocal of e3 at F-point (masked averaging of e3t)
+         DO jk = 1, jpk
+            DO jj = 1, jpjm1
+               DO ji = 1, jpim1
+                  ze3f(ji,jj,jk) = ( fse3t(ji,jj+1,jk)*tmask(ji,jj+1,jk) + fse3t(ji+1,jj+1,jk)*tmask(ji+1,jj+1,jk)   &
+                     &             + fse3t(ji,jj  ,jk)*tmask(ji,jj  ,jk) + fse3t(ji+1,jj  ,jk)*tmask(ji+1,jj  ,jk) ) * 0.25
+                  IF( ze3f(ji,jj,jk) /= 0._wp )   ze3f(ji,jj,jk) = 1._wp / ze3f(ji,jj,jk)
+               END DO
+            END DO
+         END DO
+         CALL lbc_lnk( ze3f, 'F', 1. )
+      ENDIF
+}}}}
+
 ----
 == Testing ==