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 5189 for branches/2015/dev_r5151_UKMO_ISF/NEMOGCM/NEMO/OPA_SRC/step.F90 – NEMO

Ignore:
Timestamp:
2015-03-31T19:58:23+02:00 (9 years ago)
Author:
mathiot
Message:

ISF cleaning branch: simplification and bug correction in hpg_isf, zps_hde_isf, mixed layer definition, slope, diffusion, vertical advection and top friction

File:
1 edited

Legend:

Unmodified
Added
Removed

  • branches/2015/dev_r5151_UKMO_ISF/NEMOGCM/NEMO/OPA_SRC/step.F90

    r5147 r5189  
    150150            &                                          rhd, gru , grv    )  ! of t, s, rd at the last ocean level 
    151151         IF( ln_zps .AND.       ln_isfcav)                               & 
    152             &            CALL zps_hde_isf( kstp, jpts, tsb, gtsu, gtsv,  &    ! Partial steps for top cell (ISF) 
    153             &                                          rhd, gru , grv , aru , arv , gzu , gzv , ge3ru , ge3rv ,   & 
    154             &                                   gtui, gtvi, grui, grvi, arui, arvi, gzui, gzvi, ge3rui, ge3rvi    ) ! of t, s, rd at the first ocean level 
     152            &            CALL zps_hde_isf( kstp, jpts, tsb, gtsu, gtsv, gtui, gtvi,  &  ! Partial steps for top cell (ISF) 
     153            &                                          rhd, gru , grv , grui, grvi   )  ! of t, s, rd at the first ocean level 
    155154         IF( ln_traldf_grif ) THEN                           ! before slope for Griffies operator 
    156155                         CALL ldf_slp_grif( kstp ) 
    … …  
    181180          ! is necessary to compute momentum advection for the rhs of barotropic loop: 
    182181                            CALL eos    ( tsn, rhd, rhop, fsdept_n(:,:,:) ) ! now in situ density for hpg computation 
    183             IF( ln_zps .AND. .NOT. ln_isfcav)                               & 
    184                &            CALL zps_hde    ( kstp, jpts, tsn, gtsu, gtsv,  &    ! Partial steps: before horizontal gradient 
    185                &                                          rhd, gru , grv    )  ! of t, s, rd at the last ocean level 
    186             IF( ln_zps .AND.       ln_isfcav)                               & 
    187                &            CALL zps_hde_isf( kstp, jpts, tsn, gtsu, gtsv,  &    ! Partial steps for top cell (ISF) 
    188                &                                          rhd, gru , grv , aru , arv , gzu , gzv , ge3ru , ge3rv ,   & 
    189                &                                   gtui, gtvi, grui, grvi, arui, arvi, gzui, gzvi, ge3rui, ge3rvi    ) ! of t, s, rd at the last ocean level 
     182            IF( ln_zps .AND. .NOT. ln_isfcav )                               & 
     183               &            CALL zps_hde    ( kstp, jpts, tsn, gtsu, gtsv,   &  ! Partial steps: before horizontal gradient 
     184               &                                          rhd, gru , grv     )  ! of t, s, rd at the last ocean level 
     185            IF( ln_zps .AND.       ln_isfcav )                                          & 
     186               &            CALL zps_hde_isf( kstp, jpts, tsn, gtsu, gtsv, gtui, gtvi,  &  ! Partial steps for top cell (ISF) 
     187               &                                          rhd, gru , grv , grui, grvi   )  ! of t, s, rd at the first ocean level 
    190188 
    191189                                  ua(:,:,:) = 0.e0             ! set dynamics trends to zero 
    … …  
    266264               &                                           rhd, gru , grv    )  ! of t, s, rd at the last ocean level 
    267265            IF( ln_zps .AND.       ln_isfcav)                                & 
    268                &             CALL zps_hde_isf( kstp, jpts, tsa, gtsu, gtsv,  &    ! Partial steps for top cell (ISF) 
    269                &                                           rhd, gru , grv , aru , arv , gzu , gzv , ge3ru , ge3rv ,   & 
    270                &                                    gtui, gtvi, grui, grvi, arui, arvi, gzui, gzvi, ge3rui, ge3rvi    ) ! of t, s, rd at the last ocean level 
     266               &             CALL zps_hde_isf( kstp, jpts, tsa, gtsu, gtsv, gtui, gtvi,  &  ! Partial steps for top cell (ISF) 
     267               &                                           rhd, gru , grv , grui, grvi   )  ! of t, s, rd at the first ocean level 
    271268      ELSE                                                  ! centered hpg  (eos then time stepping) 
    272269         IF ( .NOT. lk_dynspg_ts ) THEN                     ! eos already called in time-split case 
    … …  
    276273               &                                           rhd, gru , grv    )  ! of t, s, rd at the last ocean level 
    277274         IF( ln_zps .AND.       ln_isfcav)                                   &  
    278                &             CALL zps_hde_isf( kstp, jpts, tsn, gtsu, gtsv,  &    ! Partial steps for top cell (ISF) 
    279                &                                           rhd, gru , grv , aru , arv , gzu , gzv , ge3ru , ge3rv ,   & 
    280                &                                    gtui, gtvi, grui, grvi, arui, arvi, gzui, gzvi, ge3rui, ge3rvi    ) ! of t, s, rd at the last ocean level 
     275               &             CALL zps_hde_isf( kstp, jpts, tsn, gtsu, gtsv, gtui, gtvi,  &  ! Partial steps for top cell (ISF) 
     276               &                                           rhd, gru , grv , grui, grvi   )  ! of t, s, rd at the first ocean level 
    281277         ENDIF 
    282278         IF( ln_zdfnpc   )   CALL tra_npc( kstp )                ! update after fields by non-penetrative convection 
Note: See TracChangeset for help on using the changeset viewer.