Changeset 8489 for branches/NERC/dev_r5518_GO6_MEDUSA_conserv/NEMOGCM/NEMO/TOP_SRC/MEDUSA/bio_medusa_mod.F90

Timestamp:

2017-09-01T16:55:53+02:00 (7 years ago)

Author:

jpalmier

Message:

JPALM -- gmed ticket #346 : improve MEDUSA conservation -- import BBL bug fix from NEMO ticket #1932 to GO6 branch

File:

• branches/NERC/dev_r5518_GO6_MEDUSA_conserv/NEMOGCM/NEMO/TOP_SRC/MEDUSA/bio_medusa_mod.F90 (modified) (5 diffs)

branches/NERC/dev_r5518_GO6_MEDUSA_conserv/NEMOGCM/NEMO/TOP_SRC/MEDUSA/bio_medusa_mod.F90

@@ -7 +7 @@
    !! History :
    !!   -   ! 2017-04 (M. Stringer)        Code taken from trcbio_medusa.F90
+   !!   -   ! 2017-08 (A. Yool)            Slow detritus, ML-avg chl variables
    !!----------------------------------------------------------------------
 #if defined key_medusa
@@ -54 +55 @@
    !! AXY (01/03/10): add in mixed layer PP diagnostics
    REAL(wp), ALLOCATABLE, DIMENSION(:,:) :: fprn_ml,fprd_ml
+   !! AXY (16/08/17): add in mixed layer chlorophyll diagnostic
+   REAL(wp), ALLOCATABLE, DIMENSION(:,:) :: fchl_ml
    !!
    !! nutrient limiting factors
@@ -94 +97 @@
    REAL(wp), ALLOCATABLE, DIMENSION(:,:) :: fregenfastc
 # endif
-
+   !!
+   !! AXY (08/08/17): sinking of detritus moved here
+   REAL(wp), ALLOCATABLE, DIMENSION(:,:) ::    fslowsink, fslowgain, fslowloss
+# if defined key_roam
+   REAL(wp), ALLOCATABLE, DIMENSION(:,:) ::    fslowsinkc, fslowgainc, fslowlossc
+# endif
+   !!
    !! Particle flux
    REAL(wp), ALLOCATABLE, DIMENSION(:,:) :: fdep1
@@ -287 +296 @@
                fjlim_pn(jpi,jpj),fjlim_pd(jpi,jpj),                   &
                fun_T(jpi,jpj),fun_Q10(jpi,jpj),                       &
-               fprn_ml(jpi,jpj),fprd_ml(jpi,jpj),                     &
+               fprn_ml(jpi,jpj),fprd_ml(jpi,jpj),fchl_ml(jpi,jpj),    &
                fnln(jpi,jpj),ffln2(jpi,jpj),                          &
                fnld(jpi,jpj),ffld(jpi,jpj),fsld(jpi,jpj),             &
@@ -316 +325 @@
                fregenfastc(jpi,jpj),                                  &
 # endif
+          fslowsink(jpi,jpj),fslowgain(jpi,jpj),                 &
+               fslowloss(jpi,jpj),                                    &
+# if defined key_roam
+          fslowsinkc(jpi,jpj),fslowgainc(jpi,jpj),               &
+               fslowlossc(jpi,jpj),                                   &
+# endif
                fdep1(jpi,jpj),                                        &
                ftempn(jpi,jpj),ftempsi(jpi,jpj),ftempfe(jpi,jpj),     &