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 1531 for trunk/NEMO/OPA_SRC/ZDF/zdftke2.F90 – NEMO

Ignore:
Timestamp:
2009-07-24T08:48:42+02:00 (15 years ago)
Author:
ctlod
Message:

change TKE cpp key names, see ticket: #485

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMO/OPA_SRC/ZDF/zdftke2.F90

    r1492 r1531  
    1 MODULE zdftke2 
     1MODULE zdftke 
    22   !!====================================================================== 
    3    !!                       ***  MODULE  zdftke2  *** 
     3   !!                       ***  MODULE  zdftke  *** 
    44   !! Ocean physics:  vertical mixing coefficient computed from the tke  
    55   !!                 turbulent closure parameterization 
     
    2727   !!                 !                                + cleaning of the parameters + bugs correction 
    2828   !!---------------------------------------------------------------------- 
    29 #if defined key_zdftke2   ||   defined key_esopa 
     29#if defined key_zdftke   ||   defined key_esopa 
    3030   !!---------------------------------------------------------------------- 
    31    !!   'key_zdftke2'                                  TKE vertical physics 
     31   !!   'key_zdftke'                                   TKE vertical physics 
    3232   !!---------------------------------------------------------------------- 
    33    !!   zdf_tke2      : update momentum and tracer Kz from a tke scheme 
     33   !!   zdf_tke       : update momentum and tracer Kz from a tke scheme 
    3434   !!   tke_tke       : tke time stepping: update tke at now time step (en) 
    3535   !!   tke_avn       : compute mixing length scale and deduce avm and avt 
    3636   !!   tke_init      : initialization, namelist read, and parameters control 
    37    !!   tke2_rst      : read/write tke restart in ocean restart file 
     37   !!   tke_rst       : read/write tke restart in ocean restart file 
    3838   !!---------------------------------------------------------------------- 
    3939   USE oce            ! ocean dynamics and active tracers  
     
    5353   PRIVATE 
    5454 
    55    PUBLIC   zdf_tke2   ! routine called in step module 
    56    PUBLIC   tke2_rst   ! routine called in step module 
    57  
    58    LOGICAL , PUBLIC, PARAMETER              ::   lk_zdftke2 = .TRUE.  !: TKE vertical mixing flag 
     55   PUBLIC   zdf_tke    ! routine called in step module 
     56   PUBLIC   tke_rst    ! routine called in step module 
     57 
     58   LOGICAL , PUBLIC, PARAMETER              ::   lk_zdftke = .TRUE.  !: TKE vertical mixing flag 
    5959 
    6060#if defined key_c1d 
     
    9999 
    100100   REAL(wp), DIMENSION(jpi,jpj)     ::   htau      ! depth of tke penetration (nn_htau) 
    101    REAL(wp), DIMENSION(jpi,jpj)     ::   avtb_2d   ! set in tke2_init, for other modif than ice 
     101   REAL(wp), DIMENSION(jpi,jpj)     ::   avtb_2d   ! set in tke_init, for other modif than ice 
    102102   REAL(wp), DIMENSION(jpi,jpj,jpk) ::   en        ! now turbulent kinetic energy   [m2/s2] 
    103103   REAL(wp), DIMENSION(jpi,jpj,jpk) ::   dissl     ! now mixing lenght of dissipation 
     
    114114CONTAINS 
    115115 
    116    SUBROUTINE zdf_tke2( kt ) 
     116   SUBROUTINE zdf_tke( kt ) 
    117117      !!---------------------------------------------------------------------- 
    118       !!                   ***  ROUTINE zdf_tke2  *** 
     118      !!                   ***  ROUTINE zdf_tke  *** 
    119119      !! 
    120120      !! ** Purpose :   Compute the vertical eddy viscosity and diffusivity 
     
    168168                           CALL tke_avn      ! now avt, avm, avmu, avmv 
    169169      ! 
    170    END SUBROUTINE zdf_tke2 
     170   END SUBROUTINE zdf_tke 
    171171 
    172172 
     
    627627   SUBROUTINE tke_init 
    628628      !!---------------------------------------------------------------------- 
    629       !!                  ***  ROUTINE zdf_tke2_init  *** 
     629      !!                  ***  ROUTINE tke_init  *** 
    630630      !!                      
    631631      !! ** Purpose :   Initialization of the vertical eddy diffivity and  
     
    654654      IF(lwp) THEN                    !* Control print 
    655655         WRITE(numout,*) 
    656          WRITE(numout,*) 'zdf_tke2 : tke turbulent closure scheme - initialisation' 
     656         WRITE(numout,*) 'zdf_tke : tke turbulent closure scheme - initialisation' 
    657657         WRITE(numout,*) '~~~~~~~~' 
    658658         WRITE(numout,*) '          Namelist namtke : set tke mixing parameters' 
     
    738738      dissl(:,:,:) = 1.e-12 
    739739      !                               !* read or initialize all required files  
    740       CALL tke2_rst( nit000, 'READ' ) 
     740      CALL tke_rst( nit000, 'READ' ) 
    741741      ! 
    742742   END SUBROUTINE tke_init 
    743743 
    744744 
    745    SUBROUTINE tke2_rst( kt, cdrw ) 
     745   SUBROUTINE tke_rst( kt, cdrw ) 
    746746     !!--------------------------------------------------------------------- 
    747      !!                   ***  ROUTINE ts_rst  *** 
     747     !!                   ***  ROUTINE tke_rst  *** 
    748748     !!                      
    749749     !! ** Purpose :   Read or write TKE file (en) in restart file 
     
    785785              en (:,:,:) = rn_emin * tmask(:,:,:) 
    786786              CALL tke_avn                               ! recompute avt, avm, avmu, avmv and dissl (approximation) 
    787               DO jit = nit000 + 1, nit000 + 10   ;   CALL zdf_tke2( jit )   ;   END DO 
     787              DO jit = nit000 + 1, nit000 + 10   ;   CALL zdf_tke( jit )   ;   END DO 
    788788           ENDIF 
    789789        ELSE                                   !* Start from rest 
     
    799799     ELSEIF( TRIM(cdrw) == 'WRITE' ) THEN   ! Create restart file 
    800800        !                                   ! ------------------- 
    801         IF(lwp) WRITE(numout,*) '---- tke2-rst ----' 
     801        IF(lwp) WRITE(numout,*) '---- tke-rst ----' 
    802802        CALL iom_rstput( kt, nitrst, numrow, 'en'   , en     ) 
    803803        CALL iom_rstput( kt, nitrst, numrow, 'avt'  , avt    ) 
     
    809809     ENDIF 
    810810     ! 
    811    END SUBROUTINE tke2_rst 
     811   END SUBROUTINE tke_rst 
    812812 
    813813#else 
     
    815815   !!   Dummy module :                                        NO TKE scheme 
    816816   !!---------------------------------------------------------------------- 
    817    LOGICAL, PUBLIC, PARAMETER ::   lk_zdftke2 = .FALSE.   !: TKE flag 
     817   LOGICAL, PUBLIC, PARAMETER ::   lk_zdftke = .FALSE.   !: TKE flag 
    818818CONTAINS 
    819    SUBROUTINE zdf_tke2( kt )          ! Empty routine 
    820       WRITE(*,*) 'zdf_tke2: You should not have seen this print! error?', kt 
    821    END SUBROUTINE zdf_tke2 
    822    SUBROUTINE tke2_rst( kt, cdrw ) 
     819   SUBROUTINE zdf_tke( kt )          ! Empty routine 
     820      WRITE(*,*) 'zdf_tke: You should not have seen this print! error?', kt 
     821   END SUBROUTINE zdf_tke 
     822   SUBROUTINE tke_rst( kt, cdrw ) 
    823823     CHARACTER(len=*) ::   cdrw 
    824      WRITE(*,*) 'tke2_rst: You should not have seen this print! error?', kt, cdwr 
    825    END SUBROUTINE tke2_rst 
     824     WRITE(*,*) 'tke_rst: You should not have seen this print! error?', kt, cdwr 
     825   END SUBROUTINE tke_rst 
    826826#endif 
    827827 
    828828   !!====================================================================== 
    829 END MODULE zdftke2 
     829END MODULE zdftke 
Note: See TracChangeset for help on using the changeset viewer.