Changeset 14385 for NEMO/branches/2021/dev_r14383_PISCES_NEWDEV_PISCO/src/TOP/PISCES/trcini_pisces.F90
- Timestamp:
- 2021-02-03T16:03:34+01:00 (3 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2021/dev_r14383_PISCES_NEWDEV_PISCO/src/TOP/PISCES/trcini_pisces.F90
r14086 r14385 3 3 !! *** MODULE trcini_pisces *** 4 4 !! TOP : initialisation of the PISCES biochemical model 5 !! This module is for LOBSTER, PISCES and PISCES-QUOTA 5 6 !!====================================================================== 6 7 !! History : - ! 1988-07 (E. Maier-Reiner) Original code … … 37 38 !! 38 39 !! ** Purpose : Initialisation of the PISCES biochemical model 40 !! Allocation of the dynamic arrays 39 41 !!---------------------------------------------------------------------- 40 42 INTEGER, INTENT(in) :: Kmm ! time level indices 41 43 ! 44 ! Read the PISCES namelist 42 45 CALL trc_nam_pisces 43 46 ! … … 55 58 !! ** Purpose : Initialisation of the PISCES biochemical model 56 59 !!---------------------------------------------------------------------- 57 USE p4zsms ! Main P4Z routine60 USE p4zsms ! Main P4Z routine 58 61 USE p4zche ! Chemical model 59 62 USE p4zsink ! vertical flux of particulate matter due to sinking … … 72 75 USE p4zpoc ! Remineralization of organic particles 73 76 USE p4zligand ! Remineralization of organic ligands 74 USE p5zlim ! Co-limitations of differents nutrients 75 USE p5zprod ! Growth rate of the 2 phyto groups76 USE p5zmicro ! Sources and sinks of microzooplankton 77 USE p5zmeso ! Sources and sinks of mesozooplankton 78 USE p5zmort ! Mortality terms for phytoplankton 77 USE p5zlim ! Co-limitations of differents nutrients (QUOTA) 78 USE p5zprod ! Growth rate of the 3 phyto groups (QUOTA) 79 USE p5zmicro ! Sources and sinks of microzooplankton (QUOTA) 80 USE p5zmeso ! Sources and sinks of mesozooplankton (QUOTA) 81 USE p5zmort ! Mortality terms for phytoplankton (QUOTA) 79 82 ! 80 83 INTEGER, INTENT(in) :: Kmm ! time level indices … … 115 118 ierr = ierr + p4z_lim_alloc() 116 119 IF( ln_p4z ) THEN 120 ! PISCES part 117 121 ierr = ierr + p4z_prod_alloc() 122 ierr = ierr + p4z_meso_alloc() 118 123 ELSE 124 ! PISCES-QUOTA part 119 125 ierr = ierr + p5z_lim_alloc() 120 126 ierr = ierr + p5z_prod_alloc() 127 ierr = ierr + p5z_meso_alloc() 121 128 ENDIF 122 129 ierr = ierr + p4z_rem_alloc() … … 130 137 131 138 ! assign an index in trc arrays for each prognostic variables 139 ! This is based on the information read in the namelist_top 132 140 DO jn = 1, jptra 133 141 cltra = ctrcnm(jn) … … 179 187 ! Set biological ratios 180 188 ! --------------------- 181 rno3 = 16._wp / 122._wp 182 po4r = 1._wp / 122._wp 183 o2nit = 32._wp / 122._wp 184 o2ut = 133._wp / 122._wp 185 rdenit = ( ( o2ut + o2nit ) * 0.80 - rno3 - rno3 * 0.60 ) / rno3 186 rdenita = 3._wp / 5._wp 189 rno3 = 16._wp / 122._wp ! C/N 190 po4r = 1._wp / 122._wp ! C/P 191 o2nit = 32._wp / 122._wp ! O2/C for nitrification 192 o2ut = 133._wp / 122._wp ! O2/C for ammonification 193 rdenit = ( ( o2ut + o2nit ) * 0.80 - rno3 - rno3 * 0.60 ) / rno3 ! Denitrification 194 rdenita = 3._wp / 5._wp ! Denitrification 187 195 IF( ln_p5z ) THEN 188 no3rat3 = no3rat3 / rno3 189 po4rat3 = po4rat3 / po4r 196 no3rat3 = no3rat3 / rno3 ! C/N ratio in zooplankton 197 po4rat3 = po4rat3 / po4r ! C/P ratio in zooplankton 190 198 ENDIF 191 199 … … 249 257 250 258 259 ! Initialization of the different PISCES modules 260 ! Mainly corresponds to the namelist use 261 ! ---------------------------------------------- 251 262 CALL p4z_sink_init ! vertical flux of particulate organic matter 252 263 CALL p4z_opt_init ! Optic: PAR in the water column 253 264 IF( ln_p4z ) THEN 265 ! PISCES part 254 266 CALL p4z_lim_init ! co-limitations by the various nutrients 255 267 CALL p4z_prod_init ! phytoplankton growth rate over the global ocean. 256 268 ELSE 269 ! PISCES-QUOTA part 257 270 CALL p5z_lim_init ! co-limitations by the various nutrients 258 271 CALL p5z_prod_init ! phytoplankton growth rate over the global ocean. … … 265 278 & CALL p4z_ligand_init ! remineralisation of organic ligands 266 279 267 IF( ln_p4z ) THEN 280 IF( ln_p4z ) THEN ! PISCES-std 268 281 CALL p4z_mort_init ! phytoplankton mortality 269 282 CALL p4z_micro_init ! microzooplankton 270 283 CALL p4z_meso_init ! mesozooplankton 271 ELSE 284 ELSE ! PISCES-QUOTA 272 285 CALL p5z_mort_init ! phytoplankton mortality 273 286 CALL p5z_micro_init ! microzooplankton
Note: See TracChangeset
for help on using the changeset viewer.