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 3318 for branches/2012/dev_r3309_LOCEAN12_Ediag/NEMOGCM/NEMO/OPA_SRC/DYN/dynnept.F90 – NEMO

Ignore:
Timestamp:
2012-02-25T16:50:01+01:00 (12 years ago)
Author:
gm
Message:

Ediag branche: #927 split TRA/DYN trd computation

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2012/dev_r3309_LOCEAN12_Ediag/NEMOGCM/NEMO/OPA_SRC/DYN/dynnept.F90

    r3294 r3318  
    55   !!                 recoded version of simplest case (u*, v* only) 
    66   !!====================================================================== 
    7    !! History :  1.0  !  2007-06  (Michael Dunphy)  Modular form: - new namelist parameters 
    8    !!                                                             - horizontal diffusion for Neptune 
    9    !!                                                             - vertical diffusion for gm in momentum eqns 
    10    !!                                                             - option to use Neptune in Coriolis eqn 
     7   !! History :  1.0  !  2007-06  (Zeliang Wang, Michael Dunphy, BIO)  Original code 
     8   !!                             Modular form: - new namelist parameters 
     9   !!                                           - horizontal diffusion for Neptune 
     10   !!                                           - vertical diffusion for gm in momentum eqns 
     11   !!                                           - option to use Neptune in Coriolis eqn 
    1112   !!                    2011-08  (Jeff Blundell, NOCS) Simplified form for temporally invariant u*, v* 
    12    !!                                               Horizontal and vertical diffusivity formulations removed 
    13    !!                                               Dynamic allocation of storage added 
    14    !!                                               Option of ramping Neptune vel. down 
    15    !!                                               to zero added in shallow depths added 
     13   !!                             Horizontal and vertical diffusivity formulations removed 
     14   !!                             Dynamic allocation of storage added 
     15   !!                             Option of ramping Neptune vel. down to zero added in shallow depths added 
     16   !!---------------------------------------------------------------------- 
     17 
    1618   !!---------------------------------------------------------------------- 
    1719   !! dynnept_alloc        : 
     
    3032   USE phycst 
    3133   USE lbclnk 
    32    USE wrk_nemo        ! Memory Allocation 
     34   USE wrk_nemo         ! Memory Allocation 
    3335 
    3436   IMPLICIT NONE 
    3537   PRIVATE 
    3638 
    37    !! * Routine accessibility 
    38    PUBLIC dyn_nept_init      ! routine called by nemogcm.F90 
    39    PUBLIC dyn_nept_cor       ! routine called by step.F90 
    40    !! dynnept_alloc()         is called only by dyn_nept_init, within this module 
    41    !! dyn_nept_div_cur_init   is called only by dyn_nept_init, within this module 
    42    !! dyn_nept_vel            is called only by dyn_nept_cor,  within this module 
    43  
    44    !! * Shared module variables 
    45    REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:)    :: zunep, zvnep  ! Neptune u and v 
    46    REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:)  :: zhdivnep      ! hor. div for Neptune vel. 
    47    REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:)  :: zmrotnep      ! curl for Neptune vel. 
    48  
    49  
    50    !! * Namelist namdyn_nept variables 
     39   PUBLIC   dyn_nept_init   ! routine called by nemogcm.F90 
     40   PUBLIC   dyn_nept_cor    ! routine called by step.F90 
     41 
     42   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:)   :: zunep, zvnep   ! Neptune u and v 
     43   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: zhdivnep       ! hor. div for Neptune vel. 
     44   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: zmrotnep       ! curl for Neptune vel. 
     45 
     46   !                                                 !!* Namelist namdyn_nept variables 
    5147   LOGICAL, PUBLIC  ::  ln_neptsimp        = .FALSE.  ! yes/no simplified neptune 
    5248 
     
    6056   REAL(wp)         ::  rn_htrmax          =  200.0   ! max. depth of transition range 
    6157 
    62    !! * Module variables 
    63  
    64  
    6558   !! * Substitutions 
    6659#  include "vectopt_loop_substitute.h90" 
    6760#  include "domzgr_substitute.h90" 
    6861   !!---------------------------------------------------------------------- 
    69    !!   OPA 9.0 , implemented by Bedford Institute of Oceanography 
    70    !!---------------------------------------------------------------------- 
    71  
     62   !! NEMO/OPA 3.3 , NEMO Consortium (2011) 
     63   !! $Id: dynadv_cen2.F90 3316 2012-02-21 16:00:02Z gm $ 
     64   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
     65   !!---------------------------------------------------------------------- 
    7266 CONTAINS 
    7367 
     
    9084      !!                and compute the arrays zunep and zvnep 
    9185      !! 
    92       !! ** Method  :   zunep = 
    93       !!                zvnep = 
    94       !! 
    95       !! ** History :  1.0  !   07-05  (Zeliang Wang)   Original code for zunep, zvnep 
    96       !!               1.1  !   07-06  (Michael Dunphy) namelist and  initialisation 
    97       !!               2.0  ! 2011-07  (Jeff Blundell, NOCS) 
    98       !!                    ! Simplified form for temporally invariant u*, v* 
    99       !!                    ! Horizontal and vertical diffusivity formulations removed 
    100       !!                    ! Includes optional tapering-off in shallow depths 
     86      !! ** Method  :   Simplified form for temporally invariant u*, v* 
     87      !!                Horizontal and vertical diffusivity formulations removed 
     88      !!                Includes optional tapering-off in shallow depths 
    10189      !!---------------------------------------------------------------------- 
    10290      USE iom 
    103       !! 
     91      ! 
    10492      INTEGER  ::   ji, jj, jk    ! dummy loop indices 
    10593      REAL(wp) :: unemin,unemax,vnemin,vnemax   ! extrema of (u*, v*) fields 
     
    122110      ! ========================================== 
    123111 
    124 !!    WRITE(numout,*) ' start dynnept namelist' 
    125 !!    CALL FLUSH(numout) 
    126112      REWIND( numnam )                  ! Read Namelist namdyn_nept:  Simplified Neptune 
    127113      READ  ( numnam, namdyn_nept ) 
    128 !!    WRITE(numout,*) ' dynnept namelist done' 
    129 !!    CALL FLUSH(numout) 
    130114 
    131115      IF(lwp) THEN                      ! Control print 
     
    172156 
    173157      IF( .not. ln_rstart ) THEN      ! If restarting, these arrays are read from the restart file 
    174          zhdivnep(:,:,:) = 0.0_wp 
    175          zmrotnep(:,:,:) = 0.0_wp 
     158         zhdivnep(:,:,:) = 0._wp 
     159         zmrotnep(:,:,:) = 0._wp 
    176160      END IF 
    177161 
     
    361345      !! ** Action  : - compute zhdivnep, the hor. divergence of (u*, v*) 
    362346      !!              - compute zmrotnep, the rel. vorticity  of (u*, v*) 
    363       !! 
    364       !! History :  OPA  ! 1987-06  (P. Andrich, D. L Hostis)  Original code 
    365       !!            4.0  ! 1991-11  (G. Madec) 
    366       !!            6.0  ! 1993-03  (M. Guyon)  symetrical conditions 
    367       !!            7.0  ! 1996-01  (G. Madec)  s-coordinates 
    368       !!            8.0  ! 1997-06  (G. Madec)  lateral boundary cond., lbc 
    369       !!            8.1  ! 1997-08  (J.M. Molines)  Open boundaries 
    370       !!            8.2  ! 2000-03  (G. Madec)  no slip accurate 
    371       !!  NEMO      1.0  ! 2002-09  (G. Madec, E. Durand)  Free form, F90 
    372       !!             -   ! 2005-01  (J. Chanut) Unstructured open boundaries 
    373       !!             -   ! 2003-08  (G. Madec)  merged of cur and div, free form, F90 
    374       !!             -   ! 2005-01  (J. Chanut, A. Sellar) unstructured open boundaries 
    375       !!            3.3  ! 2010-09  (D.Storkey and E.O'Dea) bug fixes for BDY module 
    376       !!                 ! 2011-06  (Jeff Blundell, NOCS) Adapt code from divcur.F90 
    377       !!                 !           to compute Neptune effect fields only 
    378347      !!---------------------------------------------------------------------- 
    379348      ! 
     
    517486 
    518487   SUBROUTINE dyn_nept_smooth_vel( htold, htnew, ld_option ) 
    519  
    520488      !!---------------------------------------------------------------------- 
    521489      !!                  ***  ROUTINE dyn_nept_smooth_vel  *** 
     
    607575   END SUBROUTINE dyn_nept_smooth_vel 
    608576 
     577   !!============================================================================== 
    609578END MODULE dynnept 
Note: See TracChangeset for help on using the changeset viewer.