Changeset 6683

Timestamp:

2016-06-10T16:48:58+02:00 (8 years ago)

Author:

frrh

Message:

Changes for compatibility with UKMO/dev_r5518_GSI7_GSI8_landice_bitcomp_medusa.

File:

• branches/UKMO/dev_r5518_RH_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC/MEDUSA/trcbio_medusa.F90 (modified) (6 diffs)

branches/UKMO/dev_r5518_RH_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC/MEDUSA/trcbio_medusa.F90

@@ -58 +58 @@
    !!   trc_bio_medusa        :  
    !!----------------------------------------------------------------------
+      USE oce, ONLY: CO2Flux_out_cpl, DMS_out_cpl, PCO2a_in_cpl
       USE oce_trc
       USE trc
@@ -96 +97 @@
       PUBLIC   trc_bio_medusa    ! called in ???
 
+      REAL(wp), PARAMETER :: dms_unit_conv = 1.0e+6 ! Coversion factor to get 
+                                                    ! outgong DMS in standard 
+                                                    ! units for coupling, i.e.
+                                               ! specifically nmol/L 
+                                                    !(= umol/m3)
    !!* Substitution
 #  include "domzgr_substitute.h90"
@@ -357 +363 @@
       REAL(wp), POINTER, DIMENSION(:,:  ) :: oben_fe2d, oben_c2d, oben_si2d, oben_ca2d, sfr_ocal2d
       REAL(wp), POINTER, DIMENSION(:,:  ) :: sfr_oarg2d, lyso_ca2d 
+
 # endif
       !! 2D var for diagnostics.
@@ -1538 +1545 @@
                        ENDIF
 #   if defined key_debug_medusa
-                       IF (lwp) write (numout,*) 'trc_bio_medusa: finnish calculating dms'
+                       IF (lwp) write (numout,*) 'trc_bio_medusa: finish calculating dms'
                      CALL flush(numout)
 #   endif 
@@ -4278 +4285 @@
                  !!-----------------------------------------
                  IF (jk.eq.1) THEN
+
+#if defined key_oasis3
+                     ! Capture outgoing CO2 flux coupling field
+                     CO2Flux_out_cpl(:,:) = f_co2flux2d(:,:) 
+                     ! Ensure halos are populated in the couping field
+                     CALL lbc_lnk(CO2Flux_out_cpl(:,:),'T',1. )
+
+                     ! Capture outgoing DMS coupling field and convert
+                     ! values to standard units
+                     DMS_out_cpl(:,:) = dms_surf2d(:,:) * dms_unit_conv
+                     ! Ensure halos are populated in the couping field
+                     CALL lbc_lnk(DMS_out_cpl(:,:),'T',1. )
+#endif
+
 #   if defined key_debug_medusa
                      IF (lwp) write (numout,*) 'trc_bio_medusa: diag jk = 1'
@@ -4368 +4389 @@
                        IF( med_diag%DMS_SURF%dgsave ) THEN
                          CALL iom_put( "DMS_SURF"  , dms_surf2d )
+
+
                           CALL wrk_dealloc( jpi, jpj,   dms_surf2d   )
                        ENDIF