- Timestamp:
- 2016-08-03T16:59:29+02:00 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/CNRS/dev_r6270_PISCES_QUOTA/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zagg.F90
r6453 r6841 58 58 IF( nn_timing == 1 ) CALL timing_start('p4z_agg') 59 59 60 ! Initialization of some global variables61 ! ---------------------------------------62 prodpoc(:,:,:) = 0.63 conspoc(:,:,:) = 0.64 prodgoc(:,:,:) = 0.65 consgoc(:,:,:) = 0.66 60 ! 67 61 ! Exchange between organic matter compartments due to coagulation/disaggregation … … 100 94 zaggdoc2 = ( 3.53E3 * zfact + 0.1 * zstep ) * trb(ji,jj,jk,jpgoc) * 0.3 * trb(ji,jj,jk,jpdoc) 101 95 ! tranfer of DOC to POC due to brownian motion 102 zaggdoc3 = ( 5095. * 0. * trb(ji,jj,jk,jppoc) +114. * 0.3 * trb(ji,jj,jk,jpdoc) ) *zstep * 0.3 * trb(ji,jj,jk,jpdoc)96 zaggdoc3 = ( 114. * 0.3 * trb(ji,jj,jk,jpdoc) ) *zstep * 0.3 * trb(ji,jj,jk,jpdoc) 103 97 104 98 ! Update the trends … … 109 103 tra(ji,jj,jk,jpdoc) = tra(ji,jj,jk,jpdoc) - zaggdoc - zaggdoc2 - zaggdoc3 110 104 ! 111 prodpoc(ji,jj,jk) = prodpoc(ji,jj,jk) + zaggdoc + zaggdoc3 112 conspoc(ji,jj,jk) = conspoc(ji,jj,jk) - zagg 105 conspoc(ji,jj,jk) = conspoc(ji,jj,jk) - zagg + zaggdoc + zaggdoc3 113 106 prodgoc(ji,jj,jk) = prodgoc(ji,jj,jk) + zagg + zaggdoc2 114 107 !
Note: See TracChangeset
for help on using the changeset viewer.