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 – NEMO

Changeset 1531


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

change TKE cpp key names, see ticket: #485

Location:
trunk/NEMO
Files:
9 edited

Legend:

Unmodified
Added
Removed

  • trunk/NEMO/C1D_SRC/diawri_c1d.F90

    r1465 r1531  
    221221         CALL histdef( nid_T, "sometauy", "Wind Stress along j-axis"           , "N/m2"   ,   &  ! vtau 
    222222            &          jpi, jpj, nh_T, 1  , 1, 1  , - 99, 32, clop, zsto, zout ) 
    223 #if defined key_zdftke 
     223#if defined key_zdftke_old 
    224224         CALL histdef( nid_T, "votlsdis", " Dissipation Turbulent Lenght Scale", "m"      ,   &  ! e_dis 
    225225            &          jpi, jpj, nh_T, ipk, 1, ipk, nz_T, 32, clop, zsto, zout ) 
    … …  
    321321      CALL histwrite( nid_T, "vomecrty", it, vn            , ndim_T , ndex_T  )   ! j-current 
    322322      CALL histwrite( nid_T, "sometauy", it, vtau          , ndim_hT, ndex_hT )   ! j-wind stress 
    323 #if defined key_zdftke 
     323#if defined key_zdftke_old 
    324324      CALL histwrite( nid_T, "votlsdis", it, e_dis         , ndim_T , ndex_T )    ! Diss. Turb. lenght scale 
    325325      CALL histwrite( nid_T, "votlsmix", it, e_mix         , ndim_T , ndex_T )    ! Mixing Turb. lenght scale 

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

    r1239 r1531  
    116116 
    117117      ! Force implicit schemes 
    118       IF( lk_zdftke .OR. lk_zdftke2 .OR. lk_zdfkpp   ) nzdf = 1   ! TKE or KPP physics 
    119       IF( ln_dynldf_iso              )                 nzdf = 1   ! iso-neutral lateral physics 
    120       IF( ln_dynldf_hor .AND. ln_sco )                 nzdf = 1   ! horizontal lateral physics in s-coordinate 
     118      IF( lk_zdftke_old .OR. lk_zdftke .OR. lk_zdfkpp )  nzdf = 1   ! TKE or KPP physics 
     119      IF( ln_dynldf_iso                               )   nzdf = 1   ! iso-neutral lateral physics 
     120      IF( ln_dynldf_hor .AND. ln_sco                  )   nzdf = 1   ! horizontal lateral physics in s-coordinate 
    121121 
    122122      IF( lk_esopa )    nzdf = -1                   ! Esopa key: All schemes used 

  • trunk/NEMO/OPA_SRC/TRA/trazdf.F90

    r1438 r1531  
    139139 
    140140      ! Force implicit schemes 
    141       IF( lk_zdftke .OR. lk_zdftke2 .OR. lk_zdfkpp   ) nzdf = 1      ! TKE or KPP physics 
    142       IF( ln_traldf_iso              )                 nzdf = 1      ! iso-neutral lateral physics 
    143       IF( ln_traldf_hor .AND. ln_sco )                 nzdf = 1      ! horizontal lateral physics in s-coordinate 
     141      IF( lk_zdftke_old .OR. lk_zdftke .OR. lk_zdfkpp )  nzdf = 1      ! TKE or KPP physics 
     142      IF( ln_traldf_iso                               )   nzdf = 1      ! iso-neutral lateral physics 
     143      IF( ln_traldf_hor .AND. ln_sco                  )   nzdf = 1      ! horizontal lateral physics in s-coordinate 
    144144 
    145145      IF( ln_zdfexp .AND. nzdf == 1 )   THEN 

  • trunk/NEMO/OPA_SRC/ZDF/zdfini.F90

    r1239 r1531  
    101101         ioptio = ioptio+1 
    102102      ENDIF 
    103       IF( lk_zdftke ) THEN 
     103      IF( lk_zdftke_old ) THEN 
    104104         IF(lwp) WRITE(numout,*) '             TKE dependent eddy coefficients' 
    105105         ioptio = ioptio+1 
    106106      ENDIF 
    107       IF( lk_zdftke2 ) THEN 
     107      IF( lk_zdftke ) THEN 
    108108         IF(lwp) WRITE(numout,*) '             TKE2 dependent eddy coefficients' 
    109109         ioptio = ioptio+1 
    … …  
    128128         ioptio = ioptio+1 
    129129      ENDIF 
    130       IF( lk_zdftke .OR. lk_zdftke2 ) THEN 
     130      IF( lk_zdftke_old .OR. lk_zdftke ) THEN 
    131131         IF(lwp) WRITE(numout,*) '             use the 1.5 turbulent closure' 
    132132      ENDIF 
    … …  
    140140      IF ( ioptio > 1 .AND. .NOT. lk_esopa ) & 
    141141           CALL ctl_stop( ' chose between ln_zdfnpc', '           and ln_zdfevd' ) 
    142       IF( ioptio == 0 .AND. .NOT.( lk_zdftke .OR. lk_zdftke2 .OR. lk_zdfkpp ) ) & 
     142      IF( ioptio == 0 .AND. .NOT.( lk_zdftke_old .OR. lk_zdftke .OR. lk_zdfkpp ) ) & 
    143143         CALL ctl_stop( ' except for TKE sor KPP physics, a convection scheme is', & 
    144144         &              ' required: ln_zdfevd or ln_zdfnpc logicals' ) 

  • trunk/NEMO/OPA_SRC/ZDF/zdftke.F90

    r1438 r1531  
    1 MODULE zdftke 
     1MODULE zdftke_old 
    22   !!====================================================================== 
    3    !!                       ***  MODULE  zdftke  *** 
     3   !!                       ***  MODULE  zdftke_old  *** 
    44   !! Ocean physics:  vertical mixing coefficient computed from the tke  
    55   !!                 turbulent closure parameterization 
    … …  
    2424   !!              -   !  2008-06  (G. Madec)  style + DOCTOR name for namelist parameters 
    2525   !!---------------------------------------------------------------------- 
    26 #if defined key_zdftke   ||   defined key_esopa 
     26#if defined key_zdftke_old   ||   defined key_esopa 
    2727   !!---------------------------------------------------------------------- 
    28    !!   'key_zdftke'                                   TKE vertical physics 
     28   !!   'key_zdftke_old'                               TKE vertical physics 
    2929   !!---------------------------------------------------------------------- 
    3030   !!---------------------------------------------------------------------- 
    31    !!   zdf_tke      : update momentum and tracer Kz from a tke scheme 
     31   !!   zdf_tke_old  : update momentum and tracer Kz from a tke scheme 
    3232   !!   zdf_tke_init : initialization, namelist read, and parameters control 
    3333   !!   tke_rst      : read/write tke restart in ocean restart file 
    … …  
    4848   PRIVATE 
    4949 
    50    PUBLIC   zdf_tke    ! routine called in step module 
    51  
    52    LOGICAL , PUBLIC, PARAMETER              ::   lk_zdftke = .TRUE.  !: TKE vertical mixing flag 
     50   PUBLIC   zdf_tke_old   ! routine called in step module 
     51 
     52   LOGICAL , PUBLIC, PARAMETER              ::   lk_zdftke_old = .TRUE.  !: TKE vertical mixing flag 
    5353   REAL(wp), PUBLIC                         ::   eboost              !: multiplicative coeff of the shear product. 
    5454   REAL(wp), PUBLIC, DIMENSION(jpi,jpj,jpk) ::   en                  !: now turbulent kinetic energy 
    … …  
    102102CONTAINS 
    103103 
    104    SUBROUTINE zdf_tke( kt ) 
     104   SUBROUTINE zdf_tke_old( kt ) 
    105105      !!---------------------------------------------------------------------- 
    106       !!                   ***  ROUTINE zdf_tke  *** 
     106      !!                   ***  ROUTINE zdf_tke_old  *** 
    107107      !! 
    108108      !! ** Purpose :   Compute the vertical eddy viscosity and diffusivity 
    … …  
    704704      ENDIF 
    705705      ! 
    706    END SUBROUTINE zdf_tke 
     706   END SUBROUTINE zdf_tke_old 
    707707 
    708708 
    … …  
    754754      IF(lwp) THEN 
    755755         WRITE(numout,*) 
    756          WRITE(numout,*) 'zdf_tke_init : tke turbulent closure scheme' 
     756         WRITE(numout,*) 'zdf_tke_init : tke turbulent closure scheme (old scheme)' 
    757757         WRITE(numout,*) '~~~~~~~~~~~~' 
    758758         WRITE(numout,*) '          Namelist namtke : set tke mixing parameters' 
    … …  
    920920              en (:,:,:) = rn_emin * tmask(:,:,:) 
    921921              DO jit = 2, nn_itke + 1 
    922                  CALL zdf_tke( jit ) 
     922                 CALL zdf_tke_old( jit ) 
    923923              END DO 
    924924           ENDIF 
    … …  
    936936   !!   Dummy module :                                        NO TKE scheme 
    937937   !!---------------------------------------------------------------------- 
    938    LOGICAL, PUBLIC, PARAMETER ::   lk_zdftke = .FALSE.   !: TKE flag 
     938   LOGICAL, PUBLIC, PARAMETER ::   lk_zdftke_old = .FALSE.   !: TKE flag 
    939939CONTAINS 
    940    SUBROUTINE zdf_tke( kt )          ! Empty routine 
    941       WRITE(*,*) 'zdf_tke: You should not have seen this print! error?', kt 
    942    END SUBROUTINE zdf_tke 
     940   SUBROUTINE zdf_tke_old( kt )          ! Empty routine 
     941      WRITE(*,*) 'zdf_tke_old: You should not have seen this print! error?', kt 
     942   END SUBROUTINE zdf_tke_old 
    943943#endif 
    944944 
    945945   !!====================================================================== 
    946 END MODULE zdftke 
     946END MODULE zdftke_old 

  • 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 

  • trunk/NEMO/OPA_SRC/restart.F90

    r1483 r1531  
    150150      !! ** Purpose :   Read files for restart (format fixed by jprstlib ) 
    151151      !!  
    152       !! ** Method  :   Read the previous fields on the NetCDF/binary file  
    153       !!      the first record indicates previous characterics 
    154       !!      after control with the present run, we read : 
    155       !!      - prognostic variables on the second record 
    156       !!      - elliptic solver arrays  
    157       !!      - barotropic stream function arrays ("key_dynspg_rl" defined) 
    158       !!        or free surface arrays  
    159       !!      - tke arrays (lk_zdftke=T .OR. lk_zdftke2=T) 
    160       !!      for this last three records,  the previous characteristics  
    161       !!      could be different with those used in the present run.  
     152      !! ** Method  :   Read in restart.nc file fields which are necessary for restart 
    162153      !!---------------------------------------------------------------------- 
    163154      REAL(wp) ::   zrdt, zrdttra1 

  • trunk/NEMO/OPA_SRC/step.F90

    r1528 r1531  
    8787   USE zdftmx          ! tide-induced vertical mixing     (zdf_tmx routine) 
    8888   USE zdfbfr          ! bottom friction                  (zdf_bfr routine) 
     89   USE zdftke_old      ! old TKE vertical mixing      (zdf_tke_old routine) 
    8990   USE zdftke          ! TKE vertical mixing              (zdf_tke routine) 
    90    USE zdftke2         ! TKE2 vertical mixing             (zdf_tke2 routine) 
    9191   USE zdfkpp          ! KPP vertical mixing              (zdf_kpp routine) 
    9292   USE zdfddm          ! double diffusion mixing          (zdf_ddm routine) 
    … …  
    200200      !  VERTICAL PHYSICS    
    201201      !                                               ! Vertical eddy viscosity and diffusivity coefficients 
    202       IF( lk_zdfric  )   CALL zdf_ric( kstp )         ! Richardson number dependent Kz 
    203       IF( lk_zdftke  )   CALL zdf_tke ( kstp )        ! TKE  closure scheme for Kz 
    204       IF( lk_zdftke2 )   CALL zdf_tke2( kstp )        ! TKE2 closure scheme for Kz 
    205       IF( lk_zdfkpp  )   CALL zdf_kpp( kstp )         ! KPP closure scheme for Kz 
    206       IF( lk_zdfcst  )   THEN                         ! Constant Kz (reset avt, avm[uv] to the background value) 
     202      IF( lk_zdfric     )   CALL zdf_ric    ( kstp )  ! Richardson number dependent Kz 
     203      IF( lk_zdftke_old )   CALL zdf_tke_old( kstp )  ! TKE closure scheme for Kz (old scheme) 
     204      IF( lk_zdftke     )   CALL zdf_tke    ( kstp )  ! TKE closure scheme for Kz 
     205      IF( lk_zdfkpp     )   CALL zdf_kpp    ( kstp )  ! KPP closure scheme for Kz 
     206      IF( lk_zdfcst     )   THEN                      ! Constant Kz (reset avt, avm[uv] to the background value) 
    207207         avt (:,:,:) = avt0 * tmask(:,:,:) 
    208208         avmu(:,:,:) = avm0 * umask(:,:,:) 
    … …  
    223223                         CALL zdf_mxl( kstp )         ! mixed layer depth 
    224224 
    225                                                       ! write tke2 information in the restart file 
    226       IF( lrst_oce .AND. lk_zdftke2 )   CALL tke2_rst( kstp, 'WRITE' ) 
     225                                                      ! write tke information in the restart file 
     226      IF( lrst_oce .AND. lk_zdftke )   CALL tke_rst( kstp, 'WRITE' ) 
    227227      ! 
    228228      !  LATERAL  PHYSICS  

  • trunk/NEMO/TOP_SRC/TRP/trctrp_ctl.F90

    r1239 r1531  
    249249      ! ... vertical mixing 
    250250      ! time stepping scheme (N.B. TKE && KPP scheme => force the use of implicit scheme) 
    251 #if defined key_zdftke || defined key_zdftke2 || defined key_zdfkpp 
     251#if defined key_zdftke || defined key_zdftke_old || defined key_zdfkpp 
    252252      l_trczdf_exp = .FALSE.          ! use implicit scheme 
    253253      l_trczdf_imp = .TRUE.  
Note: See TracChangeset for help on using the changeset viewer.