Changeset 5196 for branches/2012/dev_v3_4_STABLE_2012/NEMOGCM/NEMO/OPA_SRC/DOM/istate.F90

Timestamp:

2015-04-07T10:28:07+02:00 (9 years ago)

Author:

pabouttier

Message:

Add SEABASS reference configuration for this NEMO version for now; See Ticket #1505

File:

• branches/2012/dev_v3_4_STABLE_2012/NEMOGCM/NEMO/OPA_SRC/DOM/istate.F90 (modified) (2 diffs)

branches/2012/dev_v3_4_STABLE_2012/NEMOGCM/NEMO/OPA_SRC/DOM/istate.F90

@@ -113 +113 @@
          ELSEIF( cp_cfg == 'gyre' ) THEN
             CALL istate_gyre                     ! GYRE  configuration : start from pre-defined T-S fields
+         ELSEIF( cp_cfg == 'seabass' ) THEN
+            CALL istate_seabass
          ELSEIF( ln_tsd_init      ) THEN         ! Initial T-S fields read in files
             CALL dta_tsd( nit000, tsb )                  ! read 3D T and S data at nit000
@@ -412 +414 @@
 
 
+
+   SUBROUTINE istate_seabass
+      !!----------------------------------------------------------------------
+      !!                   ***  ROUTINE istate_seabass  ***
+      !!
+      !! ** Purpose :   Initialization of the dynamics and tracers for seabass
+      !!      configuration (double gyre)
+      !!
+      !! ** Method  : - set temperature field following Chassignet and Gent, JPO
+      !!                    21, pp1290-1299, 1991, and the law
+      !!                    rho/rho0=1-2.e-4(T-T0)
+      !!              - set salinity field constant
+      !!
+      !!----------------------------------------------------------------------
+      !! * Local variables
+      INTEGER :: ji, jj, jk     ! dummy loop indices
+      REAL(wp) ::   zsal = 35.5
+      !!----------------------------------------------------------------------
+
+      IF(lwp) WRITE(numout,*)
+      IF(lwp) WRITE(numout,*) 'istate_seabass : initial analytical T and constant S profiles '
+      IF(lwp) WRITE(numout,*) '~~~~~~~~~~~'
+
+      DO jk = 1, jpk
+         DO jj = 1, jpj
+            DO ji = 1, jpi
+               tsn(ji,jj,jk,jp_tem) = ( 25.+5.9e-5*800./9.81/2.e-4*   &
+               &    (exp(-gdept_0(jk)/800.)-1.))  * tmask(ji,jj,jk)
+               tsb(ji,jj,jk,jp_tem) = tsn(ji,jj,jk, jp_tem)
+          END DO
+        END DO
+      END DO
+
+      tsn(:,:,:,jp_sal) = zsal  * tmask(:,:,:)
+      tsb(:,:,:,jp_sal) = tsn(:,:,:,jp_sal)
+
+      IF(lwp) THEN
+         WRITE(numout,*)
+         WRITE(numout,*) '              Initial temperature and salinity profiles:'
+         WRITE(numout, "(9x,' level   gdept   temperature   salinity   ')" )
+         WRITE(numout, "(10x, i4, 3f10.2)" ) ( jk, gdept_0(jk), tsn(2,2,jk,jp_tem), tsn(2,2,jk,jp_sal), jk = 1, jpk )
+      ENDIF
+
+
+   END SUBROUTINE istate_seabass
+
+
+
    SUBROUTINE istate_uvg
       !!----------------------------------------------------------------------