Changeset 7978 for branches/UKMO/dev_r5518_medusa_chg_trc_bio_medusa/NEMOGCM/NEMO/TOP_SRC/MEDUSA/bio_medusa_mod.F90

Timestamp:

2017-04-27T13:13:32+02:00 (7 years ago)

Author:

marc

Message:

Moved the iron chemistry and scavenging out of trcbio_medusa.F90

File:

• branches/UKMO/dev_r5518_medusa_chg_trc_bio_medusa/NEMOGCM/NEMO/TOP_SRC/MEDUSA/bio_medusa_mod.F90 (modified) (6 diffs)

branches/UKMO/dev_r5518_medusa_chg_trc_bio_medusa/NEMOGCM/NEMO/TOP_SRC/MEDUSA/bio_medusa_mod.F90

@@ -58 +58 @@
    REAL(wp), ALLOCATABLE, DIMENSION(:,:) :: fsin,fnsi,fprds,fsdiss
 
+   !! iron cycle; includes parameters for Parekh et al. (2005) iron scheme
+   REAL(wp), ALLOCATABLE, DIMENSION(:,:) :: ffetop,ffebot,ffescav
    !! Variable for iron-ligand system
    REAL(wp), ALLOCATABLE, DIMENSION(:,:) :: xFree
@@ -74 +76 @@
    !! mortality/Remineralisation (defunct parameter "fz" removed)
    REAL(wp), ALLOCATABLE, DIMENSION(:,:) :: fdpn,fdpd,fdpds,fdzmi,fdzme,fdd
+# if defined key_roam
+   REAL(wp), ALLOCATABLE, DIMENSION(:,:) :: fddc
+# endif
    REAL(wp), ALLOCATABLE, DIMENSION(:,:) :: fdpn2,fdpd2,fdpds2,fdzmi2,fdzme2
 
@@ -85 +90 @@
    !! Particle flux
    REAL(wp), ALLOCATABLE, DIMENSION(:,:) :: fdep1,fcaco3
-
+   REAL(wp), ALLOCATABLE, DIMENSION(:,:) :: ftempn,ftempsi,ftempfe
+   REAL(wp), ALLOCATABLE, DIMENSION(:,:) :: ftempc,ftempca
+   REAL(wp), ALLOCATABLE, DIMENSION(:,:) :: freminn,freminsi,freminfe
+   REAL(wp), ALLOCATABLE, DIMENSION(:,:) :: freminc,freminca
    REAL(wp), ALLOCATABLE, DIMENSION(:,:) :: ffastn,ffastsi,ffastfe
    REAL(wp), ALLOCATABLE, DIMENSION(:,:) :: ffastc,ffastca
+   REAL(wp), ALLOCATABLE, DIMENSION(:,:) :: fprotf
    REAL(wp), ALLOCATABLE, DIMENSION(:,:) :: fsedn,fsedsi,fsedfe,fsedc,fsedca
    REAL(wp), ALLOCATABLE, DIMENSION(:,:) :: fccd
+   REAL(wp), ALLOCATABLE, DIMENSION(:,:) :: fccd_dep
 
    !! water column nutrient and flux integrals
@@ -269 +279 @@
                fsin(jpi,jpj),fnsi(jpi,jpj),fprds(jpi,jpj),            &
                fsdiss(jpi,jpj),                                       &
+               ffetop(jpi,jpj),ffebot(jpi,jpj),ffescav(jpi,jpj),      &
                xFree(jpi,jpj),                                        &
                fmi1(jpi,jpj),fmi(jpi,jpj),fgmipn(jpi,jpj),            &
@@ -283 +294 @@
                fdpn(jpi,jpj),fdpd(jpi,jpj),fdpds(jpi,jpj),            &
                fdzmi(jpi,jpj),fdzme(jpi,jpj),fdd(jpi,jpj),            &
+               fddc(jpi,jpj),                                         &
                fdpn2(jpi,jpj),fdpd2(jpi,jpj),fdpds2(jpi,jpj),         &
                fdzmi2(jpi,jpj),fdzme2(jpi,jpj),                       &
@@ -291 +303 @@
 # endif
                fdep1(jpi,jpj),fcaco3(jpi,jpj),                        &
+               ftempn(jpi,jpj),ftempsi(jpi,jpj),ftempfe(jpi,jpj),     &
+               ftempc(jpi,jpj),ftempca(jpi,jpj),                      &
+               freminn(jpi,jpj),freminsi(jpi,jpj),freminfe(jpi,jpj),  &
+               freminc(jpi,jpj),freminca(jpi,jpj),                    &
                ffastn(jpi,jpj),ffastsi(jpi,jpj),ffastfe(jpi,jpj),     &
                ffastc(jpi,jpj),ffastca(jpi,jpj),                      &
+               fprotf(jpi,jpj),                                       &
                fsedn(jpi,jpj),fsedsi(jpi,jpj),fsedfe(jpi,jpj),        &
                fsedc(jpi,jpj),fsedca(jpi,jpj),                        &
                fccd(jpi,jpj),                                         &
+               fccd_dep(jpi,jpj),                                     &
                ftot_n(jpi,jpj),ftot_si(jpi,jpj),ftot_fe(jpi,jpj),     &
                fflx_n(jpi,jpj),fflx_si(jpi,jpj),fflx_fe(jpi,jpj),     &