Changeset 5196 for branches/2012/dev_v3_4_STABLE_2012/NEMOGCM/NEMO/OPA_SRC/SBC/sbcmod.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/SBC/sbcmod.F90 (modified) (3 diffs)

branches/2012/dev_v3_4_STABLE_2012/NEMOGCM/NEMO/OPA_SRC/SBC/sbcmod.F90

@@ -101 +101 @@
         IF( lk_cice )   nn_ice      = 4
       ENDIF
-      IF( cp_cfg == 'gyre' ) THEN            ! GYRE configuration
+      IF( cp_cfg == 'gyre' .OR. cp_cfg == 'seabass' ) THEN            ! GYRE configuration
           ln_ana      = .TRUE.   
           nn_ice      =   0
@@ -178 +178 @@
       IF( ln_blk_mfs      ) THEN   ;   nsbc =  6   ; icpt = icpt + 1   ;   ENDIF       ! MFS  bulk       formulation
       IF( ln_cpl          ) THEN   ;   nsbc =  5   ; icpt = icpt + 1   ;   ENDIF       ! Coupled         formulation
-      IF( cp_cfg == 'gyre') THEN   ;   nsbc =  0                       ;   ENDIF       ! GYRE analytical formulation
+      IF( cp_cfg == 'gyre' .OR. cp_cfg == 'seabass' ) THEN   ;   nsbc =  0                       ;   ENDIF       ! GYRE analytical formulation
       IF( lk_esopa        )            nsbc = -1                                       ! esopa test, ALL formulations
       !
@@ -256 +256 @@
       SELECT CASE( nsbc )                                ! Compute ocean surface boundary condition
       !                                                  ! (i.e. utau,vtau, qns, qsr, emp, emps)
-      CASE(  0 )   ;   CALL sbc_gyre    ( kt )                    ! analytical formulation : GYRE configuration
+      CASE(  0 )   ;                                     ! analytical formulation : GYRE or 
+                                                         ! SEABASS configurations
+         IF (cp_cfg == 'gyre') THEN
+            CALL sbc_gyre    ( kt )                    ! analytical formulation : GYRE configuration
+         ELSEIF (cp_cfg == 'seabass') THEN
+            CALL sbc_seabass     ( kt )
+         ENDIF
       CASE(  1 )   ;   CALL sbc_ana     ( kt )                    ! analytical formulation : uniform sbc
       CASE(  2 )   ;   CALL sbc_flx     ( kt )                    ! flux formulation