Changeset 10102 for NEMO/branches/2018/dev_r9759_HPC09_ESIWACE/tests/BENCH/MY_SRC

Timestamp:

2018-09-07T15:05:33+02:00 (6 years ago)

Author:

smasson

Message:

dev_r9759_HPC09_ESIWACE: update usrdef_istate.F90 to define unique value on each point

File:

• NEMO/branches/2018/dev_r9759_HPC09_ESIWACE/tests/BENCH/MY_SRC/usrdef_istate.F90 (modified) (3 diffs)

NEMO/branches/2018/dev_r9759_HPC09_ESIWACE/tests/BENCH/MY_SRC/usrdef_istate.F90

@@ -19 +19 @@
    USE in_out_manager ! I/O manager
    USE lib_mpp        ! MPP library
+   USE lbclnk          ! lateral boundary conditions - mpp exchanges
    !   
    USE usrdef_nam
@@ -50 +51 @@
       REAL(wp), DIMENSION(jpi,jpj,jpk)     , INTENT(  out) ::   pv      ! j-component of the velocity  [m/s] 
       REAL(wp), DIMENSION(jpi,jpj)         , INTENT(  out) ::   pssh    ! sea-surface height
+      !
+      REAL(wp), DIMENSION(jpi,jpj) ::   z2d   ! 2D workspace
+      INTEGER :: ji, jj, jk
       !!----------------------------------------------------------------------
       !
@@ -56 +60 @@
       IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~~~   '
       !
+      ! define unique value on each point
+      DO jj = 1, jpj
+         DO ji = 1, jpi
+            z2d(ji,jj) = ( 0.5_wp - REAL( nimpp + ji - 1 + ( njmpp + jj - 2 ) * jpiglo) / REAL ( jpiglo * jpjglo ) ) * 0.1_wp
+         ENDDO
+      ENDDO
+      !
       ! sea level:
-      pssh(:,:) = 0.
-      ! temperature:
-      pts(:,:,:,jp_tem) = 10._wp
-      ! salinity:  
-      pts(:,:,:,jp_sal) = 35._wp
-      ! velocities:
-      pu(:,:,:) = 0.
-      pv(:,:,:) = 0.
+      pssh(:,:) = z2d(:,:)
       !
-
+      DO jk = 1, jpk
+         ! temperature:
+         pts(:,:,jk,jp_tem) = 10._wp + z2d(:,:) * 2._wp
+         ! salinity:  
+         pts(:,:,jk,jp_sal) = 35._wp + z2d(:,:)
+         ! velocities:
+         pu(:,:,jk) = z2d(:,:) * 0.1_wp
+         pv(:,:,jk) = z2d(:,:) * 0.01_wp
+      ENDDO
+      !
+      CALL lbc_lnk("usr_def_istate", pssh, 'T',  1. )            ! apply boundary conditions
+      CALL lbc_lnk("usr_def_istate",  pts, 'T',  1. )            ! apply boundary conditions
+      CALL lbc_lnk("usr_def_istate",   pu, 'U', -1. )            ! apply boundary conditions
+      CALL lbc_lnk("usr_def_istate",   pv, 'V', -1. )            ! apply boundary conditions
+      
    END SUBROUTINE usr_def_istate