1 | !!---------------------------------------------------------------------- |
---|
2 | !! ldftra_arctic.h90 2D (i,j) arctic zoom viscosity coeff. |
---|
3 | !!---------------------------------------------------------------------- |
---|
4 | !!---------------------------------------------------------------------- |
---|
5 | !! OPA 9.0 , LOCEAN-IPSL (2005) |
---|
6 | !! $Id: ldfdyn_arctic.h90 1152 2008-06-26 14:11:13Z rblod $ |
---|
7 | !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt |
---|
8 | !!---------------------------------------------------------------------- |
---|
9 | |
---|
10 | ! 1.2 Modify ahm |
---|
11 | ! -------------- |
---|
12 | IF(lwp)WRITE(numout,*) ' inildf: Arctic ocean' |
---|
13 | IF(lwp)WRITE(numout,*) ' no tropics, no reduction of ahm' |
---|
14 | IF(lwp)WRITE(numout,*) ' south and west boundary increase' |
---|
15 | |
---|
16 | |
---|
17 | ahm1(:,:) = ahm0 |
---|
18 | ahm2(:,:) = ahm0 |
---|
19 | |
---|
20 | ijpt0=max(1,min(98-jpjzoom+1-njmpp+1,jpj)) |
---|
21 | ijpt1=max(0,min(98-jpjzoom+1-njmpp+1,jpj-1)) |
---|
22 | DO jj=ijpt0,ijpt1 |
---|
23 | ahm2(:,jj)=ahm0*2. |
---|
24 | ahm1(:,jj)=ahm0*2. |
---|
25 | END DO |
---|
26 | ijpt0=max(1,min(99-jpjzoom+1-njmpp+1,jpj)) |
---|
27 | ijpt1=max(0,min(99-jpjzoom+1-njmpp+1,jpj-1)) |
---|
28 | DO jj=ijpt0,ijpt1 |
---|
29 | ahm2(:,jj)=ahm0*1.9 |
---|
30 | ahm1(:,jj)=ahm0*1.75 |
---|
31 | END DO |
---|
32 | ijpt0=max(1,min(100-jpjzoom+1-njmpp+1,jpj)) |
---|
33 | ijpt1=max(0,min(100-jpjzoom+1-njmpp+1,jpj-1)) |
---|
34 | DO jj=ijpt0,ijpt1 |
---|
35 | ahm2(:,jj)=ahm0*1.5 |
---|
36 | ahm1(:,jj)=ahm0*1.25 |
---|
37 | END DO |
---|
38 | ijpt0=max(1,min(101-jpjzoom+1-njmpp+1,jpj)) |
---|
39 | ijpt1=max(0,min(101-jpjzoom+1-njmpp+1,jpj-1)) |
---|
40 | DO jj=ijpt0,ijpt1 |
---|
41 | ahm2(:,jj)=ahm0*1.1 |
---|
42 | END DO |
---|
43 | |
---|