Ignore:
Timestamp:
2010-11-13T14:08:12+01:00 (10 years ago)
Author:
gm
Message:

v3.3beta: C1D - bug correction to compile with key_c1d

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/nemo_v3_3_beta/NEMOGCM/NEMO/OPA_SRC/step.F90

    r2370 r2382  
    3232   USE trcstp           ! passive tracer time-stepping      (trc_stp routine) 
    3333#endif 
     34#if defined key_agrif 
     35   USE agrif_opa_sponge ! Momemtum and tracers sponges 
     36#endif 
    3437   USE asminc           ! assimilation increments    (tra_asm_inc, dyn_asm_inc routines) 
    35    USE stpctl           ! time stepping control            (stp_ctl routine) 
    36    USE restart          ! ocean restart                    (rst_wri routine) 
    37    USE prtctl           ! Print control                    (prt_ctl routine) 
    38  
    39 #if defined key_agrif 
    40    USE agrif_opa_sponge ! Momemtum and tracers sponges 
    41 #endif 
    4238 
    4339   IMPLICIT NONE 
     
    116112      !  VERTICAL PHYSICS    
    117113                         CALL zdf_bfr( kstp )         ! bottom friction 
     114                          
    118115      !                                               ! Vertical eddy viscosity and diffusivity coefficients 
    119       IF( lk_zdfric     )   CALL zdf_ric    ( kstp )  ! Richardson number dependent Kz 
    120       IF( lk_zdftke     )   CALL zdf_tke    ( kstp )  ! TKE closure scheme for Kz 
    121       IF( lk_zdfgls     )   CALL zdf_gls    ( kstp )  ! GLS closure scheme for Kz 
    122       IF( lk_zdfkpp     )   CALL zdf_kpp    ( kstp )  ! KPP closure scheme for Kz 
    123       IF( lk_zdfcst     )   THEN                      ! Constant Kz (reset avt, avm[uv] to the background value) 
     116      IF( lk_zdfric  )   CALL zdf_ric( kstp )            ! Richardson number dependent Kz 
     117      IF( lk_zdftke  )   CALL zdf_tke( kstp )            ! TKE closure scheme for Kz 
     118      IF( lk_zdfgls  )   CALL zdf_gls( kstp )            ! GLS closure scheme for Kz 
     119      IF( lk_zdfkpp  )   CALL zdf_kpp( kstp )            ! KPP closure scheme for Kz 
     120      IF( lk_zdfcst  )   THEN                            ! Constant Kz (reset avt, avm[uv] to the background value) 
    124121         avt (:,:,:) = rn_avt0 * tmask(:,:,:) 
    125122         avmu(:,:,:) = rn_avm0 * umask(:,:,:) 
     
    135132      IF( lk_zdfddm .AND. .NOT. lk_zdfkpp )   & 
    136133         &               CALL zdf_ddm( kstp )         ! double diffusive mixing 
     134          
    137135                         CALL zdf_mxl( kstp )         ! mixed layer depth 
    138136 
    139                                                       ! write tke information in the restart file 
     137                                                      ! write TKE or GLS information in the restart file 
    140138      IF( lrst_oce .AND. lk_zdftke )   CALL tke_rst( kstp, 'WRITE' ) 
    141                                                       ! write gls information in the restart file 
    142139      IF( lrst_oce .AND. lk_zdfgls )   CALL gls_rst( kstp, 'WRITE' ) 
    143140      ! 
     
    156153#if defined key_traldf_c2d 
    157154      IF( lk_traldf_eiv )   CALL ldf_eiv( kstp )      ! eddy induced velocity coefficient 
    158 #  endif 
     155#endif 
    159156 
    160157      !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
     
    179176      ! Active tracers                              (ua, va used as workspace) 
    180177      !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< 
    181                              tsa(:,:,:,:) = 0.e0               ! set tracer trends to zero 
     178                             tsa(:,:,:,:) = 0.e0            ! set tracer trends to zero 
    182179 
    183180      IF(  ln_asmiau .AND. & 
     
    190187                             CALL tra_adv    ( kstp )       ! horizontal & vertical advection 
    191188      IF( n_cla == 1     )   CALL tra_cla    ( kstp )       ! Cross Land Advection (Update Hor. advection) 
    192       IF( lk_zdfkpp )        CALL tra_kpp    ( kstp )       ! KPP non-local tracer fluxes 
     189      IF( lk_zdfkpp      )   CALL tra_kpp    ( kstp )       ! KPP non-local tracer fluxes 
    193190                             CALL tra_ldf    ( kstp )       ! lateral mixing 
    194191#if defined key_agrif 
     
    213210                             CALL tra_nxt    ( kstp )       ! tracer fields at next time step 
    214211      ENDIF  
    215                              CALL tra_unswap 
     212                             CALL tra_unswap                ! udate T & S 3D arrays  (to be suppressed) 
    216213 
    217214      !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
     
    222219 
    223220      IF(  ln_asmiau .AND. & 
    224          & ln_dyninc     )     CALL dyn_asm_inc( kstp )       ! apply dynamics assimilation increment 
     221         & ln_dyninc       )   CALL dyn_asm_inc( kstp )     ! apply dynamics assimilation increment 
    225222                               CALL dyn_adv( kstp )         ! advection (vector or flux form) 
    226223                               CALL dyn_vor( kstp )         ! vorticity term including Coriolis 
     
    230227#endif 
    231228                               CALL dyn_hpg( kstp )         ! horizontal gradient of Hydrostatic pressure 
    232                                CALL dyn_bfr( kstp )           ! bottom friction    
     229                               CALL dyn_bfr( kstp )         ! bottom friction    
    233230                               CALL dyn_zdf( kstp )         ! vertical diffusion 
    234231                               CALL dyn_spg( kstp, indic )  ! surface pressure gradient 
Note: See TracChangeset for help on using the changeset viewer.