Ignore:
Timestamp:
2015-02-23T17:07:54+01:00 (6 years ago)
Author:
mathiot
Message:

correction of minor bug with isf + add some extra test on ln_isfcav + restore ssu/ssv

File:
1 edited

Legend:

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

    r5098 r5104  
    146146      IF( lk_ldfslp ) THEN                            ! slope of lateral mixing 
    147147                         CALL eos( tsb, rhd, gdept_0(:,:,:) )               ! before in situ density 
    148          IF( ln_zps )    CALL zps_hde    ( kstp, jpts, tsb, gtsu, gtsv,  &  ! Partial steps: before horizontal gradient 
     148         IF( ln_zps .AND. .NOT. ln_isfcav)                               & 
     149            &            CALL zps_hde    ( kstp, jpts, tsb, gtsu, gtsv,  &  ! Partial steps: before horizontal gradient 
    149150            &                                          rhd, gru , grv    )  ! of t, s, rd at the last ocean level 
    150          IF( ln_zps .AND. ln_isfcav) & 
     151         IF( ln_zps .AND.       ln_isfcav)                              & 
    151152            &            CALL zps_hde_isf( kstp, jpts, tsb, gtsu, gtsv,  &    ! Partial steps for top cell (ISF) 
    152153            &                                          rhd, gru , grv , aru , arv , gzu , gzv , ge3ru , ge3rv ,   & 
     
    180181          ! is necessary to compute momentum advection for the rhs of barotropic loop: 
    181182                            CALL eos    ( tsn, rhd, rhop, fsdept_n(:,:,:) ) ! now in situ density for hpg computation 
    182             IF( ln_zps )    CALL zps_hde    ( kstp, jpts, tsn, gtsu, gtsv,  &    ! Partial steps: before horizontal gradient 
     183            IF( ln_zps .AND. .NOT. ln_isfcav)                               & 
     184               &            CALL zps_hde    ( kstp, jpts, tsn, gtsu, gtsv,  &    ! Partial steps: before horizontal gradient 
    183185               &                                          rhd, gru , grv    )  ! of t, s, rd at the last ocean level 
    184             IF( ln_zps .AND. ln_isfcav) & 
     186            IF( ln_zps .AND.       ln_isfcav)                              & 
    185187               &            CALL zps_hde_isf( kstp, jpts, tsn, gtsu, gtsv,  &    ! Partial steps for top cell (ISF) 
    186188               &                                          rhd, gru , grv , aru , arv , gzu , gzv , ge3ru , ge3rv ,   & 
     
    259261                             CALL tra_nxt( kstp )                ! tracer fields at next time step 
    260262                             CALL eos    ( tsa, rhd, rhop, fsdept_n(:,:,:) )  ! Time-filtered in situ density for hpg computation 
    261             IF( ln_zps )     CALL zps_hde    ( kstp, jpts, tsa, gtsu, gtsv,  &    ! Partial steps: before horizontal gradient 
     263            IF( ln_zps .AND. .NOT. ln_isfcav)                                & 
     264               &             CALL zps_hde    ( kstp, jpts, tsa, gtsu, gtsv,  &    ! Partial steps: before horizontal gradient 
    262265               &                                           rhd, gru , grv    )  ! of t, s, rd at the last ocean level 
    263             IF( ln_zps .AND. ln_isfcav) & 
     266            IF( ln_zps .AND.       ln_isfcav)                                & 
    264267               &             CALL zps_hde_isf( kstp, jpts, tsa, gtsu, gtsv,  &    ! Partial steps for top cell (ISF) 
    265268               &                                           rhd, gru , grv , aru , arv , gzu , gzv , ge3ru , ge3rv ,   & 
     
    268271         IF ( .NOT. lk_dynspg_ts ) THEN                     ! eos already called in time-split case 
    269272                             CALL eos    ( tsn, rhd, rhop, fsdept_n(:,:,:) )  ! now in situ density for hpg computation 
    270          IF( ln_zps )        CALL zps_hde    ( kstp, jpts, tsn, gtsu, gtsv,  &    ! Partial steps: before horizontal gradient 
    271             &                                              rhd, gru , grv    )  ! of t, s, rd at the last ocean level 
    272          IF( ln_zps .AND. ln_isfcav) & 
     273         IF( ln_zps .AND. .NOT. ln_isfcav)                                   & 
     274               &             CALL zps_hde    ( kstp, jpts, tsn, gtsu, gtsv,  &    ! Partial steps: before horizontal gradient 
     275               &                                           rhd, gru , grv    )  ! of t, s, rd at the last ocean level 
     276         IF( ln_zps .AND.       ln_isfcav)                                   &  
    273277               &             CALL zps_hde_isf( kstp, jpts, tsn, gtsu, gtsv,  &    ! Partial steps for top cell (ISF) 
    274278               &                                           rhd, gru , grv , aru , arv , gzu , gzv , ge3ru , ge3rv ,   & 
Note: See TracChangeset for help on using the changeset viewer.