source:
branches/UKMO/r6232_CO6_CO5_zenv_pomsdwl/NEMOGCM/NEMO/OPA_SRC/SBC/inv_bar_vel_mod.F90
@
8186
Last change on this file since 8186 was 7454, checked in by jcastill, 8 years ago | |
---|---|
File size: 1.0 KB |
Line | |
---|---|
1 | MODULE inv_bar_vel_mod |
2 | |
3 | !! Modules used |
4 | USE oce |
5 | USE dom_oce |
6 | USE sbcflx |
7 | USE sbc_oce |
8 | |
9 | ! REAL(wp), DIMENSION(jpi,jpj) :: apgu, apgv !Pres grad force u and v component |
10 | PUBLIC inv |
11 | |
12 | CONTAINS |
13 | |
14 | SUBROUTINE inv(kt) |
15 | USE oce |
16 | USE dom_oce |
17 | USE sbcflx |
18 | USE sbc_oce |
19 | |
20 | INTEGER, INTENT ( in ) :: kt ! ocean time step |
21 | |
22 | ! REAL(wp), INTENT(in, out) :: ua,va |
23 | ! REAL(wp), DIMENSION(jpi,jpj) :: apgu, apgv !Pres grad force u and v component |
24 | |
25 | ! IF ( ln_shelf_flx ) then |
26 | ! Adjust the velocity |
27 | ! write(10,*)'kt = ', kt |
28 | |
29 | Do jk = 1, jpk |
30 | Do jj = 1, jpj !should this be jpjm1 because the apgu and apgv were calcd with jpjm1? |
31 | Do ji = 1, jpi |
32 | ! Do jk = 1, jpk-1 |
33 | ! Do jj = 2, jpj-1 !should this be jpjm1 because the apgu and apgv were calcd with jpjm1? |
34 | ! Do ji = 2, jpi-1 |
35 | ua(ji,jj,jk) = ua(ji,jj,jk) + apgu(ji,jj) |
36 | va(ji,jj,jk) = va(ji,jj,jk) + apgv(ji,jj) |
37 | ! write(10,*)ua(ji,jj,jk) |
38 | ! write(11,*)va(ji,jj,jk) |
39 | End do |
40 | End do |
41 | End do |
42 | ! ENDIF |
43 | END SUBROUTINE inv |
44 | END MODULE inv_bar_vel_mod |
Note: See TracBrowser
for help on using the repository browser.