| 1 | CCC $Header$ |
|---|
| 2 | C $Id: |
|---|
| 3 | SUBROUTINE h3cadj |
|---|
| 4 | #if defined key_passivetrc |
|---|
| 5 | #if defined key_trc_hamocc3 |
|---|
| 6 | CCC--------------------------------------------------------------------- |
|---|
| 7 | CCC |
|---|
| 8 | CCC ROUTINE h3cadj |
|---|
| 9 | CCC ****************** |
|---|
| 10 | CCC |
|---|
| 11 | CC |
|---|
| 12 | CC PURPOSE. |
|---|
| 13 | CC -------- |
|---|
| 14 | CC *H3CADJ* CORRECT THE TRACERS CONCENTRATION IN THE WATER |
|---|
| 15 | CC BY THE MASS LOST IN THE SEDIMENT |
|---|
| 16 | CC |
|---|
| 17 | CC METHOD. |
|---|
| 18 | CC ------- |
|---|
| 19 | CC |
|---|
| 20 | CC EXTERNALS. |
|---|
| 21 | CC ---------- |
|---|
| 22 | CC NONE. |
|---|
| 23 | CC |
|---|
| 24 | CC REFERENCE. |
|---|
| 25 | CC ---------- |
|---|
| 26 | CC |
|---|
| 27 | CC MODIFICATIONS: |
|---|
| 28 | CC -------------- |
|---|
| 29 | CC original : 2001 O. Aumont (& EK) |
|---|
| 30 | CC |
|---|
| 31 | CC --------------------------------------------------------------------- |
|---|
| 32 | CC parameters and commons |
|---|
| 33 | CC ====================== |
|---|
| 34 | CDIR$ NOLIST |
|---|
| 35 | USE oce_trc |
|---|
| 36 | USE trp_trc |
|---|
| 37 | USE sms |
|---|
| 38 | IMPLICIT NONE |
|---|
| 39 | CDIR$ LIST |
|---|
| 40 | CC---------------------------------------------------------------------- |
|---|
| 41 | CC local declarations |
|---|
| 42 | CC ================== |
|---|
| 43 | C |
|---|
| 44 | INTEGER ji, jj, jk |
|---|
| 45 | REAL ztal, zsil, zpo4, zvol, alka0, silic1, po4 |
|---|
| 46 | CC---------------------------------------------------------------------- |
|---|
| 47 | C |
|---|
| 48 | C Total mass of tracers that are sequestered in the sediment |
|---|
| 49 | C |
|---|
| 50 | ztal=0. |
|---|
| 51 | zsil=0. |
|---|
| 52 | zpo4=0. |
|---|
| 53 | zvol=0. |
|---|
| 54 | DO jk= 1,jpk-1 |
|---|
| 55 | DO jj=2,jpj-1 |
|---|
| 56 | DO ji=2,jpi-1 |
|---|
| 57 | ztal = ztal + tmask(ji,jj,jk)* |
|---|
| 58 | . trn(ji,jj,jk,jptal)*e1t(ji,jj)*e2t(ji,jj)*e3t(jk) |
|---|
| 59 | zsil = zsil + tmask(ji,jj,jk)* |
|---|
| 60 | . trn(ji,jj,jk,jpsil)*e1t(ji,jj)*e2t(ji,jj)*e3t(jk) |
|---|
| 61 | zpo4 = zpo4 + tmask(ji,jj,jk)* |
|---|
| 62 | . trn(ji,jj,jk,jppo4)*e1t(ji,jj)*e2t(ji,jj)*e3t(jk) |
|---|
| 63 | zvol = zvol + tmask(ji,jj,jk)* |
|---|
| 64 | . e1t(ji,jj)*e2t(ji,jj)*e3t(jk) |
|---|
| 65 | ENDDO |
|---|
| 66 | ENDDO |
|---|
| 67 | ENDDO |
|---|
| 68 | C |
|---|
| 69 | C Correction to the concentration in the water column |
|---|
| 70 | C to conserve the total mass |
|---|
| 71 | C |
|---|
| 72 | alka0 = 2.45e-3 |
|---|
| 73 | silic1 = 85.e-6 |
|---|
| 74 | po4 = 2.1e-6/po4r |
|---|
| 75 | ztal = alka0/(ztal/zvol) |
|---|
| 76 | zsil = silic1/(zsil/zvol) |
|---|
| 77 | zpo4 = po4/(zpo4/zvol) |
|---|
| 78 | C |
|---|
| 79 | trn(:,:,:,jptal) = trn(:,:,:,jptal)*ztal |
|---|
| 80 | trb(:,:,:,jptal) = trb(:,:,:,jptal)*ztal |
|---|
| 81 | trn(:,:,:,jpsil) = trn(:,:,:,jpsil)*zsil |
|---|
| 82 | trb(:,:,:,jpsil) = trb(:,:,:,jpsil)*zsil |
|---|
| 83 | trn(:,:,:,jppo4) = trn(:,:,:,jppo4)*zpo4 |
|---|
| 84 | trb(:,:,:,jppo4) = trb(:,:,:,jppo4)*zpo4 |
|---|
| 85 | C |
|---|
| 86 | #endif |
|---|
| 87 | #endif |
|---|
| 88 | RETURN |
|---|
| 89 | END |
|---|
