Changeset 314 for trunk/NEMO/OPA_SRC/DYN/dynspg_rl.F90
- Timestamp:
- 2005-09-30T12:24:23+02:00 (19 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/OPA_SRC/DYN/dynspg_rl.F90
r247 r314 25 25 USE solsor ! Successive Over-relaxation solver 26 26 USE solfet ! FETI solver 27 USE solsor_e ! Successive Over-relaxation solver with MPP optimization 27 28 USE solisl ! ??? 28 29 USE obc_oce ! Lateral open boundary condition … … 106 107 !! * Local declarations 107 108 INTEGER :: ji, jj, jk ! dummy loop indices 108 REAL(wp) :: zgwgt,zbsfa, zgcx, z2dt109 REAL(wp) :: zbsfa, zgcx, z2dt 109 110 # if defined key_obc 110 111 INTEGER :: ip, ii, ij … … 203 204 END DO 204 205 END DO 206 ! applied the lateral boundary conditions 207 IF( nsolv == 4) CALL lbc_lnk_e( gcb, c_solver_pt, 1. ) 205 208 206 209 !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,synchro,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, … … 208 211 ! Relative precision (computation on one processor) 209 212 rnorme = 0.e0 210 DO jj = 1, jpj 211 DO ji = 1, jpi 212 zgwgt = gcdmat(ji,jj) * gcb(ji,jj) 213 rnorme = rnorme + gcb(ji,jj) * zgwgt 214 END DO 215 END DO 213 rnorme = SUM( gcb(1:nlci,1:nlcj) * gcdmat(1:nlci,1:nlcj) * gcb(1:nlci,1:nlcj) * bmask(:,:) ) 216 214 IF( lk_mpp ) CALL mpp_sum( rnorme ) ! sum over the global domain 217 215 … … 236 234 CASE( 3 ) ! FETI solver 237 235 CALL sol_fet( kindic ) 236 CASE( 4 ) ! successive-over-relaxation with extra outer halo 237 CALL sol_sor_e( kindic ) 238 238 CASE DEFAULT ! e r r o r in nsolv namelist parameter 239 239 IF(lwp) WRITE(numout,cform_err) 240 IF(lwp) WRITE(numout,*) ' dyn_spg_rl : e r r o r, nsolv = 1, 2 or 3'240 IF(lwp) WRITE(numout,*) ' dyn_spg_rl : e r r o r, nsolv = 1, 2 ,3 or 4' 241 241 IF(lwp) WRITE(numout,*) ' ~~~~~~~~~~ not = ', nsolv 242 242 nstop = nstop + 1 … … 248 248 ! ---------------- 249 249 250 bsfd( :,:) = gcx(:,:)250 bsfd(1:nlci,1:nlcj) = gcx(1:nlci,1:nlcj) 251 251 252 252
Note: See TracChangeset
for help on using the changeset viewer.