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 3294 for trunk/NEMOGCM/NEMO/OPA_SRC/step.F90 – NEMO

Ignore:
Timestamp:
2012-01-28T17:44:18+01:00 (12 years ago)
Author:
rblod
Message:

Merge of 3.4beta into the trunk

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMOGCM/NEMO/OPA_SRC/step.F90

    r2715 r3294  
    2323   !!            3.3  !  2010-05  (K. Mogensen, A. Weaver, M. Martin, D. Lea) Assimilation interface 
    2424   !!             -   !  2010-10  (C. Ethe, G. Madec) reorganisation of initialisation phase + merge TRC-TRA 
     25   !!            3.4  !  2011-04  (G. Madec, C. Ethe) Merge of dtatem and dtasal 
    2526   !!---------------------------------------------------------------------- 
    2627 
     
    3536   USE agrif_opa_sponge ! Momemtum and tracers sponges 
    3637#endif 
    37    USE asminc           ! assimilation increments    (tra_asm_inc, dyn_asm_inc routines) 
    3838 
    3939   IMPLICIT NONE 
     
    9494      ! Update data, open boundaries, surface boundary condition (including sea-ice) 
    9595      !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< 
    96       IF( lk_dtatem  )   CALL dta_tem( kstp )         ! update 3D temperature data 
    97       IF( lk_dtasal  )   CALL dta_sal( kstp )         ! update 3D salinity data 
    9896                         CALL sbc    ( kstp )         ! Sea Boundary Condition (including sea-ice) 
     97      IF( lk_tide    )   CALL sbc_tide( kstp ) 
    9998      IF( lk_obc     )   CALL obc_dta( kstp )         ! update dynamic and tracer data at open boundaries 
    10099      IF( lk_obc     )   CALL obc_rad( kstp )         ! compute phase velocities at open boundaries 
    101       IF( lk_bdy     )   CALL bdy_dta_frs( kstp )     ! update dynamic and tracer data for FRS conditions (BDY) 
     100      IF( lk_bdy     )   CALL bdy_dta( kstp, time_offset=+1 ) ! update dynamic and tracer data at open boundaries 
    102101 
    103102      !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
     
    107106 
    108107      !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    109       ! Ocean physics update                (ua, va, ta, sa used as workspace) 
     108      ! Ocean physics update                (ua, va, tsa used as workspace) 
    110109      !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< 
    111110                         CALL bn2( tsb, rn2b )        ! before Brunt-Vaisala frequency 
     
    158157 
    159158      !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    160       ! diagnostics and outputs             (ua, va, ta, sa used as workspace) 
     159      ! diagnostics and outputs             (ua, va, tsa used as workspace) 
    161160      !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< 
    162161      IF( lk_floats  )   CALL flo_stp( kstp )         ! drifting Floats 
     
    164163      IF( lk_diafwb  )   CALL dia_fwb( kstp )         ! Fresh water budget diagnostics 
    165164      IF( ln_diaptr  )   CALL dia_ptr( kstp )         ! Poleward TRansports diagnostics 
     165      IF( lk_diadct  )   CALL dia_dct( kstp )         ! Transports 
    166166      IF( lk_diaar5  )   CALL dia_ar5( kstp )         ! ar5 diag 
     167      IF( lk_diaharm )   CALL dia_harm( kstp )        ! Tidal harmonic analysis 
    167168                         CALL dia_wri( kstp )         ! ocean model: outputs 
    168169 
     
    185186      IF( ln_trabbc      )   CALL tra_bbc    ( kstp )       ! bottom heat flux 
    186187      IF( lk_trabbl      )   CALL tra_bbl    ( kstp )       ! advective (and/or diffusive) bottom boundary layer scheme 
    187       IF( lk_tradmp      )   CALL tra_dmp    ( kstp )       ! internal damping trends 
     188      IF( ln_tradmp      )   CALL tra_dmp    ( kstp )       ! internal damping trends 
    188189                             CALL tra_adv    ( kstp )       ! horizontal & vertical advection 
    189190      IF( lk_zdfkpp      )   CALL tra_kpp    ( kstp )       ! KPP non-local tracer fluxes 
    190191                             CALL tra_ldf    ( kstp )       ! lateral mixing 
    191192#if defined key_agrif 
    192                              CALL tra_unswap 
    193193      IF(.NOT. Agrif_Root()) CALL Agrif_Sponge_tra          ! tracers sponge 
    194                              CALL tra_swap 
    195194#endif 
    196195                             CALL tra_zdf    ( kstp )       ! vertical mixing and after tracer fields 
     
    210209                             CALL tra_nxt( kstp )                ! tracer fields at next time step 
    211210      ENDIF  
    212                              CALL tra_unswap                ! udate T & S 3D arrays  (to be suppressed) 
    213  
    214       !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    215       ! Dynamics                                    (ta, sa used as workspace) 
     211 
     212      !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
     213      ! Dynamics                                    (tsa used as workspace) 
    216214      !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< 
    217215                               ua(:,:,:) = 0.e0             ! set dynamics trends to zero 
     
    220218      IF(  ln_asmiau .AND. & 
    221219         & ln_dyninc       )   CALL dyn_asm_inc( kstp )     ! apply dynamics assimilation increment 
     220      IF( ln_neptsimp )        CALL dyn_nept_cor( kstp )    ! subtract Neptune velocities (simplified) 
    222221                               CALL dyn_adv( kstp )         ! advection (vector or flux form) 
    223222                               CALL dyn_vor( kstp )         ! vorticity term including Coriolis 
    224223                               CALL dyn_ldf( kstp )         ! lateral mixing 
     224      IF( ln_neptsimp )        CALL dyn_nept_cor( kstp )    ! add Neptune velocities (simplified) 
    225225#if defined key_agrif 
    226226      IF(.NOT. Agrif_Root())   CALL Agrif_Sponge_dyn        ! momemtum sponge 
     
    250250 
    251251      !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    252       ! Trends                              (ua, va, ta, sa used as workspace) 
     252      ! Trends                              (ua, va, tsa used as workspace) 
    253253      !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< 
    254254      IF( nstop == 0 ) THEN                          
     
    264264      IF( lk_cpl           )   CALL sbc_cpl_snd( kstp )     ! coupled mode : field exchanges 
    265265      ! 
     266      IF( nn_timing == 1 .AND.  kstp == nit000  )   CALL timing_reset 
    266267      ! 
    267268   END SUBROUTINE stp 
Note: See TracChangeset for help on using the changeset viewer.