Changeset 4389 for branches/2013/dev_r4050_NOC_WaD/NEMOGCM/NEMO/OPA_SRC
- Timestamp:
- 2014-01-31T15:55:04+01:00 (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2013/dev_r4050_NOC_WaD/NEMOGCM/NEMO/OPA_SRC/DYN/wadlmt.F90
r4385 r4389 52 52 REAL(wp) :: zzflxp, zzflxn ! local scalars 53 53 REAL(wp) :: ztmp ! local scalars 54 INTEGER, POINTER, DIMENSION(:,:) :: iwdflg !: W/D cell flag 54 INTEGER, POINTER, DIMENSION(:,:) :: iwdflg !: W/D cell flag (should become global) 55 55 REAL(wp), POINTER, DIMENSION(:,:) :: zwdlmtu, zwdlmtv !: W/D flux limiters 56 56 REAL(wp), POINTER, DIMENSION(:,:) :: zflxp, zflxn ! specific 2D workspace … … 158 158 159 159 IF(jk1 == nn_waditr) THEN 160 161 IF(zflxu1(ji, jj ) >= 0._wp) zwdlmtu(ji, jj ) = 0._wp 162 IF(zflxu1(ji-1,jj ) < 0._wp) zwdlmtu(ji-1,jj ) = 0._wp 163 IF(zflxv1(ji, jj ) >= 0._wp) zwdlmtv(ji, jj ) = 0._wp 164 IF(zflxv1(ji, jj-1) < 0._wp) zwdlmtv(ji, jj-1) = 0._wp 165 166 CALL lbc_lnk( zwdlmtu, 'U', 1. ) 167 CALL lbc_lnk( zwdlmtv, 'V', 1. ) 160 DO jj = 2, jpjm1 161 DO ji = fs_2, fs_jpim1 ! vector opt. 162 163 IF(zflxu1(ji, jj ) >= 0._wp) zwdlmtu(ji, jj ) = 0._wp 164 IF(zflxu1(ji-1,jj ) < 0._wp) zwdlmtu(ji-1,jj ) = 0._wp 165 IF(zflxv1(ji, jj ) >= 0._wp) zwdlmtv(ji, jj ) = 0._wp 166 IF(zflxv1(ji, jj-1) < 0._wp) zwdlmtv(ji, jj-1) = 0._wp 167 168 CALL lbc_lnk( zwdlmtu, 'U', 1. ) 169 CALL lbc_lnk( zwdlmtv, 'V', 1. ) 170 END DO ! ji loop 171 END DO ! jj loop 168 172 END IF 169 173
Note: See TracChangeset
for help on using the changeset viewer.