Changeset 6906 for branches/2016/dev_r6409_SIMPLIF_2_usrdef/NEMOGCM/CONFIG/LOCK_EXCHANGE/MY_SRC/usrdef_istate.F90

Timestamp:

2016-09-02T15:06:17+02:00 (8 years ago)

Author:

flavoni

Message:

#1692 and ROBUST-3 : commit zgr and some update for LOCK EXCHANGE case

File:

• branches/2016/dev_r6409_SIMPLIF_2_usrdef/NEMOGCM/CONFIG/LOCK_EXCHANGE/MY_SRC/usrdef_istate.F90 (modified) (3 diffs)

branches/2016/dev_r6409_SIMPLIF_2_usrdef/NEMOGCM/CONFIG/LOCK_EXCHANGE/MY_SRC/usrdef_istate.F90

@@ -14 +14 @@
    !!----------------------------------------------------------------------
    USE par_oce        ! ocean space and time domain
+   USE dom_oce , ONLY : glamt 
    USE phycst         ! physical constants
    !
@@ -49 +50 @@
       !
       INTEGER  ::   jk     ! dummy loop indices
+      REAL(wp) ::   zdam   ! location of dam [Km]
       !!----------------------------------------------------------------------
       !
       IF(lwp) WRITE(numout,*)
-      IF(lwp) WRITE(numout,*) 'usr_def_istate : LOCK_EXCHANGE configuration, analytical definition of initial state '
-      IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~~~   Ocean at rest, with a dam on T profile, and uniform S profile'
+      IF(lwp) WRITE(numout,*) 'usr_def_istate : LOCK_EXCHANGE configuration, analytical definition of initial state'
+      IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~~~   Ocean at rest, with a constant salinity (not used as rho=F(T) '
+      IF(lwp) WRITE(numout,*) '                 and a vertical density front with a 5 kg/m3 difference located at glam=32km'
+      IF(lwp) WRITE(numout,*) '                 (i.e. a temperature difference of 25 degrees with rn_a0 = 0.2'
+      !
+      !  rn_a0 =  2.e-1   !  thermal expension coefficient (nn_eos= 1)
+      !  rho = rau0 - rn_a0 * (T-10) 
+      !  delta_T = 25 degrees  ==>>  delta_rho = 25 * rn_a0 = 5 kg/m3
       !
       pu  (:,:,:) = 0._wp        ! ocean at rest
@@ -60 +68 @@
       !
       !                          ! T & S profiles
-      pts(:,:,:,jp_tem) = 30._wp * tmask(:,:,:)
+      zdam = 32.                      ! in kilometers
+      pts(:,:,:,jp_tem) = 30._wp * ptmask(:,:,:)
       DO jk = 1, jpkm1
-!SF   ADD rn_dx2
-         WHERE( gphit(:,:) <= zdam )   pts(:,:,jk,jp_tem) = 5._wp * tmask(:,:,jk)
+         WHERE( glamt(:,:) <= zdam )   pts(:,:,jk,jp_tem) = 5._wp * ptmask(:,:,jk)
       END DO
       !
-      pts(:,:,:,jp_sal) = 35. * tmask(:,:,:)
+      pts(:,:,:,jp_sal) = 35._wp * ptmask(:,:,:)
       !   
    END SUBROUTINE usr_def_istate