Changeset 7716 for branches/UKMO/dev_r5518_v3.4_asm_nemovar_community_bgc_ersem/NEMOGCM/NEMO/OPA_SRC/ASM/asminc.F90
- Timestamp:
- 2017-02-22T16:11:16+01:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/dev_r5518_v3.4_asm_nemovar_community_bgc_ersem/NEMOGCM/NEMO/OPA_SRC/ASM/asminc.F90
r7472 r7716 53 53 USE bdy_oce, ONLY: bdytmask 54 54 #endif 55 #if defined key_top 56 USE trc, ONLY: & 57 & trn, & 58 & trb 59 USE par_trc, ONLY: & 60 & jptra 61 #endif 55 62 #if defined key_fabm 56 63 USE asmlogchlbal_ersem, ONLY: & 57 64 & asm_logchl_bal_ersem 58 USE trc, ONLY: &59 & trn, &60 & trb61 65 USE par_fabm 62 66 #elif defined key_medusa && defined key_foam_medusa 63 67 USE asmlogchlbal_medusa, ONLY: & 64 68 & asm_logchl_bal_medusa 65 USE trc, ONLY: &66 & trn, &67 & trb68 69 USE par_medusa 69 70 #elif defined key_hadocc 70 71 USE asmlogchlbal_hadocc, ONLY: & 71 72 & asm_logchl_bal_hadocc 72 USE trc, ONLY: &73 & trn, &74 & trb75 73 USE par_hadocc 76 74 #endif … … 130 128 131 129 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 184 132 #endif 185 133 … … 544 492 #endif 545 493 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)) 551 495 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 653 499 #endif 654 500 ENDIF … … 1561 1407 1562 1408 #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 ) 1579 1411 #elif defined key_medusa && defined key_foam_medusa 1580 1412 !CALL asm_logchl_bal_medusa() … … 1611 1443 ENDIF 1612 1444 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 1614 1447 #if defined key_fabm 1615 1448 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 1730 1453 END DO 1731 1732 1454 #elif defined key_medusa && defined key_foam_medusa 1733 1455 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 1766 1460 END DO 1767 1768 1461 #elif defined key_hadocc 1769 1462 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 1784 1467 END DO 1785 1468 #endif … … 1805 1488 CALL ctl_warn( ' Doing direct initialisation of ERSEM with chlorophyll assimilation', & 1806 1489 & ' 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) 1893 1493 #elif defined key_medusa && defined key_foam_medusa 1894 1494 ! Initialize the now fields with the background + increment … … 1896 1496 CALL ctl_warn( ' Doing direct initialisation of MEDUSA with chlorophyll assimilation', & 1897 1497 & ' 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) 1930 1501 #elif defined key_hadocc 1931 1502 ! Initialize the now fields with the background + increment … … 1933 1504 CALL ctl_warn( ' Doing direct initialisation of HadOCC with chlorophyll assimilation', & 1934 1505 & ' 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) 1948 1509 #endif 1949 1510
Note: See TracChangeset
for help on using the changeset viewer.
