source:
tags/nemo_dev_x5/NEMO/OPA_SRC/LDF/ldfdyn_arctic.h90
@
4038
Last change on this file since 4038 was 139, checked in by cvs2svn, 20 years ago | |
---|---|
|
|
File size: 1.2 KB |
Line | |
---|---|
1 | !!---------------------------------------------------------------------- |
2 | !! ldftra_arctic.h90 2D (i,j) arctic zoom viscosity coeff. |
3 | !!---------------------------------------------------------------------- |
4 | |
5 | |
6 | ! 1.2 Modify ahm |
7 | ! -------------- |
8 | IF(lwp)WRITE(numout,*) ' inildf: Arctic ocean' |
9 | IF(lwp)WRITE(numout,*) ' no tropics, no reduction of ahm' |
10 | IF(lwp)WRITE(numout,*) ' south and west boundary increase' |
11 | |
12 | |
13 | ahm1(:,:) = ahm0 |
14 | ahm2(:,:) = ahm0 |
15 | |
16 | ijpt0=max(1,min(98-jpjzoom+1-njmpp+1,jpj)) |
17 | ijpt1=max(0,min(98-jpjzoom+1-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(99-jpjzoom+1-njmpp+1,jpj)) |
23 | ijpt1=max(0,min(99-jpjzoom+1-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(100-jpjzoom+1-njmpp+1,jpj)) |
29 | ijpt1=max(0,min(100-jpjzoom+1-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(101-jpjzoom+1-njmpp+1,jpj)) |
35 | ijpt1=max(0,min(101-jpjzoom+1-njmpp+1,jpj-1)) |
36 | DO jj=ijpt0,ijpt1 |
37 | ahm2(:,jj)=ahm0*1.1 |
38 | END DO |
39 |
Note: See TracBrowser
for help on using the repository browser.