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

Changeset 2329


Ignore:
Timestamp:
2010-10-28T12:21:58+02:00 (13 years ago)
Author:
gm
Message:

v3.3beta: Suppress old keys (key_diaspr, key_flx..., key_vectopt_memory) & phasing of zdfgls interface

Location:
branches/nemo_v3_3_beta/NEMOGCM/NEMO/OPA_SRC
Files:
7 edited

Legend:

Unmodified
Added
Removed

  • branches/nemo_v3_3_beta/NEMOGCM/NEMO/OPA_SRC/DIA/diawri.F90

    r2287 r2329  
    44   !! Ocean diagnostics :  write ocean output files 
    55   !!===================================================================== 
    6  
    7    !!---------------------------------------------------------------------- 
    8    !! * Modules used 
     6   !! History :  OPA  ! 1991-03  (M.-A. Foujols)  Original code 
     7   !!            4.0  ! 1991-11  (G. Madec) 
     8   !!                 ! 1992-06  (M. Imbard)  correction restart file 
     9   !!                 ! 1992-07  (M. Imbard)  split into diawri and rstwri 
     10   !!                 ! 1993-03  (M. Imbard)  suppress writibm 
     11   !!                 ! 1998-01  (C. Levy)  NETCDF format using ioipsl INTERFACE 
     12   !!                 ! 1999-02  (E. Guilyardi)  name of netCDF files + variables 
     13   !!            8.2  ! 2000-06  (M. Imbard)  Original code (diabort.F) 
     14   !!   NEMO     1.0  ! 2002-06  (A.Bozec, E. Durand)  Original code (diainit.F) 
     15   !!             -   ! 2002-09  (G. Madec)  F90: Free form and module 
     16   !!             -   ! 2002-12  (G. Madec)  merge of diabort and diainit, F90 
     17   !!                 ! 2005-11  (V. Garnier) Surface pressure gradient organization 
     18   !!            3.2  ! 2008-11  (B. Lemaire) creation from old diawri 
     19   !!---------------------------------------------------------------------- 
     20 
     21   !!---------------------------------------------------------------------- 
     22   !!   dia_wri       : create the standart output files 
     23   !!   dia_wri_state : create an output NetCDF file for a single instantaeous ocean state and forcing fields 
     24   !!---------------------------------------------------------------------- 
    925   USE oce             ! ocean dynamics and tracers  
    1026   USE dom_oce         ! ocean space and time domain 
    … …  
    3652   PRIVATE 
    3753 
    38    !! * Accessibility 
    39    PUBLIC dia_wri                 ! routines called by step.F90 
    40    PUBLIC dia_wri_state 
    41  
    42    !! * Module variables 
    43    INTEGER ::   & 
    44       nid_T, nz_T, nh_T, ndim_T, ndim_hT,      &   ! grid_T file 
    45       nid_U, nz_U, nh_U, ndim_U, ndim_hU,      &   ! grid_U file 
    46       nid_V, nz_V, nh_V, ndim_V, ndim_hV,      &   ! grid_V file 
    47       nid_W, nz_W, nh_W,                       &   ! grid_W file 
    48       ndex(1)                                      ! ??? 
    49    INTEGER, DIMENSION(jpi*jpj) ::   & 
    50       ndex_hT, ndex_hU, ndex_hV 
    51    INTEGER, DIMENSION(jpi*jpj*jpk) ::   & 
    52       ndex_T, ndex_U, ndex_V 
     54   PUBLIC   dia_wri                 ! routines called by step.F90 
     55   PUBLIC   dia_wri_state 
     56 
     57   INTEGER ::   nid_T, nz_T, nh_T, ndim_T, ndim_hT   ! grid_T file 
     58   INTEGER ::   nid_U, nz_U, nh_U, ndim_U, ndim_hU   ! grid_U file 
     59   INTEGER ::   nid_V, nz_V, nh_V, ndim_V, ndim_hV   ! grid_V file 
     60   INTEGER ::   nid_W, nz_W, nh_W                    ! grid_W file 
     61   INTEGER ::   ndex(1)                              ! ??? 
     62   INTEGER, DIMENSION(jpi*jpj)     ::   ndex_hT, ndex_hU, ndex_hV 
     63   INTEGER, DIMENSION(jpi*jpj*jpk) ::   ndex_T, ndex_U, ndex_V 
    5364 
    5465   !! * Substitutions 
    … …  
    5869   !!---------------------------------------------------------------------- 
    5970   !! NEMO/OPA 3.3 , NEMO Consortium (2010) 
    60    !! $Id$ 
    61    !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
    62    !!---------------------------------------------------------------------- 
    63  
     71   !! $Id $ 
     72   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
     73   !!---------------------------------------------------------------------- 
    6474CONTAINS 
    6575 
    6676#if defined key_dimgout 
    6777   !!---------------------------------------------------------------------- 
    68    !!   dia_wri       : create the dimg direct access output file (mpp) 
     78   !!   'key_dimgout'                                      DIMG output file 
    6979   !!---------------------------------------------------------------------- 
    7080#   include "diawri_dimg.h90" 
    … …  
    7484   !!   Default option                                   NetCDF output file 
    7585   !!---------------------------------------------------------------------- 
    76    !!   dia_wri       : create the standart NetCDF output files 
    77    !!   dia_wri_state : create an output NetCDF file for a single 
    78    !!                   instantaeous ocean state and forcing fields 
    79    !!---------------------------------------------------------------------- 
    8086# if defined key_iomput 
     87   !!---------------------------------------------------------------------- 
     88   !!   'key_iomput'                                        use IOM library 
     89   !!---------------------------------------------------------------------- 
    8190   SUBROUTINE dia_wri( kt ) 
    8291      !!--------------------------------------------------------------------- 
    … …  
    8796      !! 
    8897      !! ** Method  :  use iom_put 
    89       !! 
    90       !! History : 
    91       !!   3.2  !  05-11  (B. Lemaire) creation from old diawri 
    9298      !!---------------------------------------------------------------------- 
    9399      USE oce, ONLY :   z3d => ta   ! use ta as 3D workspace 
    … …  
    168174         CALL iom_put( "v_heattr", z2d )                  !  heat transport in i-direction 
    169175      ENDIF 
    170  
     176      ! 
    171177   END SUBROUTINE dia_wri 
    172178 
    173179#else 
     180   !!---------------------------------------------------------------------- 
     181   !!   Default option                                  use IOIPSL  library 
     182   !!---------------------------------------------------------------------- 
     183 
    174184   SUBROUTINE dia_wri( kt ) 
    175185      !!--------------------------------------------------------------------- 
    … …  
    183193      !!      At each time step call histdef to compute the mean if ncessary 
    184194      !!      Each nwrite time step, output the instantaneous or mean fields 
     195      !!---------------------------------------------------------------------- 
     196      INTEGER, INTENT( in ) ::   kt      ! ocean time-step index 
    185197      !! 
    186       !! History : 
    187       !!        !  91-03  (M.-A. Foujols)  Original code 
    188       !!        !  91-11  (G. Madec) 
    189       !!        !  92-06  (M. Imbard)  correction restart file 
    190       !!        !  92-07  (M. Imbard)  split into diawri and rstwri 
    191       !!        !  93-03  (M. Imbard)  suppress writibm 
    192       !!        !  98-01  (C. Levy)  NETCDF format using ioipsl INTERFACE 
    193       !!        !  99-02  (E. Guilyardi)  name of netCDF files + variables 
    194       !!   8.5  !  02-09  (G. Madec)  F90: Free form and module 
    195       !!   9.0  !  05-11  (V. Garnier) Surface pressure gradient organization 
    196       !!---------------------------------------------------------------------- 
    197       !! * Arguments 
    198       INTEGER, INTENT( in ) ::   kt      ! ocean time-step index 
    199  
    200       !! * Local declarations 
    201       LOGICAL ::   ll_print = .FALSE.    ! =T print and flush numout 
    202       CHARACTER (len=40) ::           & 
    203          clhstnam, clop, clmx            ! temporary names 
    204       INTEGER ::   inum = 11             ! temporary logical unit 
    205       INTEGER ::   & 
    206          iimi, iima, ipk, it, itmod,  &  ! temporary integers 
    207          ijmi, ijma                      !    "          " 
    208       REAL(wp) ::   & 
    209          zsto, zout, zmax,            &  ! temporary scalars 
    210          zjulian, zdt                    !    "         " 
    211       REAL(wp), DIMENSION(jpi,jpj) :: & 
    212          zw2d                            ! temporary workspace 
     198      LOGICAL ::   ll_print = .FALSE.                        ! =T print and flush numout 
     199      CHARACTER (len=40) ::   clhstnam, clop, clmx           ! local names 
     200      INTEGER  ::   inum = 11                                ! temporary logical unit 
     201      INTEGER  ::   iimi, iima, ipk, it, itmod, ijmi, ijma   ! local integers 
     202      REAL(wp) ::   zsto, zout, zmax, zjulian, zdt           ! local scalars 
     203      REAL(wp), DIMENSION(jpi,jpj) ::   zw2d                 ! 2D workspace 
    213204      !!---------------------------------------------------------------------- 
    214205      ! 
    … …  
    386377            &          jpi, jpj, nh_T, 1  , 1, 1  , -99 , 32, clop, zsto, zout ) 
    387378         CALL histdef( nid_T, "sosafldp", "Surface salt flux: Damping"         , "Kg/m2/s",   &  ! erp * sn 
    388             &          jpi, jpj, nh_T, 1  , 1, 1  , -99 , 32, clop, zsto, zout ) 
    389 #endif 
    390 #if defined key_diaspr 
    391          CALL histdef( nid_T, "sosurfps", "Surface Pressure"                   , "cm"     ,   &  ! sp 
    392379            &          jpi, jpj, nh_T, 1  , 1, 1  , -99 , 32, clop, zsto, zout ) 
    393380#endif 
    … …  
    524511      CALL histwrite( nid_T, "sosafldp", it, zw2d          , ndim_hT, ndex_hT )   ! salt flux damping 
    525512#endif 
    526 #if defined key_diaspr 
    527       CALL histwrite( nid_T, "sosurfps", it, gps           , ndim_hT, ndex_hT )   ! surface pressure 
    528 #endif 
    529          zw2d(:,:) = FLOAT( nmln(:,:) ) * tmask(:,:,1) 
     513      zw2d(:,:) = FLOAT( nmln(:,:) ) * tmask(:,:,1) 
    530514      CALL histwrite( nid_T, "sobowlin", it, zw2d          , ndim_hT, ndex_hT )   ! ??? 
    531515 
    … …  
    586570         CALL histclo( nid_W ) 
    587571      ENDIF 
    588  
     572      ! 
    589573   END SUBROUTINE dia_wri 
    590574# endif 
    … …  
    604588      !!      File 'output.init.nc'  is created if ninist = 1 (namelist) 
    605589      !!      File 'output.abort.nc' is created in case of abnormal job end 
    606       !! 
    607       !! History : 
    608       !!   8.2  !  00-06  (M. Imbard)  Original code (diabort.F) 
    609       !!   8.5  !  02-06  (A.Bozec, E. Durand)  Original code (diainit.F) 
    610       !!   9.0  !  02-12  (G. Madec)  merge of diabort and diainit, F90 
    611       !!    "   !  05-11  (V. Garnier) Surface pressure gradient organization 
    612590      !!---------------------------------------------------------------------- 
    613       !! * Arguments 
    614591      CHARACTER (len=* ), INTENT( in ) ::   cdfile_name      ! name of the file created 
    615592      INTEGER           , INTENT( in ) ::   kt               ! ocean time-step index 
    616  
    617       !! * Local declarations 
     593      !!  
    618594      CHARACTER (len=32) :: clname 
    619595      CHARACTER (len=40) :: clop 
    620       INTEGER  ::   & 
    621          id_i , nz_i, nh_i        
    622       INTEGER, DIMENSION(1) ::   & 
    623          idex             ! temprary workspace 
    624       REAL(wp) ::   & 
    625          zsto, zout, zmax,   & 
    626          zjulian, zdt 
     596      INTEGER  ::   id_i , nz_i, nh_i        
     597      INTEGER, DIMENSION(1) ::   idex             ! local workspace 
     598      REAL(wp) ::   zsto, zout, zmax, zjulian, zdt 
    627599      !!---------------------------------------------------------------------- 
    628600 

  • branches/nemo_v3_3_beta/NEMOGCM/NEMO/OPA_SRC/DIA/diawri_dimg.h90

    r2281 r2329  
    22  !!                        ***  diawri_dimg.h90  *** 
    33  !!---------------------------------------------------------------------- 
    4   !!   OPA 9.0 , LOCEAN-IPSL (2005)  
    5   !! $Id$ 
    6   !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt  
     4  !! NEMO/OPA 3.3 , NEMO Consortium (2010) 
     5  !! $Id $ 
     6  !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    77  !!---------------------------------------------------------------------- 
    88 
    9   SUBROUTINE dia_wri (kt) 
     9  SUBROUTINE dia_wri( kt ) 
    1010    !!---------------------------------------------------------------------- 
    1111    !!           *** routine dia_wri *** 
    … …  
    6060    !!  level 20: spgv(:,:)                the surface pressure gradient in Y direction. 
    6161    !!  
    62     !! History  
    63     !!      original  : 91-03 () 
    64     !!      additions : 91-11 (G. Madec) 
    65     !!      additions : 92-06 (M. Imbard) correction restart file 
    66     !!      additions : 92-07 (M. Imbard) split into diawri and rstwri 
    67     !!      additions : 93-03 (M. Imbard) suppress writibm 
    68     !!      additions : 94-12 (M. Imbard) acces direct files 
    69     !!      additions : 97-2002 ( Clipper Group ) dimg files 
    70     !!                  dec 2003 ( J.M. Molines) f90, mpp output for OPA9.0 
    71     !!   9.0  !  05-05  (S. Theetten) add emps fsalt move gps spgu spgv 2 lines below 
    72     !!   9.0  !  05-11  (V. Garnier) Surface pressure gradient organization 
     62    !! History:  OPA  ! 1997-02 ( Clipper Group ) dimg files 
     63    !!            -   ! 2003-12 ( J.M. Molines) f90, mpp output for OPA9.0 
     64    !!   NEMO    1.0  ! 2005-05  (S. Theetten) add emps fsalt move gps spgu spgv 2 lines below 
     65    !!            -   ! 2005-11  (V. Garnier) Surface pressure gradient organization 
    7366    !!---------------------------------------------------------------------- 
    74     !! * modules used 
    7567    USE lib_mpp 
    76  
    77     !! * Arguments 
     68    !! 
    7869    INTEGER ,INTENT(in) :: kt 
    79  
    80     !! * local declarations 
     70    !! 
    8171    INTEGER :: inbsel, jk 
    82 !!  INTEGER :: iwrite 
    8372    INTEGER :: iyear,imon,iday 
    8473    INTEGER, SAVE :: nmoyct  
    … …  
    10190    CHARACTER(LEN=80) :: clmode 
    10291    CHARACTER(LEN= 4) :: clver 
     92    !!---------------------------------------------------------------------- 
    10393    ! 
    10494    !  Initialization 
    10595    !  --------------- 
    10696    ! 
    107 #ifdef key_diaspr 
    108     inbsel = 20 
    109 #else 
    11097    inbsel = 17 
    111 #endif 
    112 #if defined key_flx_core 
    113     inbsel = 23 
    114 #endif 
    115  
    116     IF( inbsel >  jpk) THEN 
    117        IF( lwp) WRITE(numout,*)  & 
    118             ' STOP inbsel =',inbsel,' is larger than jpk=',jpk 
     98 
     99    IF( inbsel >  jpk ) THEN 
     100       IF(lwp) WRITE(numout,*)  ' STOP inbsel =',inbsel,' is larger than jpk=',jpk 
    119101       STOP 
    120102    ENDIF 
    121  
    122103 
    123104    iyear = ndastp/10000 
    … …  
    180161       !        fsel(:,:,15) = fsel(:,:,15) + fbt(:,:) 
    181162       fsel(:,:,16) = fsel(:,:,16) + ( emps(:,:)-rnf(:,:) )  
    182 #ifdef key_diaspr    
    183        fsel(:,:,18) = fsel(:,:,18) + gps(:,:)/g  
    184 #endif 
    185163       ! 
    186164       ! Output of dynamics and tracer fields and selected fields 
    … …  
    210188          ! computed from the averaged gradients. 
    211189          ! 
    212 #ifdef key_diaspr 
    213           fsel(:,:,18)= gps(:,:)/g 
    214           fsel(:,:,19)= spgu(:,:) 
    215           fsel(:,:,20)= spgv(:,:) 
    216 #endif 
    217190          ! mask mean field with tmask except utau vtau (1,2) 
    218191          DO jk=3,inbsel 
    … …  
    256229          !         fsel(:,:,15) =  fbt(:,:) 
    257230          fsel(:,:,16) = ( emps(:,:)-rnf(:,:) ) * tmask(:,:,1)  
    258 #ifdef key_diaspr            
    259           fsel(:,:,18) =      gps(:,:) /g 
    260           fsel(:,:,19) =      spgu(:,:) 
    261           fsel(:,:,20) =      spgv(:,:) 
    262 #endif 
    263231          ! 
    264232          !         qct(:,:) = 0._wp 
    … …  
    383351    ! 
    3843529000 FORMAT(a,"_",a,"_y",i4.4,"m",i2.2,"d",i2.2,".dimgproc") 
    385  
     353    ! 
    386354  END SUBROUTINE dia_wri 

  • branches/nemo_v3_3_beta/NEMOGCM/NEMO/OPA_SRC/SBC/cpl_oasis4.F90

    r2287 r2329  
    154154      IF(lwp) WRITE(numout,*) 
    155155      
    156 #if defined key_flx_bulk_monthly || defined key_flx_bulk_daily || defined key_flx_forced_daily 
    157  
    158 #if !defined key_oasis4 
    159 ! The following is not necessarily a valid peice of checking 
    160  
    161       IF(lwp) WRITE(numout,cform_err) 
    162       IF(lwp) WRITE(numout,*) ' key_coupled and key_flx_bulk_* key_flx_forced_daily are incompatible' 
    163       nstop = nstop + 1 
    164  
    165 #endif 
    166  
    167 #endif 
    168  
    169156      REWIND( numnam ) 
    170157      READ  ( numnam, nam_mpp ) 
    … …  
    308295      IF(lwp) WRITE(numout,*) 
    309296      
    310 #if defined key_flx_bulk_monthly || defined key_flx_bulk_daily || defined key_flx_forced_daily 
    311  
    312 #if !defined key_oasis4 
    313 ! The problem with the following is that it ASSUMES we're only ever coupling to an atmosphere 
    314 ! which is not necessarily the case. Prevent this test temporarily for NEMOGAM development. 
    315  
    316       IF(lwp) WRITE(numout,cform_err) 
    317       IF(lwp) WRITE(numout,*) ' key_coupled and key_flx_bulk_... are incompatible' 
    318       nstop = nstop + 1 
    319  
    320 #endif 
    321  
    322 #endif 
    323  
    324297      ! ----------------------------------------------------------------- 
    325298      ! ... Some initialisation 

  • branches/nemo_v3_3_beta/NEMOGCM/NEMO/OPA_SRC/ZDF/zdfgls.F90

    r2299 r2329  
    55   !!                 turbulent closure parameterization 
    66   !!====================================================================== 
    7    !! History :   3.0  !  2009-09  (G. Reffray) : Original code 
     7   !! History :   3.0  ! 2009-09 (G. Reffray)  Original code 
     8   !!             3.3  ! 2010-10 (C. Bricaud)  add in the reference 
    89   !!---------------------------------------------------------------------- 
    910#if defined key_zdfgls   ||   defined key_esopa 
    1011   !!---------------------------------------------------------------------- 
    1112   !!   'key_zdfgls'                 Generic Length Scale vertical physics 
    12    !!---------------------------------------------------------------------- 
    1313   !!---------------------------------------------------------------------- 
    1414   !!   zdf_gls      : update momentum and tracer Kz from a gls scheme 
    … …  
    3333   PRIVATE 
    3434 
    35    PUBLIC   zdf_gls    ! routine called in step module 
    36    PUBLIC   gls_rst    ! routine called in step module 
     35   PUBLIC   zdf_gls        ! routine called in step module 
     36   PUBLIC   zdf_gls_init   ! routine called in step module 
     37   PUBLIC   gls_rst        ! routine called in step module 
    3738 
    3839   LOGICAL , PUBLIC, PARAMETER              ::   lk_zdfgls = .TRUE.  !: TKE vertical mixing flag 
    3940   REAL(wp), PUBLIC, DIMENSION(jpi,jpj,jpk) ::   en                  !: now turbulent kinetic energy 
    4041   REAL(wp), PUBLIC, DIMENSION(jpi,jpj,jpk) ::   mxln                !: now mixing length 
    41    REAL(wp), PUBLIC, DIMENSION(jpi,jpj,jpk) ::   zwall              !: wall function 
    42    REAL(wp), PUBLIC, DIMENSION(jpi,jpj)     ::   ustars2            !: Squared surface velocity scale at T-points 
    43    REAL(wp), PUBLIC, DIMENSION(jpi,jpj)     ::   ustarb2            !: Squared bottom  velocity scale at T-points 
     42   REAL(wp), PUBLIC, DIMENSION(jpi,jpj,jpk) ::   zwall               !: wall function 
     43   REAL(wp), PUBLIC, DIMENSION(jpi,jpj)     ::   ustars2             !: Squared surface velocity scale at T-points 
     44   REAL(wp), PUBLIC, DIMENSION(jpi,jpj)     ::   ustarb2             !: Squared bottom  velocity scale at T-points 
    4445 
    4546   !                                         !!! ** Namelist  namzdf_gls  ** 
    4647   LOGICAL  ::   ln_crban      = .FALSE.      ! =T use Craig and Banner scheme 
    47    LOGICAL  ::   ln_length_lim = .FALSE.      ! use limit on the dissipation rate understable stratification (Galperin et al., 1988) 
     48   LOGICAL  ::   ln_length_lim = .FALSE.      ! use limit on the dissipation rate under stable stratif. (Galperin et al., 1988) 
    4849   LOGICAL  ::   ln_sigpsi     = .FALSE.      ! Activate Burchard (2003) modification for k-eps closure AND wave breaking mixing 
    4950   REAL(wp) ::   rn_epsmin     = 1.e-12_wp    ! minimum value of dissipation (m2/s3) 
    … …  
    113114   !!---------------------------------------------------------------------- 
    114115   !! NEMO/OPA 3.3 , NEMO Consortium (2010) 
    115    !! $Id$ 
    116    !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
     116   !! $Id $ 
     117   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    117118   !!---------------------------------------------------------------------- 
    118  
    119119CONTAINS 
    120120 
    … …  
    150150      !!-------------------------------------------------------------------- 
    151151 
    152       IF( kt == nit000  ) CALL zdf_gls_init        ! Initialization (first time-step only) 
    153  
    154       !!-------------------------------------------------------------------- 
    155152      ! Preliminary computing 
    156153 
    … …  
    905902   END SUBROUTINE zdf_gls 
    906903 
     904 
    907905   SUBROUTINE zdf_gls_init 
    908906      !!---------------------------------------------------------------------- 
    … …  
    920918      !! 
    921919      !!---------------------------------------------------------------------- 
    922       USE dynzdf_exp 
    923       USE trazdf_exp 
    924       ! 
    925 # if defined key_vectopt_memory 
    926       INTEGER ::   ji, jj, jk   ! dummy loop indices 
    927 # else 
    928       INTEGER ::           jk   ! dummy loop indices 
    929 # endif 
    930       REAL(wp)::   zcr 
     920      INTEGER ::   jk    ! dummy loop indices 
     921      REAL(wp)::   zcr   ! local scalar 
    931922      !! 
    932923      NAMELIST/namzdf_gls/rn_emin, rn_epsmin, ln_length_lim, & 
    … …  
    12321223   END SUBROUTINE zdf_gls_init 
    12331224 
     1225 
    12341226   SUBROUTINE gls_rst( kt, cdrw ) 
    12351227     !!--------------------------------------------------------------------- 
    … …  
    12441236     INTEGER         , INTENT(in) ::   kt         ! ocean time-step 
    12451237     CHARACTER(len=*), INTENT(in) ::   cdrw       ! "READ"/"WRITE" flag 
    1246      ! 
     1238     !! 
    12471239     INTEGER ::   jit, jk   ! dummy loop indices 
    12481240     INTEGER ::   id1, id2, id3, id4, id5, id6, id7, id8 

  • branches/nemo_v3_3_beta/NEMOGCM/NEMO/OPA_SRC/opa.F90

    r2325 r2329  
    8282   !!---------------------------------------------------------------------- 
    8383   !! NEMO/OPA 3.3 , NEMO Consortium (2010) 
    84    !! $Id$ 
    85    !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
     84   !! $Id $ 
     85   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    8686   !!---------------------------------------------------------------------- 
    87  
    8887CONTAINS 
    8988 
    … …  
    276275                            CALL zdf_bfr_init      ! bottom friction 
    277276      IF( lk_zdfric     )   CALL zdf_ric_init      ! Richardson number dependent Kz 
    278       IF( lk_zdftke     )   CALL zdf_tke_init      ! TKE closure scheme for Kz 
    279       IF( lk_zdfkpp     )   CALL zdf_kpp_init      ! KPP closure scheme for Kz 
     277      IF( lk_zdftke     )   CALL zdf_tke_init      ! TKE closure scheme 
     278      IF( lk_zdfgls     )   CALL zdf_gls_init      ! GLS closure scheme 
     279      IF( lk_zdfkpp     )   CALL zdf_kpp_init      ! KPP closure scheme 
    280280      IF( lk_zdftmx     )   CALL zdf_tmx_init      ! tidal vertical mixing 
    281281      IF( lk_zdfddm .AND. .NOT. lk_zdfkpp )   &  

  • branches/nemo_v3_3_beta/NEMOGCM/NEMO/OPA_SRC/step.F90

    r2325 r2329  
    3232   USE trcstp          ! passive tracer time-stepping      (trc_stp routine) 
    3333#endif 
    34  
    35    USE zdfgls          ! generic length scale case 
    3634   USE asminc          ! assimilation increments    (tra_asm_inc, dyn_asm_inc routines) 
    3735   USE stpctl          ! time stepping control            (stp_ctl routine) 

  • branches/nemo_v3_3_beta/NEMOGCM/NEMO/OPA_SRC/step_oce.F90

    r2303 r2329  
    6868   USE zdfbfr           ! bottom friction                  (zdf_bfr routine) 
    6969   USE zdftke           ! TKE vertical mixing              (zdf_tke routine) 
     70   USE zdfgls           ! GLS vertical mixing              (zdf_gls routine) 
    7071   USE zdfkpp           ! KPP vertical mixing              (zdf_kpp routine) 
    7172   USE zdfddm           ! double diffusion mixing          (zdf_ddm routine) 
Note: See TracChangeset for help on using the changeset viewer.