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 5446 for branches/UKMO/dev_r5021_nn_etau_revision/NEMOGCM/NEMO/OFF_SRC/dtadyn.F90 – NEMO

Ignore:
Timestamp:
2015-06-19T17:52:53+02:00 (9 years ago)
Author:
davestorkey
Message:

Update UKMO/dev_r5021_nn_etau_revision branch to revision 5442 of the trunk.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/UKMO/dev_r5021_nn_etau_revision/NEMOGCM/NEMO/OFF_SRC/dtadyn.F90

    r5239 r5446  
    245245      tsn(:,:,:,jp_sal) = sf_dyn(jf_sal)%fnow(:,:,:)  * tmask(:,:,:)    ! salinity 
    246246      ! 
    247       CALL eos    ( tsn, rhd, rhop, gdept_0(:,:,:) )                                       ! In any case, we need rhop 
     247      ! 
     248      CALL eos    ( tsn, rhd, rhop, gdept_0(:,:,:) ) ! In any case, we need rhop 
     249      CALL eos_rab( tsn, rab_n )       ! now    local thermal/haline expension ratio at T-points 
     250      CALL bn2    ( tsn, rab_n, rn2 ) ! before Brunt-Vaisala frequency need for zdfmxl 
     251 
     252      rn2b(:,:,:) = rn2(:,:,:)         ! need for zdfmxl 
    248253      CALL zdf_mxl( kt )                                                   ! In any case, we need mxl  
    249254      ! 
     
    259264      fr_i(:,:)        = sf_dyn(jf_ice)%fnow(:,:,1) * tmask(:,:,1)    ! Sea-ice fraction 
    260265      qsr (:,:)        = sf_dyn(jf_qsr)%fnow(:,:,1) * tmask(:,:,1)    ! solar radiation 
    261       IF ( ln_dynrnf ) & 
     266      IF( ln_dynrnf ) & 
    262267      rnf (:,:)        = sf_dyn(jf_rnf)%fnow(:,:,1) * tmask(:,:,1)    ! river runoffs  
    263268 
     
    383388 
    384389      ! 
    385       IF ( ln_dynrnf ) THEN 
     390      IF( ln_dynrnf ) THEN 
    386391                jf_rnf = jfld + 1  ;  jfld  = jf_rnf 
    387392         slf_d(jf_rnf) = sn_rnf 
     
    535540      !!--------------------------------------------------------------------- 
    536541#if defined key_ldfslp && ! defined key_c1d 
     542      CALL eos    ( pts, rhd, rhop, gdept_0(:,:,:) ) 
    537543      CALL eos_rab( pts, rab_n )       ! now local thermal/haline expension ratio at T-points 
    538544      CALL bn2    ( pts, rab_n, rn2  ) ! now    Brunt-Vaisala 
    539       IF( ln_zps )   & ! Partial steps: before Horizontal DErivative 
    540         &    CALL zps_hde( kt, jpts, pts, gtsu, gtsv,  &        ! Partial steps: before horizontal gradient 
    541         &                                        rhd, gru , grv , aru , arv , gzu , gzv , ge3ru , ge3rv ,   &             ! 
    542         &                                 gtui, gtvi, grui, grvi, arui, arvi, gzui, gzvi, ge3rui, ge3rvi    )  ! of t, s, rd at the last ocean level 
    543         ! only gtsu, gtsv, rhd, gru , grv are used  
    544  
    545  
    546          !                                                            ! of t, s, rd at the bottom ocean level 
     545 
     546      ! Partial steps: before Horizontal DErivative 
     547      IF( ln_zps  .AND. .NOT. ln_isfcav)                            & 
     548         &            CALL zps_hde    ( kt, jpts, pts, gtsu, gtsv,  &  ! Partial steps: before horizontal gradient 
     549         &                                        rhd, gru , grv    )  ! of t, s, rd at the last ocean level 
     550      IF( ln_zps .AND.        ln_isfcav)                            & 
     551         &            CALL zps_hde_isf( kt, jpts, pts, gtsu, gtsv,  &    ! Partial steps for top cell (ISF) 
     552         &                                        rhd, gru , grv , aru , arv , gzu , gzv , ge3ru , ge3rv ,   & 
     553         &                                 gtui, gtvi, grui, grvi, arui, arvi, gzui, gzvi, ge3rui, ge3rvi    ) ! of t, s, rd at the first ocean level 
     554 
     555      rn2b(:,:,:) = rn2(:,:,:)         ! need for zdfmxl 
    547556      CALL zdf_mxl( kt )            ! mixed layer depth 
    548557      CALL ldf_slp( kt, rhd, rn2 )  ! slopes 
Note: See TracChangeset for help on using the changeset viewer.