source:
tags/nemo_v2/NEMO/OPA_SRC/LDF/ldfdyn_antarctic.h90
@
587
Last change on this file since 587 was 247, checked in by opalod, 19 years ago | |
---|---|
|
|
File size: 1.4 KB |
Line | |
---|---|
1 | !!---------------------------------------------------------------------- |
2 | !! ldfdyn_antarctic.h90 2D (i,j) antarctic zoom viscosity coeff. |
3 | !!---------------------------------------------------------------------- |
4 | !!---------------------------------------------------------------------- |
5 | !! OPA 9.0 , LOCEAN-IPSL (2005) |
6 | !! $Header$ |
7 | !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt |
8 | !!---------------------------------------------------------------------- |
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.