Changeset 6256

Timestamp:

2016-01-15T10:59:16+01:00 (8 years ago)

Author:

frrh

Message:

Add some TEMPORARY code to facilitate testing MEDUSA coupling
interface fields with UM atmosphere.

File:

• branches/UKMO/nemo_v3_6_STABLE_pkg/NEMOGCM/NEMO/OPA_SRC/SBC/sbccpl.F90 (modified) (4 diffs)

branches/UKMO/nemo_v3_6_STABLE_pkg/NEMOGCM/NEMO/OPA_SRC/SBC/sbccpl.F90

@@ -53 +53 @@
 #endif
 #if defined key_medusa
-   USE trc
+! RSRH: Temporarily commented out until we get MEDUSA code included
+!   USE trc
 #endif
 
@@ -2056 +2057 @@
       REAL(wp), POINTER, DIMENSION(:,:)   ::   zfr_l, ztmp1, ztmp2, zotx1, zoty1, zotz1, zitx1, zity1, zitz1
       REAL(wp), POINTER, DIMENSION(:,:,:) ::   ztmp3, ztmp4   
+
+#if defined key_medusa 
+! RSRH Temporarily definitition of array used in medusa coupling until we actually have the medusa code included in the build 
+      REAL(wp)  ::  trc2d_TEST(jpi,jpj) 
+#endif 
       !!----------------------------------------------------------------------
       !
@@ -2272 +2278 @@
 !! May add also a coupling MED-UKCA key
 
+! RSRH set up a test field purely for use during development:
+      trc2d_TEST(:,:) = 0.0 
+
       ! RSRH. We don't want to use magic numbers in the final code (i.e. 98 and 221).
       ! These need moving to a parameter statement (as part of MEDUSA code) or specifying in a namelist
       ! so the following code MUST NOT be viewed as anything more than temporary.
       IF ( ssnd(jps_bio_co2)%laction ) THEN
-         CALL cpl_prism_snd( jps_bio_co2, isec, trc2d(:,:,98:98), info )
+!         CALL cpl_prism_snd( jps_bio_co2, isec, trc2d(:,:,98:98), info )
+! RSRH Send our test field in place of the real bio_co2
+         CALL cpl_prism_snd( jps_bio_co2, isec, trc2d_TEST(:,:), info )
       END IF
 
@@ -2282 +2293 @@
           ! We need to multiply DMS by a conversion factor to get values in the standard units expected in
           ! the coupling space.
-          ztmp1(:,: ) = trc2d(:,:,221) * dms_unit_conv
+!          ztmp1(:,: ) = trc2d(:,:,221) * dms_unit_conv
+!RSRH Use our test field in place of the real bio_dms 
+         ztmp1(:,: ) = trc2d_TEST(:,:) * dms_unit_conv
          CALL cpl_prism_snd( jps_bio_dms, isec, RESHAPE ( ztmp1, (/jpi,jpj,1/) ), info )
       ENDIF