New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 3294 for trunk/NEMOGCM/NEMO/OPA_SRC/SOL/solsor.F90 – NEMO

Ignore:
Timestamp:
2012-01-28T17:44:18+01:00 (12 years ago)
Author:
rblod
Message:

Merge of 3.4beta into the trunk

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMOGCM/NEMO/OPA_SRC/SOL/solsor.F90

    r2715 r3294  
    2222   USE lib_mpp         ! distributed memory computing 
    2323   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 
    2527 
    2628   IMPLICIT NONE 
     
    5759      !!                Beare and Stevens 1997 Ann. Geophysicae 15, 1369-1377 
    5860      !!---------------------------------------------------------------------- 
    59       USE wrk_nemo, ONLY:   wrk_in_use, wrk_not_released 
    60       USE wrk_nemo, ONLY:   ztab => wrk_2d_1    ! 2D workspace 
    6161      !! 
    6262      INTEGER, INTENT(inout) ::   kindic   ! solver indicator, < 0 if the convergence is not reached: 
     
    6666      INTEGER  ::   ishift, icount, ijmppodd, ijmppeven, ijpr2d   ! local integers 
    6767      REAL(wp) ::   ztmp, zres, zres2                             ! local scalars 
     68      REAL(wp), POINTER, DIMENSION(:,:) ::   ztab                 ! 2D workspace 
    6869      !!---------------------------------------------------------------------- 
    69        
    70       IF( wrk_in_use(2, 1) )THEN 
    71          CALL ctl_stop('sol_sor: requested workspace array is unavailable')   ;   RETURN 
    72       ENDIF 
    73  
     70      ! 
     71      IF( nn_timing == 1 )  CALL timing_start('sol_sor') 
     72      ! 
     73      CALL wrk_alloc( jpi, jpj, ztab ) 
     74      ! 
    7475      ijmppeven = MOD( nimpp+njmpp+jpr2di+jpr2dj   , 2 ) 
    7576      ijmppodd  = MOD( nimpp+njmpp+jpr2di+jpr2dj+1 , 2 ) 
     
    167168      !  ------------- 
    168169      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') 
    171174      ! 
    172175   END SUBROUTINE sol_sor 
Note: See TracChangeset for help on using the changeset viewer.