source:
trunk/NEMOGCM/NEMO/OPA_SRC/LDF/ldfdyn_antarctic.h90
@
3389
Last change on this file since 3389 was 2528, checked in by rblod, 13 years ago | |
---|---|
|
|
File size: 1.4 KB |
Rev | Line | |
---|---|---|
[3] | 1 | !!---------------------------------------------------------------------- |
2 | !! ldfdyn_antarctic.h90 2D (i,j) antarctic zoom viscosity coeff. | |
3 | !!---------------------------------------------------------------------- | |
[247] | 4 | !!---------------------------------------------------------------------- |
[2528] | 5 | !! NEMO/OPA 3.3 , NEMO Consortium (2010) |
[1152] | 6 | !! $Id$ |
[2528] | 7 | !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) |
[247] | 8 | !!---------------------------------------------------------------------- |
[3] | 9 | ! |
10 | ! 1.2 Modify ahm | |
11 | ! -------------- | |
12 | ||
13 | IF(lwp)WRITE(numout,*) ' inildf: Antarctic ocean' | |
14 | IF(lwp)WRITE(numout,*) ' no tropics, no reduction of ahm' | |
15 | IF(lwp)WRITE(numout,*) ' north boundary increase' | |
16 | ||
17 | ahm1(:,:) = ahm0 | |
18 | ahm2(:,:) = ahm0 | |
19 | ||
20 | ijpt0=max(1,min(49 -njmpp+1,jpj)) | |
21 | ijpt1=max(0,min(49-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(48 -njmpp+1,jpj)) | |
27 | ijpt1=max(0,min(48-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(47 -njmpp+1,jpj)) | |
33 | ijpt1=max(0,min(47-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(46 -njmpp+1,jpj)) | |
39 | ijpt1=max(0,min(46-njmpp+1,jpj-1)) | |
40 | DO jj=ijpt0,ijpt1 | |
41 | ahm2(:,jj)=ahm0*1.1 | |
42 | END DO | |
43 | C |
Note: See TracBrowser
for help on using the repository browser.