[260] | 1 | CCC $Header$ |
---|
| 2 | C $Id: |
---|
[186] | 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 |
---|