Changeset 3294 for trunk/NEMOGCM/NEMO/OPA_SRC/SOL/solsor.F90
- Timestamp:
- 2012-01-28T17:44:18+01:00 (12 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMOGCM/NEMO/OPA_SRC/SOL/solsor.F90
r2715 r3294 22 22 USE lib_mpp ! distributed memory computing 23 23 USE lbclnk ! ocean lateral boundary conditions (or mpp link) 24 USE lib_fortran 24 USE lib_fortran ! Fortran routines library 25 USE wrk_nemo ! Memory allocation 26 USE timing ! Timing 25 27 26 28 IMPLICIT NONE … … 57 59 !! Beare and Stevens 1997 Ann. Geophysicae 15, 1369-1377 58 60 !!---------------------------------------------------------------------- 59 USE wrk_nemo, ONLY: wrk_in_use, wrk_not_released60 USE wrk_nemo, ONLY: ztab => wrk_2d_1 ! 2D workspace61 61 !! 62 62 INTEGER, INTENT(inout) :: kindic ! solver indicator, < 0 if the convergence is not reached: … … 66 66 INTEGER :: ishift, icount, ijmppodd, ijmppeven, ijpr2d ! local integers 67 67 REAL(wp) :: ztmp, zres, zres2 ! local scalars 68 REAL(wp), POINTER, DIMENSION(:,:) :: ztab ! 2D workspace 68 69 !!---------------------------------------------------------------------- 69 70 IF( wrk_in_use(2, 1) )THEN71 CALL ctl_stop('sol_sor: requested workspace array is unavailable') ; RETURN72 ENDIF73 70 ! 71 IF( nn_timing == 1 ) CALL timing_start('sol_sor') 72 ! 73 CALL wrk_alloc( jpi, jpj, ztab ) 74 ! 74 75 ijmppeven = MOD( nimpp+njmpp+jpr2di+jpr2dj , 2 ) 75 76 ijmppodd = MOD( nimpp+njmpp+jpr2di+jpr2dj+1 , 2 ) … … 167 168 ! ------------- 168 169 CALL lbc_lnk_e( gcx, c_solver_pt, 1. ) ! boundary conditions 169 ! 170 IF( wrk_not_released(2, 1) ) CALL ctl_stop('sol_sor: failed to release workspace array') 170 ! 171 CALL wrk_dealloc( jpi, jpj, ztab ) 172 ! 173 IF( nn_timing == 1 ) CALL timing_stop('sol_sor') 171 174 ! 172 175 END SUBROUTINE sol_sor
Note: See TracChangeset
for help on using the changeset viewer.