!!----------------------------------------------------------------------
   !!  ldfdyn_antarctic.h90        2D (i,j) antarctic zoom viscosity coeff.
   !!----------------------------------------------------------------------
   !!----------------------------------------------------------------------
   !!  OPA 9.0 , LOCEAN-IPSL (2005) 
   !! $Id$ 
   !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt 
   !!----------------------------------------------------------------------
!
! 1.2 Modify ahm
! --------------
 
      IF(lwp)WRITE(numout,*) ' inildf: Antarctic ocean'
      IF(lwp)WRITE(numout,*) '         no tropics, no reduction of ahm'
      IF(lwp)WRITE(numout,*) '         north boundary increase'

      ahm1(:,:) = ahm0
      ahm2(:,:) = ahm0

      ijpt0=max(1,min(49 -njmpp+1,jpj))
      ijpt1=max(0,min(49-njmpp+1,jpj-1))
      DO jj=ijpt0,ijpt1
         ahm2(:,jj)=ahm0*2.
         ahm1(:,jj)=ahm0*2.
      END DO
      ijpt0=max(1,min(48 -njmpp+1,jpj))
      ijpt1=max(0,min(48-njmpp+1,jpj-1))
      DO jj=ijpt0,ijpt1
         ahm2(:,jj)=ahm0*1.9
         ahm1(:,jj)=ahm0*1.75
      END DO
      ijpt0=max(1,min(47 -njmpp+1,jpj))
      ijpt1=max(0,min(47-njmpp+1,jpj-1))
      DO jj=ijpt0,ijpt1
         ahm2(:,jj)=ahm0*1.5
         ahm1(:,jj)=ahm0*1.25
      END DO
      ijpt0=max(1,min(46 -njmpp+1,jpj))
      ijpt1=max(0,min(46-njmpp+1,jpj-1))
      DO jj=ijpt0,ijpt1
         ahm2(:,jj)=ahm0*1.1
      END DO
C