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

Timestamp:

2017-05-18T11:14:31+02:00 (7 years ago)

Author:

marc

Message:

Removed about 40 2d arrays from bio_medusa_mod.F90

File:

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

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

@@ -41 +41 @@
    REAL(wp) ::    b0
 
-   REAL(wp), ALLOCATABLE, DIMENSION(:,:) :: fthetan,faln,fchn1,fchn
-   REAL(wp), ALLOCATABLE, DIMENSION(:,:) :: fjln,fprn,frn
-   REAL(wp), ALLOCATABLE, DIMENSION(:,:) :: fthetad,fald,fchd1,fchd
-   REAL(wp), ALLOCATABLE, DIMENSION(:,:) :: fjld,fprd,frd
-
-   REAL(wp), ALLOCATABLE, DIMENSION(:,:) :: fjlim_pn, fjlim_pd
-   !! AXY (03/02/11): add in Liebig terms
-   REAL(wp), ALLOCATABLE, DIMENSION(:,:) :: fpnlim, fpdlim
+   REAL(wp), ALLOCATABLE, DIMENSION(:,:) :: fthetan,fprn,frn
+   REAL(wp), ALLOCATABLE, DIMENSION(:,:) :: fthetad,fprd,frd
+
+   REAL(wp), ALLOCATABLE, DIMENSION(:,:) :: fjlim_pn,fjlim_pd
    !! AXY (16/07/09): add in Eppley curve functionality
    REAL(wp), ALLOCATABLE, DIMENSION(:,:) :: fun_T,xvpnT,xvpdT
@@ -67 +63 @@
    !!
    !! silicon cycle
-   REAL(wp), ALLOCATABLE, DIMENSION(:,:) :: fsin,fnsi,fprds,fsdiss
+   REAL(wp), ALLOCATABLE, DIMENSION(:,:) :: fsin,fprds,fsdiss
 
    !! iron cycle; includes parameters for Parekh et al. (2005) iron scheme
@@ -75 +71 @@
 
    !! Microzooplankton grazing
-   REAL(wp), ALLOCATABLE, DIMENSION(:,:) :: fmi1,fmi,fgmipn,fgmid,fgmidc
-   REAL(wp), ALLOCATABLE, DIMENSION(:,:) :: finmi,ficmi,fstarmi,fmith
+   REAL(wp), ALLOCATABLE, DIMENSION(:,:) :: fgmipn,fgmid,fgmidc
+   REAL(wp), ALLOCATABLE, DIMENSION(:,:) :: finmi,ficmi
    REAL(wp), ALLOCATABLE, DIMENSION(:,:) :: fmigrow,fmiexcr,fmiresp
    !!
    !! Mesozooplankton grazing
-   REAL(wp), ALLOCATABLE, DIMENSION(:,:) :: fme1,fme,fgmepn,fgmepd
+   REAL(wp), ALLOCATABLE, DIMENSION(:,:) :: fgmepn,fgmepd
    REAL(wp), ALLOCATABLE, DIMENSION(:,:) :: fgmepds,fgmezmi,fgmed,fgmedc
-   REAL(wp), ALLOCATABLE, DIMENSION(:,:) :: finme,ficme,fstarme,fmeth
+   REAL(wp), ALLOCATABLE, DIMENSION(:,:) :: finme,ficme
    REAL(wp), ALLOCATABLE, DIMENSION(:,:) :: fmegrow,fmeexcr,fmeresp
    !!
@@ -100 +96 @@
 
    !! Particle flux
-   REAL(wp), ALLOCATABLE, DIMENSION(:,:) :: fdep1,fcaco3
+   REAL(wp), ALLOCATABLE, DIMENSION(:,:) :: fdep1
    REAL(wp), ALLOCATABLE, DIMENSION(:,:) :: ftempn,ftempsi,ftempfe
    REAL(wp), ALLOCATABLE, DIMENSION(:,:) :: ftempc,ftempca
@@ -107 +103 @@
    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
 
    !! AXY (08/07/11): fate of fast detritus reaching the seafloor
@@ -157 +151 @@
 
    !! AXY (23/06/15): additional diagnostics for MOCSY and oxygen
-   REAL(wp), ALLOCATABLE, DIMENSION(:,:) :: f_fco2w,f_BetaD,f_rhosw,f_opres
-   REAL(wp), ALLOCATABLE, DIMENSION(:,:) :: f_insitut,f_pco2atm,f_fco2atm
-   REAL(wp), ALLOCATABLE, DIMENSION(:,:) :: f_schmidtco2,f_kwco2,f_K0
-   REAL(wp), ALLOCATABLE, DIMENSION(:,:) :: f_co2starair,f_dpco2,f_kwo2
-
-   INTEGER, ALLOCATABLE, DIMENSION(:,:)  :: iters 
+   REAL(wp), ALLOCATABLE, DIMENSION(:,:) :: f_pco2atm
 
    !! Carbon, alkalinity production and consumption
@@ -292 +281 @@
                zpho(jpi,jpj),                                         &
 # endif
-               fthetan(jpi,jpj),faln(jpi,jpj),fchn1(jpi,jpj),         &
-               fchn(jpi,jpj),fjln(jpi,jpj),fprn(jpi,jpj),             &
-               frn(jpi,jpj),                                          &
-               fthetad(jpi,jpj),fald(jpi,jpj),fchd1(jpi,jpj),         &
-               fchd(jpi,jpj),fjld(jpi,jpj),fprd(jpi,jpj),             &
-               frd(jpi,jpj),                                          &
-               fjlim_pn(jpi,jpj), fjlim_pd(jpi,jpj),                  &
-               fpnlim(jpi,jpj), fpdlim(jpi,jpj),                      &
-               fun_T(jpi,jpj),xvpnT(jpi,jpj),xvpdT(jpi,jpj),          &
-               fun_Q10(jpi,jpj),                                      &
+               fthetan(jpi,jpj),fprn(jpi,jpj),frn(jpi,jpj),           &
+               fthetad(jpi,jpj),fprd(jpi,jpj),frd(jpi,jpj),           &
+               fjlim_pn(jpi,jpj),fjlim_pd(jpi,jpj),                   &
+               fun_T(jpi,jpj),fun_Q10(jpi,jpj),                       &
                fprn_ml(jpi,jpj),fprd_ml(jpi,jpj),                     &
                fnln(jpi,jpj),ffln2(jpi,jpj),                          &
                fnld(jpi,jpj),ffld(jpi,jpj),fsld(jpi,jpj),             &
                fsld2(jpi,jpj),                                        &
-               fsin(jpi,jpj),fnsi(jpi,jpj),fprds(jpi,jpj),            &
-               fsdiss(jpi,jpj),                                       &
+               fsin(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),            &
-               fgmid(jpi,jpj),fgmidc(jpi,jpj),                        &
-               finmi(jpi,jpj),ficmi(jpi,jpj),fstarmi(jpi,jpj),        &
-               fmith(jpi,jpj),fmigrow(jpi,jpj),fmiexcr(jpi,jpj),      &
-               fmiresp(jpi,jpj),                                      &
-               fme1(jpi,jpj),fme(jpi,jpj),fgmepn(jpi,jpj),            &
-               fgmepd(jpi,jpj),fgmepds(jpi,jpj),fgmezmi(jpi,jpj),     &
-               fgmed(jpi,jpj),fgmedc(jpi,jpj),                        &
-               finme(jpi,jpj),ficme(jpi,jpj),fstarme(jpi,jpj),        &
-               fmeth(jpi,jpj),fmegrow(jpi,jpj),fmeexcr(jpi,jpj),      &
-               fmeresp(jpi,jpj),                                      &
+               fgmipn(jpi,jpj),fgmid(jpi,jpj),fgmidc(jpi,jpj),        &
+               finmi(jpi,jpj),ficmi(jpi,jpj),                         &
+               fmigrow(jpi,jpj),fmiexcr(jpi,jpj),fmiresp(jpi,jpj),    &
+               fgmepn(jpi,jpj),fgmepd(jpi,jpj),                       &
+               fgmepds(jpi,jpj),fgmezmi(jpi,jpj),fgmed(jpi,jpj),      &
+               fgmedc(jpi,jpj),                                       &
+               finme(jpi,jpj),ficme(jpi,jpj),                         &
+               fmegrow(jpi,jpj),fmeexcr(jpi,jpj),fmeresp(jpi,jpj),    &
                fdpn(jpi,jpj),fdpd(jpi,jpj),fdpds(jpi,jpj),            &
                fdzmi(jpi,jpj),fdzme(jpi,jpj),fdd(jpi,jpj),            &
@@ -335 +314 @@
                fregenfastc(jpi,jpj),                                  &
 # endif
-               fdep1(jpi,jpj),fcaco3(jpi,jpj),                        &
+               fdep1(jpi,jpj),                                        &
                ftempn(jpi,jpj),ftempsi(jpi,jpj),ftempfe(jpi,jpj),     &
                ftempc(jpi,jpj),ftempca(jpi,jpj),                      &
@@ -342 +321 @@
                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),                                     &
                ffast2slown(jpi,jpj),ffast2slowc(jpi,jpj),             &
                ftot_n(jpi,jpj),ftot_si(jpi,jpj),ftot_fe(jpi,jpj),     &
@@ -377 +354 @@
                f_kw660(jpi,jpj),f_o2flux(jpi,jpj),f_o2sat(jpi,jpj),   &
                f_omcal(jpi,jpj),f_omarg(jpi,jpj),                     &
-               f_fco2w(jpi,jpj),f_BetaD(jpi,jpj),f_rhosw(jpi,jpj),    &
-               f_opres(jpi,jpj),f_insitut(jpi,jpj),                   &
-               f_pco2atm(jpi,jpj),f_fco2atm(jpi,jpj),                 &
-               f_schmidtco2(jpi,jpj),f_kwco2(jpi,jpj),f_K0(jpi,jpj),  &
-               f_co2starair(jpi,jpj),f_dpco2(jpi,jpj),                &
-               f_kwo2(jpi,jpj),                                       &
-               iters(jpi,jpj),                                        &
+               f_pco2atm(jpi,jpj),                 &
                fcomm_resp(jpi,jpj),                                   &
                fcar_prod(jpi,jpj),fcar_cons(jpi,jpj),                 &