Changeset 6220 for branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM
- Timestamp:
- 2016-01-06T17:21:02+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
r6213 r6220 300 300 fq1 = MAXVAL(qint_idtra(:,:,1)) 301 301 fq2 = SUM(qint_idtra(:,:,1)) 302 if (lwp) write (numout,'(a,3f15.5)') 'qint_ idtra', fq0, fq1, fq2302 if (lwp) write (numout,'(a,3f15.5)') 'qint_IDTRA ', fq0, fq1, fq2 303 303 #endif 304 304 ! 305 305 #if defined key_cfc 306 DO jn = jp_cfc0, jp_cfc1 306 DO jl = 1, jp_cfc 307 jn = jp_cfc0 + jl - 1 307 308 IF( iom_varid( numrtr, 'qint_'//ctrcnm(jn), ldstop = .FALSE. ) > 0 ) THEN 308 309 !! YES; in which case read them 309 310 !! 310 311 IF(lwp) WRITE(numout,*) ' CFC averaged properties present - reading in ...' 311 CALL iom_get( numrtr, jpdom_autoglo, 'qint_'//ctrcnm(jn), qint_cfc(:,:,j n) )312 CALL iom_get( numrtr, jpdom_autoglo, 'qint_'//ctrcnm(jn), qint_cfc(:,:,jl) ) 312 313 ELSE 313 314 !! NO; in which case set them to zero … … 319 320 !! calculate stats on these fields 320 321 IF(lwp) WRITE(numout,*) ' CFC averaged properties stats (min, max, sum) ...' 321 fq0 = MINVAL(qint_cfc(:,:,j n))322 fq1 = MAXVAL(qint_cfc(:,:,j n))323 fq2 = SUM(qint_cfc(:,:,j n))322 fq0 = MINVAL(qint_cfc(:,:,jl)) 323 fq1 = MAXVAL(qint_cfc(:,:,jl)) 324 fq2 = SUM(qint_cfc(:,:,jl)) 324 325 if (lwp) write (numout,'(a,3f15.5)') 'qint_'//ctrcnm(jn), fq0, fq1, fq2 325 326 END DO … … 336 337 INTEGER, INTENT( in ) :: kt ! ocean time-step index 337 338 !! 338 INTEGER :: jn 339 INTEGER :: jn, jl 339 340 REAL(wp) :: zarak0 340 341 !! AXY (05/11/13): temporary variables … … 459 460 fq1 = MAXVAL(qint_idtra(:,:,1)) 460 461 fq2 = SUM(qint_idtra(:,:,1)) 461 if (lwp) write (numout,'(a,3f15.5)') 'qint_ idtra', fq0, fq1, fq2462 if (lwp) write (numout,'(a,3f15.5)') 'qint_IDTRA ', fq0, fq1, fq2 462 463 #endif 463 464 ! 464 465 #if defined key_cfc 465 DO jn = jp_cfc0, jp_cfc1 466 DO jl = 1, jp_cfc 467 jn = jp_cfc0 + jl - 1 466 468 IF(lwp) WRITE(numout,*) ' CFC averaged properties - writing out ...' 467 CALL iom_rstput( kt, nitrst, numrtw, 'qint_'//ctrcnm(jn), qint_cfc(:,:,j n) )469 CALL iom_rstput( kt, nitrst, numrtw, 'qint_'//ctrcnm(jn), qint_cfc(:,:,jl) ) 468 470 !! 469 471 !! calculate stats on these fields 470 472 IF(lwp) WRITE(numout,*) ' CFC averaged properties stats (min, max, sum) ...' 471 fq0 = MINVAL(qint_cfc(:,:,j n))472 fq1 = MAXVAL(qint_cfc(:,:,j n))473 fq2 = SUM(qint_cfc(:,:,j n))473 fq0 = MINVAL(qint_cfc(:,:,jl)) 474 fq1 = MAXVAL(qint_cfc(:,:,jl)) 475 fq2 = SUM(qint_cfc(:,:,jl)) 474 476 if (lwp) write (numout,'(a,3f15.5)') 'qint_'//ctrcnm(jn), fq0, fq1, fq2 475 477 END DO
Note: See TracChangeset
for help on using the changeset viewer.