Changes between Version 28 and Version 29 of ticket/0863


Ignore:
Timestamp:
2013-01-08T11:43:19+01:00 (8 years ago)
Author:
mlelod
Comment:

Legend:

Unmodified
Added
Removed
Modified
  • ticket/0863

    v28 v29  
    153153 
    154154}}} 
     155 
     156== To be done == 
     157 
     158=== Compatibility with een vorticity scheme === 
     159The vertical scale factor at f points should remain constant in time if the f-point belongs to an "earth cell". 
     160{{{{ 
     161      IF( kt == nit000 .OR. lk_vvl ) THEN      ! reciprocal of e3 at F-point (masked averaging of e3t) 
     162         DO jk = 1, jpk 
     163            DO jj = 1, jpjm1 
     164               DO ji = 1, jpim1 
     165                  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)   & 
     166                     &             + fse3t(ji,jj  ,jk)*tmask(ji,jj  ,jk) + fse3t(ji+1,jj  ,jk)*tmask(ji+1,jj  ,jk) ) * 0.25 
     167                  IF( ze3f(ji,jj,jk) /= 0._wp )   ze3f(ji,jj,jk) = 1._wp / ze3f(ji,jj,jk) 
     168               END DO 
     169            END DO 
     170         END DO 
     171         CALL lbc_lnk( ze3f, 'F', 1. ) 
     172      ENDIF 
     173}}}} 
     174 
    155175---- 
    156176== Testing ==