Changeset 10141 for branches/UKMO/dev_r5518_GO6_package_asm_3d_bgc_v3/NEMOGCM/NEMO/TOP_SRC/MEDUSA/bio_medusa_fin.F90
- Timestamp:
- 2018-09-19T17:14:46+02:00 (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/dev_r5518_GO6_package_asm_3d_bgc_v3/NEMOGCM/NEMO/TOP_SRC/MEDUSA/bio_medusa_fin.F90
r10055 r10141 46 46 f3_co3, f3_h2co3, f3_hco3, & 47 47 f3_omarg, f3_omcal, f3_pH, & 48 # if defined key_foam_medusa 49 mld_max, pgrow_avg, & 48 ln_foam_medusa, mld_max, pgrow_avg, & 50 49 ploss_avg, phyt_avg, & 51 # endif52 50 za_sed_c, za_sed_ca, za_sed_fe, & 53 51 za_sed_n, za_sed_si, & … … 59 57 USE trc, ONLY: med_diag, nittrc000, trn 60 58 USE trcnam_trp, ONLY: ln_trcadv_cen2, ln_trcadv_tvd 61 # if defined key_foam_medusa62 59 USE zdfmxl, ONLY: hmld 63 # endif64 60 65 61 !! time (integer timestep) … … 120 116 # endif 121 117 122 # if defined key_foam_medusa 123 !!---------------------------------------------------------------------- 124 !! Diagnostics required for ocean colour assimilation: 125 !! Mixed layer average phytoplankton growth, loss and concentration 126 !! Maximum mixed layer depth 127 !!---------------------------------------------------------------------- 128 !! 129 DO jj = 2,jpjm1 130 DO ji = 2,jpim1 131 IF ( hmld(ji,jj) .GT. 0.0 ) THEN 132 pgrow_avg(ji,jj) = pgrow_avg(ji,jj) / hmld(ji,jj) 133 ploss_avg(ji,jj) = ploss_avg(ji,jj) / hmld(ji,jj) 134 phyt_avg(ji,jj) = phyt_avg(ji,jj) / hmld(ji,jj) 135 IF ( hmld(ji,jj) .GT. mld_max(ji,jj) ) THEN 136 mld_max(ji,jj) = hmld(ji,jj) 137 ENDIF 138 ENDIF 118 IF ( ln_foam_medusa ) THEN 119 !!---------------------------------------------------------------------- 120 !! Diagnostics required for ocean colour assimilation: 121 !! Mixed layer average phytoplankton growth, loss and concentration 122 !! Maximum mixed layer depth 123 !!---------------------------------------------------------------------- 124 !! 125 DO jj = 2,jpjm1 126 DO ji = 2,jpim1 127 IF ( hmld(ji,jj) .GT. 0.0 ) THEN 128 pgrow_avg(ji,jj) = pgrow_avg(ji,jj) / hmld(ji,jj) 129 ploss_avg(ji,jj) = ploss_avg(ji,jj) / hmld(ji,jj) 130 phyt_avg(ji,jj) = phyt_avg(ji,jj) / hmld(ji,jj) 131 IF ( hmld(ji,jj) .GT. mld_max(ji,jj) ) THEN 132 mld_max(ji,jj) = hmld(ji,jj) 133 ENDIF 134 ENDIF 135 END DO 139 136 END DO 140 END DO 141 # endif 137 ENDIF 142 138 143 139 # if defined key_debug_medusa
Note: See TracChangeset
for help on using the changeset viewer.