Changeset 7045
- Timestamp:
- 2016-10-19T17:46:53+02:00 (8 years ago)
- Location:
- branches/UKMO/dev_r5518_v3.4_asm_nemovar_community_bgc_ersem/NEMOGCM/NEMO/OPA_SRC/ASM
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/dev_r5518_v3.4_asm_nemovar_community_bgc_ersem/NEMOGCM/NEMO/OPA_SRC/ASM/asminc.F90
r6983 r7045 188 188 !: 4) hmld_tref - Temperature criterion (0.2 K change from surface) [T points] 189 189 190 INTEGER :: mld_choice_bgc = 5!: choice of mld criteria to use for physics assimilation190 INTEGER :: mld_choice_bgc = 4 !: choice of mld criteria to use for physics assimilation 191 191 !: 1) hmld - Turbocline/mixing depth [W points] 192 192 !: 2) hmlp - Density criterion (0.01 kg/m^3 change from 10m) [W points] -
branches/UKMO/dev_r5518_v3.4_asm_nemovar_community_bgc_ersem/NEMOGCM/NEMO/OPA_SRC/ASM/asmlogchlbal_ersem.F90
r6983 r7045 206 206 ! Maintain existing background ratios 207 207 ! 208 DO jk = 1, jpk209 DO jj = 1, jpj210 DO ji = 1, jpi211 IF ( trn(ji,jj,1,jp_fabm_m1+jp_fabm_chl1) > 0.0 ) THEN212 logchl_balinc_ersem_p1c(ji,jj,jk) = logchl_balinc_ersem_chl1(ji,jj,jk) * &213 & trn(ji,jj,jk,jp_fabm_m1+jp_fabm_p1c) / &214 & trn(ji,jj,jk,jp_fabm_m1+jp_fabm_chl1)215 logchl_balinc_ersem_p1n(ji,jj,jk) = logchl_balinc_ersem_chl1(ji,jj,jk) * &216 & trn(ji,jj,jk,jp_fabm_m1+jp_fabm_p1n) / &217 & trn(ji,jj,jk,jp_fabm_m1+jp_fabm_chl1)218 logchl_balinc_ersem_p1p(ji,jj,jk) = logchl_balinc_ersem_chl1(ji,jj,jk) * &219 & trn(ji,jj,jk,jp_fabm_m1+jp_fabm_p1p) / &220 & trn(ji,jj,jk,jp_fabm_m1+jp_fabm_chl1)221 logchl_balinc_ersem_p1s(ji,jj,jk) = logchl_balinc_ersem_chl1(ji,jj,jk) * &222 & trn(ji,jj,jk,jp_fabm_m1+jp_fabm_p1s) / &223 & trn(ji,jj,jk,jp_fabm_m1+jp_fabm_chl1)224 ENDIF225 IF ( trn(ji,jj,1,jp_fabm_m1+jp_fabm_chl2) > 0.0 ) THEN226 logchl_balinc_ersem_p2c(ji,jj,jk) = logchl_balinc_ersem_chl2(ji,jj,jk) * &227 & trn(ji,jj,jk,jp_fabm_m1+jp_fabm_p2c) / &228 & trn(ji,jj,jk,jp_fabm_m1+jp_fabm_chl2)229 logchl_balinc_ersem_p2n(ji,jj,jk) = logchl_balinc_ersem_chl2(ji,jj,jk) * &230 & trn(ji,jj,jk,jp_fabm_m1+jp_fabm_p2n) / &231 & trn(ji,jj,jk,jp_fabm_m1+jp_fabm_chl2)232 logchl_balinc_ersem_p2p(ji,jj,jk) = logchl_balinc_ersem_chl2(ji,jj,jk) * &233 & trn(ji,jj,jk,jp_fabm_m1+jp_fabm_p2p) / &234 & trn(ji,jj,jk,jp_fabm_m1+jp_fabm_chl2)235 ENDIF236 IF ( trn(ji,jj,1,jp_fabm_m1+jp_fabm_chl3) > 0.0 ) THEN237 logchl_balinc_ersem_p3c(ji,jj,jk) = logchl_balinc_ersem_chl3(ji,jj,jk) * &238 & trn(ji,jj,jk,jp_fabm_m1+jp_fabm_p3c) / &239 & trn(ji,jj,jk,jp_fabm_m1+jp_fabm_chl3)240 logchl_balinc_ersem_p3n(ji,jj,jk) = logchl_balinc_ersem_chl3(ji,jj,jk) * &241 & trn(ji,jj,jk,jp_fabm_m1+jp_fabm_p3n) / &242 & trn(ji,jj,jk,jp_fabm_m1+jp_fabm_chl3)243 logchl_balinc_ersem_p3p(ji,jj,jk) = logchl_balinc_ersem_chl3(ji,jj,jk) * &244 & trn(ji,jj,jk,jp_fabm_m1+jp_fabm_p3p) / &245 & trn(ji,jj,jk,jp_fabm_m1+jp_fabm_chl3)246 ENDIF247 IF ( trn(ji,jj,1,jp_fabm_m1+jp_fabm_chl4) > 0.0 ) THEN248 logchl_balinc_ersem_p4c(ji,jj,jk) = logchl_balinc_ersem_chl4(ji,jj,jk) * &249 & trn(ji,jj,jk,jp_fabm_m1+jp_fabm_p4c) / &250 & trn(ji,jj,jk,jp_fabm_m1+jp_fabm_chl4)251 logchl_balinc_ersem_p4n(ji,jj,jk) = logchl_balinc_ersem_chl4(ji,jj,jk) * &252 & trn(ji,jj,jk,jp_fabm_m1+jp_fabm_p4n) / &253 & trn(ji,jj,jk,jp_fabm_m1+jp_fabm_chl4)254 logchl_balinc_ersem_p4p(ji,jj,jk) = logchl_balinc_ersem_chl4(ji,jj,jk) * &255 & trn(ji,jj,jk,jp_fabm_m1+jp_fabm_p4p) / &256 & trn(ji,jj,jk,jp_fabm_m1+jp_fabm_chl4)257 ENDIF258 END DO259 END DO260 END DO208 !DO jk = 1, jpk 209 ! DO jj = 1, jpj 210 ! DO ji = 1, jpi 211 ! IF ( trn(ji,jj,1,jp_fabm_m1+jp_fabm_chl1) > 0.0 ) THEN 212 ! logchl_balinc_ersem_p1c(ji,jj,jk) = logchl_balinc_ersem_chl1(ji,jj,jk) * & 213 ! & trn(ji,jj,jk,jp_fabm_m1+jp_fabm_p1c) / & 214 ! & trn(ji,jj,jk,jp_fabm_m1+jp_fabm_chl1) 215 ! logchl_balinc_ersem_p1n(ji,jj,jk) = logchl_balinc_ersem_chl1(ji,jj,jk) * & 216 ! & trn(ji,jj,jk,jp_fabm_m1+jp_fabm_p1n) / & 217 ! & trn(ji,jj,jk,jp_fabm_m1+jp_fabm_chl1) 218 ! logchl_balinc_ersem_p1p(ji,jj,jk) = logchl_balinc_ersem_chl1(ji,jj,jk) * & 219 ! & trn(ji,jj,jk,jp_fabm_m1+jp_fabm_p1p) / & 220 ! & trn(ji,jj,jk,jp_fabm_m1+jp_fabm_chl1) 221 ! logchl_balinc_ersem_p1s(ji,jj,jk) = logchl_balinc_ersem_chl1(ji,jj,jk) * & 222 ! & trn(ji,jj,jk,jp_fabm_m1+jp_fabm_p1s) / & 223 ! & trn(ji,jj,jk,jp_fabm_m1+jp_fabm_chl1) 224 ! ENDIF 225 ! IF ( trn(ji,jj,1,jp_fabm_m1+jp_fabm_chl2) > 0.0 ) THEN 226 ! logchl_balinc_ersem_p2c(ji,jj,jk) = logchl_balinc_ersem_chl2(ji,jj,jk) * & 227 ! & trn(ji,jj,jk,jp_fabm_m1+jp_fabm_p2c) / & 228 ! & trn(ji,jj,jk,jp_fabm_m1+jp_fabm_chl2) 229 ! logchl_balinc_ersem_p2n(ji,jj,jk) = logchl_balinc_ersem_chl2(ji,jj,jk) * & 230 ! & trn(ji,jj,jk,jp_fabm_m1+jp_fabm_p2n) / & 231 ! & trn(ji,jj,jk,jp_fabm_m1+jp_fabm_chl2) 232 ! logchl_balinc_ersem_p2p(ji,jj,jk) = logchl_balinc_ersem_chl2(ji,jj,jk) * & 233 ! & trn(ji,jj,jk,jp_fabm_m1+jp_fabm_p2p) / & 234 ! & trn(ji,jj,jk,jp_fabm_m1+jp_fabm_chl2) 235 ! ENDIF 236 ! IF ( trn(ji,jj,1,jp_fabm_m1+jp_fabm_chl3) > 0.0 ) THEN 237 ! logchl_balinc_ersem_p3c(ji,jj,jk) = logchl_balinc_ersem_chl3(ji,jj,jk) * & 238 ! & trn(ji,jj,jk,jp_fabm_m1+jp_fabm_p3c) / & 239 ! & trn(ji,jj,jk,jp_fabm_m1+jp_fabm_chl3) 240 ! logchl_balinc_ersem_p3n(ji,jj,jk) = logchl_balinc_ersem_chl3(ji,jj,jk) * & 241 ! & trn(ji,jj,jk,jp_fabm_m1+jp_fabm_p3n) / & 242 ! & trn(ji,jj,jk,jp_fabm_m1+jp_fabm_chl3) 243 ! logchl_balinc_ersem_p3p(ji,jj,jk) = logchl_balinc_ersem_chl3(ji,jj,jk) * & 244 ! & trn(ji,jj,jk,jp_fabm_m1+jp_fabm_p3p) / & 245 ! & trn(ji,jj,jk,jp_fabm_m1+jp_fabm_chl3) 246 ! ENDIF 247 ! IF ( trn(ji,jj,1,jp_fabm_m1+jp_fabm_chl4) > 0.0 ) THEN 248 ! logchl_balinc_ersem_p4c(ji,jj,jk) = logchl_balinc_ersem_chl4(ji,jj,jk) * & 249 ! & trn(ji,jj,jk,jp_fabm_m1+jp_fabm_p4c) / & 250 ! & trn(ji,jj,jk,jp_fabm_m1+jp_fabm_chl4) 251 ! logchl_balinc_ersem_p4n(ji,jj,jk) = logchl_balinc_ersem_chl4(ji,jj,jk) * & 252 ! & trn(ji,jj,jk,jp_fabm_m1+jp_fabm_p4n) / & 253 ! & trn(ji,jj,jk,jp_fabm_m1+jp_fabm_chl4) 254 ! logchl_balinc_ersem_p4p(ji,jj,jk) = logchl_balinc_ersem_chl4(ji,jj,jk) * & 255 ! & trn(ji,jj,jk,jp_fabm_m1+jp_fabm_p4p) / & 256 ! & trn(ji,jj,jk,jp_fabm_m1+jp_fabm_chl4) 257 ! ENDIF 258 ! END DO 259 ! END DO 260 !END DO 261 261 262 262 END SUBROUTINE asm_logchl_bal_ersem
Note: See TracChangeset
for help on using the changeset viewer.