Changeset 8648 for branches/UKMO/dev_r5518_GO6_package_asm_surf_bgc_v2/NEMOGCM/NEMO/OPA_SRC/ASM/asminc.F90

Timestamp:

2017-10-20T15:31:59+02:00 (7 years ago)

Author:

dford

Message:

Change source of MLD for direct initialisation of logchl and pCO2 to work.

File:

• branches/UKMO/dev_r5518_GO6_package_asm_surf_bgc_v2/NEMOGCM/NEMO/OPA_SRC/ASM/asminc.F90 (modified) (7 diffs)

branches/UKMO/dev_r5518_GO6_package_asm_surf_bgc_v2/NEMOGCM/NEMO/OPA_SRC/ASM/asminc.F90

@@ -750 +750 @@
             CALL iom_get( inum, jpdom_autoglo, 'medusa_alk', tracer_bkg(:,:,:,jpalk) )
 #endif
+            CALL iom_get( inum, jpdom_autoglo, 'mld_max',   mld_max_bkg   )
+            mld_max_bkg(:,:) = mld_max_bkg(:,:) * tmask(:,:,1)
          ENDIF
 
@@ -761 +763 @@
 
          ALLOCATE( tracer_bkg(jpi,jpj,jpk,jptra) )
+         ALLOCATE( mld_max_bkg(jpi,jpj)          )
          tracer_bkg(:,:,:,:) = 0.0
+         mld_max_bkg(:,:)    = 0.0
 
          CALL iom_open( c_asmbkg, inum )
@@ -772 +776 @@
          CALL iom_get( inum, jpdom_autoglo, 'medusa_alk', tracer_bkg(:,:,:,jpalk) )
 #endif
+         CALL iom_get( inum, jpdom_autoglo, 'mld_max',   mld_max_bkg   )
 
          CALL iom_close( inum )
@@ -778 +783 @@
             tracer_bkg(:,:,:,jt) = tracer_bkg(:,:,:,jt) * tmask(:,:,:)
          END DO
+         mld_max_bkg(:,:) = mld_max_bkg(:,:) * tmask(:,:,1)
       
       ENDIF
@@ -1428 +1434 @@
 #if defined key_medusa && defined key_foam_medusa
          CALL asm_logchl_bal_medusa( logchl_bkginc, zincper, mld_choice_bgc, &
-            &                        rn_maxchlinc, ln_logchlbal,             &
+            &                        rn_maxchlinc, ln_logchlbal, ln_asmdin,  &
             &                        pgrow_avg_bkg, ploss_avg_bkg,           &
             &                        phyt_avg_bkg, mld_max_bkg,              &
@@ -1434 +1440 @@
 #elif defined key_hadocc
          CALL asm_logchl_bal_hadocc( logchl_bkginc, zincper, mld_choice_bgc, &
-            &                        rn_maxchlinc, ln_logchlbal,             &
+            &                        rn_maxchlinc, ln_logchlbal, ln_asmdin,  &
             &                        pgrow_avg_bkg, ploss_avg_bkg,           &
             &                        phyt_avg_bkg, mld_max_bkg,              &
@@ -1637 +1643 @@
 
          ! Select mixed layer
-         SELECT CASE( mld_choice_bgc )
-         CASE ( 1 )                   ! Turbocline/mixing depth [W points]
-            zmld(:,:) = hmld(:,:)
-         CASE ( 2 )                   ! Density criterion (0.01 kg/m^3 change from 10m) [W points]
-            zmld(:,:) = hmlp(:,:)
-         CASE ( 3 )                   ! Kara MLD [Interpolated]
+         IF ( ln_asmdin ) THEN
+            CALL ctl_warn( ' Doing direct initialisation with pCO2 assimilation', &
+               &           ' Mixed layer depth taken to be background maximum mld_max_bkg' )
+            zmld(:,:) = mld_max_bkg(:,:)
+         ELSE
+            SELECT CASE( mld_choice_bgc )
+            CASE ( 1 )                   ! Turbocline/mixing depth [W points]
+               zmld(:,:) = hmld(:,:)
+            CASE ( 2 )                   ! Density criterion (0.01 kg/m^3 change from 10m) [W points]
+               zmld(:,:) = hmlp(:,:)
+            CASE ( 3 )                   ! Kara MLD [Interpolated]
 #if defined key_karaml
-            IF ( ln_kara ) THEN
-               zmld(:,:) = hmld_kara(:,:)
-            ELSE
+               IF ( ln_kara ) THEN
+                  zmld(:,:) = hmld_kara(:,:)
+               ELSE
+                  CALL ctl_stop( ' Kara mixed layer requested for pCO2 assimilation,', &
+                     &           ' but ln_kara=.false.' )
+               ENDIF
+#else
                CALL ctl_stop( ' Kara mixed layer requested for pCO2 assimilation,', &
-                  &           ' but ln_kara=.false.' )
-            ENDIF
-#else
-            CALL ctl_stop( ' Kara mixed layer requested for pCO2 assimilation,', &
-               &           ' but is not defined' )
-#endif
-         CASE ( 4 )                   ! Temperature criterion (0.2 K change from surface) [T points]
-            !zmld(:,:) = hmld_tref(:,:)
-            CALL ctl_stop( ' hmld_tref mixed layer requested for pCO2 assimilation,', &
-               &           ' but is not available in this version' )
-         CASE ( 5 )                   ! Density criterion (0.01 kg/m^3 change from 10m) [T points]
-            zmld(:,:) = hmlpt(:,:)
-         END SELECT
+                  &           ' but is not defined' )
+#endif
+            CASE ( 4 )                   ! Temperature criterion (0.2 K change from surface) [T points]
+               !zmld(:,:) = hmld_tref(:,:)
+               CALL ctl_stop( ' hmld_tref mixed layer requested for pCO2 assimilation,', &
+                  &           ' but is not available in this version' )
+            CASE ( 5 )                   ! Density criterion (0.01 kg/m^3 change from 10m) [T points]
+               zmld(:,:) = hmlpt(:,:)
+            END SELECT
+         ENDIF
 
          ! Propagate through mixed layer