source: trunk/NEMOGCM/NEMO/OPA_SRC/LDF/ldfdyn_antarctic.h90 @ 2528

   !!----------------------------------------------------------------------
   !!  ldfdyn_antarctic.h90        2D (i,j) antarctic zoom viscosity coeff.
   !!----------------------------------------------------------------------
   !!----------------------------------------------------------------------
   !! NEMO/OPA 3.3 , NEMO Consortium (2010)
   !! $Id$ 
   !! Software governed by the CeCILL licence (NEMOGCM/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