source: lmdz_wrf/WRFV3/share/module_interp_fcn.F @ 1

Last change on this file since 1 was 1, checked in by lfita, 10 years ago
  • -- --- Opening of the WRF+LMDZ coupling repository --- -- -

WRF: version v3.3
LMDZ: version v1818

More details in:

File size: 2.0 KB
Line 
1!WRF:MEDIATION_LAYER:INTERPOLATIONFUNCTION
2!
3
4   SUBROUTINE interp_fcn ( cfld,                                 &  ! CD field
5                           cids, cide, ckds, ckde, cjds, cjde,   &
6                           cims, cime, ckms, ckme, cjms, cjme,   &
7                           cits, cite, ckts, ckte, cjts, cjte,   &
8                           nfld,                                 &  ! ND field
9                           nids, nide, nkds, nkde, njds, njde,   &
10                           nims, nime, nkms, nkme, njms, njme,   &
11                           nits, nite, nkts, nkte, njts, njte,   &
12                           nri, nrj                             )   ! nest ratios
13     IMPLICIT NONE
14
15
16     INTEGER, INTENT(IN) :: cids, cide, ckds, ckde, cjds, cjde,   &
17                            cims, cime, ckms, ckme, cjms, cjme,   &
18                            cits, cite, ckts, ckte, cjts, cjte,   &
19                            nids, nide, nkds, nkde, njds, njde,   &
20                            nims, nime, nkms, nkme, njms, njme,   &
21                            nits, nite, nkts, nkte, njts, njte,   &
22                            nri, nrj
23
24     REAL, DIMENSION ( cims:cime, ckms:ckme, cjms:cjme ) :: cfld
25     REAL, DIMENSION ( nims:nime, nkms:nkme, njms:njme ) :: nfld
26
27     ! Local
28
29     INTEGER ci, cj, ck, ni, nj, nk, ip, jp
30
31     ! Iterate over the ND tile and comtute the values
32     ! from the CD tile.
33
34     DO nj = njts, njte
35        cj = cjts + nj / nrj          ! coord of CD point
36        jp = cjts + mod ( nj , nrj )  ! coord of ND w/i CD point
37        DO nk = nkts, nkte
38           ck = nk
39           DO ni = nits, nite
40              ci = cits + ni / nri          ! coord of CD point
41              ip = cits + mod ( ni , nri )  ! coord of ND w/i CD point
42
43              ! This is a trivial implementation of the interp_fcn; just copies
44              ! the values from the CD into the ND
45
46              nfld( ni, nk, nj ) = cfld( ci , ck , cj )
47
48           ENDDO
49        ENDDO
50     ENDDO
51
52     RETURN
53
54   END SUBROUTINE interp_fcn
55
56
Note: See TracBrowser for help on using the repository browser.