Changeset 13443 for NEMO/branches/UKMO/NEMO_4.0.1_NGMS_couple_stage2/src/OCE/SBC/sbcmod.F90

Timestamp:

2020-08-28T17:01:04+02:00 (4 years ago)

Author:

frrh

Message:

Commit changes which allow first working 2 way coupling exchanges between GO8
and LFRic aquaplanet.

File:

• NEMO/branches/UKMO/NEMO_4.0.1_NGMS_couple_stage2/src/OCE/SBC/sbcmod.F90 (modified) (7 diffs)

NEMO/branches/UKMO/NEMO_4.0.1_NGMS_couple_stage2/src/OCE/SBC/sbcmod.F90

@@ -98 +98 @@
          &             ln_rnf   , nn_fwb   , ln_ssr   , ln_isf    , ln_apr_dyn ,     &
          &             ln_wave  , ln_cdgw  , ln_sdw   , ln_tauwoc  , ln_stcor   ,     &
-         &             ln_tauw  , nn_lsm, nn_sdrift
+         &             ln_tauw  , nn_lsm, nn_sdrift, ln_couple_test
       !!----------------------------------------------------------------------
       !
@@ -139 +139 @@
          WRITE(numout,*) '      Type of coupling (Ocean/Ice/Atmosphere) : '
          WRITE(numout,*) '         ocean-atmosphere coupled formulation       ln_cpl        = ', ln_cpl
+         WRITE(numout,*) '         ocean-atmos test coupling ONLY             ln_couple_test   = ', ln_couple_test
          WRITE(numout,*) '         mixed forced-coupled     formulation       ln_mixcpl     = ', ln_mixcpl
 !!gm  lk_oasis is controlled by key_oasis3  ===>>>  It shoud be removed from the namelist 
@@ -208 +209 @@
       END SELECT
       !                             !* coupled options
-      IF( ln_cpl ) THEN
+      IF( ln_cpl .OR. ln_couple_test ) THEN
          IF( .NOT. lk_oasis )   CALL ctl_stop( 'sbc_init : coupled mode with an atmosphere model (ln_cpl=T)',   &
             &                                  '           required to defined key_oasis3' )
@@ -275 +276 @@
       IF( icpt /= 1 )    CALL ctl_stop( 'sbc_init : choose ONE and only ONE sbc option' )
       !
-      IF(lwp) THEN                     !- print the choice of surface flux formulation
+!RSRH      IF(lwp) THEN                     !- print the choice of surface flux formulation
          WRITE(numout,*)
          SELECT CASE( nsbc )
@@ -287 +288 @@
          END SELECT
          IF( ll_not_nemo  )       WRITE(numout,*) '               + OASIS coupled SAS'
-      ENDIF
+!RSRH      ENDIF
       !
       !                             !* OASIS initialization
@@ -421 +422 @@
       !                                            !==  sbc formulation  ==!
       !                                                   
+WRITE(numout,*) "RSRH nsbc setting=", nsbc ;flush(numout)
+      
+
+
       SELECT CASE( nsbc )                                ! Compute ocean surface boundary condition
       !                                                  ! (i.e. utau,vtau, qns, qsr, emp, sfx)
@@ -432 +437 @@
       CASE( jp_none    )
          IF( ll_opa    )       CALL sbc_cpl_rcv   ( kt, nn_fsbc, nn_ice )   ! OPA-SAS coupling: OPA receiving fields from SAS
+
       END SELECT
       !
       IF( ln_mixcpl )          CALL sbc_cpl_rcv   ( kt, nn_fsbc, nn_ice )   ! forced-coupled mixed formulation after forcing
+
+      ! RSRH for test/dev only  
+      IF( ln_couple_test )      CALL sbc_cpl_rcv   ( kt, nn_fsbc, nn_ice )   ! Purely for testing/development of
+                                                                             ! coupled models if none of the above
+                                                                             ! cases are true. In this case we are likely
+                                                                             ! be running a forced configuration and
+                                                                             ! simply exchanging dummy fields which
+                                                                             ! have no bearing on model evolution.         
       !
       IF ( ln_wave .AND. (ln_tauwoc .OR. ln_tauw) ) CALL sbc_wstress( )      ! Wind stress provided by waves