Changeset 9336 for branches/UKMO
- Timestamp:
- 2018-02-19T18:37:36+01:00 (6 years ago)
- Location:
- branches/UKMO/CO6_KD490_amm7_oper_fabm_hem08/NEMOGCM/NEMO/TOP_SRC/FABM
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/CO6_KD490_amm7_oper_fabm_hem08/NEMOGCM/NEMO/TOP_SRC/FABM/par_fabm.F90
r9330 r9336 36 36 37 37 INTEGER, PUBLIC :: jp_fabmdia_chltot, jp_fabmdia_netpp, jp_fabmdia_phytot 38 INTEGER, PUBLIC :: jp_fabmdia_pgrow, jp_fabmdia_ploss 38 39 39 40 #if defined key_fabm -
branches/UKMO/CO6_KD490_amm7_oper_fabm_hem08/NEMOGCM/NEMO/TOP_SRC/FABM/trcini_fabm.F90
r9330 r9336 121 121 jp_fabm_xeps = fabm_diag_index( 'light_xEPS' ) 122 122 jp_fabmdia_phytot = fabm_diag_index( 'total_phytoplankton_result' ) 123 jp_fabmdia_pgrow = fabm_diag_index( 'p_grow_sum_result' ) 124 jp_fabmdia_ploss = fabm_diag_index( 'p_loss_sum_result' ) 125 123 126 124 127 IF(lwp) WRITE(numout,*) 'DAF: jp_fabm_n5s = ', jp_fabm_n5s -
branches/UKMO/CO6_KD490_amm7_oper_fabm_hem08/NEMOGCM/NEMO/TOP_SRC/FABM/trcsms_fabm.F90
r9330 r9336 136 136 SUBROUTINE asmdiags_fabm 137 137 INTEGER :: ji,jj,jk,jkmax 138 REAL(wp), DIMENSION(jpi,jpj,jpk) :: pgrow_3d, ploss_3d 138 139 139 PGROW_AVG(:,:) = 1.0140 PLOSS_AVG(:,:) = 1.0140 PGROW_AVG(:,:) = 0.0 141 PLOSS_AVG(:,:) = 0.0 141 142 PHYT_AVG(:,:) = 0.0 143 144 pgrow_3d(:,:,:) = fabm_get_bulk_diagnostic_data(model, jp_fabmdia_pgrow) 145 ploss_3d(:,:,:) = fabm_get_bulk_diagnostic_data(model, jp_fabmdia_ploss) 142 146 143 147 DO jj = 1, jpj … … 159 163 & trn(ji,jj,jk,jp_fabm_p3n) + & 160 164 & trn(ji,jj,jk,jp_fabm_p4n) 165 IF ( pgrow_3d(ji,jj,jk) .GT. 0.0 ) THEN 166 PGROW_AVG(ji,jj) = PGROW_AVG(ji,jj) + & 167 & pgrow_3d(ji,jj,jk) 168 ENDIF 169 IF ( ploss_3d(ji,jj,jk) .GT. 0.0 ) THEN 170 PLOSS_AVG(ji,jj) = PLOSS_AVG(ji,jj) + & 171 & ploss_3d(ji,jj,jk) 172 ENDIF 161 173 END DO 162 174 163 PHYT_AVG(ji,jj) = PHYT_AVG(ji,jj) / REAL(jkmax) 175 PHYT_AVG(ji,jj) = PHYT_AVG(ji,jj) / REAL(jkmax) 176 PGROW_AVG(ji,jj) = PGROW_AVG(ji,jj) / REAL(jkmax) 177 PLOSS_AVG(ji,jj) = PLOSS_AVG(ji,jj) / REAL(jkmax) 164 178 165 179 IF ( hmld_tref(ji,jj) .GT. MLD_MAX(ji,jj) ) THEN
Note: See TracChangeset
for help on using the changeset viewer.