- Timestamp:
- 2010-09-17T14:35:46+02:00 (14 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/DEV_r2006_merge_TRA_TRC/NEMO/TOP_SRC/PISCES/p4zlim.F90
r2082 r2104 23 23 24 24 PUBLIC p4z_lim 25 PUBLIC p4z_lim_init 25 26 26 27 !! * Shared module variables … … 50 51 CONTAINS 51 52 52 SUBROUTINE p4z_lim( kt , jnt)53 SUBROUTINE p4z_lim( kt ) 53 54 !!--------------------------------------------------------------------- 54 55 !! *** ROUTINE p4z_lim *** … … 59 60 !! ** Method : - ??? 60 61 !!--------------------------------------------------------------------- 61 INTEGER, INTENT(in) :: kt, jnt ! ocean time step62 INTEGER, INTENT(in) :: kt 62 63 INTEGER :: ji, jj, jk 63 64 REAL(wp) :: zlim1, zlim2, zlim3, zlim4, zno3, zferlim … … 67 68 68 69 69 IF( ( kt * jnt ) == nittrc000 ) CALL p4z_lim_init ! Initialization (first time-step only) 70 71 72 ! Tuning of the iron concentration to a minimum 73 ! level that is set to the detection limit 74 ! ------------------------------------- 70 ! Tuning of the iron concentration to a minimum 71 ! level that is set to the detection limit 72 ! ------------------------------------- 75 73 76 74 DO jk = 1, jpkm1 … … 85 83 END DO 86 84 87 ! Computation of a variable Ks for iron on diatoms 88 ! taking into account that increasing biomass is 89 ! made of generally bigger cells 90 ! ------------------------------------------------ 85 ! Computation of a variable Ks for iron on diatoms taking into account 86 ! that increasing biomass is made of generally bigger cells 87 ! ------------------------------------------------ 91 88 92 89 DO jk = 1, jpkm1 … … 107 104 END DO 108 105 109 DO jk = 1, jpkm1 110 DO jj = 1, jpj 111 DO ji = 1, jpi 112 113 ! Michaelis-Menten Limitation term for nutrients 114 ! Small flagellates 115 ! ----------------------------------------------- 106 ! Michaelis-Menten Limitation term for nutrients Small flagellates 107 ! ----------------------------------------------- 108 DO jk = 1, jpkm1 109 DO jj = 1, jpj 110 DO ji = 1, jpi 116 111 zdenom = 1. / & 117 112 & ( conc0 * concnnh4 + concnnh4 * trn(ji,jj,jk,jpno3) + conc0 * trn(ji,jj,jk,jpnh4) ) … … 132 127 END DO 133 128 134 DO jk = 1, jpkm1 135 DO jj = 1, jpj 136 DO ji = 1, jpi 137 138 ! Michaelis-Menten Limitation term for nutrients Diatoms 139 ! ---------------------------------------------- 129 ! Michaelis-Menten Limitation term for nutrients Diatoms 130 ! ---------------------------------------------- 131 DO jk = 1, jpkm1 132 DO jj = 1, jpj 133 DO ji = 1, jpi 140 134 zdenom = 1. / & 141 135 & ( conc1 * concdnh4 + concdnh4 * trn(ji,jj,jk,jpno3) + conc1 * trn(ji,jj,jk,jpnh4) ) … … 181 175 !! 182 176 !! ** Method : Read the nampislim namelist and check the parameters 183 !! called at the first timestep (nit trc000)177 !! called at the first timestep (nit000) 184 178 !! 185 179 !! ** input : Namelist nampislim
Note: See TracChangeset
for help on using the changeset viewer.