/[MITgcm]/MITgcm/pkg/diagnostics/diagstats_calc.F
ViewVC logotype

Diff of /MITgcm/pkg/diagnostics/diagstats_calc.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.4 by mlosch, Mon Dec 12 15:28:16 2011 UTC revision 1.5 by mlosch, Fri Jan 20 14:24:08 2012 UTC
# Line 467  C     to clean this fields first by copy Line 467  C     to clean this fields first by copy
467        tmpFld = 0. _d 0        tmpFld = 0. _d 0
468        DO j = 1,jRun        DO j = 1,jRun
469         DO i = 1,iRun         DO i = 1,iRun
470          tmpFld(i,j) = inpArr(i,j)          tmpFld(i,j) = inpArr(i,j)*tmpFac
471         ENDDO         ENDDO
472        ENDDO        ENDDO
473        IF ( power.EQ.2) THEN        IF ( power.EQ.2) THEN
# Line 478  C     sum up the volume Line 478  C     sum up the volume
478        statArr(0)  = SUM(tmpVol)*arrDr        statArr(0)  = SUM(tmpVol)*arrDr
479  C     compute and sum up volume*field  C     compute and sum up volume*field
480        tmpVol = tmpVol*tmpFld        tmpVol = tmpVol*tmpFld
481        statArr(1)  = SUM(tmpVol)*tmpFac*arrDr        statArr(1)  = SUM(tmpVol)*arrDr
482  C     compute and sum up volume*field**2  C     compute and sum up volume*field**2
483        tmpVol = tmpVol*tmpFld        tmpVol = tmpVol*tmpFld
484        statArr(2)  = SUM(tmpVol)*tmpFac*tmpFac*arrDr        statArr(2)  = SUM(tmpVol)*arrDr
485        statArr(im) = MINVAL(tmpFld, MASK = arrMaskL>0.)*tmpFac        statArr(im) = MINVAL(tmpFld, MASK = arrMaskL>0.)
486        statArr(ix) = MAXVAL(tmpFld, MASK = arrMaskL>0.)*tmpFac        statArr(ix) = MAXVAL(tmpFld, MASK = arrMaskL>0.)
487    
488  #endif /* TARGET_NEC_SX */  #endif /* TARGET_NEC_SX */
489    
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.22