The restart file produced by the REBUILD TOOL do not account for the filling with some value of the neglected land subdomains, when the simulation is carried out only for the ocean subdomains (inimpp2).

In this case the rebuilt restart file presents NaNs? in all land areas, which are obviously not accounted in the loop over the single output files.

In turns out that using a single restart file produced by the rebuild to initialize NEMO with a different number of processors the set up of halos is not correct. This causes the failure of the model, e.g., in the schemes for computing the ssh.

A possible solution is to initialize to zero all the values in the ouput netcdf file, or a more refined one, which accounts for the domain decomposition in the rebuild tool and fill with zero only the land domain for which no output files are generated.

The computation of avmu and avmv in zdfgls and zdftke is not correct in 1D case. The values are systematically the half of the expected one (not the case for 3D configurations).

The following lines:

avmu(ji,jj,jk) = 0.5 * ( avm(ji,jj,jk) + avm(ji+1,jj ,jk) ) * umask(ji,jj,jk) avmv(ji,jj,jk) = 0.5 * ( avm(ji,jj,jk) + avm(ji ,jj+1,jk) ) * vmask(ji,jj,jk)

should be replaced by: avmu(ji,jj,jk) = ( avm(ji,jj,jk) + avm(ji+1,jj ,jk) ) / MAX( 1.e0 , umask(ji,jj,jk) + umask(ji+1,jj ,jk) ) avmv(ji,jj,jk) = ( avm(ji,jj,jk) + avm(ji ,jj+1,jk) ) / MAX( 1.e0 , vmask(ji,jj,jk) + vmask(ji ,jj+1,jk) )

From Jerome Chanut:
local variable named ztab are not seen as local variable by the conv...
for example the following piece of code:

      IF (MOD(nbcline,nbclineupdate) == 0) THEN
         CALL Agrif_Update_Variable(ztab,tsn_id, procname=updateTS)
      ELSE
         CALL Agrif_Update_Variable(ztab,tsn_id,locupdate=(/0,2/), procname=updateTS)
      ENDIF

is converted into

      IF (MOD(nbcline,nbclineupdate) == 0) THEN
         CALL Agrif_Update_Variable( 72, 171, procname=updateTS)
      ELSE
         CALL Agrif_Update_Variable( 72, 171,locupdate=(/0,2/), procname=updateTS)
      ENDIF

instead of

      IF (MOD(nbcline,nbclineupdate) == 0) THEN
         CALL Agrif_Update_Variable(ztab, 171, procname=updateTS)
      ELSE
         CALL Agrif_Update_Variable(ztab, 171,locupdate=(/0,2/), procname=updateTS)
      ENDIF

Solution (before a bufix in the conv): rename ztab (for example into ztab4d)

      IF (MOD(nbcline,nbclineupdate) == 0) THEN
         CALL Agrif_Update_Variable(ztab4d, 171, procname=updateTS)
      ELSE
         CALL Agrif_Update_Variable(ztab4d, 171,locupdate=(/0,2/), procname=updateTS)
      ENDIF