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 15540 for NEMO/branches/2021/dev_r14116_HPC-10_mcastril_Mixed_Precision_implementation/src/OCE/DYN/dynzdf.F90 – NEMO

Ignore:
Timestamp:
2021-11-26T12:27:56+01:00 (3 years ago)
Author:
sparonuz
Message:

Mixed precision version, tested up to 30 years on ORCA2.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2021/dev_r14116_HPC-10_mcastril_Mixed_Precision_implementation/src/OCE/DYN/dynzdf.F90

    r14986 r15540  
    3838   PUBLIC   dyn_zdf   !  routine called by step.F90 
    3939 
    40    REAL(wp) ::  r_vvl     ! non-linear free surface indicator: =0 if ln_linssh=T, =1 otherwise  
     40   REAL(dp) ::  r_vvl     ! non-linear free surface indicator: =0 if ln_linssh=T, =1 otherwise  
    4141 
    4242   !! * Substitutions 
     
    7979      REAL(wp) ::   zzwi, ze3ua, zdt   ! local scalars 
    8080      REAL(wp) ::   zzws, ze3va        !   -      - 
    81       REAL(wp) ::   z1_e3ua, z1_e3va   !   -      - 
    82       REAL(wp) ::   zWu , zWv          !   -      - 
     81      REAL(dp) ::   z1_e3ua, z1_e3va   !   -      - 
     82      REAL(dp) ::   zWu , zWv          !   -      - 
    8383      REAL(wp) ::   zWui, zWvi         !   -      - 
    8484      REAL(wp) ::   zWus, zWvs         !   -      - 
     
    445445      ENDIF 
    446446      !                                          ! print mean trends (used for debugging) 
    447       IF(sn_cfctl%l_prtctl)   CALL prt_ctl( tab3d_1=CASTWP(puu(:,:,:,Kaa)), clinfo1=' zdf  - Ua: ', mask1=umask,               & 
    448          &                                  tab3d_2=CASTWP(pvv(:,:,:,Kaa)), clinfo2=       ' Va: ', mask2=vmask, clinfo3='dyn' ) 
     447      IF(sn_cfctl%l_prtctl)   CALL prt_ctl( tab3d_1=puu(:,:,:,Kaa), clinfo1=' zdf  - Ua: ', mask1=umask,               & 
     448         &                                  tab3d_2=pvv(:,:,:,Kaa), clinfo2=       ' Va: ', mask2=vmask, clinfo3='dyn' ) 
    449449         ! 
    450450      IF( ln_timing )   CALL timing_stop('dyn_zdf') 
Note: See TracChangeset for help on using the changeset viewer.