New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
ldfdyn_antarctic.h90 in trunk/NEMO/OPA_SRC/LDF – NEMO

source: trunk/NEMO/OPA_SRC/LDF/ldfdyn_antarctic.h90 @ 3

Last change on this file since 3 was 3, checked in by opalod, 20 years ago

Initial revision

  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 1.2 KB
Line 
1   !!----------------------------------------------------------------------
2   !!  ldfdyn_antarctic.h90        2D (i,j) antarctic zoom viscosity coeff.
3   !!----------------------------------------------------------------------
4
5!
6! 1.2 Modify ahm
7! --------------
8 
9      IF(lwp)WRITE(numout,*) ' inildf: Antarctic ocean'
10      IF(lwp)WRITE(numout,*) '         no tropics, no reduction of ahm'
11      IF(lwp)WRITE(numout,*) '         north boundary increase'
12
13      ahm1(:,:) = ahm0
14      ahm2(:,:) = ahm0
15
16      ijpt0=max(1,min(49 -njmpp+1,jpj))
17      ijpt1=max(0,min(49-njmpp+1,jpj-1))
18      DO jj=ijpt0,ijpt1
19         ahm2(:,jj)=ahm0*2.
20         ahm1(:,jj)=ahm0*2.
21      END DO
22      ijpt0=max(1,min(48 -njmpp+1,jpj))
23      ijpt1=max(0,min(48-njmpp+1,jpj-1))
24      DO jj=ijpt0,ijpt1
25         ahm2(:,jj)=ahm0*1.9
26         ahm1(:,jj)=ahm0*1.75
27      END DO
28      ijpt0=max(1,min(47 -njmpp+1,jpj))
29      ijpt1=max(0,min(47-njmpp+1,jpj-1))
30      DO jj=ijpt0,ijpt1
31         ahm2(:,jj)=ahm0*1.5
32         ahm1(:,jj)=ahm0*1.25
33      END DO
34      ijpt0=max(1,min(46 -njmpp+1,jpj))
35      ijpt1=max(0,min(46-njmpp+1,jpj-1))
36      DO jj=ijpt0,ijpt1
37         ahm2(:,jj)=ahm0*1.1
38      END DO
39C
Note: See TracBrowser for help on using the repository browser.