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_arctic.h90 in trunk/NEMOGCM/NEMO/OPA_SRC/LDF – NEMO

source: trunk/NEMOGCM/NEMO/OPA_SRC/LDF/ldfdyn_arctic.h90 @ 2528

Last change on this file since 2528 was 2528, checked in by rblod, 13 years ago

Update NEMOGCM from branch nemo_v3_3_beta

  • Property svn:keywords set to Id
File size: 1.5 KB
Line 
1   !!----------------------------------------------------------------------
2   !!  ldftra_arctic.h90              2D (i,j) arctic zoom viscosity coeff.
3   !!----------------------------------------------------------------------
4   !!----------------------------------------------------------------------
5   !! NEMO/OPA 3.3 , NEMO Consortium (2010)
6   !! $Id$
7   !! Software governed by the CeCILL licence (NEMOGCM/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
Note: See TracBrowser for help on using the repository browser.