Changeset 7716


Ignore:
Timestamp:
2017-02-22T16:11:16+01:00 (4 years ago)
Author:
dford
Message:

Tidy up and simplify use of arrays for logchl assimilation balancing.

Location:
branches/UKMO/dev_r5518_v3.4_asm_nemovar_community_bgc_ersem/NEMOGCM/NEMO/OPA_SRC/ASM
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • branches/UKMO/dev_r5518_v3.4_asm_nemovar_community_bgc_ersem/NEMOGCM/NEMO/OPA_SRC/ASM/asmbal.F90

    r6983 r7716  
    7070         IF ( ln_logchltotinc .OR. ln_logchlpftinc ) THEN 
    7171#if defined key_fabm 
    72             CALL iom_rstput( kt, kt, inum, 'logchl_balinc_ersem_chl1', logchl_balinc_ersem_chl1 ) 
    73             CALL iom_rstput( kt, kt, inum, 'logchl_balinc_ersem_chl2', logchl_balinc_ersem_chl2 ) 
    74             CALL iom_rstput( kt, kt, inum, 'logchl_balinc_ersem_chl3', logchl_balinc_ersem_chl3 ) 
    75             CALL iom_rstput( kt, kt, inum, 'logchl_balinc_ersem_chl4', logchl_balinc_ersem_chl4 ) 
    76             CALL iom_rstput( kt, kt, inum, 'logchl_balinc_ersem_p1c',  logchl_balinc_ersem_p1c  ) 
    77             CALL iom_rstput( kt, kt, inum, 'logchl_balinc_ersem_p1n',  logchl_balinc_ersem_p1n  ) 
    78             CALL iom_rstput( kt, kt, inum, 'logchl_balinc_ersem_p1p',  logchl_balinc_ersem_p1p  ) 
    79             CALL iom_rstput( kt, kt, inum, 'logchl_balinc_ersem_p1s',  logchl_balinc_ersem_p1s  ) 
    80             CALL iom_rstput( kt, kt, inum, 'logchl_balinc_ersem_p2c',  logchl_balinc_ersem_p2c  ) 
    81             CALL iom_rstput( kt, kt, inum, 'logchl_balinc_ersem_p2n',  logchl_balinc_ersem_p2n  ) 
    82             CALL iom_rstput( kt, kt, inum, 'logchl_balinc_ersem_p2p',  logchl_balinc_ersem_p2p  ) 
    83             CALL iom_rstput( kt, kt, inum, 'logchl_balinc_ersem_p3c',  logchl_balinc_ersem_p3c  ) 
    84             CALL iom_rstput( kt, kt, inum, 'logchl_balinc_ersem_p3n',  logchl_balinc_ersem_p3n  ) 
    85             CALL iom_rstput( kt, kt, inum, 'logchl_balinc_ersem_p3p',  logchl_balinc_ersem_p3p  ) 
    86             CALL iom_rstput( kt, kt, inum, 'logchl_balinc_ersem_p4c',  logchl_balinc_ersem_p4c  ) 
    87             CALL iom_rstput( kt, kt, inum, 'logchl_balinc_ersem_p4n',  logchl_balinc_ersem_p4n  ) 
    88             CALL iom_rstput( kt, kt, inum, 'logchl_balinc_ersem_p4p',  logchl_balinc_ersem_p4p  ) 
    89             CALL iom_rstput( kt, kt, inum, 'logchl_balinc_ersem_z4c',  logchl_balinc_ersem_z4c  ) 
    90             CALL iom_rstput( kt, kt, inum, 'logchl_balinc_ersem_z5c',  logchl_balinc_ersem_z5c  ) 
    91             CALL iom_rstput( kt, kt, inum, 'logchl_balinc_ersem_z5n',  logchl_balinc_ersem_z5n  ) 
    92             CALL iom_rstput( kt, kt, inum, 'logchl_balinc_ersem_z5p',  logchl_balinc_ersem_z5p  ) 
    93             CALL iom_rstput( kt, kt, inum, 'logchl_balinc_ersem_z6c',  logchl_balinc_ersem_z6c  ) 
    94             CALL iom_rstput( kt, kt, inum, 'logchl_balinc_ersem_z6n',  logchl_balinc_ersem_z6n  ) 
    95             CALL iom_rstput( kt, kt, inum, 'logchl_balinc_ersem_z6p',  logchl_balinc_ersem_z6p  ) 
    96             CALL iom_rstput( kt, kt, inum, 'logchl_balinc_ersem_n1p',  logchl_balinc_ersem_n1p  ) 
    97             CALL iom_rstput( kt, kt, inum, 'logchl_balinc_ersem_n3n',  logchl_balinc_ersem_n3n  ) 
    98             CALL iom_rstput( kt, kt, inum, 'logchl_balinc_ersem_n4n',  logchl_balinc_ersem_n4n  ) 
    99             CALL iom_rstput( kt, kt, inum, 'logchl_balinc_ersem_n5s',  logchl_balinc_ersem_n5s  ) 
     72            CALL iom_rstput( kt, kt, inum, 'logchl_balinc_chl1', logchl_balinc(:,:,:,jp_fabm_m1+jp_fabm_chl1) ) 
     73            CALL iom_rstput( kt, kt, inum, 'logchl_balinc_chl2', logchl_balinc(:,:,:,jp_fabm_m1+jp_fabm_chl2) ) 
     74            CALL iom_rstput( kt, kt, inum, 'logchl_balinc_chl3', logchl_balinc(:,:,:,jp_fabm_m1+jp_fabm_chl3) ) 
     75            CALL iom_rstput( kt, kt, inum, 'logchl_balinc_chl4', logchl_balinc(:,:,:,jp_fabm_m1+jp_fabm_chl4) ) 
    10076#elif defined key_medusa && defined key_foam_medusa 
    101             CALL iom_rstput( kt, kt, inum, 'logchl_balinc_medusa_chn', logchl_balinc_medusa_chn ) 
    102             CALL iom_rstput( kt, kt, inum, 'logchl_balinc_medusa_chd', logchl_balinc_medusa_chd ) 
    103             CALL iom_rstput( kt, kt, inum, 'logchl_balinc_medusa_phn', logchl_balinc_medusa_phn ) 
    104             CALL iom_rstput( kt, kt, inum, 'logchl_balinc_medusa_phd', logchl_balinc_medusa_phd ) 
    105             CALL iom_rstput( kt, kt, inum, 'logchl_balinc_medusa_pds', logchl_balinc_medusa_pds ) 
    106             CALL iom_rstput( kt, kt, inum, 'logchl_balinc_medusa_zmi', logchl_balinc_medusa_zmi ) 
    107             CALL iom_rstput( kt, kt, inum, 'logchl_balinc_medusa_zme', logchl_balinc_medusa_zme ) 
    108             CALL iom_rstput( kt, kt, inum, 'logchl_balinc_medusa_din', logchl_balinc_medusa_din ) 
    109             CALL iom_rstput( kt, kt, inum, 'logchl_balinc_medusa_sil', logchl_balinc_medusa_sil ) 
    110             CALL iom_rstput( kt, kt, inum, 'logchl_balinc_medusa_fer', logchl_balinc_medusa_fer ) 
    111             CALL iom_rstput( kt, kt, inum, 'logchl_balinc_medusa_det', logchl_balinc_medusa_det ) 
    112             CALL iom_rstput( kt, kt, inum, 'logchl_balinc_medusa_dtc', logchl_balinc_medusa_dtc ) 
    113             CALL iom_rstput( kt, kt, inum, 'logchl_balinc_medusa_dic', logchl_balinc_medusa_dic ) 
    114             CALL iom_rstput( kt, kt, inum, 'logchl_balinc_medusa_alk', logchl_balinc_medusa_alk ) 
    115             CALL iom_rstput( kt, kt, inum, 'logchl_balinc_medusa_oxy', logchl_balinc_medusa_oxy ) 
     77            CALL iom_rstput( kt, kt, inum, 'logchl_balinc_chn', logchl_balinc(:,:,:,jpchn) ) 
     78            CALL iom_rstput( kt, kt, inum, 'logchl_balinc_chd', logchl_balinc(:,:,:,jpchd) ) 
     79            CALL iom_rstput( kt, kt, inum, 'logchl_balinc_phn', logchl_balinc(:,:,:,jpphn) ) 
     80            CALL iom_rstput( kt, kt, inum, 'logchl_balinc_phd', logchl_balinc(:,:,:,jpphd) ) 
     81            CALL iom_rstput( kt, kt, inum, 'logchl_balinc_pds', logchl_balinc(:,:,:,jppds) ) 
     82            CALL iom_rstput( kt, kt, inum, 'logchl_balinc_zmi', logchl_balinc(:,:,:,jpzmi) ) 
     83            CALL iom_rstput( kt, kt, inum, 'logchl_balinc_zme', logchl_balinc(:,:,:,jpzme) ) 
     84            CALL iom_rstput( kt, kt, inum, 'logchl_balinc_din', logchl_balinc(:,:,:,jpdin) ) 
     85            CALL iom_rstput( kt, kt, inum, 'logchl_balinc_sil', logchl_balinc(:,:,:,jpsil) ) 
     86            CALL iom_rstput( kt, kt, inum, 'logchl_balinc_fer', logchl_balinc(:,:,:,jpfer) ) 
     87            CALL iom_rstput( kt, kt, inum, 'logchl_balinc_det', logchl_balinc(:,:,:,jpdet) ) 
     88            CALL iom_rstput( kt, kt, inum, 'logchl_balinc_dtc', logchl_balinc(:,:,:,jpdtc) ) 
     89            CALL iom_rstput( kt, kt, inum, 'logchl_balinc_dic', logchl_balinc(:,:,:,jpdic) ) 
     90            CALL iom_rstput( kt, kt, inum, 'logchl_balinc_alk', logchl_balinc(:,:,:,jpalk) ) 
     91            CALL iom_rstput( kt, kt, inum, 'logchl_balinc_oxy', logchl_balinc(:,:,:,jpoxy) ) 
    11692#elif defined key_hadocc 
    117             CALL iom_rstput( kt, kt, inum, 'logchl_balinc_hadocc_nut', logchl_balinc_hadocc_nut ) 
    118             CALL iom_rstput( kt, kt, inum, 'logchl_balinc_hadocc_phy', logchl_balinc_hadocc_phy ) 
    119             CALL iom_rstput( kt, kt, inum, 'logchl_balinc_hadocc_zoo', logchl_balinc_hadocc_zoo ) 
    120             CALL iom_rstput( kt, kt, inum, 'logchl_balinc_hadocc_det', logchl_balinc_hadocc_det ) 
    121             CALL iom_rstput( kt, kt, inum, 'logchl_balinc_hadocc_dic', logchl_balinc_hadocc_dic ) 
    122             CALL iom_rstput( kt, kt, inum, 'logchl_balinc_hadocc_alk', logchl_balinc_hadocc_alk ) 
     93            CALL iom_rstput( kt, kt, inum, 'logchl_balinc_nut', logchl_balinc(:,:,:,jp_had_nut) ) 
     94            CALL iom_rstput( kt, kt, inum, 'logchl_balinc_phy', logchl_balinc(:,:,:,jp_had_phy) ) 
     95            CALL iom_rstput( kt, kt, inum, 'logchl_balinc_zoo', logchl_balinc(:,:,:,jp_had_zoo) ) 
     96            CALL iom_rstput( kt, kt, inum, 'logchl_balinc_det', logchl_balinc(:,:,:,jp_had_det) ) 
     97            CALL iom_rstput( kt, kt, inum, 'logchl_balinc_dic', logchl_balinc(:,:,:,jp_had_dic) ) 
     98            CALL iom_rstput( kt, kt, inum, 'logchl_balinc_alk', logchl_balinc(:,:,:,jp_had_alk) ) 
    12399#endif 
    124100         ENDIF 
  • branches/UKMO/dev_r5518_v3.4_asm_nemovar_community_bgc_ersem/NEMOGCM/NEMO/OPA_SRC/ASM/asminc.F90

    r7472 r7716  
    5353   USE bdy_oce, ONLY: bdytmask   
    5454#endif   
     55#if defined key_top 
     56   USE trc, ONLY: & 
     57      & trn,      & 
     58      & trb 
     59   USE par_trc, ONLY: & 
     60      & jptra 
     61#endif 
    5562#if defined key_fabm 
    5663   USE asmlogchlbal_ersem, ONLY: & 
    5764      & asm_logchl_bal_ersem 
    58    USE trc, ONLY: & 
    59       & trn,      & 
    60       & trb 
    6165   USE par_fabm 
    6266#elif defined key_medusa && defined key_foam_medusa 
    6367   USE asmlogchlbal_medusa, ONLY: & 
    6468      & asm_logchl_bal_medusa 
    65    USE trc, ONLY: & 
    66       & trn,      & 
    67       & trb 
    6869   USE par_medusa 
    6970#elif defined key_hadocc 
    7071   USE asmlogchlbal_hadocc, ONLY: & 
    7172      & asm_logchl_bal_hadocc 
    72    USE trc, ONLY: & 
    73       & trn,      & 
    74       & trb 
    7573   USE par_hadocc 
    7674#endif 
     
    130128 
    131129   REAL(wp), PUBLIC, DIMENSION(:,:,:), ALLOCATABLE :: logchl_bkginc !: Increment to background logchl 
    132 #if defined key_fabm 
    133    REAL(wp), PUBLIC, DIMENSION(:,:,:), ALLOCATABLE :: logchl_balinc_ersem_chl1  !: Increment to ERSEM diatom chl   from logchl 
    134    REAL(wp), PUBLIC, DIMENSION(:,:,:), ALLOCATABLE :: logchl_balinc_ersem_chl2  !: Increment to ERSEM nanophy chl  from logchl 
    135    REAL(wp), PUBLIC, DIMENSION(:,:,:), ALLOCATABLE :: logchl_balinc_ersem_chl3  !: Increment to ERSEM picophy chl  from logchl 
    136    REAL(wp), PUBLIC, DIMENSION(:,:,:), ALLOCATABLE :: logchl_balinc_ersem_chl4  !: Increment to ERSEM microphy chl from logchl 
    137    REAL(wp), PUBLIC, DIMENSION(:,:,:), ALLOCATABLE :: logchl_balinc_ersem_p1c   !: Increment to ERSEM diatom c     from logchl 
    138    REAL(wp), PUBLIC, DIMENSION(:,:,:), ALLOCATABLE :: logchl_balinc_ersem_p1n   !: Increment to ERSEM diatom n     from logchl 
    139    REAL(wp), PUBLIC, DIMENSION(:,:,:), ALLOCATABLE :: logchl_balinc_ersem_p1p   !: Increment to ERSEM diatom p     from logchl 
    140    REAL(wp), PUBLIC, DIMENSION(:,:,:), ALLOCATABLE :: logchl_balinc_ersem_p1s   !: Increment to ERSEM diatom s     from logchl 
    141    REAL(wp), PUBLIC, DIMENSION(:,:,:), ALLOCATABLE :: logchl_balinc_ersem_p2c   !: Increment to ERSEM nanophy c    from logchl 
    142    REAL(wp), PUBLIC, DIMENSION(:,:,:), ALLOCATABLE :: logchl_balinc_ersem_p2n   !: Increment to ERSEM nanophy n    from logchl 
    143    REAL(wp), PUBLIC, DIMENSION(:,:,:), ALLOCATABLE :: logchl_balinc_ersem_p2p   !: Increment to ERSEM nanophy p    from logchl 
    144    REAL(wp), PUBLIC, DIMENSION(:,:,:), ALLOCATABLE :: logchl_balinc_ersem_p3c   !: Increment to ERSEM picophy c    from logchl 
    145    REAL(wp), PUBLIC, DIMENSION(:,:,:), ALLOCATABLE :: logchl_balinc_ersem_p3n   !: Increment to ERSEM picophy n    from logchl 
    146    REAL(wp), PUBLIC, DIMENSION(:,:,:), ALLOCATABLE :: logchl_balinc_ersem_p3p   !: Increment to ERSEM picophy p    from logchl 
    147    REAL(wp), PUBLIC, DIMENSION(:,:,:), ALLOCATABLE :: logchl_balinc_ersem_p4c   !: Increment to ERSEM microphy c   from logchl 
    148    REAL(wp), PUBLIC, DIMENSION(:,:,:), ALLOCATABLE :: logchl_balinc_ersem_p4n   !: Increment to ERSEM microphy n   from logchl 
    149    REAL(wp), PUBLIC, DIMENSION(:,:,:), ALLOCATABLE :: logchl_balinc_ersem_p4p   !: Increment to ERSEM microphy p   from logchl 
    150    REAL(wp), PUBLIC, DIMENSION(:,:,:), ALLOCATABLE :: logchl_balinc_ersem_z4c   !: Increment to ERSEM mesozoo c    from logchl 
    151    REAL(wp), PUBLIC, DIMENSION(:,:,:), ALLOCATABLE :: logchl_balinc_ersem_z5c   !: Increment to ERSEM microzoo c   from logchl 
    152    REAL(wp), PUBLIC, DIMENSION(:,:,:), ALLOCATABLE :: logchl_balinc_ersem_z5n   !: Increment to ERSEM microzoo n   from logchl 
    153    REAL(wp), PUBLIC, DIMENSION(:,:,:), ALLOCATABLE :: logchl_balinc_ersem_z5p   !: Increment to ERSEM microzoo p   from logchl 
    154    REAL(wp), PUBLIC, DIMENSION(:,:,:), ALLOCATABLE :: logchl_balinc_ersem_z6c   !: Increment to ERSEM het flag c   from logchl 
    155    REAL(wp), PUBLIC, DIMENSION(:,:,:), ALLOCATABLE :: logchl_balinc_ersem_z6n   !: Increment to ERSEM het flag n   from logchl 
    156    REAL(wp), PUBLIC, DIMENSION(:,:,:), ALLOCATABLE :: logchl_balinc_ersem_z6p   !: Increment to ERSEM het flag p   from logchl 
    157    REAL(wp), PUBLIC, DIMENSION(:,:,:), ALLOCATABLE :: logchl_balinc_ersem_n1p   !: Increment to ERSEM phosphate    from logchl 
    158    REAL(wp), PUBLIC, DIMENSION(:,:,:), ALLOCATABLE :: logchl_balinc_ersem_n3n   !: Increment to ERSEM nitrate      from logchl 
    159    REAL(wp), PUBLIC, DIMENSION(:,:,:), ALLOCATABLE :: logchl_balinc_ersem_n4n   !: Increment to ERSEM ammonium     from logchl 
    160    REAL(wp), PUBLIC, DIMENSION(:,:,:), ALLOCATABLE :: logchl_balinc_ersem_n5s   !: Increment to ERSEM silicate     from logchl 
    161 #elif defined key_medusa && defined key_foam_medusa 
    162    REAL(wp), PUBLIC, DIMENSION(:,:,:), ALLOCATABLE :: logchl_balinc_medusa_chn  !: Increment to MEDUSA nondiatom chl from logchl 
    163    REAL(wp), PUBLIC, DIMENSION(:,:,:), ALLOCATABLE :: logchl_balinc_medusa_chd  !: Increment to MEDUSA diatom chl    from logchl 
    164    REAL(wp), PUBLIC, DIMENSION(:,:,:), ALLOCATABLE :: logchl_balinc_medusa_phn  !: Increment to MEDUSA nondiatom n   from logchl 
    165    REAL(wp), PUBLIC, DIMENSION(:,:,:), ALLOCATABLE :: logchl_balinc_medusa_phd  !: Increment to MEDUSA diatom n      from logchl 
    166    REAL(wp), PUBLIC, DIMENSION(:,:,:), ALLOCATABLE :: logchl_balinc_medusa_pds  !: Increment to MEDUSA diatom s      from logchl 
    167    REAL(wp), PUBLIC, DIMENSION(:,:,:), ALLOCATABLE :: logchl_balinc_medusa_zmi  !: Increment to MEDUSA microzoop n   from logchl 
    168    REAL(wp), PUBLIC, DIMENSION(:,:,:), ALLOCATABLE :: logchl_balinc_medusa_zme  !: Increment to MEDUSA mesozoop n    from logchl 
    169    REAL(wp), PUBLIC, DIMENSION(:,:,:), ALLOCATABLE :: logchl_balinc_medusa_din  !: Increment to MEDUSA nitrate       from logchl 
    170    REAL(wp), PUBLIC, DIMENSION(:,:,:), ALLOCATABLE :: logchl_balinc_medusa_sil  !: Increment to MEDUSA silicate      from logchl 
    171    REAL(wp), PUBLIC, DIMENSION(:,:,:), ALLOCATABLE :: logchl_balinc_medusa_fer  !: Increment to MEDUSA iron          from logchl 
    172    REAL(wp), PUBLIC, DIMENSION(:,:,:), ALLOCATABLE :: logchl_balinc_medusa_det  !: Increment to MEDUSA detritus n    from logchl 
    173    REAL(wp), PUBLIC, DIMENSION(:,:,:), ALLOCATABLE :: logchl_balinc_medusa_dtc  !: Increment to MEDUSA detritus c    from logchl 
    174    REAL(wp), PUBLIC, DIMENSION(:,:,:), ALLOCATABLE :: logchl_balinc_medusa_dic  !: Increment to MEDUSA dic           from logchl 
    175    REAL(wp), PUBLIC, DIMENSION(:,:,:), ALLOCATABLE :: logchl_balinc_medusa_alk  !: Increment to MEDUSA alkalinity    from logchl 
    176    REAL(wp), PUBLIC, DIMENSION(:,:,:), ALLOCATABLE :: logchl_balinc_medusa_oxy  !: Increment to MEDUSA oxygen        from logchl 
    177 #elif defined key_hadocc 
    178    REAL(wp), PUBLIC, DIMENSION(:,:,:), ALLOCATABLE :: logchl_balinc_hadocc_nut  !: Increment to HadOCC nutrient      from logchl 
    179    REAL(wp), PUBLIC, DIMENSION(:,:,:), ALLOCATABLE :: logchl_balinc_hadocc_phy  !: Increment to HadOCC phytoplankton from logchl 
    180    REAL(wp), PUBLIC, DIMENSION(:,:,:), ALLOCATABLE :: logchl_balinc_hadocc_zoo  !: Increment to HadOCC zooplankton   from logchl 
    181    REAL(wp), PUBLIC, DIMENSION(:,:,:), ALLOCATABLE :: logchl_balinc_hadocc_det  !: Increment to HadOCC detritus      from logchl 
    182    REAL(wp), PUBLIC, DIMENSION(:,:,:), ALLOCATABLE :: logchl_balinc_hadocc_dic  !: Increment to HadOCC DIC           from logchl 
    183    REAL(wp), PUBLIC, DIMENSION(:,:,:), ALLOCATABLE :: logchl_balinc_hadocc_alk  !: Increment to HadOCC alkalinity    from logchl 
     130#if defined key_top 
     131   REAL(wp), PUBLIC, DIMENSION(:,:,:,:), ALLOCATABLE :: logchl_balinc  !: Increment to BGC variables from logchl assim 
    184132#endif 
    185133 
     
    544492#endif 
    545493      IF ( ( ln_logchltotinc ).OR.( ln_logchlpftinc ) ) THEN 
    546          IF ( ln_logchltotinc ) THEN 
    547             ALLOCATE( logchl_bkginc(jpi,jpj,1)) 
    548          ELSE IF ( ln_logchlpftinc ) THEN 
    549             ALLOCATE( logchl_bkginc(jpi,jpj,nn_asmpfts)) 
    550          ENDIF 
     494         ALLOCATE( logchl_bkginc(jpi,jpj,nn_asmpfts)) 
    551495         logchl_bkginc(:,:,:) = 0.0 
    552 #if defined key_fabm 
    553          ALLOCATE( logchl_balinc_ersem_chl1(jpi,jpj,jpk) ) 
    554          ALLOCATE( logchl_balinc_ersem_chl2(jpi,jpj,jpk) ) 
    555          ALLOCATE( logchl_balinc_ersem_chl3(jpi,jpj,jpk) ) 
    556          ALLOCATE( logchl_balinc_ersem_chl4(jpi,jpj,jpk) ) 
    557          ALLOCATE( logchl_balinc_ersem_p1c(jpi,jpj,jpk)  ) 
    558          ALLOCATE( logchl_balinc_ersem_p1n(jpi,jpj,jpk)  ) 
    559          ALLOCATE( logchl_balinc_ersem_p1p(jpi,jpj,jpk)  ) 
    560          ALLOCATE( logchl_balinc_ersem_p1s(jpi,jpj,jpk)  ) 
    561          ALLOCATE( logchl_balinc_ersem_p2c(jpi,jpj,jpk)  ) 
    562          ALLOCATE( logchl_balinc_ersem_p2n(jpi,jpj,jpk)  ) 
    563          ALLOCATE( logchl_balinc_ersem_p2p(jpi,jpj,jpk)  ) 
    564          ALLOCATE( logchl_balinc_ersem_p3c(jpi,jpj,jpk)  ) 
    565          ALLOCATE( logchl_balinc_ersem_p3n(jpi,jpj,jpk)  ) 
    566          ALLOCATE( logchl_balinc_ersem_p3p(jpi,jpj,jpk)  ) 
    567          ALLOCATE( logchl_balinc_ersem_p4c(jpi,jpj,jpk)  ) 
    568          ALLOCATE( logchl_balinc_ersem_p4n(jpi,jpj,jpk)  ) 
    569          ALLOCATE( logchl_balinc_ersem_p4p(jpi,jpj,jpk)  ) 
    570          ALLOCATE( logchl_balinc_ersem_z4c(jpi,jpj,jpk)  ) 
    571          ALLOCATE( logchl_balinc_ersem_z5c(jpi,jpj,jpk)  ) 
    572          ALLOCATE( logchl_balinc_ersem_z5n(jpi,jpj,jpk)  ) 
    573          ALLOCATE( logchl_balinc_ersem_z5p(jpi,jpj,jpk)  ) 
    574          ALLOCATE( logchl_balinc_ersem_z6c(jpi,jpj,jpk)  ) 
    575          ALLOCATE( logchl_balinc_ersem_z6n(jpi,jpj,jpk)  ) 
    576          ALLOCATE( logchl_balinc_ersem_z6p(jpi,jpj,jpk)  ) 
    577          ALLOCATE( logchl_balinc_ersem_n1p(jpi,jpj,jpk)  ) 
    578          ALLOCATE( logchl_balinc_ersem_n3n(jpi,jpj,jpk)  ) 
    579          ALLOCATE( logchl_balinc_ersem_n4n(jpi,jpj,jpk)  ) 
    580          ALLOCATE( logchl_balinc_ersem_n5s(jpi,jpj,jpk)  ) 
    581          logchl_balinc_ersem_chl1(:,:,:) = 0.0 
    582          logchl_balinc_ersem_chl2(:,:,:) = 0.0 
    583          logchl_balinc_ersem_chl3(:,:,:) = 0.0 
    584          logchl_balinc_ersem_chl4(:,:,:) = 0.0 
    585          logchl_balinc_ersem_p1c(:,:,:)  = 0.0 
    586          logchl_balinc_ersem_p1n(:,:,:)  = 0.0 
    587          logchl_balinc_ersem_p1p(:,:,:)  = 0.0 
    588          logchl_balinc_ersem_p1s(:,:,:)  = 0.0 
    589          logchl_balinc_ersem_p2c(:,:,:)  = 0.0 
    590          logchl_balinc_ersem_p2n(:,:,:)  = 0.0 
    591          logchl_balinc_ersem_p2p(:,:,:)  = 0.0 
    592          logchl_balinc_ersem_p3c(:,:,:)  = 0.0 
    593          logchl_balinc_ersem_p3n(:,:,:)  = 0.0 
    594          logchl_balinc_ersem_p3p(:,:,:)  = 0.0 
    595          logchl_balinc_ersem_p4c(:,:,:)  = 0.0 
    596          logchl_balinc_ersem_p4n(:,:,:)  = 0.0 
    597          logchl_balinc_ersem_p4p(:,:,:)  = 0.0 
    598          logchl_balinc_ersem_z4c(:,:,:)  = 0.0 
    599          logchl_balinc_ersem_z5c(:,:,:)  = 0.0 
    600          logchl_balinc_ersem_z5n(:,:,:)  = 0.0 
    601          logchl_balinc_ersem_z5p(:,:,:)  = 0.0 
    602          logchl_balinc_ersem_z6c(:,:,:)  = 0.0 
    603          logchl_balinc_ersem_z6n(:,:,:)  = 0.0 
    604          logchl_balinc_ersem_z6p(:,:,:)  = 0.0 
    605          logchl_balinc_ersem_n1p(:,:,:)  = 0.0 
    606          logchl_balinc_ersem_n3n(:,:,:)  = 0.0 
    607          logchl_balinc_ersem_n4n(:,:,:)  = 0.0 
    608          logchl_balinc_ersem_n5s(:,:,:)  = 0.0 
    609 #elif defined key_medusa && defined key_foam_medusa 
    610          ALLOCATE( logchl_balinc_medusa_chn(jpi,jpj,jpk) ) 
    611          ALLOCATE( logchl_balinc_medusa_chd(jpi,jpj,jpk) ) 
    612          ALLOCATE( logchl_balinc_medusa_phn(jpi,jpj,jpk) ) 
    613          ALLOCATE( logchl_balinc_medusa_phd(jpi,jpj,jpk) ) 
    614          ALLOCATE( logchl_balinc_medusa_pds(jpi,jpj,jpk) ) 
    615          ALLOCATE( logchl_balinc_medusa_zmi(jpi,jpj,jpk) ) 
    616          ALLOCATE( logchl_balinc_medusa_zme(jpi,jpj,jpk) ) 
    617          ALLOCATE( logchl_balinc_medusa_din(jpi,jpj,jpk) ) 
    618          ALLOCATE( logchl_balinc_medusa_sil(jpi,jpj,jpk) ) 
    619          ALLOCATE( logchl_balinc_medusa_fer(jpi,jpj,jpk) ) 
    620          ALLOCATE( logchl_balinc_medusa_det(jpi,jpj,jpk) ) 
    621          ALLOCATE( logchl_balinc_medusa_dtc(jpi,jpj,jpk) ) 
    622          ALLOCATE( logchl_balinc_medusa_dic(jpi,jpj,jpk) ) 
    623          ALLOCATE( logchl_balinc_medusa_alk(jpi,jpj,jpk) ) 
    624          ALLOCATE( logchl_balinc_medusa_oxy(jpi,jpj,jpk) ) 
    625          logchl_balinc_medusa_chn(:,:,:) = 0.0 
    626          logchl_balinc_medusa_chd(:,:,:) = 0.0 
    627          logchl_balinc_medusa_phn(:,:,:) = 0.0 
    628          logchl_balinc_medusa_phd(:,:,:) = 0.0 
    629          logchl_balinc_medusa_pds(:,:,:) = 0.0 
    630          logchl_balinc_medusa_zmi(:,:,:) = 0.0 
    631          logchl_balinc_medusa_zme(:,:,:) = 0.0 
    632          logchl_balinc_medusa_din(:,:,:) = 0.0 
    633          logchl_balinc_medusa_sil(:,:,:) = 0.0 
    634          logchl_balinc_medusa_fer(:,:,:) = 0.0 
    635          logchl_balinc_medusa_det(:,:,:) = 0.0 
    636          logchl_balinc_medusa_dtc(:,:,:) = 0.0 
    637          logchl_balinc_medusa_dic(:,:,:) = 0.0 
    638          logchl_balinc_medusa_alk(:,:,:) = 0.0 
    639          logchl_balinc_medusa_oxy(:,:,:) = 0.0 
    640 #elif defined key_hadocc 
    641          ALLOCATE( logchl_balinc_hadocc_nut(jpi,jpj,jpk) ) 
    642          ALLOCATE( logchl_balinc_hadocc_phy(jpi,jpj,jpk) ) 
    643          ALLOCATE( logchl_balinc_hadocc_zoo(jpi,jpj,jpk) ) 
    644          ALLOCATE( logchl_balinc_hadocc_det(jpi,jpj,jpk) ) 
    645          ALLOCATE( logchl_balinc_hadocc_dic(jpi,jpj,jpk) ) 
    646          ALLOCATE( logchl_balinc_hadocc_alk(jpi,jpj,jpk) ) 
    647          logchl_balinc_hadocc_nut(:,:,:) = 0.0 
    648          logchl_balinc_hadocc_phy(:,:,:) = 0.0 
    649          logchl_balinc_hadocc_zoo(:,:,:) = 0.0 
    650          logchl_balinc_hadocc_det(:,:,:) = 0.0 
    651          logchl_balinc_hadocc_dic(:,:,:) = 0.0 
    652          logchl_balinc_hadocc_alk(:,:,:) = 0.0 
     496#if defined key_top 
     497         ALLOCATE( logchl_balinc(jpi,jpj,jpk,jptra) ) 
     498         logchl_balinc(:,:,:,:) = 0.0 
    653499#endif 
    654500      ENDIF 
     
    15611407 
    15621408#if defined key_fabm 
    1563          CALL asm_logchl_bal_ersem( ln_logchlpftinc,          nn_asmpfts,               & 
    1564             &                       mld_choice_bgc,           logchl_bkginc,            & 
    1565             &                       logchl_balinc_ersem_chl1, logchl_balinc_ersem_chl2, & 
    1566             &                       logchl_balinc_ersem_chl3, logchl_balinc_ersem_chl4, & 
    1567             &                       logchl_balinc_ersem_p1c,  logchl_balinc_ersem_p1n,  & 
    1568             &                       logchl_balinc_ersem_p1p,  logchl_balinc_ersem_p1s,  & 
    1569             &                       logchl_balinc_ersem_p2c,  logchl_balinc_ersem_p2n,  & 
    1570             &                       logchl_balinc_ersem_p2p,  logchl_balinc_ersem_p3c,  & 
    1571             &                       logchl_balinc_ersem_p3n,  logchl_balinc_ersem_p3p,  & 
    1572             &                       logchl_balinc_ersem_p4c,  logchl_balinc_ersem_p4n,  & 
    1573             &                       logchl_balinc_ersem_p4p,  logchl_balinc_ersem_z4c,  & 
    1574             &                       logchl_balinc_ersem_z5c,  logchl_balinc_ersem_z5n,  & 
    1575             &                       logchl_balinc_ersem_z5p,  logchl_balinc_ersem_z6c,  & 
    1576             &                       logchl_balinc_ersem_z6n,  logchl_balinc_ersem_z6p,  & 
    1577             &                       logchl_balinc_ersem_n1p,  logchl_balinc_ersem_n3n,  & 
    1578             &                       logchl_balinc_ersem_n4n,  logchl_balinc_ersem_n5s ) 
     1409         CALL asm_logchl_bal_ersem( ln_logchlpftinc, nn_asmpfts, mld_choice_bgc, & 
     1410            &                       logchl_bkginc, logchl_balinc ) 
    15791411#elif defined key_medusa && defined key_foam_medusa 
    15801412         !CALL asm_logchl_bal_medusa() 
     
    16111443            ENDIF 
    16121444 
    1613             ! Update the biogeochemical tendencies 
     1445            ! Update the biogeochemical variables 
     1446            ! Add directly to trn and trb, rather than to tra, as not a tendency 
    16141447#if defined key_fabm 
    16151448            DO jk = 1, jpkm1 
    1616                ! Add directly to trn and trb, rather than to tra, as not a tendency 
    1617                trn(:,:,jk,jp_fabm_m1+jp_fabm_chl1) = trn(:,:,jk,jp_fabm_m1+jp_fabm_chl1) + & 
    1618                   &                                  logchl_balinc_ersem_chl1(:,:,jk) * zincwgt 
    1619                trn(:,:,jk,jp_fabm_m1+jp_fabm_chl2) = trn(:,:,jk,jp_fabm_m1+jp_fabm_chl2) + & 
    1620                   &                                  logchl_balinc_ersem_chl2(:,:,jk) * zincwgt 
    1621                trn(:,:,jk,jp_fabm_m1+jp_fabm_chl3) = trn(:,:,jk,jp_fabm_m1+jp_fabm_chl3) + & 
    1622                   &                                  logchl_balinc_ersem_chl3(:,:,jk) * zincwgt 
    1623                trn(:,:,jk,jp_fabm_m1+jp_fabm_chl4) = trn(:,:,jk,jp_fabm_m1+jp_fabm_chl4) + & 
    1624                   &                                  logchl_balinc_ersem_chl4(:,:,jk) * zincwgt 
    1625                trn(:,:,jk,jp_fabm_m1+jp_fabm_p1c)  = trn(:,:,jk,jp_fabm_m1+jp_fabm_p1c) + & 
    1626                   &                                  logchl_balinc_ersem_p1c(:,:,jk)  * zincwgt 
    1627                trn(:,:,jk,jp_fabm_m1+jp_fabm_p1n)  = trn(:,:,jk,jp_fabm_m1+jp_fabm_p1n) + & 
    1628                   &                                  logchl_balinc_ersem_p1n(:,:,jk)  * zincwgt 
    1629                trn(:,:,jk,jp_fabm_m1+jp_fabm_p1p)  = trn(:,:,jk,jp_fabm_m1+jp_fabm_p1p) + & 
    1630                   &                                  logchl_balinc_ersem_p1p(:,:,jk)  * zincwgt 
    1631                trn(:,:,jk,jp_fabm_m1+jp_fabm_p1s)  = trn(:,:,jk,jp_fabm_m1+jp_fabm_p1s) + & 
    1632                   &                                  logchl_balinc_ersem_p1s(:,:,jk)  * zincwgt 
    1633                trn(:,:,jk,jp_fabm_m1+jp_fabm_p2c)  = trn(:,:,jk,jp_fabm_m1+jp_fabm_p2c) + & 
    1634                   &                                  logchl_balinc_ersem_p2c(:,:,jk)  * zincwgt 
    1635                trn(:,:,jk,jp_fabm_m1+jp_fabm_p2n)  = trn(:,:,jk,jp_fabm_m1+jp_fabm_p2n) + & 
    1636                   &                                  logchl_balinc_ersem_p2n(:,:,jk)  * zincwgt 
    1637                trn(:,:,jk,jp_fabm_m1+jp_fabm_p2p)  = trn(:,:,jk,jp_fabm_m1+jp_fabm_p2p) + & 
    1638                   &                                  logchl_balinc_ersem_p2p(:,:,jk)  * zincwgt 
    1639                trn(:,:,jk,jp_fabm_m1+jp_fabm_p3c)  = trn(:,:,jk,jp_fabm_m1+jp_fabm_p3c) + & 
    1640                   &                                  logchl_balinc_ersem_p3c(:,:,jk)  * zincwgt 
    1641                trn(:,:,jk,jp_fabm_m1+jp_fabm_p3n)  = trn(:,:,jk,jp_fabm_m1+jp_fabm_p3n) + & 
    1642                   &                                  logchl_balinc_ersem_p3n(:,:,jk)  * zincwgt 
    1643                trn(:,:,jk,jp_fabm_m1+jp_fabm_p3p)  = trn(:,:,jk,jp_fabm_m1+jp_fabm_p3p) + & 
    1644                   &                                  logchl_balinc_ersem_p3p(:,:,jk)  * zincwgt 
    1645                trn(:,:,jk,jp_fabm_m1+jp_fabm_p4c)  = trn(:,:,jk,jp_fabm_m1+jp_fabm_p4c) + & 
    1646                   &                                  logchl_balinc_ersem_p4c(:,:,jk)  * zincwgt 
    1647                trn(:,:,jk,jp_fabm_m1+jp_fabm_p4n)  = trn(:,:,jk,jp_fabm_m1+jp_fabm_p4n) + & 
    1648                   &                                  logchl_balinc_ersem_p4n(:,:,jk)  * zincwgt 
    1649                trn(:,:,jk,jp_fabm_m1+jp_fabm_p4p)  = trn(:,:,jk,jp_fabm_m1+jp_fabm_p4p) + & 
    1650                   &                                  logchl_balinc_ersem_p4p(:,:,jk)  * zincwgt 
    1651                trn(:,:,jk,jp_fabm_m1+jp_fabm_z4c)  = trn(:,:,jk,jp_fabm_m1+jp_fabm_z4c) + & 
    1652                   &                                  logchl_balinc_ersem_z4c(:,:,jk)  * zincwgt 
    1653                trn(:,:,jk,jp_fabm_m1+jp_fabm_z5c)  = trn(:,:,jk,jp_fabm_m1+jp_fabm_z5c) + & 
    1654                   &                                  logchl_balinc_ersem_z5c(:,:,jk)  * zincwgt 
    1655                trn(:,:,jk,jp_fabm_m1+jp_fabm_z5n)  = trn(:,:,jk,jp_fabm_m1+jp_fabm_z5n) + & 
    1656                   &                                  logchl_balinc_ersem_z5n(:,:,jk)  * zincwgt 
    1657                trn(:,:,jk,jp_fabm_m1+jp_fabm_z5p)  = trn(:,:,jk,jp_fabm_m1+jp_fabm_z5p) + & 
    1658                   &                                  logchl_balinc_ersem_z5p(:,:,jk)  * zincwgt 
    1659                trn(:,:,jk,jp_fabm_m1+jp_fabm_z6c)  = trn(:,:,jk,jp_fabm_m1+jp_fabm_z6c) + & 
    1660                   &                                  logchl_balinc_ersem_z6c(:,:,jk)  * zincwgt 
    1661                trn(:,:,jk,jp_fabm_m1+jp_fabm_z6n)  = trn(:,:,jk,jp_fabm_m1+jp_fabm_z6n) + & 
    1662                   &                                  logchl_balinc_ersem_z6n(:,:,jk)  * zincwgt 
    1663                trn(:,:,jk,jp_fabm_m1+jp_fabm_z6p)  = trn(:,:,jk,jp_fabm_m1+jp_fabm_z6p) + & 
    1664                   &                                  logchl_balinc_ersem_z6p(:,:,jk)  * zincwgt 
    1665                trn(:,:,jk,jp_fabm_m1+jp_fabm_n1p)  = trn(:,:,jk,jp_fabm_m1+jp_fabm_n1p) + & 
    1666                   &                                  logchl_balinc_ersem_n1p(:,:,jk)  * zincwgt 
    1667                trn(:,:,jk,jp_fabm_m1+jp_fabm_n3n)  = trn(:,:,jk,jp_fabm_m1+jp_fabm_n3n) + & 
    1668                   &                                  logchl_balinc_ersem_n3n(:,:,jk)  * zincwgt 
    1669                trn(:,:,jk,jp_fabm_m1+jp_fabm_n4n)  = trn(:,:,jk,jp_fabm_m1+jp_fabm_n4n) + & 
    1670                   &                                  logchl_balinc_ersem_n4n(:,:,jk)  * zincwgt 
    1671                trn(:,:,jk,jp_fabm_m1+jp_fabm_n5s)  = trn(:,:,jk,jp_fabm_m1+jp_fabm_n5s) + & 
    1672                   &                                  logchl_balinc_ersem_n5s(:,:,jk)  * zincwgt 
    1673  
    1674                trb(:,:,jk,jp_fabm_m1+jp_fabm_chl1) = trb(:,:,jk,jp_fabm_m1+jp_fabm_chl1) + & 
    1675                   &                                  logchl_balinc_ersem_chl1(:,:,jk) * zincwgt 
    1676                trb(:,:,jk,jp_fabm_m1+jp_fabm_chl2) = trb(:,:,jk,jp_fabm_m1+jp_fabm_chl2) + & 
    1677                   &                                  logchl_balinc_ersem_chl2(:,:,jk) * zincwgt 
    1678                trb(:,:,jk,jp_fabm_m1+jp_fabm_chl3) = trb(:,:,jk,jp_fabm_m1+jp_fabm_chl3) + & 
    1679                   &                                  logchl_balinc_ersem_chl3(:,:,jk) * zincwgt 
    1680                trb(:,:,jk,jp_fabm_m1+jp_fabm_chl4) = trb(:,:,jk,jp_fabm_m1+jp_fabm_chl4) + & 
    1681                   &                                  logchl_balinc_ersem_chl4(:,:,jk) * zincwgt 
    1682                trb(:,:,jk,jp_fabm_m1+jp_fabm_p1c)  = trb(:,:,jk,jp_fabm_m1+jp_fabm_p1c) + & 
    1683                   &                                  logchl_balinc_ersem_p1c(:,:,jk)  * zincwgt 
    1684                trb(:,:,jk,jp_fabm_m1+jp_fabm_p1n)  = trb(:,:,jk,jp_fabm_m1+jp_fabm_p1n) + & 
    1685                   &                                  logchl_balinc_ersem_p1n(:,:,jk)  * zincwgt 
    1686                trb(:,:,jk,jp_fabm_m1+jp_fabm_p1p)  = trb(:,:,jk,jp_fabm_m1+jp_fabm_p1p) + & 
    1687                   &                                  logchl_balinc_ersem_p1p(:,:,jk)  * zincwgt 
    1688                trb(:,:,jk,jp_fabm_m1+jp_fabm_p1s)  = trb(:,:,jk,jp_fabm_m1+jp_fabm_p1s) + & 
    1689                   &                                  logchl_balinc_ersem_p1s(:,:,jk)  * zincwgt 
    1690                trb(:,:,jk,jp_fabm_m1+jp_fabm_p2c)  = trb(:,:,jk,jp_fabm_m1+jp_fabm_p2c) + & 
    1691                   &                                  logchl_balinc_ersem_p2c(:,:,jk)  * zincwgt 
    1692                trb(:,:,jk,jp_fabm_m1+jp_fabm_p2n)  = trb(:,:,jk,jp_fabm_m1+jp_fabm_p2n) + & 
    1693                   &                                  logchl_balinc_ersem_p2n(:,:,jk)  * zincwgt 
    1694                trb(:,:,jk,jp_fabm_m1+jp_fabm_p2p)  = trb(:,:,jk,jp_fabm_m1+jp_fabm_p2p) + & 
    1695                   &                                  logchl_balinc_ersem_p2p(:,:,jk)  * zincwgt 
    1696                trb(:,:,jk,jp_fabm_m1+jp_fabm_p3c)  = trb(:,:,jk,jp_fabm_m1+jp_fabm_p3c) + & 
    1697                   &                                  logchl_balinc_ersem_p3c(:,:,jk)  * zincwgt 
    1698                trb(:,:,jk,jp_fabm_m1+jp_fabm_p3n)  = trb(:,:,jk,jp_fabm_m1+jp_fabm_p3n) + & 
    1699                   &                                  logchl_balinc_ersem_p3n(:,:,jk)  * zincwgt 
    1700                trb(:,:,jk,jp_fabm_m1+jp_fabm_p3p)  = trb(:,:,jk,jp_fabm_m1+jp_fabm_p3p) + & 
    1701                   &                                  logchl_balinc_ersem_p3p(:,:,jk)  * zincwgt 
    1702                trb(:,:,jk,jp_fabm_m1+jp_fabm_p4c)  = trb(:,:,jk,jp_fabm_m1+jp_fabm_p4c) + & 
    1703                   &                                  logchl_balinc_ersem_p4c(:,:,jk)  * zincwgt 
    1704                trb(:,:,jk,jp_fabm_m1+jp_fabm_p4n)  = trb(:,:,jk,jp_fabm_m1+jp_fabm_p4n) + & 
    1705                   &                                  logchl_balinc_ersem_p4n(:,:,jk)  * zincwgt 
    1706                trb(:,:,jk,jp_fabm_m1+jp_fabm_p4p)  = trb(:,:,jk,jp_fabm_m1+jp_fabm_p4p) + & 
    1707                   &                                  logchl_balinc_ersem_p4p(:,:,jk)  * zincwgt 
    1708                trb(:,:,jk,jp_fabm_m1+jp_fabm_z4c)  = trb(:,:,jk,jp_fabm_m1+jp_fabm_z4c) + & 
    1709                   &                                  logchl_balinc_ersem_z4c(:,:,jk)  * zincwgt 
    1710                trb(:,:,jk,jp_fabm_m1+jp_fabm_z5c)  = trb(:,:,jk,jp_fabm_m1+jp_fabm_z5c) + & 
    1711                   &                                  logchl_balinc_ersem_z5c(:,:,jk)  * zincwgt 
    1712                trb(:,:,jk,jp_fabm_m1+jp_fabm_z5n)  = trb(:,:,jk,jp_fabm_m1+jp_fabm_z5n) + & 
    1713                   &                                  logchl_balinc_ersem_z5n(:,:,jk)  * zincwgt 
    1714                trb(:,:,jk,jp_fabm_m1+jp_fabm_z5p)  = trb(:,:,jk,jp_fabm_m1+jp_fabm_z5p) + & 
    1715                   &                                  logchl_balinc_ersem_z5p(:,:,jk)  * zincwgt 
    1716                trb(:,:,jk,jp_fabm_m1+jp_fabm_z6c)  = trb(:,:,jk,jp_fabm_m1+jp_fabm_z6c) + & 
    1717                   &                                  logchl_balinc_ersem_z6c(:,:,jk)  * zincwgt 
    1718                trb(:,:,jk,jp_fabm_m1+jp_fabm_z6n)  = trb(:,:,jk,jp_fabm_m1+jp_fabm_z6n) + & 
    1719                   &                                  logchl_balinc_ersem_z6n(:,:,jk)  * zincwgt 
    1720                trb(:,:,jk,jp_fabm_m1+jp_fabm_z6p)  = trb(:,:,jk,jp_fabm_m1+jp_fabm_z6p) + & 
    1721                   &                                  logchl_balinc_ersem_z6p(:,:,jk)  * zincwgt 
    1722                trb(:,:,jk,jp_fabm_m1+jp_fabm_n1p)  = trb(:,:,jk,jp_fabm_m1+jp_fabm_n1p) + & 
    1723                   &                                  logchl_balinc_ersem_n1p(:,:,jk)  * zincwgt 
    1724                trb(:,:,jk,jp_fabm_m1+jp_fabm_n3n)  = trb(:,:,jk,jp_fabm_m1+jp_fabm_n3n) + & 
    1725                   &                                  logchl_balinc_ersem_n3n(:,:,jk)  * zincwgt 
    1726                trb(:,:,jk,jp_fabm_m1+jp_fabm_n4n)  = trb(:,:,jk,jp_fabm_m1+jp_fabm_n4n) + & 
    1727                   &                                  logchl_balinc_ersem_n4n(:,:,jk)  * zincwgt 
    1728                trb(:,:,jk,jp_fabm_m1+jp_fabm_n5s)  = trb(:,:,jk,jp_fabm_m1+jp_fabm_n5s) + & 
    1729                   &                                  logchl_balinc_ersem_n5s(:,:,jk)  * zincwgt 
     1449               trn(:,:,jk,jp_fabm0:jp_fabm1) = trn(:,:,jk,jp_fabm0:jp_fabm1) + & 
     1450                  &                            logchl_balinc(:,:,jk,jp_fabm0:jp_fabm1) * zincwgt 
     1451               trb(:,:,jk,jp_fabm0:jp_fabm1) = trb(:,:,jk,jp_fabm0:jp_fabm1) + & 
     1452                  &                            logchl_balinc(:,:,jk,jp_fabm0:jp_fabm1) * zincwgt 
    17301453            END DO 
    1731  
    17321454#elif defined key_medusa && defined key_foam_medusa 
    17331455            DO jk = 1, jpkm1 
    1734                ! Add directly to trn and trb, rather than to tra, as not a tendency 
    1735                trn(:,:,jk,jpchn) = trn(:,:,jk,jpchn) + logchl_balinc_medusa_chn(:,:,jk) * zincwgt 
    1736                trn(:,:,jk,jpchd) = trn(:,:,jk,jpchd) + logchl_balinc_medusa_chd(:,:,jk) * zincwgt 
    1737                trn(:,:,jk,jpphn) = trn(:,:,jk,jpphn) + logchl_balinc_medusa_phn(:,:,jk) * zincwgt 
    1738                trn(:,:,jk,jpphd) = trn(:,:,jk,jpphd) + logchl_balinc_medusa_phd(:,:,jk) * zincwgt 
    1739                trn(:,:,jk,jpzmi) = trn(:,:,jk,jpzmi) + logchl_balinc_medusa_zmi(:,:,jk) * zincwgt 
    1740                trn(:,:,jk,jpzme) = trn(:,:,jk,jpzme) + logchl_balinc_medusa_zme(:,:,jk) * zincwgt 
    1741                trn(:,:,jk,jpdin) = trn(:,:,jk,jpdin) + logchl_balinc_medusa_din(:,:,jk) * zincwgt 
    1742                trn(:,:,jk,jpsil) = trn(:,:,jk,jpsil) + logchl_balinc_medusa_sil(:,:,jk) * zincwgt 
    1743                trn(:,:,jk,jpfer) = trn(:,:,jk,jpfer) + logchl_balinc_medusa_fer(:,:,jk) * zincwgt 
    1744                trn(:,:,jk,jpdet) = trn(:,:,jk,jpdet) + logchl_balinc_medusa_det(:,:,jk) * zincwgt 
    1745                trn(:,:,jk,jppds) = trn(:,:,jk,jppds) + logchl_balinc_medusa_pds(:,:,jk) * zincwgt 
    1746                trn(:,:,jk,jpdtc) = trn(:,:,jk,jpdtc) + logchl_balinc_medusa_dtc(:,:,jk) * zincwgt 
    1747                trn(:,:,jk,jpdic) = trn(:,:,jk,jpdic) + logchl_balinc_medusa_dic(:,:,jk) * zincwgt 
    1748                trn(:,:,jk,jpalk) = trn(:,:,jk,jpalk) + logchl_balinc_medusa_alk(:,:,jk) * zincwgt 
    1749                trn(:,:,jk,jpoxy) = trn(:,:,jk,jpoxy) + logchl_balinc_medusa_oxy(:,:,jk) * zincwgt 
    1750  
    1751                trb(:,:,jk,jpchn) = trb(:,:,jk,jpchn) + logchl_balinc_medusa_chn(:,:,jk) * zincwgt 
    1752                trb(:,:,jk,jpchd) = trb(:,:,jk,jpchd) + logchl_balinc_medusa_chd(:,:,jk) * zincwgt 
    1753                trb(:,:,jk,jpphn) = trb(:,:,jk,jpphn) + logchl_balinc_medusa_phn(:,:,jk) * zincwgt 
    1754                trb(:,:,jk,jpphd) = trb(:,:,jk,jpphd) + logchl_balinc_medusa_phd(:,:,jk) * zincwgt 
    1755                trb(:,:,jk,jpzmi) = trb(:,:,jk,jpzmi) + logchl_balinc_medusa_zmi(:,:,jk) * zincwgt 
    1756                trb(:,:,jk,jpzme) = trb(:,:,jk,jpzme) + logchl_balinc_medusa_zme(:,:,jk) * zincwgt 
    1757                trb(:,:,jk,jpdin) = trb(:,:,jk,jpdin) + logchl_balinc_medusa_din(:,:,jk) * zincwgt 
    1758                trb(:,:,jk,jpsil) = trb(:,:,jk,jpsil) + logchl_balinc_medusa_sil(:,:,jk) * zincwgt 
    1759                trb(:,:,jk,jpfer) = trb(:,:,jk,jpfer) + logchl_balinc_medusa_fer(:,:,jk) * zincwgt 
    1760                trb(:,:,jk,jpdet) = trb(:,:,jk,jpdet) + logchl_balinc_medusa_det(:,:,jk) * zincwgt 
    1761                trb(:,:,jk,jppds) = trb(:,:,jk,jppds) + logchl_balinc_medusa_pds(:,:,jk) * zincwgt 
    1762                trb(:,:,jk,jpdtc) = trb(:,:,jk,jpdtc) + logchl_balinc_medusa_dtc(:,:,jk) * zincwgt 
    1763                trb(:,:,jk,jpdic) = trb(:,:,jk,jpdic) + logchl_balinc_medusa_dic(:,:,jk) * zincwgt 
    1764                trb(:,:,jk,jpalk) = trb(:,:,jk,jpalk) + logchl_balinc_medusa_alk(:,:,jk) * zincwgt 
    1765                trb(:,:,jk,jpoxy) = trb(:,:,jk,jpoxy) + logchl_balinc_medusa_oxy(:,:,jk) * zincwgt 
     1456               trn(:,:,jk,jp_msa0:jp_msa1) = trn(:,:,jk,jp_msa0:jp_msa1) + & 
     1457                  &                          logchl_balinc(:,:,jk,jp_msa0:jp_msa1) * zincwgt 
     1458               trb(:,:,jk,jp_msa0:jp_msa1) = trb(:,:,jk,jp_msa0:jp_msa1) + & 
     1459                  &                          logchl_balinc(:,:,jk,jp_msa0:jp_msa1) * zincwgt 
    17661460            END DO 
    1767  
    17681461#elif defined key_hadocc 
    17691462            DO jk = 1, jpkm1 
    1770                ! Add directly to trn and trb, rather than to tra, as not a tendency 
    1771                trn(:,:,jk,jp_had_nut) = trn(:,:,jk,jp_had_nut) + logchl_balinc_hadocc_nut(:,:,jk) * zincwgt 
    1772                trn(:,:,jk,jp_had_phy) = trn(:,:,jk,jp_had_phy) + logchl_balinc_hadocc_phy(:,:,jk) * zincwgt 
    1773                trn(:,:,jk,jp_had_zoo) = trn(:,:,jk,jp_had_zoo) + logchl_balinc_hadocc_zoo(:,:,jk) * zincwgt 
    1774                trn(:,:,jk,jp_had_pdn) = trn(:,:,jk,jp_had_pdn) + logchl_balinc_hadocc_det(:,:,jk) * zincwgt 
    1775                trn(:,:,jk,jp_had_dic) = trn(:,:,jk,jp_had_dic) + logchl_balinc_hadocc_dic(:,:,jk) * zincwgt 
    1776                trn(:,:,jk,jp_had_alk) = trn(:,:,jk,jp_had_alk) + logchl_balinc_hadocc_alk(:,:,jk) * zincwgt 
    1777              
    1778                trb(:,:,jk,jp_had_nut) = trb(:,:,jk,jp_had_nut) + logchl_balinc_hadocc_nut(:,:,jk) * zincwgt 
    1779                trb(:,:,jk,jp_had_phy) = trb(:,:,jk,jp_had_phy) + logchl_balinc_hadocc_phy(:,:,jk) * zincwgt 
    1780                trb(:,:,jk,jp_had_zoo) = trb(:,:,jk,jp_had_zoo) + logchl_balinc_hadocc_zoo(:,:,jk) * zincwgt 
    1781                trb(:,:,jk,jp_had_pdn) = trb(:,:,jk,jp_had_pdn) + logchl_balinc_hadocc_det(:,:,jk) * zincwgt 
    1782                trb(:,:,jk,jp_had_dic) = trb(:,:,jk,jp_had_dic) + logchl_balinc_hadocc_dic(:,:,jk) * zincwgt 
    1783                trb(:,:,jk,jp_had_alk) = trb(:,:,jk,jp_had_alk) + logchl_balinc_hadocc_alk(:,:,jk) * zincwgt 
     1463               trn(:,:,jk,jp_had0:jp_had1) = trn(:,:,jk,jp_had0:jp_had1) + & 
     1464                  &                          logchl_balinc(:,:,jk,jp_had0:jp_had1) * zincwgt 
     1465               trb(:,:,jk,jp_had0:jp_had1) = trb(:,:,jk,jp_had0:jp_had1) + & 
     1466                  &                          logchl_balinc(:,:,jk,jp_had0:jp_had1) * zincwgt 
    17841467            END DO 
    17851468#endif 
     
    18051488            CALL ctl_warn( ' Doing direct initialisation of ERSEM with chlorophyll assimilation', & 
    18061489               &           ' Background state is taken from model rather than background file' ) 
    1807             trn(:,:,:,jp_fabm_m1+jp_fabm_chl1) = trn(:,:,:,jp_fabm_m1+jp_fabm_chl1) + & 
    1808                &                                 logchl_balinc_ersem_chl1(:,:,:) 
    1809             trn(:,:,:,jp_fabm_m1+jp_fabm_chl2) = trn(:,:,:,jp_fabm_m1+jp_fabm_chl2) + & 
    1810                &                                 logchl_balinc_ersem_chl2(:,:,:) 
    1811             trn(:,:,:,jp_fabm_m1+jp_fabm_chl3) = trn(:,:,:,jp_fabm_m1+jp_fabm_chl3) + & 
    1812                &                                 logchl_balinc_ersem_chl3(:,:,:) 
    1813             trn(:,:,:,jp_fabm_m1+jp_fabm_chl4) = trn(:,:,:,jp_fabm_m1+jp_fabm_chl4) + & 
    1814                &                                 logchl_balinc_ersem_chl4(:,:,:) 
    1815             trn(:,:,:,jp_fabm_m1+jp_fabm_p1c)  = trn(:,:,:,jp_fabm_m1+jp_fabm_p1c)  + & 
    1816                &                                 logchl_balinc_ersem_p1c(:,:,:) 
    1817             trn(:,:,:,jp_fabm_m1+jp_fabm_p1n)  = trn(:,:,:,jp_fabm_m1+jp_fabm_p1n)  + & 
    1818                &                                 logchl_balinc_ersem_p1n(:,:,:) 
    1819             trn(:,:,:,jp_fabm_m1+jp_fabm_p1p)  = trn(:,:,:,jp_fabm_m1+jp_fabm_p1p)  + & 
    1820                &                                 logchl_balinc_ersem_p1p(:,:,:) 
    1821             trn(:,:,:,jp_fabm_m1+jp_fabm_p1s)  = trn(:,:,:,jp_fabm_m1+jp_fabm_p1s)  + & 
    1822                &                                 logchl_balinc_ersem_p1s(:,:,:) 
    1823             trn(:,:,:,jp_fabm_m1+jp_fabm_p2c)  = trn(:,:,:,jp_fabm_m1+jp_fabm_p2c)  + & 
    1824                &                                 logchl_balinc_ersem_p2c(:,:,:) 
    1825             trn(:,:,:,jp_fabm_m1+jp_fabm_p2n)  = trn(:,:,:,jp_fabm_m1+jp_fabm_p2n)  + & 
    1826                &                                 logchl_balinc_ersem_p2n(:,:,:) 
    1827             trn(:,:,:,jp_fabm_m1+jp_fabm_p2p)  = trn(:,:,:,jp_fabm_m1+jp_fabm_p2p)  + & 
    1828                &                                 logchl_balinc_ersem_p2p(:,:,:) 
    1829             trn(:,:,:,jp_fabm_m1+jp_fabm_p3c)  = trn(:,:,:,jp_fabm_m1+jp_fabm_p3c)  + & 
    1830                &                                 logchl_balinc_ersem_p3c(:,:,:) 
    1831             trn(:,:,:,jp_fabm_m1+jp_fabm_p3n)  = trn(:,:,:,jp_fabm_m1+jp_fabm_p3n)  + & 
    1832                &                                 logchl_balinc_ersem_p3n(:,:,:) 
    1833             trn(:,:,:,jp_fabm_m1+jp_fabm_p3p)  = trn(:,:,:,jp_fabm_m1+jp_fabm_p3p)  + & 
    1834                &                                 logchl_balinc_ersem_p3p(:,:,:) 
    1835             trn(:,:,:,jp_fabm_m1+jp_fabm_p4c)  = trn(:,:,:,jp_fabm_m1+jp_fabm_p4c)  + & 
    1836                &                                 logchl_balinc_ersem_p4c(:,:,:) 
    1837             trn(:,:,:,jp_fabm_m1+jp_fabm_p4n)  = trn(:,:,:,jp_fabm_m1+jp_fabm_p4n)  + & 
    1838                &                                 logchl_balinc_ersem_p4n(:,:,:) 
    1839             trn(:,:,:,jp_fabm_m1+jp_fabm_p4p)  = trn(:,:,:,jp_fabm_m1+jp_fabm_p4p)  + & 
    1840                &                                 logchl_balinc_ersem_p4p(:,:,:) 
    1841             trn(:,:,:,jp_fabm_m1+jp_fabm_z4c)  = trn(:,:,:,jp_fabm_m1+jp_fabm_z4c)  + & 
    1842                &                                 logchl_balinc_ersem_z4c(:,:,:) 
    1843             trn(:,:,:,jp_fabm_m1+jp_fabm_z5c)  = trn(:,:,:,jp_fabm_m1+jp_fabm_z5c)  + & 
    1844                &                                 logchl_balinc_ersem_z5c(:,:,:) 
    1845             trn(:,:,:,jp_fabm_m1+jp_fabm_z5n)  = trn(:,:,:,jp_fabm_m1+jp_fabm_z5n)  + & 
    1846                &                                 logchl_balinc_ersem_z5n(:,:,:) 
    1847             trn(:,:,:,jp_fabm_m1+jp_fabm_z5p)  = trn(:,:,:,jp_fabm_m1+jp_fabm_z5p)  + & 
    1848                &                                 logchl_balinc_ersem_z5p(:,:,:) 
    1849             trn(:,:,:,jp_fabm_m1+jp_fabm_z6c)  = trn(:,:,:,jp_fabm_m1+jp_fabm_z6c)  + & 
    1850                &                                 logchl_balinc_ersem_z6c(:,:,:) 
    1851             trn(:,:,:,jp_fabm_m1+jp_fabm_z6n)  = trn(:,:,:,jp_fabm_m1+jp_fabm_z6n)  + & 
    1852                &                                 logchl_balinc_ersem_z6n(:,:,:) 
    1853             trn(:,:,:,jp_fabm_m1+jp_fabm_z6p)  = trn(:,:,:,jp_fabm_m1+jp_fabm_z6p)  + & 
    1854                &                                 logchl_balinc_ersem_z6p(:,:,:) 
    1855             trn(:,:,:,jp_fabm_m1+jp_fabm_n1p)  = trn(:,:,:,jp_fabm_m1+jp_fabm_n1p)  + & 
    1856                &                                 logchl_balinc_ersem_n1p(:,:,:) 
    1857             trn(:,:,:,jp_fabm_m1+jp_fabm_n3n)  = trn(:,:,:,jp_fabm_m1+jp_fabm_n3n)  + & 
    1858                &                                 logchl_balinc_ersem_n3n(:,:,:) 
    1859             trn(:,:,:,jp_fabm_m1+jp_fabm_n4n)  = trn(:,:,:,jp_fabm_m1+jp_fabm_n4n)  + & 
    1860                &                                 logchl_balinc_ersem_n4n(:,:,:) 
    1861             trn(:,:,:,jp_fabm_m1+jp_fabm_n5s)  = trn(:,:,:,jp_fabm_m1+jp_fabm_n5s)  + & 
    1862                &                                 logchl_balinc_ersem_n5s(:,:,:) 
    1863  
    1864             trb(:,:,:,jp_fabm_m1+jp_fabm_chl1) = trn(:,:,:,jp_fabm_m1+jp_fabm_chl1) 
    1865             trb(:,:,:,jp_fabm_m1+jp_fabm_chl2) = trn(:,:,:,jp_fabm_m1+jp_fabm_chl2) 
    1866             trb(:,:,:,jp_fabm_m1+jp_fabm_chl3) = trn(:,:,:,jp_fabm_m1+jp_fabm_chl3) 
    1867             trb(:,:,:,jp_fabm_m1+jp_fabm_chl4) = trn(:,:,:,jp_fabm_m1+jp_fabm_chl4) 
    1868             trb(:,:,:,jp_fabm_m1+jp_fabm_p1c)  = trn(:,:,:,jp_fabm_m1+jp_fabm_p1c) 
    1869             trb(:,:,:,jp_fabm_m1+jp_fabm_p1n)  = trn(:,:,:,jp_fabm_m1+jp_fabm_p1n) 
    1870             trb(:,:,:,jp_fabm_m1+jp_fabm_p1p)  = trn(:,:,:,jp_fabm_m1+jp_fabm_p1p) 
    1871             trb(:,:,:,jp_fabm_m1+jp_fabm_p1s)  = trn(:,:,:,jp_fabm_m1+jp_fabm_p1s) 
    1872             trb(:,:,:,jp_fabm_m1+jp_fabm_p2c)  = trn(:,:,:,jp_fabm_m1+jp_fabm_p2c) 
    1873             trb(:,:,:,jp_fabm_m1+jp_fabm_p2n)  = trn(:,:,:,jp_fabm_m1+jp_fabm_p2n) 
    1874             trb(:,:,:,jp_fabm_m1+jp_fabm_p2p)  = trn(:,:,:,jp_fabm_m1+jp_fabm_p2p) 
    1875             trb(:,:,:,jp_fabm_m1+jp_fabm_p3c)  = trn(:,:,:,jp_fabm_m1+jp_fabm_p3c) 
    1876             trb(:,:,:,jp_fabm_m1+jp_fabm_p3n)  = trn(:,:,:,jp_fabm_m1+jp_fabm_p3n) 
    1877             trb(:,:,:,jp_fabm_m1+jp_fabm_p3p)  = trn(:,:,:,jp_fabm_m1+jp_fabm_p3p) 
    1878             trb(:,:,:,jp_fabm_m1+jp_fabm_p4c)  = trn(:,:,:,jp_fabm_m1+jp_fabm_p4c) 
    1879             trb(:,:,:,jp_fabm_m1+jp_fabm_p4n)  = trn(:,:,:,jp_fabm_m1+jp_fabm_p4n) 
    1880             trb(:,:,:,jp_fabm_m1+jp_fabm_p4p)  = trn(:,:,:,jp_fabm_m1+jp_fabm_p4p) 
    1881             trb(:,:,:,jp_fabm_m1+jp_fabm_z4c)  = trn(:,:,:,jp_fabm_m1+jp_fabm_z4c) 
    1882             trb(:,:,:,jp_fabm_m1+jp_fabm_z5c)  = trn(:,:,:,jp_fabm_m1+jp_fabm_z5c) 
    1883             trb(:,:,:,jp_fabm_m1+jp_fabm_z5n)  = trn(:,:,:,jp_fabm_m1+jp_fabm_z5n) 
    1884             trb(:,:,:,jp_fabm_m1+jp_fabm_z5p)  = trn(:,:,:,jp_fabm_m1+jp_fabm_z5p) 
    1885             trb(:,:,:,jp_fabm_m1+jp_fabm_z6c)  = trn(:,:,:,jp_fabm_m1+jp_fabm_z6c) 
    1886             trb(:,:,:,jp_fabm_m1+jp_fabm_z6n)  = trn(:,:,:,jp_fabm_m1+jp_fabm_z6n) 
    1887             trb(:,:,:,jp_fabm_m1+jp_fabm_z6p)  = trn(:,:,:,jp_fabm_m1+jp_fabm_z6p) 
    1888             trb(:,:,:,jp_fabm_m1+jp_fabm_n1p)  = trn(:,:,:,jp_fabm_m1+jp_fabm_n1p) 
    1889             trb(:,:,:,jp_fabm_m1+jp_fabm_n3n)  = trn(:,:,:,jp_fabm_m1+jp_fabm_n3n) 
    1890             trb(:,:,:,jp_fabm_m1+jp_fabm_n4n)  = trn(:,:,:,jp_fabm_m1+jp_fabm_n4n) 
    1891             trb(:,:,:,jp_fabm_m1+jp_fabm_n5s)  = trn(:,:,:,jp_fabm_m1+jp_fabm_n5s) 
    1892  
     1490            trn(:,:,:,jp_fabm0:jp_fabm1) = trn(:,:,:,jp_fabm0:jp_fabm1) + & 
     1491               &                           logchl_balinc(:,:,:,jp_fabm0:jp_fabm1) 
     1492            trb(:,:,:,jp_fabm0:jp_fabm1) = trn(:,:,:,jp_fabm0:jp_fabm1) 
    18931493#elif defined key_medusa && defined key_foam_medusa 
    18941494            ! Initialize the now fields with the background + increment 
     
    18961496            CALL ctl_warn( ' Doing direct initialisation of MEDUSA with chlorophyll assimilation', & 
    18971497               &           ' Background state is taken from model rather than background file' ) 
    1898             trn(:,:,:,jpchn) = trn(:,:,:,jpchn) + logchl_balinc_medusa_chn(:,:,:) 
    1899             trn(:,:,:,jpchd) = trn(:,:,:,jpchd) + logchl_balinc_medusa_chd(:,:,:) 
    1900             trn(:,:,:,jpphn) = trn(:,:,:,jpphn) + logchl_balinc_medusa_phn(:,:,:) 
    1901             trn(:,:,:,jpphd) = trn(:,:,:,jpphd) + logchl_balinc_medusa_phd(:,:,:) 
    1902             trn(:,:,:,jpzmi) = trn(:,:,:,jpzmi) + logchl_balinc_medusa_zmi(:,:,:) 
    1903             trn(:,:,:,jpzme) = trn(:,:,:,jpzme) + logchl_balinc_medusa_zme(:,:,:) 
    1904             trn(:,:,:,jpdin) = trn(:,:,:,jpdin) + logchl_balinc_medusa_din(:,:,:) 
    1905             trn(:,:,:,jpsil) = trn(:,:,:,jpsil) + logchl_balinc_medusa_sil(:,:,:) 
    1906             trn(:,:,:,jpfer) = trn(:,:,:,jpfer) + logchl_balinc_medusa_fer(:,:,:) 
    1907             trn(:,:,:,jpdet) = trn(:,:,:,jpdet) + logchl_balinc_medusa_det(:,:,:) 
    1908             trn(:,:,:,jppds) = trn(:,:,:,jppds) + logchl_balinc_medusa_pds(:,:,:) 
    1909             trn(:,:,:,jpdtc) = trn(:,:,:,jpdtc) + logchl_balinc_medusa_dtc(:,:,:) 
    1910             trn(:,:,:,jpdic) = trn(:,:,:,jpdic) + logchl_balinc_medusa_dic(:,:,:) 
    1911             trn(:,:,:,jpalk) = trn(:,:,:,jpalk) + logchl_balinc_medusa_alk(:,:,:) 
    1912             trn(:,:,:,jpoxy) = trn(:,:,:,jpoxy) + logchl_balinc_medusa_oxy(:,:,:) 
    1913  
    1914             trb(:,:,:,jpchn) = trn(:,:,:,jpchn) 
    1915             trb(:,:,:,jpchd) = trn(:,:,:,jpchd) 
    1916             trb(:,:,:,jpphn) = trn(:,:,:,jpphn) 
    1917             trb(:,:,:,jpphd) = trn(:,:,:,jpphd) 
    1918             trb(:,:,:,jpzmi) = trn(:,:,:,jpzmi) 
    1919             trb(:,:,:,jpzme) = trn(:,:,:,jpzme) 
    1920             trb(:,:,:,jpdin) = trn(:,:,:,jpdin) 
    1921             trb(:,:,:,jpsil) = trn(:,:,:,jpsil) 
    1922             trb(:,:,:,jpfer) = trn(:,:,:,jpfer) 
    1923             trb(:,:,:,jpdet) = trn(:,:,:,jpdet) 
    1924             trb(:,:,:,jppds) = trn(:,:,:,jppds) 
    1925             trb(:,:,:,jpdtc) = trn(:,:,:,jpdtc) 
    1926             trb(:,:,:,jpdic) = trn(:,:,:,jpdic) 
    1927             trb(:,:,:,jpalk) = trn(:,:,:,jpalk) 
    1928             trb(:,:,:,jpoxy) = trn(:,:,:,jpoxy) 
    1929  
     1498            trn(:,:,:,jp_msa0:jp_msa1) = trn(:,:,:,jp_msa0:jp_msa1) + & 
     1499               &                         logchl_balinc(:,:,:,jp_msa0:jp_msa1) 
     1500            trb(:,:,:,jp_msa0:jp_msa1) = trn(:,:,:,jp_msa0:jp_msa1) 
    19301501#elif defined key_hadocc 
    19311502            ! Initialize the now fields with the background + increment 
     
    19331504            CALL ctl_warn( ' Doing direct initialisation of HadOCC with chlorophyll assimilation', & 
    19341505               &           ' Background state is taken from model rather than background file' ) 
    1935             trn(:,:,:,jp_had_nut) = trn(:,:,:,jp_had_nut) + logchl_balinc_hadocc_nut(:,:,:) 
    1936             trn(:,:,:,jp_had_phy) = trn(:,:,:,jp_had_phy) + logchl_balinc_hadocc_phy(:,:,:) 
    1937             trn(:,:,:,jp_had_zoo) = trn(:,:,:,jp_had_zoo) + logchl_balinc_hadocc_zoo(:,:,:) 
    1938             trn(:,:,:,jp_had_pdn) = trn(:,:,:,jp_had_pdn) + logchl_balinc_hadocc_det(:,:,:) 
    1939             trn(:,:,:,jp_had_dic) = trn(:,:,:,jp_had_dic) + logchl_balinc_hadocc_dic(:,:,:) 
    1940             trn(:,:,:,jp_had_alk) = trn(:,:,:,jp_had_alk) + logchl_balinc_hadocc_alk(:,:,:) 
    1941  
    1942             trb(:,:,:,jp_had_nut) = trn(:,:,:,jp_had_nut) 
    1943             trb(:,:,:,jp_had_phy) = trn(:,:,:,jp_had_phy) 
    1944             trb(:,:,:,jp_had_zoo) = trn(:,:,:,jp_had_zoo) 
    1945             trb(:,:,:,jp_had_pdn) = trn(:,:,:,jp_had_pdn) 
    1946             trb(:,:,:,jp_had_dic) = trn(:,:,:,jp_had_dic) 
    1947             trb(:,:,:,jp_had_alk) = trn(:,:,:,jp_had_alk) 
     1506            trn(:,:,:,jp_had0:jp_had1) = trn(:,:,:,jp_had0:jp_had1) + & 
     1507               &                         logchl_balinc(:,:,:,jp_had0:jp_had1) 
     1508            trb(:,:,:,jp_had0:jp_had1) = trn(:,:,:,jp_had0:jp_had1) 
    19481509#endif 
    19491510  
  • branches/UKMO/dev_r5518_v3.4_asm_nemovar_community_bgc_ersem/NEMOGCM/NEMO/OPA_SRC/ASM/asmlogchlbal_ersem.F90

    r7045 r7716  
    2020   USE trc,         ONLY: trn            ! ERSEM state variables 
    2121   USE par_fabm                          ! FABM parameters 
     22   USE par_trc,     ONLY: jptra          ! Tracer parameters 
    2223 
    2324   IMPLICIT NONE 
     
    3031CONTAINS 
    3132 
    32    SUBROUTINE asm_logchl_bal_ersem( ln_logchlpftinc,          npfts,                    & 
    33       &                             mld_choice_bgc,           logchl_bkginc,            &       
    34       &                             logchl_balinc_ersem_chl1, logchl_balinc_ersem_chl2, & 
    35       &                             logchl_balinc_ersem_chl3, logchl_balinc_ersem_chl4, & 
    36       &                             logchl_balinc_ersem_p1c,  logchl_balinc_ersem_p1n,  & 
    37       &                             logchl_balinc_ersem_p1p,  logchl_balinc_ersem_p1s,  & 
    38       &                             logchl_balinc_ersem_p2c,  logchl_balinc_ersem_p2n,  & 
    39       &                             logchl_balinc_ersem_p2p,  logchl_balinc_ersem_p3c,  & 
    40       &                             logchl_balinc_ersem_p3n,  logchl_balinc_ersem_p3p,  & 
    41       &                             logchl_balinc_ersem_p4c,  logchl_balinc_ersem_p4n,  & 
    42       &                             logchl_balinc_ersem_p4p,  logchl_balinc_ersem_z4c,  & 
    43       &                             logchl_balinc_ersem_z5c,  logchl_balinc_ersem_z5n,  & 
    44       &                             logchl_balinc_ersem_z5p,  logchl_balinc_ersem_z6c,  & 
    45       &                             logchl_balinc_ersem_z6n,  logchl_balinc_ersem_z6p,  & 
    46       &                             logchl_balinc_ersem_n1p,  logchl_balinc_ersem_n3n,  & 
    47       &                             logchl_balinc_ersem_n4n,  logchl_balinc_ersem_n5s ) 
     33   SUBROUTINE asm_logchl_bal_ersem( ln_logchlpftinc, npfts, mld_choice_bgc, & 
     34      &                             logchl_bkginc, logchl_balinc ) 
    4835      !!--------------------------------------------------------------------------- 
    4936      !!                    ***  ROUTINE asm_logchl_bal_ersem  *** 
     
    5845      !!--------------------------------------------------------------------------- 
    5946      !! 
    60       LOGICAL,  INTENT(in   )                           :: ln_logchlpftinc 
    61       INTEGER,  INTENT(in   )                           :: npfts 
    62       INTEGER,  INTENT(in   )                           :: mld_choice_bgc 
    63       REAL(wp), INTENT(in   ), DIMENSION(jpi,jpj,npfts) :: logchl_bkginc 
    64       REAL(wp), INTENT(  out), DIMENSION(jpi,jpj,jpk)   :: logchl_balinc_ersem_chl1 
    65       REAL(wp), INTENT(  out), DIMENSION(jpi,jpj,jpk)   :: logchl_balinc_ersem_chl2 
    66       REAL(wp), INTENT(  out), DIMENSION(jpi,jpj,jpk)   :: logchl_balinc_ersem_chl3 
    67       REAL(wp), INTENT(  out), DIMENSION(jpi,jpj,jpk)   :: logchl_balinc_ersem_chl4 
    68       REAL(wp), INTENT(  out), DIMENSION(jpi,jpj,jpk)   :: logchl_balinc_ersem_p1c 
    69       REAL(wp), INTENT(  out), DIMENSION(jpi,jpj,jpk)   :: logchl_balinc_ersem_p1n 
    70       REAL(wp), INTENT(  out), DIMENSION(jpi,jpj,jpk)   :: logchl_balinc_ersem_p1p 
    71       REAL(wp), INTENT(  out), DIMENSION(jpi,jpj,jpk)   :: logchl_balinc_ersem_p1s 
    72       REAL(wp), INTENT(  out), DIMENSION(jpi,jpj,jpk)   :: logchl_balinc_ersem_p2c 
    73       REAL(wp), INTENT(  out), DIMENSION(jpi,jpj,jpk)   :: logchl_balinc_ersem_p2n 
    74       REAL(wp), INTENT(  out), DIMENSION(jpi,jpj,jpk)   :: logchl_balinc_ersem_p2p 
    75       REAL(wp), INTENT(  out), DIMENSION(jpi,jpj,jpk)   :: logchl_balinc_ersem_p3c 
    76       REAL(wp), INTENT(  out), DIMENSION(jpi,jpj,jpk)   :: logchl_balinc_ersem_p3n 
    77       REAL(wp), INTENT(  out), DIMENSION(jpi,jpj,jpk)   :: logchl_balinc_ersem_p3p 
    78       REAL(wp), INTENT(  out), DIMENSION(jpi,jpj,jpk)   :: logchl_balinc_ersem_p4c 
    79       REAL(wp), INTENT(  out), DIMENSION(jpi,jpj,jpk)   :: logchl_balinc_ersem_p4n 
    80       REAL(wp), INTENT(  out), DIMENSION(jpi,jpj,jpk)   :: logchl_balinc_ersem_p4p 
    81       REAL(wp), INTENT(  out), DIMENSION(jpi,jpj,jpk)   :: logchl_balinc_ersem_z4c 
    82       REAL(wp), INTENT(  out), DIMENSION(jpi,jpj,jpk)   :: logchl_balinc_ersem_z5c 
    83       REAL(wp), INTENT(  out), DIMENSION(jpi,jpj,jpk)   :: logchl_balinc_ersem_z5n 
    84       REAL(wp), INTENT(  out), DIMENSION(jpi,jpj,jpk)   :: logchl_balinc_ersem_z5p 
    85       REAL(wp), INTENT(  out), DIMENSION(jpi,jpj,jpk)   :: logchl_balinc_ersem_z6c 
    86       REAL(wp), INTENT(  out), DIMENSION(jpi,jpj,jpk)   :: logchl_balinc_ersem_z6n 
    87       REAL(wp), INTENT(  out), DIMENSION(jpi,jpj,jpk)   :: logchl_balinc_ersem_z6p 
    88       REAL(wp), INTENT(  out), DIMENSION(jpi,jpj,jpk)   :: logchl_balinc_ersem_n1p 
    89       REAL(wp), INTENT(  out), DIMENSION(jpi,jpj,jpk)   :: logchl_balinc_ersem_n3n 
    90       REAL(wp), INTENT(  out), DIMENSION(jpi,jpj,jpk)   :: logchl_balinc_ersem_n4n 
    91       REAL(wp), INTENT(  out), DIMENSION(jpi,jpj,jpk)   :: logchl_balinc_ersem_n5s 
     47      LOGICAL,  INTENT(in   )                               :: ln_logchlpftinc 
     48      INTEGER,  INTENT(in   )                               :: npfts 
     49      INTEGER,  INTENT(in   )                               :: mld_choice_bgc 
     50      REAL(wp), INTENT(in   ), DIMENSION(jpi,jpj,npfts)     :: logchl_bkginc 
     51      REAL(wp), INTENT(  out), DIMENSION(jpi,jpj,jpk,jptra) :: logchl_balinc 
    9252      !! 
    9353      INTEGER                      :: ji, jj, jk 
     
    9858 
    9959      ! Split surface logchl incs into surface Chl1-4 incs 
     60      ! 
     61      ! In order to transform logchl incs to chl incs, need to account for the background, 
     62      ! cannot simply do 10^logchl_bkginc. Need to: 
     63      ! 1) Add logchl inc to log10(background) to get log10(analysis) 
     64      ! 2) Take 10^log10(analysis) to get analysis 
     65      ! 3) Subtract background from analysis to get chl incs 
    10066      ! 
    10167      IF ( ln_logchlpftinc ) THEN 
     
    10874         DO jj = 1, jpj 
    10975            DO ji = 1, jpi 
    110                IF ( trn(ji,jj,1,jp_fabm_m1+jp_fabm_chl1) > 0.0 ) THEN 
    111                   logchl_balinc_ersem_chl1(ji,jj,1) = 10**( LOG10( trn(ji,jj,1,jp_fabm_m1+jp_fabm_chl1) + & 
    112                      &                                             tiny ) + logchl_bkginc(ji,jj,1) ) -    & 
    113                      &                                trn(ji,jj,1,jp_fabm_m1+jp_fabm_chl1) - tiny 
    114                ENDIF 
    115                IF ( trn(ji,jj,1,jp_fabm_m1+jp_fabm_chl2) > 0.0 ) THEN 
    116                   logchl_balinc_ersem_chl2(ji,jj,1) = 10**( LOG10( trn(ji,jj,1,jp_fabm_m1+jp_fabm_chl2) + & 
    117                      &                                             tiny ) + logchl_bkginc(ji,jj,2) ) -    & 
    118                      &                                trn(ji,jj,1,jp_fabm_m1+jp_fabm_chl2) - tiny 
    119                ENDIF 
    120                IF ( trn(ji,jj,1,jp_fabm_m1+jp_fabm_chl3) > 0.0 ) THEN 
    121                   logchl_balinc_ersem_chl3(ji,jj,1) = 10**( LOG10( trn(ji,jj,1,jp_fabm_m1+jp_fabm_chl3) + & 
    122                      &                                             tiny ) + logchl_bkginc(ji,jj,3) ) -    & 
    123                      &                                trn(ji,jj,1,jp_fabm_m1+jp_fabm_chl3) - tiny 
    124                ENDIF 
    125                IF ( trn(ji,jj,1,jp_fabm_m1+jp_fabm_chl4) > 0.0 ) THEN 
    126                   logchl_balinc_ersem_chl4(ji,jj,1) = 10**( LOG10( trn(ji,jj,1,jp_fabm_m1+jp_fabm_chl4) + & 
    127                      &                                             tiny ) + logchl_bkginc(ji,jj,4) ) -    & 
    128                      &                                trn(ji,jj,1,jp_fabm_m1+jp_fabm_chl4) - tiny 
     76               IF ( ( trn(ji,jj,1,jp_fabm_m1+jp_fabm_chl1) > 0.0 ) .AND. & 
     77                  & ( trn(ji,jj,1,jp_fabm_m1+jp_fabm_chl2) > 0.0 ) .AND. & 
     78                  & ( trn(ji,jj,1,jp_fabm_m1+jp_fabm_chl3) > 0.0 ) .AND. & 
     79                  & ( trn(ji,jj,1,jp_fabm_m1+jp_fabm_chl4) > 0.0 ) ) THEN 
     80                  logchl_balinc(ji,jj,1,jp_fabm_m1+jp_fabm_chl1) = 10**( LOG10( trn(ji,jj,1,jp_fabm_m1+jp_fabm_chl1) ) + & 
     81                     &                                                   logchl_bkginc(ji,jj,1) )                      - & 
     82                     &                                             trn(ji,jj,1,jp_fabm_m1+jp_fabm_chl1) 
     83                  logchl_balinc(ji,jj,1,jp_fabm_m1+jp_fabm_chl2) = 10**( LOG10( trn(ji,jj,1,jp_fabm_m1+jp_fabm_chl2) ) + & 
     84                     &                                                   logchl_bkginc(ji,jj,2) )                      - & 
     85                     &                                             trn(ji,jj,1,jp_fabm_m1+jp_fabm_chl2) 
     86                  logchl_balinc(ji,jj,1,jp_fabm_m1+jp_fabm_chl3) = 10**( LOG10( trn(ji,jj,1,jp_fabm_m1+jp_fabm_chl3) ) + & 
     87                     &                                                   logchl_bkginc(ji,jj,3) )                      - & 
     88                     &                                             trn(ji,jj,1,jp_fabm_m1+jp_fabm_chl3) 
     89                  logchl_balinc(ji,jj,1,jp_fabm_m1+jp_fabm_chl4) = 10**( LOG10( trn(ji,jj,1,jp_fabm_m1+jp_fabm_chl4) ) + & 
     90                     &                                                   logchl_bkginc(ji,jj,4) )                      - & 
     91                     &                                             trn(ji,jj,1,jp_fabm_m1+jp_fabm_chl4) 
    12992               ENDIF 
    13093            END DO 
     
    13598         ! and split between PFTs according to the existing background ratios 
    13699         ! 
     100         IF ( npfts /= 1 ) THEN 
     101            CALL ctl_stop( 'If assimilating total chlorophyll, nn_asmpfts must be 1' ) 
     102         ENDIF 
    137103         DO jj = 1, jpj 
    138104            DO ji = 1, jpi 
     
    143109                  chl_tot = trn(ji,jj,1,jp_fabm_m1+jp_fabm_chl1) + trn(ji,jj,1,jp_fabm_m1+jp_fabm_chl2) + & 
    144110                     &      trn(ji,jj,1,jp_fabm_m1+jp_fabm_chl3) + trn(ji,jj,1,jp_fabm_m1+jp_fabm_chl4) 
    145                   chl_inc = 10**( LOG10( chl_tot + tiny ) + logchl_bkginc(ji,jj,1) ) - chl_tot - tiny 
    146                   logchl_balinc_ersem_chl1(ji,jj,1) = chl_inc * trn(ji,jj,1,jp_fabm_m1+jp_fabm_chl1) / chl_tot 
    147                   logchl_balinc_ersem_chl2(ji,jj,1) = chl_inc * trn(ji,jj,1,jp_fabm_m1+jp_fabm_chl2) / chl_tot 
    148                   logchl_balinc_ersem_chl3(ji,jj,1) = chl_inc * trn(ji,jj,1,jp_fabm_m1+jp_fabm_chl3) / chl_tot 
    149                   logchl_balinc_ersem_chl4(ji,jj,1) = chl_inc * trn(ji,jj,1,jp_fabm_m1+jp_fabm_chl4) / chl_tot 
     111                  chl_inc = 10**( LOG10( chl_tot ) + logchl_bkginc(ji,jj,1) ) - chl_tot 
     112                  logchl_balinc(ji,jj,1,jp_fabm_m1+jp_fabm_chl1) = chl_inc * trn(ji,jj,1,jp_fabm_m1+jp_fabm_chl1) / chl_tot 
     113                  logchl_balinc(ji,jj,1,jp_fabm_m1+jp_fabm_chl2) = chl_inc * trn(ji,jj,1,jp_fabm_m1+jp_fabm_chl2) / chl_tot 
     114                  logchl_balinc(ji,jj,1,jp_fabm_m1+jp_fabm_chl3) = chl_inc * trn(ji,jj,1,jp_fabm_m1+jp_fabm_chl3) / chl_tot 
     115                  logchl_balinc(ji,jj,1,jp_fabm_m1+jp_fabm_chl4) = chl_inc * trn(ji,jj,1,jp_fabm_m1+jp_fabm_chl4) / chl_tot 
    150116               ENDIF 
    151117            END DO 
     
    194160            ! 
    195161            DO jk = 2, jkmax 
    196                logchl_balinc_ersem_chl1(ji,jj,jk) = logchl_balinc_ersem_chl1(ji,jj,1) 
    197                logchl_balinc_ersem_chl2(ji,jj,jk) = logchl_balinc_ersem_chl2(ji,jj,1) 
    198                logchl_balinc_ersem_chl3(ji,jj,jk) = logchl_balinc_ersem_chl3(ji,jj,1) 
    199                logchl_balinc_ersem_chl4(ji,jj,jk) = logchl_balinc_ersem_chl4(ji,jj,1) 
     162               logchl_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_chl1) = logchl_balinc(ji,jj,1,jp_fabm_m1+jp_fabm_chl1) 
     163               logchl_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_chl2) = logchl_balinc(ji,jj,1,jp_fabm_m1+jp_fabm_chl2) 
     164               logchl_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_chl3) = logchl_balinc(ji,jj,1,jp_fabm_m1+jp_fabm_chl3) 
     165               logchl_balinc(ji,jj,jk,jp_fabm_m1+jp_fabm_chl4) = logchl_balinc(ji,jj,1,jp_fabm_m1+jp_fabm_chl4) 
    200166            END DO 
    201167            ! 
     
    203169      END DO 
    204170 
    205       ! Balance P[1-4][c,n,p,s] 
    206       ! Maintain existing background ratios 
    207       ! 
    208       !DO jk = 1, jpk 
    209       !   DO jj = 1, jpj 
    210       !      DO ji = 1, jpi 
    211       !         IF ( trn(ji,jj,1,jp_fabm_m1+jp_fabm_chl1) > 0.0 ) THEN 
    212       !            logchl_balinc_ersem_p1c(ji,jj,jk) = logchl_balinc_ersem_chl1(ji,jj,jk)   * & 
    213       !               &                                trn(ji,jj,jk,jp_fabm_m1+jp_fabm_p1c) / & 
    214       !               &                                trn(ji,jj,jk,jp_fabm_m1+jp_fabm_chl1) 
    215       !            logchl_balinc_ersem_p1n(ji,jj,jk) = logchl_balinc_ersem_chl1(ji,jj,jk)   * & 
    216       !               &                                trn(ji,jj,jk,jp_fabm_m1+jp_fabm_p1n) / & 
    217       !               &                                trn(ji,jj,jk,jp_fabm_m1+jp_fabm_chl1) 
    218       !            logchl_balinc_ersem_p1p(ji,jj,jk) = logchl_balinc_ersem_chl1(ji,jj,jk)   * & 
    219       !               &                                trn(ji,jj,jk,jp_fabm_m1+jp_fabm_p1p) / & 
    220       !               &                                trn(ji,jj,jk,jp_fabm_m1+jp_fabm_chl1) 
    221       !            logchl_balinc_ersem_p1s(ji,jj,jk) = logchl_balinc_ersem_chl1(ji,jj,jk)   * & 
    222       !               &                                trn(ji,jj,jk,jp_fabm_m1+jp_fabm_p1s) / & 
    223       !               &                                trn(ji,jj,jk,jp_fabm_m1+jp_fabm_chl1) 
    224       !         ENDIF 
    225       !         IF ( trn(ji,jj,1,jp_fabm_m1+jp_fabm_chl2) > 0.0 ) THEN 
    226       !            logchl_balinc_ersem_p2c(ji,jj,jk) = logchl_balinc_ersem_chl2(ji,jj,jk)   * & 
    227       !               &                                trn(ji,jj,jk,jp_fabm_m1+jp_fabm_p2c) / & 
    228       !               &                                trn(ji,jj,jk,jp_fabm_m1+jp_fabm_chl2) 
    229       !            logchl_balinc_ersem_p2n(ji,jj,jk) = logchl_balinc_ersem_chl2(ji,jj,jk)   * & 
    230       !               &                                trn(ji,jj,jk,jp_fabm_m1+jp_fabm_p2n) / & 
    231       !               &                                trn(ji,jj,jk,jp_fabm_m1+jp_fabm_chl2) 
    232       !            logchl_balinc_ersem_p2p(ji,jj,jk) = logchl_balinc_ersem_chl2(ji,jj,jk)   * & 
    233       !               &                                trn(ji,jj,jk,jp_fabm_m1+jp_fabm_p2p) / & 
    234       !               &                                trn(ji,jj,jk,jp_fabm_m1+jp_fabm_chl2) 
    235       !         ENDIF 
    236       !         IF ( trn(ji,jj,1,jp_fabm_m1+jp_fabm_chl3) > 0.0 ) THEN 
    237       !            logchl_balinc_ersem_p3c(ji,jj,jk) = logchl_balinc_ersem_chl3(ji,jj,jk)   * & 
    238       !               &                                trn(ji,jj,jk,jp_fabm_m1+jp_fabm_p3c) / & 
    239       !               &                                trn(ji,jj,jk,jp_fabm_m1+jp_fabm_chl3) 
    240       !            logchl_balinc_ersem_p3n(ji,jj,jk) = logchl_balinc_ersem_chl3(ji,jj,jk)   * & 
    241       !               &                                trn(ji,jj,jk,jp_fabm_m1+jp_fabm_p3n) / & 
    242       !               &                                trn(ji,jj,jk,jp_fabm_m1+jp_fabm_chl3) 
    243       !            logchl_balinc_ersem_p3p(ji,jj,jk) = logchl_balinc_ersem_chl3(ji,jj,jk)   * & 
    244       !               &                                trn(ji,jj,jk,jp_fabm_m1+jp_fabm_p3p) / & 
    245       !               &                                trn(ji,jj,jk,jp_fabm_m1+jp_fabm_chl3) 
    246       !         ENDIF 
    247       !         IF ( trn(ji,jj,1,jp_fabm_m1+jp_fabm_chl4) > 0.0 ) THEN 
    248       !            logchl_balinc_ersem_p4c(ji,jj,jk) = logchl_balinc_ersem_chl4(ji,jj,jk)   * & 
    249       !               &                                trn(ji,jj,jk,jp_fabm_m1+jp_fabm_p4c) / & 
    250       !               &                                trn(ji,jj,jk,jp_fabm_m1+jp_fabm_chl4) 
    251       !            logchl_balinc_ersem_p4n(ji,jj,jk) = logchl_balinc_ersem_chl4(ji,jj,jk)   * & 
    252       !               &                                trn(ji,jj,jk,jp_fabm_m1+jp_fabm_p4n) / & 
    253       !               &                                trn(ji,jj,jk,jp_fabm_m1+jp_fabm_chl4) 
    254       !            logchl_balinc_ersem_p4p(ji,jj,jk) = logchl_balinc_ersem_chl4(ji,jj,jk)   * & 
    255       !               &                                trn(ji,jj,jk,jp_fabm_m1+jp_fabm_p4p) / & 
    256       !               &                                trn(ji,jj,jk,jp_fabm_m1+jp_fabm_chl4) 
    257       !         ENDIF 
    258       !      END DO 
    259       !   END DO 
    260       !END DO 
     171      ! Multivariate balancing forthcoming... 
    261172 
    262173   END SUBROUTINE asm_logchl_bal_ersem 
     
    267178   !!---------------------------------------------------------------------- 
    268179CONTAINS 
    269    SUBROUTINE asm_logchl_bal_ersem( ln_logchlpftinc,          npfts,                    & 
    270       &                             mld_choice_bgc,           logchl_bkginc,            &       
    271       &                             logchl_balinc_ersem_chl1, logchl_balinc_ersem_chl2, & 
    272       &                             logchl_balinc_ersem_chl3, logchl_balinc_ersem_chl4, & 
    273       &                             logchl_balinc_ersem_p1c,  logchl_balinc_ersem_p1n,  & 
    274       &                             logchl_balinc_ersem_p1p,  logchl_balinc_ersem_p1s,  & 
    275       &                             logchl_balinc_ersem_p2c,  logchl_balinc_ersem_p2n,  & 
    276       &                             logchl_balinc_ersem_p2p,  logchl_balinc_ersem_p3c,  & 
    277       &                             logchl_balinc_ersem_p3n,  logchl_balinc_ersem_p3p,  & 
    278       &                             logchl_balinc_ersem_p4c,  logchl_balinc_ersem_p4n,  & 
    279       &                             logchl_balinc_ersem_p4p,  logchl_balinc_ersem_z4c,  & 
    280       &                             logchl_balinc_ersem_z5c,  logchl_balinc_ersem_z5n,  & 
    281       &                             logchl_balinc_ersem_z5p,  logchl_balinc_ersem_z6c,  & 
    282       &                             logchl_balinc_ersem_z6n,  logchl_balinc_ersem_z6p,  & 
    283       &                             logchl_balinc_ersem_n1p,  logchl_balinc_ersem_n3n,  & 
    284       &                             logchl_balinc_ersem_n4n,  logchl_balinc_ersem_n5s ) 
     180   SUBROUTINE asm_logchl_bal_ersem( ln_logchlpftinc, npfts, mld_choice_bgc, & 
     181      &                             logchl_bkginc, logchl_balinc ) 
    285182      LOGICAL :: ln_logchlpftinc 
    286183      INTEGER :: npfts 
    287184      INTEGER :: mld_choice_bgc 
    288185      REAL    :: logchl_bkginc(:,:,:) 
    289       REAL    :: logchl_balinc_ersem_chl1(:,:,:), logchl_balinc_ersem_chl2(:,:,:) 
    290       REAL    :: logchl_balinc_ersem_chl3(:,:,:), logchl_balinc_ersem_chl4(:,:,:) 
    291       REAL    :: logchl_balinc_ersem_p1c(:,:,:),  logchl_balinc_ersem_p1n(:,:,:) 
    292       REAL    :: logchl_balinc_ersem_p1p(:,:,:),  logchl_balinc_ersem_p1s(:,:,:) 
    293       REAL    :: logchl_balinc_ersem_p2c(:,:,:),  logchl_balinc_ersem_p2n(:,:,:) 
    294       REAL    :: logchl_balinc_ersem_p2p(:,:,:),  logchl_balinc_ersem_p3c(:,:,:) 
    295       REAL    :: logchl_balinc_ersem_p3n(:,:,:),  logchl_balinc_ersem_p3p(:,:,:) 
    296       REAL    :: logchl_balinc_ersem_p4c(:,:,:),  logchl_balinc_ersem_p4n(:,:,:) 
    297       REAL    :: logchl_balinc_ersem_p4p(:,:,:),  logchl_balinc_ersem_z4c(:,:,:) 
    298       REAL    :: logchl_balinc_ersem_z5c(:,:,:),  logchl_balinc_ersem_z5n(:,:,:) 
    299       REAL    :: logchl_balinc_ersem_z5p(:,:,:),  logchl_balinc_ersem_z6c(:,:,:) 
    300       REAL    :: logchl_balinc_ersem_z6n(:,:,:),  logchl_balinc_ersem_z6p(:,:,:) 
    301       REAL    :: logchl_balinc_ersem_n1p(:,:,:),  logchl_balinc_ersem_n3n(:,:,:) 
    302       REAL    :: logchl_balinc_ersem_n4n(:,:,:),  logchl_balinc_ersem_n5s ) 
     186      REAL    :: logchl_balinc(:,:,:,:) 
    303187      WRITE(*,*) 'asm_logchl_bal_ersem: You should not have seen this print! error?' 
    304188   END SUBROUTINE asm_logchl_bal_ersem 
Note: See TracChangeset for help on using the changeset viewer.