Changeset 15267 for NEMO/trunk/src/OCE/LBC/lbc_nfd_generic.h90
- Timestamp:
- 2021-09-17T11:04:34+02:00 (3 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/trunk/src/OCE/LBC/lbc_nfd_generic.h90
r14433 r15267 29 29 SELECT CASE ( cd_nat(jf) ) 30 30 CASE ( 'T' , 'W' ) ! T-, W-point 31 DO jl = 1, ipl ;DO jk = 1, ipk31 DO jl = 1, ipl ; DO jk = 1, ipk 32 32 ! 33 33 ! last khls lines (from ipj to ipj-khls+1) : full … … 53 53 DO ji = 1, 1 ! point ipi - khls + 1 54 54 ii1 = ipi - khls + ji 55 ii2 = 55 ii2 = khls + ji 56 56 ptab(jf)%pt4d(ii1,ij1,jk,jl) = psgn(jf) * ptab(jf)%pt4d(ii2,ij2,jk,jl) 57 57 END DO … … 84 84 END DO; END DO 85 85 CASE ( 'U' ) ! U-point 86 DO jl = 1, ipl ;DO jk = 1, ipk86 DO jl = 1, ipl ; DO jk = 1, ipk 87 87 ! 88 88 ! last khls lines (from ipj to ipj-khls+1) : full … … 103 103 DO ji = 1, khls ! last khls points 104 104 ii1 = ipi - khls + ji ! ends at: ipi - khls + khls = ipi 105 ii2 = ipi - khls + 1 - ji! ends at: ipi - khls + 1 - khls = ipi - 2*khls + 1105 ii2 = ipi - khls - ji + 1 ! ends at: ipi - khls + 1 - khls = ipi - 2*khls + 1 106 106 ptab(jf)%pt4d(ii1,ij1,jk,jl) = psgn(jf) * ptab(jf)%pt4d(ii2,ij2,jk,jl) 107 107 END DO … … 129 129 END DO; END DO 130 130 CASE ( 'V' ) ! V-point 131 DO jl = 1, ipl ;DO jk = 1, ipk131 DO jl = 1, ipl ; DO jk = 1, ipk 132 132 ! 133 133 ! last khls+1 lines (from ipj to ipj-khls) : full … … 153 153 DO ji = 1, 1 ! point ipi - khls + 1 154 154 ii1 = ipi - khls + ji 155 ii2 = 155 ii2 = khls + ji 156 156 ptab(jf)%pt4d(ii1,ij1,jk,jl) = psgn(jf) * ptab(jf)%pt4d(ii2,ij2,jk,jl) 157 157 END DO … … 165 165 END DO; END DO 166 166 CASE ( 'F' ) ! F-point 167 DO jl = 1, ipl ;DO jk = 1, ipk167 DO jl = 1, ipl ; DO jk = 1, ipk 168 168 ! 169 169 ! last khls+1 lines (from ipj to ipj-khls) : full … … 184 184 DO ji = 1, khls ! last khls points 185 185 ii1 = ipi - khls + ji ! ends at: ipi - khls + khls = ipi 186 ii2 = ipi - khls + 1 - ji! ends at: ipi - khls + 1 - khls = ipi - 2*khls + 1186 ii2 = ipi - khls - ji + 1 ! ends at: ipi - khls + 1 - khls = ipi - 2*khls + 1 187 187 ptab(jf)%pt4d(ii1,ij1,jk,jl) = psgn(jf) * ptab(jf)%pt4d(ii2,ij2,jk,jl) 188 188 END DO … … 198 198 SELECT CASE ( cd_nat(jf) ) 199 199 CASE ( 'T' , 'W' ) ! T-, W-point 200 DO jl = 1, ipl ;DO jk = 1, ipk200 DO jl = 1, ipl ; DO jk = 1, ipk 201 201 ! 202 202 ! first: line number ipj-khls : 3 points … … 212 212 DO ji = 1, 1 ! points ipi - khls 213 213 ii1 = ipi - khls + ji - 1 214 ii2 = 214 ii2 = khls + ji 215 215 ptab(jf)%pt4d(ii1,ij1,jk,jl) = ptab(jf)%pt4d(ii2,ij2,jk,jl) ! Warning: pb with sign... 216 216 END DO … … 240 240 DO ji = 1, khls ! last khls points 241 241 ii1 = ipi - khls + ji ! ends at: ipi - khls + khls = ipi 242 ii2 = ipi - khls + 1 - ji! ends at: ipi - khls + 1 - khls = ipi - 2*khls + 1242 ii2 = ipi - khls - ji + 1 ! ends at: ipi - khls + 1 - khls = ipi - 2*khls + 1 243 243 ptab(jf)%pt4d(ii1,ij1,jk,jl) = psgn(jf) * ptab(jf)%pt4d(ii2,ij2,jk,jl) 244 244 END DO … … 247 247 END DO; END DO 248 248 CASE ( 'U' ) ! U-point 249 DO jl = 1, ipl ;DO jk = 1, ipk249 DO jl = 1, ipl ; DO jk = 1, ipk 250 250 ! 251 251 ! last khls lines (from ipj to ipj-khls+1) : full … … 283 283 END DO; END DO 284 284 CASE ( 'V' ) ! V-point 285 DO jl = 1, ipl ;DO jk = 1, ipk285 DO jl = 1, ipl ; DO jk = 1, ipk 286 286 ! 287 287 ! last khls lines (from ipj to ipj-khls+1) : full … … 302 302 DO ji = 1, khls ! last khls points 303 303 ii1 = ipi - khls + ji ! ends at: ipi - khls + khls = ipi 304 ii2 = ipi - khls + 1 - ji! ends at: ipi - khls + 1 - khls = ipi - 2*khls + 1304 ii2 = ipi - khls - ji + 1 ! ends at: ipi - khls + 1 - khls = ipi - 2*khls + 1 305 305 ptab(jf)%pt4d(ii1,ij1,jk,jl) = psgn(jf) * ptab(jf)%pt4d(ii2,ij2,jk,jl) 306 306 END DO … … 328 328 END DO; END DO 329 329 CASE ( 'F' ) ! F-point 330 DO jl = 1, ipl ;DO jk = 1, ipk330 DO jl = 1, ipl ; DO jk = 1, ipk 331 331 ! 332 332 ! last khls lines (from ipj to ipj-khls+1) : full
Note: See TracChangeset
for help on using the changeset viewer.