Changeset 8671 for branches/NERC/dev_r5518_GO6_ScalingCoupledChl/NEMOGCM/NEMO/TOP_SRC/MEDUSA

Timestamp:

2017-11-03T13:41:36+01:00 (7 years ago)

Author:

jpalmier

Message:

Update MLD Chl - create a coupled var output group in field_def_bgc

Location:

branches/NERC/dev_r5518_GO6_ScalingCoupledChl/NEMOGCM/NEMO/TOP_SRC/MEDUSA

Files:

• bio_medusa_diag_slice.F90 (modified) (3 diffs)
• bio_medusa_fin.F90 (modified) (5 diffs)
• trcnam_medusa.F90 (modified) (1 diff)

branches/NERC/dev_r5518_GO6_ScalingCoupledChl/NEMOGCM/NEMO/TOP_SRC/MEDUSA/bio_medusa_diag_slice.F90

@@ -39 +39 @@
 # endif
       USE lbclnk,            ONLY: lbc_lnk
-      USE trc,               ONLY: trn
-      USE oce,               ONLY: CO2Flux_out_cpl, DMS_out_cpl, chloro_out_cpl
+      USE oce,               ONLY: CO2Flux_out_cpl, DMS_out_cpl
       USE par_oce,           ONLY: jpi, jpj
       USE sbc_oce,           ONLY: lk_oasis, qsr, wndm
@@ -48 +47 @@
                                    jdms, ocal_ccd, xpar, xze,              &
                                    zb_co2_flx, zb_dms_srf,                 &
-                                   zn_co2_flx, zn_dms_srf, zn_chl_srf,     &
-                                   scl_chl, chl_out
+                                   zn_co2_flx, zn_dms_srf
       USE trc,               ONLY: med_diag
 
@@ -67 +65 @@
       !!
       IF (jk.eq.1) THEN
-         !! JPALM -- 02-06-2017 --
-         !! add Chl surf coupling
-         !! no need to output, just pass to cpl var
-         IF (lk_oasis) THEN
-            IF (chl_out.eq.1) THEN
-               !! export and scale surface chl
-               zn_chl_srf(:,:) = MAX( 0.0, (trn(:,:,1,jpchd) + trn(:,:,1,jpchn)) * scl_chl * 1.0E-6 )
-                                 !! surf Chl in Kg-chl/m3 as needed for cpl
-            ELSEIF (chl_out.eq.2) THEN
-               !! export and scale mld chl
-               zn_chl_srf(:,:) = MAX( 0.0, fchl_ml * scl_chl * 1.0E-6 )
-                                 !! mld Chl in Kg-chl/m3 as needed for cpl
-            ENDIF
-            chloro_out_cpl(:,:) = zn_chl_srf(:,:)        !! Coupling Chl
-         END IF
          IF( med_diag%MED_QSR%dgsave ) THEN
             CALL iom_put( "MED_QSR"  , qsr ) !

branches/NERC/dev_r5518_GO6_ScalingCoupledChl/NEMOGCM/NEMO/TOP_SRC/MEDUSA/bio_medusa_fin.F90

@@ -39 +39 @@
 # endif
       USE lbclnk,            ONLY: lbc_lnk
+      USE oce,               ONLY: chloro_out_cpl 
       USE par_medusa,        ONLY: jp_medusa_2d, jp_medusa_3d,          &
-                                   jp_medusa_trd
+                                   jp_medusa_trd, jpchd, jpchn
       USE par_oce,           ONLY: jpi, jpim1, jpj, jpjm1, jpk
       USE phycst,            ONLY: rsmall
+      USE sbc_oce,           ONLY: lk_oasis
       USE sms_medusa,        ONLY: jinorgben, jorgben,                  &
                                    f3_co3, f3_h2co3, f3_hco3,           &
@@ -51 +53 @@
                                    zb_sed_n, zb_sed_si,                 &
                                    zn_sed_c, zn_sed_ca, zn_sed_fe,      &
-                                   zn_sed_n, zn_sed_si, zn_chl_srf
-      USE trc,               ONLY: med_diag, nittrc000 
+                                   zn_sed_n, zn_sed_si, zn_chl_srf,     &
+                                   scl_chl, chl_out
+      USE trc,               ONLY: med_diag, nittrc000, trn 
       USE trcnam_trp,        ONLY: ln_trcadv_cen2, ln_trcadv_tvd
  
@@ -231 +234 @@
             ENDDO
          ENDDO
+
+         !!!---------------------------------------------------------------
+         !! Calculates Chl diag for UM coupling 
+         !!!---------------------------------------------------------------
+         !! JPALM -- 02-06-2017 --
+         !! add Chl surf coupling
+         !! no need to output, just pass to cpl var
+         IF (lk_oasis) THEN
+            IF (chl_out.eq.1) THEN
+               !! export and scale surface chl
+               zn_chl_srf(:,:) = MAX( 0.0, (trn(:,:,1,jpchd) + trn(:,:,1,jpchn)) * scl_chl * 1.0E-6 )
+                                 !! surf Chl in Kg-chl/m3 as needed for cpl
+            ELSEIF (chl_out.eq.2) THEN
+               !! export and scale mld chl
+               zn_chl_srf(:,:) = MAX( 0.0, fchl_ml(:,:) * scl_chl * 1.0E-6 )
+                                 !! mld Chl in Kg-chl/m3 as needed for cpl
+            ENDIF
+            chloro_out_cpl(:,:) = zn_chl_srf(:,:)        !! Coupling Chl
+         END IF
+
          !!----------------------------------------------------------------
          !! Add in XML diagnostics stuff
@@ -241 +264 @@
 #   endif
          IF ( med_diag%INVTN%dgsave ) THEN
-            !! HACK - Coupled Chl sent to UM    
-            CALL iom_put( "INVTN"  , zn_chl_srf )
-            !! CALL iom_put( "INVTN"  , ftot_n )
+            CALL iom_put( "INVTN"  , ftot_n )
          ENDIF
          IF ( med_diag%INVTSI%dgsave ) THEN
@@ -262 +283 @@
          IF ( med_diag%CHL_MLD%dgsave ) THEN
             CALL iom_put( "CHL_MLD"  , fchl_ml )
+         ENDIF
+         IF ( med_diag%CHL_CPL%dgsave ) THEN
+            CALL iom_put( "CHL_CPL"  , zn_chl_srf )
          ENDIF
          IF ( med_diag%PN_JLIM%dgsave ) THEN

branches/NERC/dev_r5518_GO6_ScalingCoupledChl/NEMOGCM/NEMO/TOP_SRC/MEDUSA/trcnam_medusa.F90

@@ -2079 +2079 @@
           med_diag%CHL_MLD%dgsave = .FALSE.
       ENDIF
+      IF  (iom_use("CHL_CPL")) THEN 
+          med_diag%CHL_CPL%dgsave = .TRUE.
+      ELSE 
+          med_diag%CHL_CPL%dgsave = .FALSE.
+      ENDIF
       !! 3D
       IF  (iom_use("TPP3")) THEN