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 7869 for branches/UKMO/dev_r5518_medusa_fix_restart/NEMOGCM/NEMO/OPA_SRC/step.F90 – NEMO

Ignore:
Timestamp:
2017-04-04T12:21:15+02:00 (7 years ago)
Author:
marc
Message:

Moved GTRU & GTRV into dump and only call ZPS_HDE in TRC_INIT at start-up.

File:
1 edited

Legend:

Unmodified
Added
Removed

  • branches/UKMO/dev_r5518_medusa_fix_restart/NEMOGCM/NEMO/OPA_SRC/step.F90

    r7865 r7869  
    3131   !!   stp             : OPA system time-stepping 
    3232   !!---------------------------------------------------------------------- 
    33    USE dom_oce, ONLY: e3w_n, gdep3w_n, gdept_n 
    3433   USE step_oce         ! time stepping definition modules 
    3534   USE iom 
    … …  
    156155                         CALL eos( tsb, rhd, gdept_0(:,:,:) )               ! before in situ density 
    157156         IF( ln_zps .AND. .NOT. ln_isfcav)                               & 
    158             &            CALL zps_hde    ( kstp, jpts, tsb,              & 
    159             &                              e3w_n, gdept_n, gtsu, gtsv,   &  ! Partial steps: before horizontal gradient 
     157            &            CALL zps_hde    ( kstp, jpts, tsb, gtsu, gtsv,  &  ! Partial steps: before horizontal gradient 
    160158            &                                          rhd, gru , grv    )  ! of t, s, rd at the last ocean level 
    161159         IF( ln_zps .AND.       ln_isfcav)                               & 
    162             &            CALL zps_hde_isf( kstp, jpts, tsb,              & 
    163             &                              e3w_n, gdept_n, gdep3w_n,     & 
    164             &                              gtsu, gtsv,                   &    ! Partial steps for top cell (ISF) 
     160            &            CALL zps_hde_isf( kstp, jpts, tsb, gtsu, gtsv,  &    ! Partial steps for top cell (ISF) 
    165161            &                                          rhd, gru , grv , aru , arv , gzu , gzv , ge3ru , ge3rv ,   & 
    166162            &                                   gtui, gtvi, grui, grvi, arui, arvi, gzui, gzvi, ge3rui, ge3rvi    ) ! of t, s, rd at the first ocean level 
    … …  
    195191                            CALL eos    ( tsn, rhd, rhop, fsdept_n(:,:,:) ) ! now in situ density for hpg computation 
    196192            IF( ln_zps .AND. .NOT. ln_isfcav)                               & 
    197                &            CALL zps_hde    ( kstp, jpts, tsn,              & 
    198                &                              e3w_n, gdept_n, gtsu, gtsv,   &    ! Partial steps: before horizontal gradient 
     193               &            CALL zps_hde    ( kstp, jpts, tsn, gtsu, gtsv,  &    ! Partial steps: before horizontal gradient 
    199194               &                                          rhd, gru , grv    )  ! of t, s, rd at the last ocean level 
    200195            IF( ln_zps .AND.       ln_isfcav)                               & 
    201                &            CALL zps_hde_isf( kstp, jpts, tsn,              & 
    202                &                              e3w_n, gdept_n, gdep3w_n,     & 
    203                &                              gtsu, gtsv,                   &    ! Partial steps for top cell (ISF) 
     196               &            CALL zps_hde_isf( kstp, jpts, tsn, gtsu, gtsv,  &    ! Partial steps for top cell (ISF) 
    204197               &                                          rhd, gru , grv , aru , arv , gzu , gzv , ge3ru , ge3rv ,   & 
    205198               &                                   gtui, gtvi, grui, grvi, arui, arvi, gzui, gzvi, ge3rui, ge3rvi    ) ! of t, s, rd at the last ocean level 
    … …  
    280273                             CALL eos    ( tsa, rhd, rhop, fsdept_n(:,:,:) )  ! Time-filtered in situ density for hpg computation 
    281274            IF( ln_zps .AND. .NOT. ln_isfcav)                                & 
    282                &             CALL zps_hde    ( kstp, jpts, tsa,              & 
    283                &                               e3w_n, gdept_n, gtsu, gtsv,   &    ! Partial steps: before horizontal gradient 
     275               &             CALL zps_hde    ( kstp, jpts, tsa, gtsu, gtsv,  &    ! Partial steps: before horizontal gradient 
    284276               &                                           rhd, gru , grv    )  ! of t, s, rd at the last ocean level 
    285277            IF( ln_zps .AND.       ln_isfcav)                                & 
    286                &             CALL zps_hde_isf( kstp, jpts, tsa,              & 
    287                &                               e3w_n, gdept_n, gdep3w_n,     & 
    288                &                               gtsu, gtsv,                   &    ! Partial steps for top cell (ISF) 
     278               &             CALL zps_hde_isf( kstp, jpts, tsa, gtsu, gtsv,  &    ! Partial steps for top cell (ISF) 
    289279               &                                           rhd, gru , grv , aru , arv , gzu , gzv , ge3ru , ge3rv ,   & 
    290280               &                                    gtui, gtvi, grui, grvi, arui, arvi, gzui, gzvi, ge3rui, ge3rvi    ) ! of t, s, rd at the last ocean level 
    … …  
    294284                             CALL eos    ( tsn, rhd, rhop, fsdept_n(:,:,:) )  ! now in situ density for hpg computation 
    295285         IF( ln_zps .AND. .NOT. ln_isfcav)                                   & 
    296                &             CALL zps_hde    ( kstp, jpts, tsn,              & 
    297                &                               e3w_n, gdept_n, gtsu, gtsv,   &    ! Partial steps: before horizontal gradient 
     286               &             CALL zps_hde    ( kstp, jpts, tsn, gtsu, gtsv,  &    ! Partial steps: before horizontal gradient 
    298287               &                                           rhd, gru , grv    )  ! of t, s, rd at the last ocean level 
    299288         IF( ln_zps .AND.       ln_isfcav)                                   &  
    300                &             CALL zps_hde_isf( kstp, jpts, tsn,              & 
    301                &                               e3w_n, gdept_n, gdep3w_n,     & 
    302                &                               gtsu, gtsv,                   &    ! Partial steps for top cell (ISF) 
     289               &             CALL zps_hde_isf( kstp, jpts, tsn, gtsu, gtsv,  &    ! Partial steps for top cell (ISF) 
    303290               &                                           rhd, gru , grv , aru , arv , gzu , gzv , ge3ru , ge3rv ,   & 
    304291               &                                    gtui, gtvi, grui, grvi, arui, arvi, gzui, gzvi, ge3rui, ge3rvi    ) ! of t, s, rd at the last ocean level 
Note: See TracChangeset for help on using the changeset viewer.