Context

try to run ORCA2_ICE_PISCES or C1D_PAPA with key_mpi_off cpp kay

Analysis

both config will crash in lbc_lnk_pt2pt routine, fillmodeoption doesn't seem to work correctly without MPI

Recommendation

find a lnb_lnk expert to have a look at it

Context

try C1D_PAPA config with post-merge trunk

Analysis

multiple problems

Recommendation

fix them all

Context

heat content diag computation in "dia_hth_htc" routine

Analysis

2 possible small errors in htc computation:

• vertical htc integration index should start at "1" instead of "2"
• missing piece of vertical integration between the last considered full cell and htc depth (when htc depth is located between two vertical levels)

Fix

in dia_hth_htc, replace:

      ilevel(:,:) = 1
      DO_3D( 1, 1, 1, 1, 2, jpkm1 )
         IF( ( gdept(ji,jj,jk,Kmm) < pdep ) .AND. ( tmask(ji,jj,jk) == 1 ) ) THEN
             ilevel(ji,jj) = jk
             zthick(ji,jj) = zthick(ji,jj) + e3t(ji,jj,jk,Kmm)
             phtc  (ji,jj) = phtc  (ji,jj) + e3t(ji,jj,jk,Kmm) * pt(ji,jj,jk)
         ENDIF
      END_3D
      !
      DO_2D( 1, 1, 1, 1 )
         ik = ilevel(ji,jj)
         zthick(ji,jj) = pdep - zthick(ji,jj)   !   remaining thickness to reach depht pdep
         phtc(ji,jj)   = phtc(ji,jj)    &
            &           + pt (ji,jj,ik+1) * MIN( e3t(ji,jj,ik+1,Kmm), zthick(ji,jj) ) &
                                                       * tmask(ji,jj,ik+1)
      END_2D

by

      ilevel(:,:) = 1
      DO_3D( 1, 1, 1, 1, 1, jpkm1 )
         IF( ( gdepw(ji,jj,jk+1,Kmm) < pdep ) .AND. ( tmask(ji,jj,jk) == 1 ) ) THEN
             ilevel(ji,jj) = jk+1
             zthick(ji,jj) = zthick(ji,jj) + e3t(ji,jj,jk,Kmm)
             phtc  (ji,jj) = phtc  (ji,jj) + e3t(ji,jj,jk,Kmm) * pt(ji,jj,jk)
         ENDIF
      END_3D
      !
      DO_2D( 1, 1, 1, 1 )
         ik = ilevel(ji,jj)
         IF( tmask(ji,jj,ik) == 1 ) THEN
            zthick(ji,jj) = MIN ( gdepw(ji,jj,ik+1,Kmm), pdep ) - zthick(ji,jj)   ! remaining thickness to reach depht pdep
            phtc(ji,jj)   = phtc(ji,jj) + ptn(ji,jj,ik) * zthick(ji,jj)
         ENDIF
      END_2D