Changeset 10975 for NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps/src/TOP/PISCES/P4Z/p4zbio.F90
- Timestamp:
- 2019-05-13T18:34:33+02:00 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps/src/TOP/PISCES/P4Z/p4zbio.F90
r10966 r10975 45 45 CONTAINS 46 46 47 SUBROUTINE p4z_bio ( kt, knt, Kbb, Kmm )47 SUBROUTINE p4z_bio ( kt, knt, Kbb, Kmm, Krhs ) 48 48 !!--------------------------------------------------------------------- 49 49 !! *** ROUTINE p4z_bio *** … … 56 56 !!--------------------------------------------------------------------- 57 57 INTEGER, INTENT(in) :: kt, knt 58 INTEGER, INTENT(in) :: Kbb, Kmm ! time level indices58 INTEGER, INTENT(in) :: Kbb, Kmm, Krhs ! time level indices 59 59 ! 60 60 INTEGER :: ji, jj, jk, jn … … 73 73 DO ji = 1, jpi 74 74 !!gm : use nmln and test on jk ... less memory acces 75 IF( gdepw _n(ji,jj,jk+1) > hmld(ji,jj) ) xdiss(ji,jj,jk) = 0.0175 IF( gdepw(ji,jj,jk+1,Kmm) > hmld(ji,jj) ) xdiss(ji,jj,jk) = 0.01 76 76 END DO 77 77 END DO 78 78 END DO 79 79 80 CALL p4z_opt ( kt, knt ) ! Optic: PAR in the water column81 CALL p4z_sink ( kt, knt, Kbb, Kmm ) ! vertical flux of particulate organic matter82 CALL p4z_fechem ( kt, knt ) ! Iron chemistry/scavenging80 CALL p4z_opt ( kt, knt, Kbb, Krhs ) ! Optic: PAR in the water column 81 CALL p4z_sink ( kt, knt, Kbb, Kmm, Krhs ) ! vertical flux of particulate organic matter 82 CALL p4z_fechem ( kt, knt, Kbb, Kmm, Krhs ) ! Iron chemistry/scavenging 83 83 ! 84 84 IF( ln_p4z ) THEN 85 CALL p4z_lim ( kt, knt ) ! co-limitations by the various nutrients86 CALL p4z_prod ( kt, knt ) ! phytoplankton growth rate over the global ocean.87 ! ! (for each element : C, Si, Fe, Chl )88 CALL p4z_mort ( kt 89 ! ! zooplankton sources/sinks routines90 CALL p4z_micro( kt, knt )! microzooplankton91 CALL p4z_meso ( kt, knt )! mesozooplankton85 CALL p4z_lim ( kt, knt, Kbb, Kmm ) ! co-limitations by the various nutrients 86 CALL p4z_prod ( kt, knt, Kbb, Kmm, Krhs ) ! phytoplankton growth rate over the global ocean. 87 ! ! (for each element : C, Si, Fe, Chl ) 88 CALL p4z_mort ( kt, Kbb, Krhs ) ! phytoplankton mortality 89 ! ! zooplankton sources/sinks routines 90 CALL p4z_micro( kt, knt, Kbb, Krhs ) ! microzooplankton 91 CALL p4z_meso ( kt, knt, Kbb, Krhs ) ! mesozooplankton 92 92 ELSE 93 CALL p5z_lim ( kt, knt ) ! co-limitations by the various nutrients94 CALL p5z_prod ( kt, knt ) ! phytoplankton growth rate over the global ocean.95 ! ! (for each element : C, Si, Fe, Chl )96 CALL p5z_mort ( kt ) ! phytoplankton mortality97 ! ! zooplankton sources/sinks routines98 CALL p5z_micro( kt, knt ) ! microzooplankton99 CALL p5z_meso ( kt, knt ) ! mesozooplankton93 CALL p5z_lim ( kt, knt, Kbb, Kmm ) ! co-limitations by the various nutrients 94 CALL p5z_prod ( kt, knt, Kbb, Kmm, Krhs ) ! phytoplankton growth rate over the global ocean. 95 ! ! (for each element : C, Si, Fe, Chl ) 96 CALL p5z_mort ( kt, Kbb, Krhs ) ! phytoplankton mortality 97 ! ! zooplankton sources/sinks routines 98 CALL p5z_micro( kt, knt, Kbb, Krhs ) ! microzooplankton 99 CALL p5z_meso ( kt, knt, Kbb, Krhs ) ! mesozooplankton 100 100 ENDIF 101 101 ! 102 CALL p4z_agg ( kt, knt ) ! Aggregation of particles103 CALL p4z_rem ( kt, knt ) ! remineralization terms of organic matter+scavenging of Fe104 CALL p4z_poc ( kt, knt ) ! Remineralization of organic particles102 CALL p4z_agg ( kt, knt, Kbb, Krhs ) ! Aggregation of particles 103 CALL p4z_rem ( kt, knt, Kbb, Kmm, Krhs ) ! remineralization terms of organic matter+scavenging of Fe 104 CALL p4z_poc ( kt, knt, Kbb, Kmm, Krhs ) ! Remineralization of organic particles 105 105 ! 106 106 IF( ln_ligand ) & 107 & CALL p4z_ligand( kt, knt )107 & CALL p4z_ligand( kt, knt, Kbb, Krhs ) 108 108 ! ! 109 109 IF(ln_ctl) THEN ! print mean trends (used for debugging) 110 110 WRITE(charout, FMT="('bio ')") 111 111 CALL prt_ctl_trc_info(charout) 112 CALL prt_ctl_trc(tab4d=tr a, mask=tmask, clinfo=ctrcnm)112 CALL prt_ctl_trc(tab4d=tr(:,:,:,:,Krhs), mask=tmask, clinfo=ctrcnm) 113 113 ENDIF 114 114 !
Note: See TracChangeset
for help on using the changeset viewer.