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 789 for trunk/NEMO/OPA_SRC/DYN/wzvmod.F90 – NEMO

Ignore:
Timestamp:
2008-01-11T19:04:56+01:00 (16 years ago)
Author:
rblod
Message:

Suppress jki routines and associated key_mpp_omp

File:
1 edited

Legend:

Unmodified
Added
Removed

  • trunk/NEMO/OPA_SRC/DYN/wzvmod.F90

    r719 r789  
    3636 
    3737CONTAINS 
    38  
    39 #if defined key_mpp_omp 
    40    !!---------------------------------------------------------------------- 
    41    !!   'key_mpp_omp'                                   j-k-i loop (j-slab) 
    42    !!---------------------------------------------------------------------- 
    43  
    44    SUBROUTINE wzv( kt ) 
    45       !!---------------------------------------------------------------------- 
    46       !!                    ***  ROUTINE wzv  *** 
    47       !!                      
    48       !! ** Purpose :   Compute the now vertical velocity after the array swap 
    49       !! 
    50       !! ** Method  :   Using the incompressibility hypothesis, the vertical 
    51       !!     velocity is computed by integrating the horizontal divergence  
    52       !!     from the bottom to the surface. 
    53       !!     The boundary conditions are w=0 at the bottom (no flux) and, 
    54       !!     in rigid-lid case, w=0 at the sea surface. 
    55       !! 
    56       !! ** action  :    wn array : the now vertical velocity 
    57       !!---------------------------------------------------------------------- 
    58       !! * Arguments 
    59       INTEGER, INTENT( in ) ::   kt      ! ocean time-step index 
    60  
    61       !! * Local declarations 
    62       INTEGER ::   jj, jk      ! dummy loop indices 
    63       !!---------------------------------------------------------------------- 
    64  
    65       IF( kt == nit000 ) THEN 
    66          IF(lwp) WRITE(numout,*) 
    67          IF(lwp) WRITE(numout,*) 'wzv     : vertical velocity from continuity eq.' 
    68          IF(lwp) WRITE(numout,*) '~~~~~~~                     j-k-i loops' 
    69  
    70          ! bottom boundary condition: w=0 (set once for all) 
    71          wn(:,:,jpk) = 0.e0 
    72       ENDIF 
    73  
    74       !                                                ! =============== 
    75       DO jj = 1, jpj                                   !  Vertical slab 
    76          !                                             ! =============== 
    77          ! Computation from the bottom 
    78          DO jk = jpkm1, 1, -1 
    79             wn(:,jj,jk) = wn(:,jj,jk+1) - fse3t(:,jj,jk) * hdivn(:,jj,jk) 
    80          END DO 
    81          !                                             ! =============== 
    82       END DO                                           !   End of slab 
    83       !                                                ! =============== 
    84  
    85       IF(ln_ctl)   CALL prt_ctl(tab3d_1=wn, clinfo1=' w**2 -   : ', mask1=wn) 
    86  
    87    END SUBROUTINE wzv 
    88  
    89 #else 
    90    !!---------------------------------------------------------------------- 
    91    !!   Default option                                           k-j-i loop 
    92    !!---------------------------------------------------------------------- 
    9338 
    9439   SUBROUTINE wzv( kt ) 
    … …  
    189134 
    190135   END SUBROUTINE wzv 
    191 #endif 
    192136 
    193137   !!====================================================================== 
Note: See TracChangeset for help on using the changeset viewer.