Support has been added for the use of key_dynldf_c3d with the ORCA1 configuration. This entails simple changes to ldfdyn_c3d.h90 and includes a generalisation of the vertical scaling coefficient from a function of jk to a function of depth so that it can be used with different numbers of vertical levels. The function (thanks to Ian Culverwell, UKMO) provides exactly the original values when using the standard 31 ORCA levels.

ldfslp.F90 is missing a dummy version of the ldf_slp_grif routine which causes a linking failure if key_ldfslp is not active. Simple fix submitted at revision 2605

Hi, With s-coordinates on multiple processors (4x4) ORCA2_LIM2 was failing with undefined e3

===>>> : E R R O R

===========

zgr_sco : e3w or e3t =< 0 at point (i,j,k)= 1 39

..etc

Gurvan suggested the following fix, which cured the problem. However might there be a better way of solving the problem involving adding a 4th argument to the preceding call to lbc_lnk CALL lbc_lnk( zenv, 'T', 1._wp ) ?

$ svn diff NEMOGCM/NEMO/OPA_SRC/DOM/domzgr.F90

Index: NEMOGCM/NEMO/OPA_SRC/DOM/domzgr.F90 =================================================================== --- NEMOGCM/NEMO/OPA_SRC/DOM/domzgr.F90 (revision 2769) +++ NEMOGCM/NEMO/OPA_SRC/DOM/domzgr.F90 (working copy) @@ -1262,6 +1262,14 @@

! ! ================ ! ! ! ! envelop bathymetry saved in hbatt

+ DO ji = nlci+1, jpi ! Fill ghost line + zenv(ji,1:nlcj) = zenv(nlci,1:nlcj) + END DO + + DO jj = nlcj+1, jpj + zenv(:,jj) = zenv(:,nlcj) + END DO +

hbatt(:,:) = zenv(:,:) IF( MINVAL( gphit(:,:) ) * MAXVAL( gphit(:,:) ) <= 0._wp ) THEN

CALL ctl_warn( ' s-coordinates are tapered in vicinity of the Equator' )