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.
limtab_2.F90 in trunk/NEMOGCM/NEMO/LIM_SRC_2 – NEMO

source: trunk/NEMOGCM/NEMO/LIM_SRC_2/limtab_2.F90 @ 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: 2.3 KB
Line 
1MODULE limtab_2
2   !!======================================================================
3   !!                       ***  MODULE limtab_2   ***
4   !!              transform 1D (2D) array to a 2D (1D) table
5   !!======================================================================
6#if defined key_lim2
7   !!----------------------------------------------------------------------
8   !!   tab_2d_1d  : 2-D to 1-D
9   !!   tab_1d_2d  : 1-D to 2-D
10   !!----------------------------------------------------------------------
11   !! * Modules used
12   USE par_kind
13
14   IMPLICIT NONE
15   PRIVATE
16
17   !! * Routine accessibility
18   PUBLIC tab_2d_1d_2  ! called by lim_ther
19   PUBLIC tab_1d_2d_2  ! called by lim_ther
20
21   !!----------------------------------------------------------------------
22   !! NEMO/LIM2 3.3 , UCL - NEMO Consortium (2010)
23   !! $Id$
24   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt)
25   !!----------------------------------------------------------------------
26CONTAINS
27
28   SUBROUTINE tab_2d_1d_2 ( ndim1d, tab1d, tab2d, ndim2d_x, ndim2d_y, tab_ind )
29
30      INTEGER, INTENT(in) :: &
31         ndim1d, ndim2d_x, ndim2d_y
32
33      REAL(wp), DIMENSION (ndim2d_x, ndim2d_y), INTENT(in) ::  &
34         tab2d
35
36      INTEGER, DIMENSION ( ndim1d), INTENT ( in) :: &
37         tab_ind
38
39      REAL(wp), DIMENSION(ndim1d), INTENT ( out) ::  & 
40         tab1d
41
42      INTEGER ::  &
43         jn , jid, jjd
44       
45      DO jn = 1, ndim1d
46         jid        = MOD( tab_ind(jn) - 1, ndim2d_x ) + 1
47         jjd        = ( tab_ind(jn) - 1 ) / ndim2d_x + 1
48         tab1d( jn) = tab2d( jid, jjd)
49      END DO
50
51   END SUBROUTINE tab_2d_1d_2
52
53
54   SUBROUTINE tab_1d_2d_2 ( ndim1d, tab2d, tab_ind, tab1d, ndim2d_x, ndim2d_y )
55
56      INTEGER, INTENT ( in) :: &
57         ndim1d, ndim2d_x, ndim2d_y
58
59      INTEGER, DIMENSION (ndim1d) , INTENT (in) :: &
60         tab_ind
61
62      REAL(wp), DIMENSION(ndim1d), INTENT (in) ::  &
63         tab1d 
64
65      REAL(wp), DIMENSION (ndim2d_x, ndim2d_y), INTENT ( out) :: &
66         tab2d
67
68      INTEGER :: &
69         jn, jid, jjd
70
71      DO jn = 1, ndim1d
72         jid             = MOD( tab_ind(jn) - 1, ndim2d_x) + 1
73         jjd             =    ( tab_ind(jn) - 1 ) / ndim2d_x  + 1
74         tab2d(jid, jjd) = tab1d( jn)
75      END DO
76
77   END SUBROUTINE tab_1d_2d_2
78
79#endif
80END MODULE limtab_2
Note: See TracBrowser for help on using the repository browser.