Changeset 8436 for branches/UKMO/dev_r5518_GO6_package_asm_surf_bgc/NEMOGCM/NEMO/OPA_SRC/ASM/asminc.F90

Timestamp:

2017-08-14T15:22:09+02:00 (7 years ago)

Author:

dford

Message:

Implement initial version of surface chlorophyll assimilation for MEDUSA.

File:

• branches/UKMO/dev_r5518_GO6_package_asm_surf_bgc/NEMOGCM/NEMO/OPA_SRC/ASM/asminc.F90 (modified) (3 diffs)

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

@@ -126 +126 @@
    REAL(wp), PUBLIC, DIMENSION(:,:,:,:), ALLOCATABLE :: logchl_balinc  !: Increment to BGC variables from logchl assim
 #endif
+#if defined key_hadocc || (defined key_medusa && defined key_foam_medusa)
+   REAL(wp), DIMENSION(:,:), ALLOCATABLE :: pgrow_avg_bkg  !: Background phyto growth for logchl balancing
+   REAL(wp), DIMENSION(:,:), ALLOCATABLE :: ploss_avg_bkg  !: Background phyto loss for logchl balancing
+   REAL(wp), DIMENSION(:,:), ALLOCATABLE :: phyt_avg_bkg   !: Background phyto for logchl balancing
+   REAL(wp), DIMENSION(:,:), ALLOCATABLE :: mld_max_bkg    !: Background max MLD for logchl balancing
+#endif
    REAL(wp) :: rn_maxchlinc = -999.0  !: maximum absolute non-log chlorophyll increment from logchl assimilation
                                       !: <= 0 implies no maximum applied (switch turned off)
@@ -620 +626 @@
 
       ENDIF
+
+#if defined key_hadocc || (defined key_medusa && defined key_foam_medusa)
+      IF ( ln_logchlinc ) THEN
+
+         ALLOCATE( pgrow_avg_bkg(jpi,jpj) )
+         ALLOCATE( ploss_avg_bkg(jpi,jpj) )
+         ALLOCATE( phyt_avg_bkg(jpi,jpj)  )
+         ALLOCATE( mld_max_bkg(jpi,jpj)   )
+         pgrow_avg_bkg(:,:) = 0.0
+         ploss_avg_bkg(:,:) = 0.0
+         phyt_avg_bkg(:,:)  = 0.0
+         mld_max_bkg(:,:)   = 0.0
+         
+         IF ( ln_logchlbal ) THEN
+         
+            !--------------------------------------------------------------------
+            ! Read background variables for logchl balancing
+            !--------------------------------------------------------------------
+
+            CALL iom_open( c_asmbkg, inum )
+
+            CALL iom_get( inum, jpdom_autoglo, 'pgrow_avg', pgrow_avg_bkg )
+            CALL iom_get( inum, jpdom_autoglo, 'ploss_avg', ploss_avg_bkg )
+            CALL iom_get( inum, jpdom_autoglo, 'phyt_avg',  phyt_avg_bkg  )
+            CALL iom_get( inum, jpdom_autoglo, 'mld_max',   mld_max_bkg   )
+            pgrow_avg_bkg(:,:) = pgrow_avg_bkg(:,:) * tmask(:,:,1)
+            ploss_avg_bkg(:,:) = ploss_avg_bkg(:,:) * tmask(:,:,1)
+            phyt_avg_bkg(:,:)  = phyt_avg_bkg(:,:)  * tmask(:,:,1)
+            mld_max_bkg(:,:)   = mld_max_bkg(:,:)   * tmask(:,:,1)
+
+            CALL iom_close( inum )
+         
+         ENDIF
+      
+      ENDIF
+#endif
       !
    END SUBROUTINE asm_inc_init
@@ -1266 +1308 @@
 
 #if defined key_medusa && defined key_foam_medusa
-         !CALL asm_logchl_bal_medusa()
-         CALL ctl_stop( 'Attempting to assimilate logchl into MEDUSA, ', &
-            &           'but not fully implemented yet' )
+         CALL asm_logchl_bal_medusa( logchl_bkginc, zincper, mld_choice_bgc, &
+            &                        rn_maxchlinc, ln_logchlbal,             &
+            &                        pgrow_avg_bkg, ploss_avg_bkg,           &
+            &                        phyt_avg_bkg, mld_max_bkg,              &
+            &                        logchl_balinc )
 #elif defined key_hadocc
          CALL asm_logchl_bal_hadocc( logchl_bkginc, zincper, mld_choice_bgc, &
-            &                        rn_maxchlinc, ln_logchlbal, logchl_balinc )
+            &                        rn_maxchlinc, ln_logchlbal,             &
+            &                        pgrow_avg_bkg, ploss_avg_bkg,           &
+            &                        phyt_avg_bkg, mld_max_bkg,              &
+            &                        logchl_balinc )
 #else
          CALL ctl_stop( 'Attempting to assimilate logchl, ', &