- Timestamp:
- 2011-03-01T20:04:06+01:00 (13 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/dev_r2586_dynamic_mem/NEMOGCM/NEMO/OPA_SRC/SOL/solsor.F90
r2633 r2636 32 32 !! NEMO/OPA 3.3 , NEMO Consortium (2010) 33 33 !! $Id$ 34 !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt)34 !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 35 35 !!---------------------------------------------------------------------- 36 37 36 CONTAINS 38 37 … … 58 57 !! Beare and Stevens 1997 Ann. Geophysicae 15, 1369-1377 59 58 !!---------------------------------------------------------------------- 60 USE wrk_nemo, ONLY: wrk_in_use, wrk_not_released61 USE wrk_nemo, ONLY: ztab => wrk_2d_159 USE wrk_nemo, ONLY: wrk_in_use, wrk_not_released 60 USE wrk_nemo, ONLY: ztab => wrk_2d_1 ! 2D workspace 62 61 !! 63 62 INTEGER, INTENT(inout) :: kindic ! solver indicator, < 0 if the convergence is not reached: 64 63 ! ! the model is stopped in step (set to zero before the call of solsor) 65 64 !! 66 INTEGER :: ji, jj, jn ! dummy loop indices 67 INTEGER :: ishift, icount 68 INTEGER :: ijmppodd, ijmppeven, ijpr2d 69 REAL(wp) :: ztmp, zres, zres2 65 INTEGER :: ji, jj, jn ! dummy loop indices 66 INTEGER :: ishift, icount, ijmppodd, ijmppeven, ijpr2d ! local integers 67 REAL(wp) :: ztmp, zres, zres2 ! local scalars 70 68 !!---------------------------------------------------------------------- 71 69 72 70 IF( wrk_in_use(2, 1) )THEN 73 CALL ctl_stop('sol_sor: requested workspace array is unavailable') 74 RETURN 75 END IF 71 CALL ctl_stop('sol_sor: requested workspace array is unavailable') ; RETURN 72 ENDIF 76 73 77 74 ijmppeven = MOD( nimpp+njmpp+jpr2di+jpr2dj , 2 ) … … 171 168 CALL lbc_lnk_e( gcx, c_solver_pt, 1. ) ! boundary conditions 172 169 ! 173 IF( wrk_not_released(2, 1) )THEN 174 CALL ctl_stop('sol_sor: failed to release workspace array') 175 END IF 170 IF( wrk_not_released(2, 1) ) CALL ctl_stop('sol_sor: failed to release workspace array') 176 171 ! 177 172 END SUBROUTINE sol_sor
Note: See TracChangeset
for help on using the changeset viewer.