Changeset 6106 for branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM
- Timestamp:
- 2015-12-17T18:09:28+01:00 (9 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC/trcrst.F90
r6029 r6106 200 200 if (lwp) write (numout,'(a,3f15.5)') 'Sediment Ca ', & 201 201 & fq0, fq1, fq2 202 !! 203 !! AXY (07/07/15): read in temporally averaged fields for DMS 204 !! calculations 205 !! 206 IF( iom_varid( knum, 'B_DMS_CHN', ldstop = .FALSE. ) > 0 ) THEN 207 !! YES; in which case read them 208 !! 209 IF(lwp) WRITE(numout,*) ' MEDUSA averaged properties for DMS present - reading in ...' 210 CALL iom_get( knum, jpdom_autoglo, 'B_DMS_CHN', zb_dms_chn(:,:) ) 211 CALL iom_get( knum, jpdom_autoglo, 'N_DMS_CHN', zn_dms_chn(:,:) ) 212 CALL iom_get( knum, jpdom_autoglo, 'B_DMS_CHD', zb_dms_chd(:,:) ) 213 CALL iom_get( knum, jpdom_autoglo, 'N_DMS_CHD', zn_dms_chd(:,:) ) 214 CALL iom_get( knum, jpdom_autoglo, 'B_DMS_MLD', zb_dms_mld(:,:) ) 215 CALL iom_get( knum, jpdom_autoglo, 'N_DMS_MLD', zn_dms_mld(:,:) ) 216 CALL iom_get( knum, jpdom_autoglo, 'B_DMS_QSR', zb_dms_qsr(:,:) ) 217 CALL iom_get( knum, jpdom_autoglo, 'N_DMS_QSR', zn_dms_qsr(:,:) ) 218 CALL iom_get( knum, jpdom_autoglo, 'B_DMS_DIN', zb_dms_din(:,:) ) 219 CALL iom_get( knum, jpdom_autoglo, 'N_DMS_DIN', zn_dms_din(:,:) ) 220 ELSE 221 !! NO; in which case set them to zero 222 !! 223 IF(lwp) WRITE(numout,*) ' MEDUSA averaged properties for DMS absent - setting to zero ...' 224 zb_dms_chn(:,:) = 0.0 !! CHN 225 zn_dms_chn(:,:) = 0.0 226 zb_dms_chd(:,:) = 0.0 !! CHD 227 zn_dms_chd(:,:) = 0.0 228 zb_dms_mld(:,:) = 0.0 !! MLD 229 zn_dms_mld(:,:) = 0.0 230 zb_dms_qsr(:,:) = 0.0 !! QSR 231 zn_dms_qsr(:,:) = 0.0 232 zb_dms_din(:,:) = 0.0 !! DIN 233 zn_dms_din(:,:) = 0.0 234 ENDIF 235 !! 236 !! calculate stats on these fields 237 IF(lwp) WRITE(numout,*) ' MEDUSA averaged properties for DMS stats (min, max, sum) ...' 238 fq0 = MINVAL(zn_dms_chn(:,:)) 239 fq1 = MAXVAL(zn_dms_chn(:,:)) 240 fq2 = SUM(zn_dms_chn(:,:)) 241 if (lwp) write (numout,'(a,3f15.5)') 'DMS, CHN ', fq0, fq1, fq2 242 fq0 = MINVAL(zn_dms_chd(:,:)) 243 fq1 = MAXVAL(zn_dms_chd(:,:)) 244 fq2 = SUM(zn_dms_chd(:,:)) 245 if (lwp) write (numout,'(a,3f15.5)') 'DMS, CHD ', fq0, fq1, fq2 246 fq0 = MINVAL(zn_dms_mld(:,:)) 247 fq1 = MAXVAL(zn_dms_mld(:,:)) 248 fq2 = SUM(zn_dms_mld(:,:)) 249 if (lwp) write (numout,'(a,3f15.5)') 'DMS, MLD ', fq0, fq1, fq2 250 fq0 = MINVAL(zn_dms_qsr(:,:)) 251 fq1 = MAXVAL(zn_dms_qsr(:,:)) 252 fq2 = SUM(zn_dms_qsr(:,:)) 253 if (lwp) write (numout,'(a,3f15.5)') 'DMS, QSR ', fq0, fq1, fq2 254 fq0 = MINVAL(zn_dms_din(:,:)) 255 fq1 = MAXVAL(zn_dms_din(:,:)) 256 fq2 = SUM(zn_dms_din(:,:)) 257 if (lwp) write (numout,'(a,3f15.5)') 'DMS, DIN ', fq0, fq1, fq2 258 202 259 #endif 203 260 … … 281 338 if (lwp) write (numout,'(a,3f15.5)') 'Sediment Ca ', & 282 339 & fq0, fq1, fq2 340 !! 341 !! AXY (07/07/15): write out temporally averaged fields for DMS 342 !! calculations 343 !! 344 IF(lwp) WRITE(numout,*) ' MEDUSA averaged properties for DMS - writing out ...' 345 CALL iom_rstput( kt, kitrst, knum, 'B_DMS_CHN', zb_dms_chn(:,:) ) 346 CALL iom_rstput( kt, kitrst, knum, 'N_DMS_CHN', zn_dms_chn(:,:) ) 347 CALL iom_rstput( kt, kitrst, knum, 'B_DMS_CHD', zb_dms_chd(:,:) ) 348 CALL iom_rstput( kt, kitrst, knum, 'N_DMS_CHD', zn_dms_chd(:,:) ) 349 CALL iom_rstput( kt, kitrst, knum, 'B_DMS_MLD', zb_dms_mld(:,:) ) 350 CALL iom_rstput( kt, kitrst, knum, 'N_DMS_MLD', zn_dms_mld(:,:) ) 351 CALL iom_rstput( kt, kitrst, knum, 'B_DMS_QSR', zb_dms_qsr(:,:) ) 352 CALL iom_rstput( kt, kitrst, knum, 'N_DMS_QSR', zn_dms_qsr(:,:) ) 353 CALL iom_rstput( kt, kitrst, knum, 'B_DMS_DIN', zb_dms_din(:,:) ) 354 CALL iom_rstput( kt, kitrst, knum, 'N_DMS_DIN', zn_dms_din(:,:) ) 355 !! 356 !! calculate stats on these fields 357 IF(lwp) WRITE(numout,*) ' MEDUSA averaged properties for DMS stats (min, max, sum) ...' 358 fq0 = MINVAL(zn_dms_chn(:,:)) 359 fq1 = MAXVAL(zn_dms_chn(:,:)) 360 fq2 = SUM(zn_dms_chn(:,:)) 361 if (lwp) write (numout,'(a,3f15.5)') 'DMS, CHN ', fq0, fq1, fq2 362 fq0 = MINVAL(zn_dms_chd(:,:)) 363 fq1 = MAXVAL(zn_dms_chd(:,:)) 364 fq2 = SUM(zn_dms_chd(:,:)) 365 if (lwp) write (numout,'(a,3f15.5)') 'DMS, CHD ', fq0, fq1, fq2 366 fq0 = MINVAL(zn_dms_mld(:,:)) 367 fq1 = MAXVAL(zn_dms_mld(:,:)) 368 fq2 = SUM(zn_dms_mld(:,:)) 369 if (lwp) write (numout,'(a,3f15.5)') 'DMS, MLD ', fq0, fq1, fq2 370 fq0 = MINVAL(zn_dms_qsr(:,:)) 371 fq1 = MAXVAL(zn_dms_qsr(:,:)) 372 fq2 = SUM(zn_dms_qsr(:,:)) 373 if (lwp) write (numout,'(a,3f15.5)') 'DMS, QSR ', fq0, fq1, fq2 374 fq0 = MINVAL(zn_dms_din(:,:)) 375 fq1 = MAXVAL(zn_dms_din(:,:)) 376 fq2 = SUM(zn_dms_din(:,:)) 377 if (lwp) write (numout,'(a,3f15.5)') 'DMS, DIN ', fq0, fq1, fq2 378 !! 283 379 #endif 284 380
Note: See TracChangeset
for help on using the changeset viewer.