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 14200 for NEMO/branches/2020/dev_r14116_HPC-04_mcastril_Mixed_Precision_implementation_final/src/OCE/ISF/isfcpl.F90 – NEMO

Ignore:
Timestamp:
2020-12-17T15:36:44+01:00 (4 years ago)
Author:
mcastril
Message:

Merging r14117 through r14199 into dev_r14116_HPC-04_mcastril_Mixed_Precision_implementation_final

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2020/dev_r14116_HPC-04_mcastril_Mixed_Precision_implementation_final/src/OCE/ISF/isfcpl.F90

    r14072 r14200  
    1515#if defined key_qco 
    1616   USE domqco  , ONLY : dom_qco_zgr      ! vertical scale factor interpolation 
     17#elif defined key_linssh 
     18   !                                     ! fix in time coordinate 
    1719#else 
    1820   USE domvvl  , ONLY : dom_vvl_zgr      ! vertical scale factor interpolation 
     
    117119      vv   (:,:,:,Kbb)   = vv   (:,:,:,Kmm) 
    118120      ssh (:,:,Kbb)     = ssh (:,:,Kmm) 
    119 #if ! defined key_qco 
     121#if ! defined key_qco   &&   ! defined key_linssh 
    120122      e3t(:,:,:,Kbb)   = e3t(:,:,:,Kmm) 
    121123#endif 
     
    217219      IF(lwp) write(numout,*) 'isfcpl_ssh : recompute scale factor from ssh (new wet cell,Kmm)' 
    218220      IF(lwp) write(numout,*) '~~~~~~~~~~~' 
    219 #if ! defined key_qco 
     221#if defined key_qco 
     222      CALL dom_qco_zgr(Kbb, Kmm) 
     223#elif defined key_linssh 
     224      ! linear ssh : fix in time coord. 
     225#else 
    220226      DO jk = 1, jpk 
    221227         e3t(:,:,jk,Kmm) = e3t_0(:,:,jk) * ( 1._wp + (ht_0(:,:) + ssh(:,:,Kmm)) * r1_ht_0(:,:) ) 
     
    223229      e3t(:,:,:,Kbb) = e3t(:,:,:,Kmm) 
    224230      CALL dom_vvl_zgr(Kbb, Kmm, Kaa) 
    225 #else 
    226       CALL dom_qco_zgr(Kbb, Kmm) 
    227231#endif 
    228232      ! 
Note: See TracChangeset for help on using the changeset viewer.