Changeset 6229 for branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM
- Timestamp:
- 2016-01-11T13:11:30+01:00 (9 years ago)
- Location:
- branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC/trcini.F90
r6164 r6229 255 255 # if defined key_debug_medusa 256 256 WRITE(numout,*) ' litle check : ', ctrcnm(1) 257 CALL trc_rst_stat 257 258 CALL flush(numout) 258 259 # endif -
branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC/trcrst.F90
r6227 r6229 196 196 !! calculate stats on these fields 197 197 IF(lwp) WRITE(numout,*) ' MEDUSA sediment field stats (min, max, sum) ...' 198 fq0 = MINVAL(zn_sed_n(:,:)) 199 fq1 = MAXVAL(zn_sed_n(:,:)) 200 fq2 = SUM(zn_sed_n(:,:)) 201 if (lwp) write (numout,'(a,3f15.5)') 'Sediment N ', & 202 & fq0, fq1, fq2 203 fq0 = MINVAL(zn_sed_fe(:,:)) 204 fq1 = MAXVAL(zn_sed_fe(:,:)) 205 fq2 = SUM(zn_sed_fe(:,:)) 206 if (lwp) write (numout,'(a,3f15.5)') 'Sediment Fe ', & 207 & fq0, fq1, fq2 208 fq0 = MINVAL(zn_sed_si(:,:)) 209 fq1 = MAXVAL(zn_sed_si(:,:)) 210 fq2 = SUM(zn_sed_si(:,:)) 211 if (lwp) write (numout,'(a,3f15.5)') 'Sediment Si ', & 212 & fq0, fq1, fq2 213 fq0 = MINVAL(zn_sed_c(:,:)) 214 fq1 = MAXVAL(zn_sed_c(:,:)) 215 fq2 = SUM(zn_sed_c(:,:)) 216 if (lwp) write (numout,'(a,3f15.5)') 'Sediment C ', & 217 & fq0, fq1, fq2 218 fq0 = MINVAL(zn_sed_ca(:,:)) 219 fq1 = MAXVAL(zn_sed_ca(:,:)) 220 fq2 = SUM(zn_sed_ca(:,:)) 221 if (lwp) write (numout,'(a,3f15.5)') 'Sediment Ca ', & 222 & fq0, fq1, fq2 198 call trc_rst_dia_stat(zn_sed_n(:,:), 'Sediment N') 199 call trc_rst_dia_stat(zn_sed_fe(:,:), 'Sediment Fe') 200 call trc_rst_dia_stat(zn_sed_si(:,:), 'Sediment Si') 201 call trc_rst_dia_stat(zn_sed_c(:,:), 'Sediment C') 202 call trc_rst_dia_stat(zn_sed_ca(:,:), 'Sediment Ca') 223 203 !! 224 204 !! AXY (07/07/15): read in temporally averaged fields for DMS … … 257 237 !! calculate stats on these fields 258 238 IF(lwp) WRITE(numout,*) ' MEDUSA averaged properties for DMS stats (min, max, sum) ...' 259 fq0 = MINVAL(zn_dms_chn(:,:)) 260 fq1 = MAXVAL(zn_dms_chn(:,:)) 261 fq2 = SUM(zn_dms_chn(:,:)) 262 if (lwp) write (numout,'(a,3f15.5)') 'DMS, CHN ', fq0, fq1, fq2 263 fq0 = MINVAL(zn_dms_chd(:,:)) 264 fq1 = MAXVAL(zn_dms_chd(:,:)) 265 fq2 = SUM(zn_dms_chd(:,:)) 266 if (lwp) write (numout,'(a,3f15.5)') 'DMS, CHD ', fq0, fq1, fq2 267 fq0 = MINVAL(zn_dms_mld(:,:)) 268 fq1 = MAXVAL(zn_dms_mld(:,:)) 269 fq2 = SUM(zn_dms_mld(:,:)) 270 if (lwp) write (numout,'(a,3f15.5)') 'DMS, MLD ', fq0, fq1, fq2 271 fq0 = MINVAL(zn_dms_qsr(:,:)) 272 fq1 = MAXVAL(zn_dms_qsr(:,:)) 273 fq2 = SUM(zn_dms_qsr(:,:)) 274 if (lwp) write (numout,'(a,3f15.5)') 'DMS, QSR ', fq0, fq1, fq2 275 fq0 = MINVAL(zn_dms_din(:,:)) 276 fq1 = MAXVAL(zn_dms_din(:,:)) 277 fq2 = SUM(zn_dms_din(:,:)) 278 if (lwp) write (numout,'(a,3f15.5)') 'DMS, DIN ', fq0, fq1, fq2 239 call trc_rst_dia_stat(zn_dms_chn(:,:), 'DMS, CHN') 240 call trc_rst_dia_stat(zn_dms_chd(:,:), 'DMS, CHD') 241 call trc_rst_dia_stat(zn_dms_mld(:,:), 'DMS, MLD') 242 call trc_rst_dia_stat(zn_dms_qsr(:,:), 'DMS, QSR') 243 call trc_rst_dia_stat(zn_dms_din(:,:), 'DMS, DIN') 279 244 #endif 280 245 ! … … 298 263 !! 299 264 !! calculate stats on these fields 265 IF(lwp) WRITE(numout,*) ' IDTRA averaged properties stats (min, max, sum) ...' 300 266 call trc_rst_dia_stat(qint_idtra(:,:,1), 'qint_IDTRA') 301 IF(lwp) WRITE(numout,*) ' IDTRA averaged properties stats (min, max, sum) ...'302 fq0 = MINVAL(qint_idtra(:,:,1))303 fq1 = MAXVAL(qint_idtra(:,:,1))304 fq2 = SUM(qint_idtra(:,:,1))305 if (lwp) write (numout,'(a,3f15.5)') 'qint_IDTRA ', fq0, fq1, fq2306 267 #endif 307 268 ! … … 322 283 !! 323 284 !! calculate stats on these fields 285 IF(lwp) WRITE(numout,*) ' CFC averaged properties stats (min, max, sum) ...' 324 286 call trc_rst_dia_stat(qint_cfc(:,:,jl), 'qint_'//ctrcnm(jn)) 325 IF(lwp) WRITE(numout,*) ' CFC averaged properties stats (min, max, sum) ...'326 fq0 = MINVAL(qint_cfc(:,:,jl))327 fq1 = MAXVAL(qint_cfc(:,:,jl))328 fq2 = SUM(qint_cfc(:,:,jl))329 if (lwp) write (numout,'(a,3f15.5)') 'qint_'//ctrcnm(jn), fq0, fq1, fq2330 287 END DO 331 288 #endif … … 384 341 !! calculate stats on these fields 385 342 IF(lwp) WRITE(numout,*) ' MEDUSA sediment field stats (min, max, sum) ...' 386 fq0 = MINVAL(zn_sed_n(:,:)) 387 fq1 = MAXVAL(zn_sed_n(:,:)) 388 fq2 = SUM(zn_sed_n(:,:)) 389 if (lwp) write (numout,'(a,3f15.5)') 'Sediment N ', & 390 & fq0, fq1, fq2 391 fq0 = MINVAL(zn_sed_fe(:,:)) 392 fq1 = MAXVAL(zn_sed_fe(:,:)) 393 fq2 = SUM(zn_sed_fe(:,:)) 394 if (lwp) write (numout,'(a,3f15.5)') 'Sediment Fe ', & 395 & fq0, fq1, fq2 396 fq0 = MINVAL(zn_sed_si(:,:)) 397 fq1 = MAXVAL(zn_sed_si(:,:)) 398 fq2 = SUM(zn_sed_si(:,:)) 399 if (lwp) write (numout,'(a,3f15.5)') 'Sediment Si ', & 400 & fq0, fq1, fq2 401 fq0 = MINVAL(zn_sed_c(:,:)) 402 fq1 = MAXVAL(zn_sed_c(:,:)) 403 fq2 = SUM(zn_sed_c(:,:)) 404 if (lwp) write (numout,'(a,3f15.5)') 'Sediment C ', & 405 & fq0, fq1, fq2 406 fq0 = MINVAL(zn_sed_ca(:,:)) 407 fq1 = MAXVAL(zn_sed_ca(:,:)) 408 fq2 = SUM(zn_sed_ca(:,:)) 409 if (lwp) write (numout,'(a,3f15.5)') 'Sediment Ca ', & 410 & fq0, fq1, fq2 343 call trc_rst_dia_stat(zn_sed_n(:,:), 'Sediment N') 344 call trc_rst_dia_stat(zn_sed_fe(:,:), 'Sediment Fe') 345 call trc_rst_dia_stat(zn_sed_si(:,:), 'Sediment Si') 346 call trc_rst_dia_stat(zn_sed_c(:,:), 'Sediment C') 347 call trc_rst_dia_stat(zn_sed_ca(:,:), 'Sediment Ca') 411 348 !! 412 349 !! AXY (07/07/15): write out temporally averaged fields for DMS … … 427 364 !! calculate stats on these fields 428 365 IF(lwp) WRITE(numout,*) ' MEDUSA averaged properties for DMS stats (min, max, sum) ...' 429 fq0 = MINVAL(zn_dms_chn(:,:)) 430 fq1 = MAXVAL(zn_dms_chn(:,:)) 431 fq2 = SUM(zn_dms_chn(:,:)) 432 if (lwp) write (numout,'(a,3f15.5)') 'DMS, CHN ', fq0, fq1, fq2 433 fq0 = MINVAL(zn_dms_chd(:,:)) 434 fq1 = MAXVAL(zn_dms_chd(:,:)) 435 fq2 = SUM(zn_dms_chd(:,:)) 436 if (lwp) write (numout,'(a,3f15.5)') 'DMS, CHD ', fq0, fq1, fq2 437 fq0 = MINVAL(zn_dms_mld(:,:)) 438 fq1 = MAXVAL(zn_dms_mld(:,:)) 439 fq2 = SUM(zn_dms_mld(:,:)) 440 if (lwp) write (numout,'(a,3f15.5)') 'DMS, MLD ', fq0, fq1, fq2 441 fq0 = MINVAL(zn_dms_qsr(:,:)) 442 fq1 = MAXVAL(zn_dms_qsr(:,:)) 443 fq2 = SUM(zn_dms_qsr(:,:)) 444 if (lwp) write (numout,'(a,3f15.5)') 'DMS, QSR ', fq0, fq1, fq2 445 fq0 = MINVAL(zn_dms_din(:,:)) 446 fq1 = MAXVAL(zn_dms_din(:,:)) 447 fq2 = SUM(zn_dms_din(:,:)) 448 if (lwp) write (numout,'(a,3f15.5)') 'DMS, DIN ', fq0, fq1, fq2 366 call trc_rst_dia_stat(zn_dms_chn(:,:), 'DMS, CHN') 367 call trc_rst_dia_stat(zn_dms_chd(:,:), 'DMS, CHD') 368 call trc_rst_dia_stat(zn_dms_mld(:,:), 'DMS, MLD') 369 call trc_rst_dia_stat(zn_dms_qsr(:,:), 'DMS, QSR') 370 call trc_rst_dia_stat(zn_dms_din(:,:), 'DMS, DIN') 449 371 !! 450 372 #endif … … 462 384 IF(lwp) WRITE(numout,*) ' IDTRA averaged properties stats (min, max, sum) ...' 463 385 call trc_rst_dia_stat(qint_idtra(:,:,1), 'qint_IDTRA') 464 fq0 = MINVAL(qint_idtra(:,:,1))465 fq1 = MAXVAL(qint_idtra(:,:,1))466 fq2 = SUM(qint_idtra(:,:,1))467 if (lwp) write (numout,'(a,3f15.5)') 'qint_IDTRA ', fq0, fq1, fq2468 386 #endif 469 387 ! … … 475 393 !! 476 394 !! calculate stats on these fields 395 IF(lwp) WRITE(numout,*) ' CFC averaged properties stats (min, max, sum) ...' 477 396 call trc_rst_dia_stat(qint_cfc(:,:,jl), 'qint_'//ctrcnm(jn)) 478 IF(lwp) WRITE(numout,*) ' CFC averaged properties stats (min, max, sum) ...'479 fq0 = MINVAL(qint_cfc(:,:,jl))480 fq1 = MAXVAL(qint_cfc(:,:,jl))481 fq2 = SUM(qint_cfc(:,:,jl))482 if (lwp) write (numout,'(a,3f15.5)') 'qint_'//ctrcnm(jn), fq0, fq1, fq2483 397 END DO 484 398 #endif
Note: See TracChangeset
for help on using the changeset viewer.