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 15062 – NEMO

Changeset 15062


Ignore:
Timestamp:
2021-06-28T13:19:48+02:00 (3 years ago)
Author:
jchanut
Message:

Suppress time varying scale factors and depths declarations with key_qco and key_linssh. Remove spaces that preclude from correct replacement of some scale factor arrays during preprocessing stage (at least with Apple clang version 11.0.3, this is problem).

Location:
NEMO/trunk/src
Files:
18 edited

Legend:

Unmodified
Added
Removed
  • NEMO/trunk/src/ICE/icedyn_rhg_eap.F90

    r15014 r15062  
    6666   !! * Substitutions 
    6767#  include "do_loop_substitute.h90" 
     68#  include "domzgr_substitute.h90" 
    6869   !!---------------------------------------------------------------------- 
    6970   !! NEMO/ICE 4.0 , NEMO Consortium (2018) 
  • NEMO/trunk/src/ICE/icedyn_rhg_evp.F90

    r15049 r15062  
    5555   !! * Substitutions 
    5656#  include "do_loop_substitute.h90" 
     57#  include "domzgr_substitute.h90" 
    5758   !!---------------------------------------------------------------------- 
    5859   !! NEMO/ICE 4.0 , NEMO Consortium (2018) 
  • NEMO/trunk/src/OCE/C1D/dtauvd.F90

    r15023 r15062  
    3333   !! * Substitutions 
    3434#  include "do_loop_substitute.h90" 
     35#  include "domzgr_substitute.h90" 
    3536   !!---------------------------------------------------------------------- 
    3637   !! NEMO/OCE 4.0 , NEMO Consortium (2018) 
     
    169170                  zvp(jk) =  pvd(ji,jj,jpkm1) 
    170171               ELSE                                      ! inbetween : vertical interpolation between jkk & jkk+1 
    171                   DO jkk = 1, jpkm1                      ! when  gdept(jkk) < zl < gdept(jkk+1) 
     172                  DO jkk = 1, jpkm1                      ! when  dept(jkk) < zl < dept(jkk+1) 
    172173                     IF( (zl-gdept_1d(jkk)) * (zl-gdept_1d(jkk+1)) <= 0._wp ) THEN 
    173174                        zi = ( zl - gdept_1d(jkk) ) / (gdept_1d(jkk+1)-gdept_1d(jkk)) 
  • NEMO/trunk/src/OCE/C1D/dyndmp.F90

    r15023 r15062  
    4444   !! * Substitutions 
    4545#  include "do_loop_substitute.h90" 
     46#  include "domzgr_substitute.h90" 
    4647   !!---------------------------------------------------------------------- 
    4748   !! NEMO/OCE 4.0 , NEMO Consortium (2018) 
  • NEMO/trunk/src/OCE/DIA/diahsb.F90

    r15048 r15062  
    186186      ! 
    187187      DO jk = 1, jpkm1           ! volume 
    188          ztmpk(:,:,jk,1) =   surf    (:,:) * e3t    (:,:,jk,Kmm)*tmask    (:,:,jk)   & 
     188         ztmpk(:,:,jk,1) =   surf    (:,:) * e3t(:,:,jk,Kmm)*tmask(:,:,jk)   & 
    189189            &              - surf_ini(:,:) * e3t_ini(:,:,jk    )*tmask_ini(:,:,jk) 
    190190      END DO 
  • NEMO/trunk/src/OCE/DOM/dom_oce.F90

    r15023 r15062  
    140140 
    141141   !                                                        !  time-dependent scale factors     (domvvl) 
     142#if defined key_qco || defined key_linssh 
     143#else 
    142144   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:,:) ::   e3t, e3u, e3v, e3w, e3uw, e3vw  !: vert. scale factor [m] 
    143145   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:)   ::   e3f                             !: F-point vert. scale factor [m] 
    144  
     146#endif 
    145147   !                                                        !  time-dependent ratio ssh / h_0   (domqco) 
    146148   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:)   ::   r3t, r3u, r3v                   !: time-dependent    ratio at t-, u- and v-point [-] 
     
    154156    
    155157   !                                                        !  time-dependent depths of cells   (domvvl) 
     158#if defined key_qco || defined key_linssh 
     159#else 
    156160   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:,:) ::  gdept, gdepw 
    157161   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:)   ::  gde3w 
    158  
     162#endif 
    159163   !                                                        !  reference heights of ocean water column and its inverse 
    160164   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:)     ::   ht_0, r1_ht_0   !: t-depth        [m] and [1/m] 
     
    164168    
    165169   !                                                        ! time-dependent heights of ocean water column   (domvvl) 
     170#if defined key_qco || defined key_linssh  
     171#else 
    166172   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:)     ::   ht          !: t-points           [m] 
    167173   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:)   ::   hu, r1_hu   !: u-depth            [m] and [1/m] 
    168174   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:)   ::   hv, r1_hv   !: v-depth            [m] and [1/m] 
    169  
     175#endif 
    170176   INTEGER, PUBLIC ::   nla10              !: deepest    W level Above  ~10m (nlb10 - 1) 
    171177   INTEGER, PUBLIC ::   nlb10              !: shallowest W level Bellow ~10m (nla10 + 1) 
  • NEMO/trunk/src/OCE/FLO/flodom.F90

    r14229 r15062  
    3333   REAL(wp), ALLOCATABLE, DIMENSION(:) ::   zgifl, zgjfl,  zgkfl      ! distances in indexes 
    3434 
     35   !! * Substitutions 
     36#  include "domzgr_substitute.h90" 
    3537   !!---------------------------------------------------------------------- 
    3638   !! NEMO/OCE 4.0 , NEMO Consortium (2018) 
  • NEMO/trunk/src/OCE/FLO/flowri.F90

    r14239 r15062  
    3333   REAL(wp), ALLOCATABLE, DIMENSION(:) ::   ztem , zsal, zrho   ! 2D workspace 
    3434 
     35   !! * Substitutions 
     36#  include "domzgr_substitute.h90" 
    3537   !!---------------------------------------------------------------------- 
    3638   !! NEMO/OCE 4.0 , NEMO Consortium (2018) 
  • NEMO/trunk/src/OCE/ISF/isfcpl.F90

    r15053 r15062  
    546546 
    547547               ! volume diff 
    548                zdvol =   e3t  (ji,jj,jk,Kmm) *  tmask  (ji,jj,jk)   & 
     548               zdvol =   e3t(ji,jj,jk,Kmm) *  tmask  (ji,jj,jk)   & 
    549549                  &   - ze3t_b(ji,jj,jk    ) * ztmask_b(ji,jj,jk) 
    550550 
  • NEMO/trunk/src/OCE/ISF/isfparmlt.F90

    r15004 r15062  
    2828 
    2929   PUBLIC  isfpar_mlt  
    30     
     30 
     31 
     32   !! * Substitutions 
     33#  include "domzgr_substitute.h90"    
    3134   !!---------------------------------------------------------------------- 
    3235   !! NEMO/OCE 4.0 , NEMO Consortium (2018) 
  • NEMO/trunk/src/OCE/LDF/ldfslp.F90

    r14834 r15062  
    208208         zfj = MAX( omlmask(ji,jj,jk), omlmask(ji,jj+1,jk) ) 
    209209         ! thickness of water column between surface and level k at u/v point 
    210          zdepu = 0.5_wp * ( ( gdept (ji,jj,jk,Kmm) + gdept (ji+1,jj,jk,Kmm) )                            & 
     210         zdepu = 0.5_wp * ( ( gdept(ji,jj,jk,Kmm) + gdept(ji+1,jj,jk,Kmm) )                            & 
    211211            &              - 2 * MAX( risfdep(ji,jj), risfdep(ji+1,jj) )        & 
    212212            &              - e3u(ji,jj,miku(ji,jj),Kmm)   ) 
    213          zdepv = 0.5_wp * ( ( gdept (ji,jj,jk,Kmm) + gdept (ji,jj+1,jk,Kmm) )                            & 
     213         zdepv = 0.5_wp * ( ( gdept(ji,jj,jk,Kmm) + gdept(ji,jj+1,jk,Kmm) )                            & 
    214214            &              - 2 * MAX( risfdep(ji,jj), risfdep(ji,jj+1) )        & 
    215215            &              - e3v(ji,jj,mikv(ji,jj),Kmm)   ) 
  • NEMO/trunk/src/OCE/OBS/diaobs.F90

    r14982 r15062  
    604604      !! ** Action  : 
    605605      !!---------------------------------------------------------------------- 
    606       USE dom_oce, ONLY : gdept, gdept_1d     ! Ocean space domain variables (Kmm time-level only) 
    607606      USE phycst , ONLY : rday                ! Physical constants 
    608607      USE oce    , ONLY : ts, uu, vv, ssh     ! Ocean dynamics and tracers variables (Kmm time-level only) 
  • NEMO/trunk/src/OCE/OBS/obs_prep.F90

    r14982 r15062  
    10721072      !!---------------------------------------------------------------------- 
    10731073      !! * Modules used 
    1074       USE dom_oce, ONLY : &       ! Geographical information  
    1075          & gdepw_1d,      & 
    1076          & gdepw_0,       &                        
    1077          & gdepw, r3t,    & 
    1078          & gdept,         & 
    1079          & ln_zco,        & 
    1080          & ln_zps              
     1074      USE dom_oce   ! Geographical information  
    10811075 
    10821076      !! * Arguments 
  • NEMO/trunk/src/OCE/SBC/sbcrnf.F90

    r15058 r15062  
    231231      ELSE                       !==   runoff put only at the surface   ==! 
    232232         DO_2D_OVR( nn_hls, nn_hls, nn_hls, nn_hls ) 
    233             h_rnf (ji,jj)   = e3t (ji,jj,1,Kmm)        ! update h_rnf to be depth of top box 
     233            h_rnf (ji,jj)   = e3t(ji,jj,1,Kmm)        ! update h_rnf to be depth of top box 
    234234         END_2D 
    235235         DO_2D_OVR( nn_hls-1, nn_hls, nn_hls-1, nn_hls ) 
  • NEMO/trunk/src/OCE/SBC/sbcssm.F90

    r14072 r15062  
    246246         sss_m(:,:) = ts  (:,:,1,jp_sal,Kmm) 
    247247         ssh_m(:,:) = ssh(:,:,Kmm) 
    248          e3t_m(:,:) = e3t (:,:,1,Kmm) 
     248         e3t_m(:,:) = e3t(:,:,1,Kmm) 
    249249         frq_m(:,:) = 1._wp 
    250250         ! 
  • NEMO/trunk/src/OCE/TRA/traldf_triad.F90

    r14834 r15062  
    210210                     &   * (  akz     (ji,jj,jk)   & 
    211211                     &      + ah_wslp2(ji,jj,jk)   & 
    212                      &        / ( e3w (ji,jj,jk,Kmm) * e3w(ji,jj,jk,Kmm) )  ) 
     212                     &        / ( e3w(ji,jj,jk,Kmm) * e3w(ji,jj,jk,Kmm) )  ) 
    213213               END_3D 
    214214            ELSEIF( ln_traldf_lap ) THEN              ! laplacian operator 
  • NEMO/trunk/src/OCE/ZDF/zdfswm.F90

    r14834 r15062  
    2929   !! * Substitutions 
    3030#  include "do_loop_substitute.h90" 
     31#  include "domzgr_substitute.h90" 
    3132   !!---------------------------------------------------------------------- 
    3233   !! NEMO/OCE 4.0 , NEMO Consortium (2018) 
  • NEMO/trunk/src/TOP/PISCES/SED/oce_sed.F90

    r14086 r15062  
    1515   USE dom_oce , ONLY :   gphit     =>   gphit          !: latitude  of t-point (degre) 
    1616  
    17 #if ! defined key_qco 
     17#if defined key_qco || defined key_linssh 
     18#else 
    1819   USE dom_oce , ONLY :   e3t       =>   e3t            !: latitude  of t-point (degre) 
    1920#endif 
Note: See TracChangeset for help on using the changeset viewer.