Changeset 10045 for branches/NERC/dev_r5518_GO6_under_ice_relax/NEMOGCM/NEMO/TOP_SRC/MEDUSA/bio_medusa_init.F90

Timestamp:

2018-08-07T14:19:25+02:00 (6 years ago)

Author:

jpalmier

Message:

Andrew's changes to add the OMIP double_DIC (activated with key_omip_dic)

File:

• branches/NERC/dev_r5518_GO6_under_ice_relax/NEMOGCM/NEMO/TOP_SRC/MEDUSA/bio_medusa_init.F90 (modified) (7 diffs)

branches/NERC/dev_r5518_GO6_under_ice_relax/NEMOGCM/NEMO/TOP_SRC/MEDUSA/bio_medusa_init.F90

@@ -7 +7 @@
    !!   -   ! 2017-04 (M. Stringer)        Code taken from trcbio_medusa.F90
    !!   -   ! 2017-08 (A. Yool)            Add slow-sinking detrius variables
+   !!   -   ! 2018-08 (A. Yool)            add OMIP preindustrial DIC
    !!----------------------------------------------------------------------
 #if defined key_medusa
@@ -186 +187 @@
       f_omarg(:,:)    = 0.0
       f_omcal(:,:)    = 0.0
+# if defined key_omip_dic
+      zomd(:,:)          = 0.0
+      f_pi_xco2a(:,:)    = 0.0
+      f_pi_pco2w(:,:)    = 0.0
+      f_pi_ph(:,:)       = 0.0
+      f_pi_co2flux(:,:)  = 0.0 
+      f_pi_pco2atm(:,:)  = 0.0
+      f_pi_h2co3(:,:)    = 0.0
+      f_pi_hco3(:,:)     = 0.0
+      f_pi_co3(:,:)      = 0.0
+      f_pi_omarg(:,:)    = 0.0
+      f_pi_omcal(:,:)    = 0.0
+# endif      
       !!
       !! AXY (08/08/17): zero slow detritus fluxes
@@ -382 +396 @@
             f_pco2a2d(:,:)      = 0.0 !!
          ENDIF
+#  if defined key_omip_dic         
+         IF( med_diag%PI_ATM_PCO2%dgsave ) THEN
+            ALLOCATE( f_pi_pco2a2d(1:jpi, 1:jpj) )
+            f_pi_pco2a2d(:,:)      = 0.0 !!
+         ENDIF
+#  endif
 !!
 !! skip OCN_PH diagnostic here
@@ -389 +409 @@
             f_pco2w2d(:,:)      = 0.0 !!
          ENDIF
+#  if defined key_omip_dic         
+         IF( med_diag%PI_OCN_PCO2%dgsave ) THEN
+            ALLOCATE( f_pi_pco2w2d(1:jpi, 1:jpj) )
+            f_pi_pco2w2d(:,:)   = 0.0 !!
+         ENDIF
+#  endif
 !!
 !! skip OCNH2CO3, OCN_HCO3, OCN_CO3 diagnostics here
@@ -396 +422 @@
             f_co2flux2d(:,:)      = 0.0 !!
          ENDIF
+#  if defined key_omip_dic         
+         IF( med_diag%PI_CO2FLUX%dgsave ) THEN
+            ALLOCATE( f_pi_co2flux2d(1:jpi, 1:jpj) )
+            f_pi_co2flux2d(:,:)   = 0.0 !!
+         ENDIF
+#  endif
 !!
 !! skip OM_XXX diagnostics here
@@ -403 +435 @@
             f_TDIC2d(:,:)      = 0.0 !!
          ENDIF
+#  if defined key_omip_dic         
+         IF( med_diag%PI_TCO2%dgsave ) THEN
+            ALLOCATE( f_pi_TDIC2d(1:jpi, 1:jpj) )
+            f_pi_TDIC2d(:,:)      = 0.0 !!
+         ENDIF
+#  endif
          IF( med_diag%TALK%dgsave ) THEN
             ALLOCATE( f_TALK2d(1:jpi, 1:jpj) )
@@ -662 +700 @@
             f_ocndpco2_2d(:,:)      = 0.0 !!
          ENDIF
+#  if defined key_omip_dic         
+         IF( med_diag%PI_ATM_XCO2%dgsave ) THEN
+            ALLOCATE( f_pi_xco2a_2d(1:jpi, 1:jpj) )
+            f_pi_xco2a_2d(:,:)      = 0.0 !!
+         ENDIF
+         IF( med_diag%PI_OCN_FCO2%dgsave ) THEN
+            ALLOCATE( f_pi_fco2w_2d(1:jpi, 1:jpj) )
+            f_pi_fco2w_2d(:,:)      = 0.0 !!
+         ENDIF
+         IF( med_diag%PI_ATM_FCO2%dgsave ) THEN
+            ALLOCATE( f_pi_fco2a_2d(1:jpi, 1:jpj) )
+            f_pi_fco2a_2d(:,:)      = 0.0 !!
+         ENDIF
+         IF( med_diag%PI_CO2STARAIR%dgsave ) THEN
+            ALLOCATE( f_pi_co2starair_2d(1:jpi, 1:jpj) )
+            f_pi_co2starair_2d(:,:)      = 0.0 !!
+         ENDIF
+         IF( med_diag%OCN_PI_DPCO2%dgsave ) THEN
+            ALLOCATE( f_pi_ocndpco2_2d(1:jpi, 1:jpj) )
+            f_pi_ocndpco2_2d(:,:)      = 0.0 !!
+         ENDIF
+         IF( med_diag%PI_FGCO2%dgsave ) THEN
+            ALLOCATE( pi_fgco2(1:jpi, 1:jpj) )
+            pi_fgco2(:,:)      = 0.0 !!
+         ENDIF
+#  endif
 # endif  
          IF( med_diag%TPP3%dgsave ) THEN