# Changeset 12740

Ignore:
Timestamp:
2020-04-12T11:03:06+02:00 (8 months ago)
Message:

trunk: update/debug of tests and C1D, see #2442

Location:
NEMO/trunk
Files:
1 deleted
31 edited

Unmodified
Removed
• ## NEMO/trunk/cfgs/C1D_PAPA/MY_SRC/usrdef_zgr.F90

 r12377 PUBLIC   usr_def_zgr        ! called by domzgr.F90 !! * Substitutions #  include "do_loop_substitute.h90" !!---------------------------------------------------------------------- !! NEMO/OCE 4.0 , NEMO Consortium (2018) pe3vw(:,:,jk) = pe3w_1d (jk) END DO DO jj = 1, jpj                      ! bottom scale factors and depth at T- and W-points DO ji = 1, jpi ik = k_bot(ji,jj) pdepw(ji,jj,ik+1) = MIN( zht(ji,jj) , pdepw_1d(ik+1) ) pe3t (ji,jj,ik  ) = pdepw(ji,jj,ik+1) - pdepw(ji,jj,ik) pe3t (ji,jj,ik+1) = pe3t (ji,jj,ik  ) ! pdept(ji,jj,ik  ) = pdepw(ji,jj,ik  ) + pe3t (ji,jj,ik  ) * 0.5_wp pdept(ji,jj,ik+1) = pdepw(ji,jj,ik+1) + pe3t (ji,jj,ik+1) * 0.5_wp pe3w (ji,jj,ik+1) = pdept(ji,jj,ik+1) - pdept(ji,jj,ik)              ! = pe3t (ji,jj,ik  ) END DO END DO ! bottom scale factors and depth at T- and W-points DO_2D_11_11 ik = k_bot(ji,jj) pdepw(ji,jj,ik+1) = MIN( zht(ji,jj) , pdepw_1d(ik+1) ) pe3t (ji,jj,ik  ) = pdepw(ji,jj,ik+1) - pdepw(ji,jj,ik) pe3t (ji,jj,ik+1) = pe3t (ji,jj,ik  ) ! pdept(ji,jj,ik  ) = pdepw(ji,jj,ik  ) + pe3t (ji,jj,ik  ) * 0.5_wp pdept(ji,jj,ik+1) = pdepw(ji,jj,ik+1) + pe3t (ji,jj,ik+1) * 0.5_wp pe3w (ji,jj,ik+1) = pdept(ji,jj,ik+1) - pdept(ji,jj,ik)              ! = pe3t (ji,jj,ik  ) END_2D !                                   ! bottom scale factors and depth at  U-, V-, UW and VW-points !                                   ! usually Computed as the minimum of neighbooring scale factors
• ## NEMO/trunk/src/OCE/C1D/step_c1d.F90

 r12377 PRIVATE PUBLIC stp_c1d      ! called by opa.F90 PUBLIC stp_c1d      ! called by nemogcm.F90 !!---------------------------------------------------------------------- !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< CALL dia_wri( kstp, Nnn )  ! ocean model: outputs IF( lk_diahth  )   CALL dia_hth( kstp, Nnn )  ! Thermocline depth (20°C) CALL dia_hth( kstp, Nnn )  ! Thermocline depth (20°C) CALL eos( ts(:,:,:,:,Nnn), rhd, rhop, gdept_0(:,:,:) )  ! now potential density for zdfmxl IF( ln_zdfnpc )   CALL tra_npc( kstp,      Nnn, Nrhs, ts, Naa   )         ! applied non penetrative convective adjustment on (t,s) CALL tra_atf( kstp, Nbb, Nnn, Nrhs,     Naa, ts   )     ! time filtering of "now" tracer fields CALL tra_atf( kstp, Nbb, Nnn, Naa, ts )                 ! time filtering of "now" tracer arrays !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> ! Control and restarts !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< CALL stp_ctl( kstp, Nnn, indic ) CALL stp_ctl( kstp, Nbb, Nnn, indic ) IF( kstp == nit000 )   CALL iom_close( numror )          ! close input  ocean restart file IF( lrst_oce       )   CALL rst_write( kstp, Nbb, Nnn )  ! write output ocean restart file
• ## NEMO/trunk/src/OCE/DOM/domvvl.F90

 r12489 e3t(:,:,:,Kbb) = e3t(:,:,:,Kmm) DO ji = 1, jpi DO jj = 1, jpj IF ( ht_0(ji,jj) .LE. 0.0 .AND. NINT( ssmask(ji,jj) ) .EQ. 1) THEN CALL ctl_stop( 'dom_vvl_rst: ht_0 must be positive at potentially wet points' ) ENDIF END DO END DO DO_2D_11_11 IF ( ht_0(ji,jj) .LE. 0.0 .AND. NINT( ssmask(ji,jj) ) .EQ. 1) THEN CALL ctl_stop( 'dom_vvl_rst: ht_0 must be positive at potentially wet points' ) ENDIF END_2D ! ELSE
• ## NEMO/trunk/src/OCE/USR/usrdef_zgr.F90

 r12377 CALL lbc_lnk( 'usrdef_zgr', z2d, 'T', 1. )           ! set surrounding land to zero (here jperio=0 ==>> closed) ! k_bot(:,:) = NINT( z2d(:,:) )           ! =jpkm1 over the ocean point, =0 elsewhere k_bot(:,:) = NINT( z2d(:,:) )          ! =jpkm1 over the ocean point, =0 elsewhere ! k_top(:,:) = MIN( 1 , k_bot(:,:) )     ! = 1    over the ocean point, =0 elsewhere
• ## NEMO/trunk/tests/BENCH/MY_SRC/usrdef_hgr.F90

 r9762 PUBLIC   usr_def_hgr   ! called by domhgr.F90 !! * Substitutions #  include "do_loop_substitute.h90" !!---------------------------------------------------------------------- !! NEMO/OPA 4.0, NEMO Consortium (2016) ! ! Position coordinates (in grid points) !                          ========== DO jj = 1, jpj DO ji = 1, jpi zti = REAL( ji - 1 + nimpp - 1, wp )          ;  ztj = REAL( jj - 1 + njmpp - 1, wp ) zui = REAL( ji - 1 + nimpp - 1, wp ) + 0.5_wp ;  zvj = REAL( jj - 1 + njmpp - 1, wp ) + 0.5_wp !                          ========== DO_2D_11_11 zti = REAL( ji - 1 + nimpp - 1, wp )          ;  ztj = REAL( jj - 1 + njmpp - 1, wp ) zui = REAL( ji - 1 + nimpp - 1, wp ) + 0.5_wp ;  zvj = REAL( jj - 1 + njmpp - 1, wp ) + 0.5_wp plamt(ji,jj) = zti plamu(ji,jj) = zui plamv(ji,jj) = zti plamf(ji,jj) = zui pphit(ji,jj) = ztj pphiv(ji,jj) = zvj pphiu(ji,jj) = ztj pphif(ji,jj) = zvj plamt(ji,jj) = zti plamu(ji,jj) = zui plamv(ji,jj) = zti plamf(ji,jj) = zui pphit(ji,jj) = ztj pphiv(ji,jj) = zvj pphiu(ji,jj) = ztj pphif(ji,jj) = zvj END DO END DO END_2D ! ! Horizontal scale factors (in meters)
• ## NEMO/trunk/tests/BENCH/MY_SRC/usrdef_istate.F90

 r11536 PUBLIC   usr_def_istate   ! called by istate.F90 !! * Substitutions #  include "do_loop_substitute.h90" !!---------------------------------------------------------------------- !! NEMO/OPA 4.0 , NEMO Consortium (2016) ! ! define unique value on each point. z2d ranging from 0.05 to -0.05 DO jj = 1, jpj DO ji = 1, jpi z2d(ji,jj) = 0.1 * ( 0.5 - REAL( mig(ji) + mjg(jj) * jpiglo, wp ) / REAL( jpiglo * jpjglo, wp ) ) ENDDO ENDDO DO_2D_11_11 z2d(ji,jj) = 0.1 * ( 0.5 - REAL( mig(ji) + mjg(jj) * jpiglo, wp ) / REAL( jpiglo * jpjglo, wp ) ) END_2D ! ! sea level:
• ## NEMO/trunk/tests/BENCH/MY_SRC/usrdef_sbc.F90

 r12377 PUBLIC   usrdef_sbc_ice_flx  ! routine called by sbcice_lim.F90 for ice thermo !! * Substitutions #  include "do_loop_substitute.h90" !!---------------------------------------------------------------------- !! NEMO/OPA 4.0 , NEMO Consortium (2016) ! ! define unique value on each point. z2d ranging from 0.05 to -0.05 DO jj = 1, jpj DO ji = 1, jpi z2d(ji,jj) = 0.1 * ( 0.5 - REAL( nimpp + ji - 1 + ( njmpp + jj - 2 ) * jpiglo, wp ) / REAL( jpiglo * jpjglo, wp ) ) ENDDO ENDDO DO_2D_11_11 z2d(ji,jj) = 0.1 * ( 0.5 - REAL( nimpp + ji - 1 + ( njmpp + jj - 2 ) * jpiglo, wp ) / REAL( jpiglo * jpjglo, wp ) ) END_2D utau_ice(:,:) = 0.1_wp +  z2d(:,:) vtau_ice(:,:) = 0.1_wp +  z2d(:,:)

• ## NEMO/trunk/tests/CANAL/MY_SRC/usrdef_hgr.F90

 r10074 PUBLIC   usr_def_hgr   ! called by domhgr.F90 !! * Substitutions #  include "do_loop_substitute.h90" !!---------------------------------------------------------------------- !! NEMO/OCE 4.0 , NEMO Consortium (2018) #endif DO jj = 1, jpj DO ji = 1, jpi zti = FLOAT( ji - 1 + nimpp - 1 )          ;  ztj = FLOAT( jj - 1 + njmpp - 1 ) zui = FLOAT( ji - 1 + nimpp - 1 ) + 0.5_wp ;  zvj = FLOAT( jj - 1 + njmpp - 1 ) + 0.5_wp plamt(ji,jj) = zlam0 + rn_dx * zti plamu(ji,jj) = zlam0 + rn_dx * zui plamv(ji,jj) = plamt(ji,jj) plamf(ji,jj) = plamu(ji,jj) pphit(ji,jj) = zphi0 + rn_dy * ztj pphiv(ji,jj) = zphi0 + rn_dy * zvj pphiu(ji,jj) = pphit(ji,jj) pphif(ji,jj) = pphiv(ji,jj) END DO END DO DO_2D_11_11 zti = FLOAT( ji - 1 + nimpp - 1 )          ;  ztj = FLOAT( jj - 1 + njmpp - 1 ) zui = FLOAT( ji - 1 + nimpp - 1 ) + 0.5_wp ;  zvj = FLOAT( jj - 1 + njmpp - 1 ) + 0.5_wp plamt(ji,jj) = zlam0 + rn_dx * zti plamu(ji,jj) = zlam0 + rn_dx * zui plamv(ji,jj) = plamt(ji,jj) plamf(ji,jj) = plamu(ji,jj) pphit(ji,jj) = zphi0 + rn_dy * ztj pphiv(ji,jj) = zphi0 + rn_dy * zvj pphiu(ji,jj) = pphit(ji,jj) pphif(ji,jj) = pphiv(ji,jj) END_2D ! ! Horizontal scale factors (in meters)

• ## NEMO/trunk/tests/CANAL/MY_SRC/usrdef_sbc.F90

 r12377 CONTAINS SUBROUTINE usrdef_sbc_oce( kt, Kmm, Kbb ) SUBROUTINE usrdef_sbc_oce( kt, Kbb ) !!--------------------------------------------------------------------- !!                    ***  ROUTINE usr_def_sbc  *** !!---------------------------------------------------------------------- INTEGER, INTENT(in) ::   kt        ! ocean time step INTEGER, INTENT(in) ::   Kbb, Kmm  ! ocean time index INTEGER, INTENT(in) ::   Kbb       ! ocean time index INTEGER  ::   ji, jj               ! dummy loop indices REAL(wp) :: zrhoair = 1.22     ! approximate air density [Kg/m3] WHERE( ABS(gphit) <= rn_windszy/2. ) zwndrel(:,:) = rn_u10 - rn_uofac * uu(:,:,1,Kmm) zwndrel(:,:) = rn_u10 - rn_uofac * uu(:,:,1,Kbb) ELSEWHERE zwndrel(:,:) =        - rn_uofac * uu(:,:,1,Kmm) zwndrel(:,:) =        - rn_uofac * uu(:,:,1,Kbb) END WHERE utau(:,:) = zrhocd * zwndrel(:,:) * zwndrel(:,:) zwndrel(:,:) = - rn_uofac * vv(:,:,1,Kmm) zwndrel(:,:) = - rn_uofac * vv(:,:,1,Kbb) vtau(:,:) = zrhocd * zwndrel(:,:) * zwndrel(:,:)
• ## NEMO/trunk/tests/CANAL/MY_SRC/usrdef_zgr.F90

 r12377 CALL lbc_lnk( 'usrdef_zgr', z2d, 'T', 1. )           ! set surrounding land to zero (here jperio=0 ==>> closed) ! k_bot(:,:) = INT( z2d(:,:) )           ! =jpkm1 over the ocean point, =0 elsewhere k_bot(:,:) = NINT( z2d(:,:) )          ! =jpkm1 over the ocean point, =0 elsewhere ! k_top(:,:) = MIN( 1 , k_bot(:,:) )     ! = 1    over the ocean point, =0 elsewhere

 r10513 PUBLIC   usr_def_hgr   ! called by domhgr.F90 !! * Substitutions #  include "do_loop_substitute.h90" !!---------------------------------------------------------------------- !! NEMO/OCE 4.0 , NEMO Consortium (2018) zphi0 = -(jpjglo-1)/2 * 1.e-3 * rn_dy DO jj = 1, jpj DO ji = 1, jpi zti = FLOAT( ji - 1 + nimpp - 1 )          ;  ztj = FLOAT( jj - 1 + njmpp - 1 ) zui = FLOAT( ji - 1 + nimpp - 1 ) + 0.5_wp ;  zvj = FLOAT( jj - 1 + njmpp - 1 ) + 0.5_wp plamt(ji,jj) = zlam0 + rn_dx * 1.e-3 * zti plamu(ji,jj) = zlam0 + rn_dx * 1.e-3 * zui plamv(ji,jj) = plamt(ji,jj) plamf(ji,jj) = plamu(ji,jj) pphit(ji,jj) = zphi0 + rn_dy * 1.e-3 * ztj pphiv(ji,jj) = zphi0 + rn_dy * 1.e-3 * zvj pphiu(ji,jj) = pphit(ji,jj) pphif(ji,jj) = pphiv(ji,jj) END DO END DO DO_2D_11_11 zti = FLOAT( ji - 1 + nimpp - 1 )          ;  ztj = FLOAT( jj - 1 + njmpp - 1 ) zui = FLOAT( ji - 1 + nimpp - 1 ) + 0.5_wp ;  zvj = FLOAT( jj - 1 + njmpp - 1 ) + 0.5_wp plamt(ji,jj) = zlam0 + rn_dx * 1.e-3 * zti plamu(ji,jj) = zlam0 + rn_dx * 1.e-3 * zui plamv(ji,jj) = plamt(ji,jj) plamf(ji,jj) = plamu(ji,jj) pphit(ji,jj) = zphi0 + rn_dy * 1.e-3 * ztj pphiv(ji,jj) = zphi0 + rn_dy * 1.e-3 * zvj pphiu(ji,jj) = pphit(ji,jj) pphif(ji,jj) = pphiv(ji,jj) END_2D ! constant scale factors

 r10515 PUBLIC   usr_def_hgr   ! called by domhgr.F90 !! * Substitutions #  include "do_loop_substitute.h90" !!---------------------------------------------------------------------- !! NEMO/OCE 4.0 , NEMO Consortium (2018) #endif DO jj = 1, jpj DO ji = 1, jpi zti = FLOAT( ji - 1 + nimpp - 1 )          ;  ztj = FLOAT( jj - 1 + njmpp - 1 ) zui = FLOAT( ji - 1 + nimpp - 1 ) + 0.5_wp ;  zvj = FLOAT( jj - 1 + njmpp - 1 ) + 0.5_wp plamt(ji,jj) = zlam0 + rn_dx * 1.e-3 * zti plamu(ji,jj) = zlam0 + rn_dx * 1.e-3 * zui plamv(ji,jj) = plamt(ji,jj) plamf(ji,jj) = plamu(ji,jj) pphit(ji,jj) = zphi0 + rn_dy * 1.e-3 * ztj pphiv(ji,jj) = zphi0 + rn_dy * 1.e-3 * zvj pphiu(ji,jj) = pphit(ji,jj) pphif(ji,jj) = pphiv(ji,jj) END DO END DO DO_2D_11_11 zti = FLOAT( ji - 1 + nimpp - 1 )          ;  ztj = FLOAT( jj - 1 + njmpp - 1 ) zui = FLOAT( ji - 1 + nimpp - 1 ) + 0.5_wp ;  zvj = FLOAT( jj - 1 + njmpp - 1 ) + 0.5_wp plamt(ji,jj) = zlam0 + rn_dx * 1.e-3 * zti plamu(ji,jj) = zlam0 + rn_dx * 1.e-3 * zui plamv(ji,jj) = plamt(ji,jj) plamf(ji,jj) = plamu(ji,jj) pphit(ji,jj) = zphi0 + rn_dy * 1.e-3 * ztj pphiv(ji,jj) = zphi0 + rn_dy * 1.e-3 * zvj pphiu(ji,jj) = pphit(ji,jj) pphif(ji,jj) = pphiv(ji,jj) END_2D ! Horizontal scale factors (in meters)

 r12377 !!   usr_def_hgr   : initialize the horizontal mesh !!---------------------------------------------------------------------- USE dom_oce  , ONLY: nimpp , njmpp            ! i- & j-indices of the local domain USE dom_oce  , ONLY: nimpp , njmpp, Agrif_Root            ! i- & j-indices of the local domain USE par_oce        ! ocean space and time domain USE phycst         ! physical constants
• ## NEMO/trunk/tests/ICE_AGRIF/MY_SRC/usrdef_hgr.F90

 r10516 PUBLIC   usr_def_hgr   ! called by domhgr.F90 !! * Substitutions #  include "do_loop_substitute.h90" !!---------------------------------------------------------------------- !! NEMO/OCE 4.0 , NEMO Consortium (2018) #endif DO jj = 1, jpj DO ji = 1, jpi zti = FLOAT( ji - 1 + nimpp - 1 )          ;  ztj = FLOAT( jj - 1 + njmpp - 1 ) zui = FLOAT( ji - 1 + nimpp - 1 ) + 0.5_wp ;  zvj = FLOAT( jj - 1 + njmpp - 1 ) + 0.5_wp plamt(ji,jj) = zlam0 + rn_dx * 1.e-3 * zti plamu(ji,jj) = zlam0 + rn_dx * 1.e-3 * zui plamv(ji,jj) = plamt(ji,jj) plamf(ji,jj) = plamu(ji,jj) pphit(ji,jj) = zphi0 + rn_dy * 1.e-3 * ztj pphiv(ji,jj) = zphi0 + rn_dy * 1.e-3 * zvj pphiu(ji,jj) = pphit(ji,jj) pphif(ji,jj) = pphiv(ji,jj) END DO END DO DO_2D_11_11 zti = FLOAT( ji - 1 + nimpp - 1 )          ;  ztj = FLOAT( jj - 1 + njmpp - 1 ) zui = FLOAT( ji - 1 + nimpp - 1 ) + 0.5_wp ;  zvj = FLOAT( jj - 1 + njmpp - 1 ) + 0.5_wp plamt(ji,jj) = zlam0 + rn_dx * 1.e-3 * zti plamu(ji,jj) = zlam0 + rn_dx * 1.e-3 * zui plamv(ji,jj) = plamt(ji,jj) plamf(ji,jj) = plamu(ji,jj) pphit(ji,jj) = zphi0 + rn_dy * 1.e-3 * ztj pphiv(ji,jj) = zphi0 + rn_dy * 1.e-3 * zvj pphiu(ji,jj) = pphit(ji,jj) pphif(ji,jj) = pphiv(ji,jj) END_2D ! Horizontal scale factors (in meters)

• ## NEMO/trunk/tests/ISOMIP/MY_SRC/usrdef_zgr.F90

 r12377 PUBLIC   usr_def_zgr   ! called by domzgr.F90 !! * Substitutions #  include "do_loop_substitute.h90" !!---------------------------------------------------------------------- !! NEMO/OCE 4.0 , NEMO Consortium (2018) pe3vw(:,:,jk) = pe3w_1d (jk) END DO DO jj = 1, jpj                      ! top scale factors and depth at T- and W-points DO ji = 1, jpi ik = k_top(ji,jj) IF ( ik > 2 ) THEN ! pdeptw at the interface pdepw(ji,jj,ik  ) = MAX( zhisf(ji,jj) , pdepw(ji,jj,ik) ) ! e3t in both side of the interface pe3t (ji,jj,ik  ) = pdepw(ji,jj,ik+1) - pdepw(ji,jj,ik) ! pdept in both side of the interface (from previous e3t) pdept(ji,jj,ik  ) = pdepw(ji,jj,ik  ) + pe3t (ji,jj,ik  ) * 0.5_wp pdept(ji,jj,ik-1) = pdepw(ji,jj,ik  ) - pe3t (ji,jj,ik  ) * 0.5_wp ! pe3w on both side of the interface pe3w (ji,jj,ik+1) = pdept(ji,jj,ik+1) - pdept(ji,jj,ik  ) pe3w (ji,jj,ik  ) = pdept(ji,jj,ik  ) - pdept(ji,jj,ik-1) ! e3t into the ice shelf pe3t (ji,jj,ik-1) = pdepw(ji,jj,ik  ) - pdepw(ji,jj,ik-1) pe3w (ji,jj,ik-1) = pdept(ji,jj,ik-1) - pdept(ji,jj,ik-2) END IF END DO END DO DO jj = 1, jpj                      ! bottom scale factors and depth at T- and W-points DO ji = 1, jpi ik = k_bot(ji,jj) pdepw(ji,jj,ik+1) = MIN( zht(ji,jj) , pdepw_1d(ik+1) ) ! top scale factors and depth at T- and W-points DO_2D_11_11 ik = k_top(ji,jj) IF ( ik > 2 ) THEN ! pdeptw at the interface pdepw(ji,jj,ik  ) = MAX( zhisf(ji,jj) , pdepw(ji,jj,ik) ) ! e3t in both side of the interface pe3t (ji,jj,ik  ) = pdepw(ji,jj,ik+1) - pdepw(ji,jj,ik) pe3t (ji,jj,ik+1) = pe3t (ji,jj,ik  ) ! ! pdept in both side of the interface (from previous e3t) pdept(ji,jj,ik  ) = pdepw(ji,jj,ik  ) + pe3t (ji,jj,ik  ) * 0.5_wp pdept(ji,jj,ik+1) = pdepw(ji,jj,ik+1) + pe3t (ji,jj,ik+1) * 0.5_wp pe3w (ji,jj,ik+1) = pdept(ji,jj,ik+1) - pdept(ji,jj,ik) END DO END DO pdept(ji,jj,ik-1) = pdepw(ji,jj,ik  ) - pe3t (ji,jj,ik  ) * 0.5_wp ! pe3w on both side of the interface pe3w (ji,jj,ik+1) = pdept(ji,jj,ik+1) - pdept(ji,jj,ik  ) pe3w (ji,jj,ik  ) = pdept(ji,jj,ik  ) - pdept(ji,jj,ik-1) ! e3t into the ice shelf pe3t (ji,jj,ik-1) = pdepw(ji,jj,ik  ) - pdepw(ji,jj,ik-1) pe3w (ji,jj,ik-1) = pdept(ji,jj,ik-1) - pdept(ji,jj,ik-2) END IF END_2D ! bottom scale factors and depth at T- and W-points DO_2D_11_11 ik = k_bot(ji,jj) pdepw(ji,jj,ik+1) = MIN( zht(ji,jj) , pdepw_1d(ik+1) ) pe3t (ji,jj,ik  ) = pdepw(ji,jj,ik+1) - pdepw(ji,jj,ik) pe3t (ji,jj,ik+1) = pe3t (ji,jj,ik  ) ! pdept(ji,jj,ik  ) = pdepw(ji,jj,ik  ) + pe3t (ji,jj,ik  ) * 0.5_wp pdept(ji,jj,ik+1) = pdepw(ji,jj,ik+1) + pe3t (ji,jj,ik+1) * 0.5_wp pe3w (ji,jj,ik+1) = pdept(ji,jj,ik+1) - pdept(ji,jj,ik) END_2D !                                   ! bottom scale factors and depth at  U-, V-, UW and VW-points pe3u (:,:,:) = pe3t(:,:,:) pe3uw(:,:,:) = pe3w(:,:,:) DO jk = 1, jpk                      ! Computed as the minimum of neighbooring scale factors DO jj = 1, jpjm1 DO ji = 1, jpi pe3v (ji,jj,jk) = MIN( pe3t(ji,jj,jk), pe3t(ji,jj+1,jk) ) pe3vw(ji,jj,jk) = MIN( pe3w(ji,jj,jk), pe3w(ji,jj+1,jk) ) pe3f (ji,jj,jk) = pe3v(ji,jj,jk) END DO END DO END DO DO_3D_00_00( 1, jpk ) !                                   ! Computed as the minimum of neighbooring scale factors pe3v (ji,jj,jk) = MIN( pe3t(ji,jj,jk), pe3t(ji,jj+1,jk) ) pe3vw(ji,jj,jk) = MIN( pe3w(ji,jj,jk), pe3w(ji,jj+1,jk) ) pe3f (ji,jj,jk) = pe3v(ji,jj,jk) END_3D CALL lbc_lnk( 'usrdef_zgr', pe3v , 'V', 1._wp )   ;   CALL lbc_lnk( 'usrdef_zgr', pe3vw, 'V', 1._wp ) CALL lbc_lnk( 'usrdef_zgr', pe3f , 'F', 1._wp )
• ## NEMO/trunk/tests/LOCK_EXCHANGE/MY_SRC/usrdef_hgr.F90

 r10074 PUBLIC   usr_def_hgr   ! called by domhgr.F90 !! * Substitutions #  include "do_loop_substitute.h90" !!---------------------------------------------------------------------- !! NEMO/OCE 4.0 , NEMO Consortium (2018) !                       !==  grid point position  ==!   (in kilometers) zfact = rn_dx * 1.e-3         ! conversion in km DO jj = 1, jpj DO ji = 1, jpi             ! longitude plamt(ji,jj) = zfact * (  - 0.5 + REAL( ji-1 + nimpp-1 , wp )  ) plamu(ji,jj) = zfact * (          REAL( ji-1 + nimpp-1 , wp )  ) plamv(ji,jj) = plamt(ji,jj) plamf(ji,jj) = plamu(ji,jj) !                       ! latitude pphit(ji,jj) = zfact * (  - 0.5 + REAL( jj-1 + njmpp-1 , wp )  ) pphiu(ji,jj) = pphit(ji,jj) pphiv(ji,jj) = zfact * (          REAL( jj-1 + njmpp-1 , wp )  ) pphif(ji,jj) = pphiv(ji,jj) END DO END DO DO_2D_11_11 !                       ! longitude plamt(ji,jj) = zfact * (  - 0.5 + REAL( ji-1 + nimpp-1 , wp )  ) plamu(ji,jj) = zfact * (          REAL( ji-1 + nimpp-1 , wp )  ) plamv(ji,jj) = plamt(ji,jj) plamf(ji,jj) = plamu(ji,jj) !                       ! latitude pphit(ji,jj) = zfact * (  - 0.5 + REAL( jj-1 + njmpp-1 , wp )  ) pphiu(ji,jj) = pphit(ji,jj) pphiv(ji,jj) = zfact * (          REAL( jj-1 + njmpp-1 , wp )  ) pphif(ji,jj) = pphiv(ji,jj) END_2D ! !                       !==  Horizontal scale factors  ==!   (in meters)
• ## NEMO/trunk/tests/OVERFLOW/MY_SRC/usrdef_hgr.F90

 r10074 PUBLIC   usr_def_hgr   ! called by domhgr.F90 !! * Substitutions #  include "do_loop_substitute.h90" !!---------------------------------------------------------------------- !! NEMO/OCE 4.0 , NEMO Consortium (2018) !                       !==  grid point position  ==!   (in kilometers) zfact = rn_dx * 1.e-3         ! conversion in km DO jj = 1, jpj DO ji = 1, jpi             ! longitude plamt(ji,jj) = zfact * (  - 0.5 + REAL( ji-1 + nimpp-1 , wp )  ) plamu(ji,jj) = zfact * (          REAL( ji-1 + nimpp-1 , wp )  ) plamv(ji,jj) = plamt(ji,jj) plamf(ji,jj) = plamu(ji,jj) !                       ! latitude pphit(ji,jj) = zfact * (  - 0.5 + REAL( jj-1 + njmpp-1 , wp )  ) pphiu(ji,jj) = pphit(ji,jj) pphiv(ji,jj) = zfact * (          REAL( jj-1 + njmpp-1 , wp )  ) pphif(ji,jj) = pphiv(ji,jj) END DO END DO DO_2D_11_11 !                       ! longitude plamt(ji,jj) = zfact * (  - 0.5 + REAL( ji-1 + nimpp-1 , wp )  ) plamu(ji,jj) = zfact * (          REAL( ji-1 + nimpp-1 , wp )  ) plamv(ji,jj) = plamt(ji,jj) plamf(ji,jj) = plamu(ji,jj) !                       ! latitude pphit(ji,jj) = zfact * (  - 0.5 + REAL( jj-1 + njmpp-1 , wp )  ) pphiu(ji,jj) = pphit(ji,jj) pphiv(ji,jj) = zfact * (          REAL( jj-1 + njmpp-1 , wp )  ) pphif(ji,jj) = pphiv(ji,jj) END_2D ! !                       !==  Horizontal scale factors  ==!   (in meters)
• ## NEMO/trunk/tests/OVERFLOW/MY_SRC/usrdef_zgr.F90

 r12377 PUBLIC   usr_def_zgr   ! called by domzgr.F90 !! * Substitutions #  include "do_loop_substitute.h90" !!---------------------------------------------------------------------- !! NEMO/OCE 4.0 , NEMO Consortium (2018) pe3vw(:,:,jk) = pe3w_1d (jk) END DO DO jj = 1, jpj                      ! bottom scale factors and depth at T- and W-points DO ji = 1, jpi ik = k_bot(ji,jj) pdepw(ji,jj,ik+1) = MIN( zht(ji,jj) , pdepw_1d(ik+1) ) pe3t (ji,jj,ik  ) = pdepw(ji,jj,ik+1) - pdepw(ji,jj,ik) pe3t (ji,jj,ik+1) = pe3t (ji,jj,ik  ) ! pdept(ji,jj,ik  ) = pdepw(ji,jj,ik  ) + pe3t (ji,jj,ik  ) * 0.5_wp pdept(ji,jj,ik+1) = pdepw(ji,jj,ik+1) + pe3t (ji,jj,ik+1) * 0.5_wp pe3w (ji,jj,ik+1) = pdept(ji,jj,ik+1) - pdept(ji,jj,ik)              ! = pe3t (ji,jj,ik  ) END DO END DO DO_2D_11_11 ik = k_bot(ji,jj) pdepw(ji,jj,ik+1) = MIN( zht(ji,jj) , pdepw_1d(ik+1) ) pe3t (ji,jj,ik  ) = pdepw(ji,jj,ik+1) - pdepw(ji,jj,ik) pe3t (ji,jj,ik+1) = pe3t (ji,jj,ik  ) ! pdept(ji,jj,ik  ) = pdepw(ji,jj,ik  ) + pe3t (ji,jj,ik  ) * 0.5_wp pdept(ji,jj,ik+1) = pdepw(ji,jj,ik+1) + pe3t (ji,jj,ik+1) * 0.5_wp pe3w (ji,jj,ik+1) = pdept(ji,jj,ik+1) - pdept(ji,jj,ik)              ! = pe3t (ji,jj,ik  ) END_2D !                                   ! bottom scale factors and depth at  U-, V-, UW and VW-points !                                   ! usually Computed as the minimum of neighbooring scale factors

• ## NEMO/trunk/tests/VORTEX/MY_SRC/usrdef_hgr.F90

 r10074 PUBLIC   usr_def_hgr   ! called by domhgr.F90 !! * Substitutions #  include "do_loop_substitute.h90" !!---------------------------------------------------------------------- !! NEMO/OCE 4.0 , NEMO Consortium (2018) #endif DO jj = 1, jpj DO ji = 1, jpi zti = FLOAT( ji - 1 + nimpp - 1 )          ;  ztj = FLOAT( jj - 1 + njmpp - 1 ) zui = FLOAT( ji - 1 + nimpp - 1 ) + 0.5_wp ;  zvj = FLOAT( jj - 1 + njmpp - 1 ) + 0.5_wp plamt(ji,jj) = zlam0 + rn_dx * 1.e-3 * zti plamu(ji,jj) = zlam0 + rn_dx * 1.e-3 * zui plamv(ji,jj) = plamt(ji,jj) plamf(ji,jj) = plamu(ji,jj) pphit(ji,jj) = zphi0 + rn_dy * 1.e-3 * ztj pphiv(ji,jj) = zphi0 + rn_dy * 1.e-3 * zvj pphiu(ji,jj) = pphit(ji,jj) pphif(ji,jj) = pphiv(ji,jj) END DO END DO DO_2D_11_11 zti = FLOAT( ji - 1 + nimpp - 1 )          ;  ztj = FLOAT( jj - 1 + njmpp - 1 ) zui = FLOAT( ji - 1 + nimpp - 1 ) + 0.5_wp ;  zvj = FLOAT( jj - 1 + njmpp - 1 ) + 0.5_wp plamt(ji,jj) = zlam0 + rn_dx * 1.e-3 * zti plamu(ji,jj) = zlam0 + rn_dx * 1.e-3 * zui plamv(ji,jj) = plamt(ji,jj) plamf(ji,jj) = plamu(ji,jj) pphit(ji,jj) = zphi0 + rn_dy * 1.e-3 * ztj pphiv(ji,jj) = zphi0 + rn_dy * 1.e-3 * zvj pphiu(ji,jj) = pphit(ji,jj) pphif(ji,jj) = pphiv(ji,jj) END_2D ! ! Horizontal scale factors (in meters)