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

Timestamp:

2018-08-08T17:12:35+02:00 (6 years ago)

Author:

jpalmier

Message:

merge with GO6_package_branch 9385-10020 ; plus debug OMIP_DIC

File:

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

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

@@ -38 +38 @@
       USE sms_medusa,        ONLY: jdms
       USE trc,               ONLY: ln_diatrc, med_diag, nittrc000 
-      USE in_out_manager,    ONLY: lwp
-
-# if defined key_iomput
-      USE iom,               ONLY: lk_iomput, numout
+      USE in_out_manager,    ONLY: lwp, numout
+
+      USE iom,               ONLY: lk_iomput
       USE trcnam_medusa,     ONLY: trc_nam_iom_medusa
-# endif
 
       !! time (integer timestep)
@@ -179 +177 @@
       zdic(:,:)  = 0.0
       zsil(:,:)  = 0.0
+# if defined key_mocsy
+      ! zpho is only defined if key_mocsy
+      ! is active, so we must protect this
+      ! initialisation accordingly. 
       zpho(:,:)  = 0.0
+# endif
       f_co2flux(:,:)  = 0.0 
       f_pco2atm(:,:)  = 0.0
@@ -398 +401 @@
 #  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 !!
+            ALLOCATE( f_pi_pco2a_2d(1:jpi, 1:jpj) )
+            f_pi_pco2a_2d(:,:)      = 0.0 !!
          ENDIF
 #  endif
@@ -405 +408 @@
 !! skip OCN_PH diagnostic here
 !!
+#  if defined key_omip_dic         
+         IF( med_diag%PI_OCN_PH%dgsave ) THEN
+            ALLOCATE( f_pi_ph_2d(1:jpi, 1:jpj) )
+            f_pi_ph_2d(:,:)   = 0.0 !!
+         ENDIF
+#  endif
          IF( med_diag%OCN_PCO2%dgsave ) THEN
             ALLOCATE( f_pco2w2d(1:jpi, 1:jpj) )
@@ -411 +420 @@
 #  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 !!
+            ALLOCATE( f_pi_pco2w_2d(1:jpi, 1:jpj) )
+            f_pi_pco2w_2d(:,:)   = 0.0 !!
          ENDIF
 #  endif
@@ -418 +427 @@
 !! skip OCNH2CO3, OCN_HCO3, OCN_CO3 diagnostics here
 !!
+#  if defined key_omip_dic         
+         IF( med_diag%PI_OCNH2CO3%dgsave ) THEN
+            ALLOCATE( f_pi_h2co3_2d(1:jpi, 1:jpj) )
+            f_pi_h2co3_2d(:,:)   = 0.0 !!
+         ENDIF
+         IF( med_diag%PI_OCN_HCO3%dgsave ) THEN
+            ALLOCATE( f_pi_hco3_2d(1:jpi, 1:jpj) )
+            f_pi_hco3_2d(:,:)   = 0.0 !!
+         ENDIF
+         IF( med_diag%PI_OCN_CO3%dgsave ) THEN
+            ALLOCATE( f_pi_co3_2d(1:jpi, 1:jpj) )
+            f_pi_co3_2d(:,:)   = 0.0 !!
+         ENDIF
+#  endif
+
          IF( med_diag%CO2FLUX%dgsave ) THEN
             ALLOCATE( f_co2flux2d(1:jpi, 1:jpj) )
@@ -431 +455 @@
 !! skip OM_XXX diagnostics here
 !!
+#  if defined key_omip_dic         
+         IF( med_diag%PI_OM_CAL%dgsave ) THEN
+            ALLOCATE( f_pi_omcal_2d(1:jpi, 1:jpj) )
+            f_pi_omcal_2d(:,:)   = 0.0 !!
+         ENDIF
+         IF( med_diag%PI_OM_ARG%dgsave ) THEN
+            ALLOCATE( f_pi_omarg_2d(1:jpi, 1:jpj) )
+            f_pi_omarg_2d(:,:)   = 0.0 !!
+         ENDIF
+#  endif
          IF( med_diag%TCO2%dgsave ) THEN
             ALLOCATE( f_TDIC2d(1:jpi, 1:jpj) )
@@ -928 +962 @@
 CONTAINS
    SUBROUTINE bio_medusa_init( )                   ! Empty routine
+      IMPLICIT NONE
       WRITE(*,*) 'bio_medusa_init: You should not have seen this print! error?'
    END SUBROUTINE bio_medusa_init