[186] | 1 | C $Id$ |
---|
| 2 | CDIR$ LIST |
---|
| 3 | SUBROUTINE h3crem |
---|
| 4 | #if defined key_passivetrc && defined key_trc_hamocc3 |
---|
| 5 | CCC--------------------------------------------------------------------- |
---|
| 6 | CCC |
---|
| 7 | CCC ROUTINE h3crem |
---|
| 8 | CCC ****************** |
---|
| 9 | CCC |
---|
| 10 | CCC |
---|
| 11 | CCC |
---|
| 12 | CCC |
---|
| 13 | CCC PURPOSE. |
---|
| 14 | CCC -------- |
---|
| 15 | CCC *H3CREM* CALCULATES OXYGEN CONSUMPTION DURING REMINERALIZATION |
---|
| 16 | CCC AND LOSS OF POC TO THE ORGANIC SEDIMENT POOL |
---|
| 17 | CCC |
---|
| 18 | CC |
---|
| 19 | CC METHOD. |
---|
| 20 | CC ------- |
---|
| 21 | CC THE OXYGEN AMOUNT AVAILABLE FOR CONSUMPTION DUE TO |
---|
| 22 | CC REMINERALIZATION OF POC IS DETERMINED. |
---|
| 23 | CC THE CORRESPONDING DECREASE OF POC AND THE INCREASES OF |
---|
| 24 | CC INORGANIC CARBON AND PHOSPHATE ARE CALCULATED. ALKALINITY |
---|
| 25 | CC IS REDUCED BY THE NITRATE SET FREE DURING REMINERALIZATION. |
---|
| 26 | CC POC SETTLES DOWN TO THE SEDIMENT FROM THE WATER LAYER |
---|
| 27 | CC DIRECTLY ABOVE THE OCEAN FLOOR (DEPENDS ON BOTTOM DEPTH) |
---|
| 28 | CC AT A RATE TO BE DEFINED (''OCSED''). |
---|
| 29 | CC |
---|
| 30 | CC EXTERNALS. |
---|
| 31 | CC ---------- |
---|
| 32 | CC NONE. |
---|
| 33 | CC |
---|
| 34 | CC REFERENCE. |
---|
| 35 | CC ---------- |
---|
| 36 | CC -. |
---|
| 37 | CC |
---|
| 38 | CC* VARIABLE TYPE PURPOSE. |
---|
| 39 | CC -------- ---- -------- |
---|
| 40 | CC |
---|
| 41 | CC *ZREM* REAL MOLES O2 AVAILABLE FOR REMINERALIZATION |
---|
| 42 | CC (LIMITING INGREDIENT FOR REMINERALIZATION) |
---|
| 43 | CC *ZREMIK* REAL REMINERALIZATION RATE (FRACTION OF POC, |
---|
| 44 | CC THAT IS REMINERALISED DURING ONE TIMESTEP) |
---|
| 45 | CC DUMMY VARIABLE |
---|
| 46 | CC |
---|
| 47 | CC MODIFICATIONS: |
---|
| 48 | CC -------------- |
---|
| 49 | CC original : 1988-07 E. MAIER-REIMER MPI HAMBURG |
---|
| 50 | CC additions : 1998 O. Aumont |
---|
| 51 | CC modifications : 1999 C. Le Quere |
---|
| 52 | CC --------------------------------------------------------------------- |
---|
| 53 | CC parameters and commons |
---|
| 54 | CC ====================== |
---|
| 55 | CDIR$ NOLIST |
---|
| 56 | USE oce_trc |
---|
| 57 | USE trp_trc |
---|
| 58 | USE sms |
---|
| 59 | IMPLICIT NONE |
---|
| 60 | CDIR$ LIST |
---|
| 61 | CC---------------------------------------------------------------------- |
---|
| 62 | CC local declarations |
---|
| 63 | CC ================== |
---|
| 64 | C |
---|
| 65 | INTEGER ji, jj, jk |
---|
| 66 | REAL zremik, zrem |
---|
| 67 | C |
---|
| 68 | C* 3. REMINERALIZATION OF POC AND O2 UTILIZATION |
---|
| 69 | C ---------------------------------------------- |
---|
| 70 | C |
---|
| 71 | C* 3.1 LAYER LOOP |
---|
| 72 | C ---------------- |
---|
| 73 | C |
---|
| 74 | DO jk = 1,jpkm1 |
---|
| 75 | DO jj = 1,jpj |
---|
| 76 | C |
---|
| 77 | C* 3.2 SET FRACTION POC TO BE REMINERALISED IN LAYER K |
---|
| 78 | C ----------------------------------------------------- |
---|
| 79 | C |
---|
| 80 | C* 3.3 MERIDIONAL/ZONAL LOOPS |
---|
| 81 | C ---------------------------- |
---|
| 82 | C |
---|
| 83 | DO ji = 1,jpi |
---|
| 84 | C |
---|
| 85 | zrem = 0. |
---|
| 86 | C |
---|
| 87 | zremik = 50.*spocri*(1.-tmask(ji,jj,jk+1)) |
---|
| 88 | $ +0.002*(tn(ji,jj,jk)+2.) |
---|
| 89 | zremik = 0.002*(tn(ji,jj,jk)+2.) |
---|
| 90 | zremik = zremik/(30.5*3600.*24.)*rfact |
---|
| 91 | # if defined key_off_degrad |
---|
| 92 | & *facvol(ji,jj,jk) |
---|
| 93 | # endif |
---|
| 94 | C |
---|
| 95 | C* 3.5 AMOUNT OF O2 AVAILABLE FOR REMINERALIZATION OF POC |
---|
| 96 | C -------------------------------------------------------- |
---|
| 97 | C |
---|
| 98 | zrem = zremik* min( |
---|
| 99 | & (trn(ji,jj,jk,jpoxy)-oxymin)/o2ut,trn(ji,jj,jk,jppoc)) |
---|
| 100 | C |
---|
| 101 | C |
---|
| 102 | C* 3.6 INCREASE OF PO4 AND SUM((12C)O2) DUE TO RELEASE TO WATER |
---|
| 103 | C COLUMN DURING REMINERALIZATION |
---|
| 104 | C -------------------------------------------------------------- |
---|
| 105 | C |
---|
| 106 | tra(ji,jj,jk,jppo4) = tra(ji,jj,jk,jppo4)+ zrem*rfactr |
---|
| 107 | tra(ji,jj,jk,jpdic) = tra(ji,jj,jk,jpdic)+ zrem*rfactr |
---|
| 108 | C |
---|
| 109 | C* 3.7 DECREASE OF ALKALINITY AND O2 DUE TO REMINERALIZATION |
---|
| 110 | C (ALKALINITY DECREASES DUE TO NITRATE RELEASE) |
---|
| 111 | C ----------------------------------------------------------- |
---|
| 112 | C |
---|
| 113 | tra(ji,jj,jk,jptal) = tra(ji,jj,jk,jptal)- |
---|
| 114 | & zrem*rno3*rfactr |
---|
| 115 | tra(ji,jj,jk,jpoxy) = tra(ji,jj,jk,jpoxy)- |
---|
| 116 | & rfactr*o2ut*zrem |
---|
| 117 | C |
---|
| 118 | C* 3.8 CHANGE IN RARE ISOTOPES OF INORGANIC CARBON AND POC |
---|
| 119 | C AFTER REMINERALIZATION AND DECREASE OF POC 12C |
---|
| 120 | C IN THE WATER COLUMN DUE TO REMINERALIZATION |
---|
| 121 | C ---------------------------------------------------------------------- |
---|
| 122 | C |
---|
| 123 | # if defined key_trc_biohamocc13 |
---|
| 124 | tra(ji,jj,jk,jp13c) = tra(ji,jj,jk,jp13c)+ |
---|
| 125 | & zrem*pdb*rfactr*plafr13 |
---|
| 126 | # endif |
---|
| 127 | tra(ji,jj,jk,jppoc) = tra(ji,jj,jk,jppoc)- |
---|
| 128 | & zrem*rfactr |
---|
| 129 | C |
---|
| 130 | ENDDO |
---|
| 131 | ENDDO |
---|
| 132 | C |
---|
| 133 | C END of slab |
---|
| 134 | C =========== |
---|
| 135 | C |
---|
| 136 | ENDDO |
---|
| 137 | # endif |
---|
| 138 | RETURN |
---|
| 139 | END |
---|
| 140 | |
---|