C $Header: /u/gcmpack/MITgcm/pkg/seaice/seaice_do_diags.F,v 1.14 2005/05/25 04:03:10 edhill Exp $
C $Name: $
#include "SEAICE_OPTIONS.h"
SUBROUTINE SEAICE_DO_DIAGS( myTime, myIter, myThid )
C /==========================================================\
C | SUBROUTINE SEAICE_DO_DIAGS |
C | o Do SEAICE diagnostic output. |
C \==========================================================/
IMPLICIT NONE
C === Global variables ===
#include "SIZE.h"
#include "EEPARAMS.h"
#include "PARAMS.h"
#include "FFIELDS.h"
#include "SEAICE_DIAGS.h"
#include "SEAICE_PARAMS.h"
#include "SEAICE_FFIELDS.h"
#include "SEAICE.h"
C == Routine arguments ==
C myTime - Current time of simulation ( s )
C myIter - Iteration number
C myThid - Number of this instance of SEAICE_DO_DIAGS
_RL myTime
INTEGER myIter
INTEGER myThid
C == Local variables ==
CHARACTER*(MAX_LEN_MBUF) suff
LOGICAL DIFFERENT_MULTIPLE
EXTERNAL
INTEGER i, j, k, bi, bj
_RS arr(1-oLx:sNx+oLx,1-oLy:sNy+oLy,nSx,nSy)
INTEGER thisdate(4), prevdate(4)
LOGICAL dumpFiles
IF (SEAICEwriteState) THEN
IF ( DIFFERENT_MULTIPLE(SEAICE_dumpFreq,myTime,deltaTClock)
& ) THEN
#ifdef ALLOW_MNC
IF (useMNC .AND. SEAICE_dump_mnc) THEN
CALL MNC_CW_SET_UDIM('sice', -1, myThid)
CALL MNC_CW_RL_W_S('D','sice',0,0,'T', myTime, myThid)
CALL MNC_CW_SET_UDIM('sice', 0, myThid)
CALL MNC_CW_I_W_S('I','sice',0,0,'iter', myIter, myThid)
CALL MNC_CW_RL_W_S('D','sice',0,0,'model_time',
& myTime,myThid)
CALL MNC_CW_RL_W('D','sice',0,0,'UWIND',uwind,myThid)
CALL MNC_CW_RL_W('D','sice',0,0,'VWIND',vwind,myThid)
CALL MNC_CW_RL_W('D','sice',0,0,'FU',fu,myThid)
CALL MNC_CW_RL_W('D','sice',0,0,'FV',fv,myThid)
CALL MNC_CW_RL_W('D','sice',0,0,'EmPmR',EmPmR,myThid)
CALL MNC_CW_RL_W('D','sice',0,0,'Qnet',Qnet,myThid)
CALL MNC_CW_RL_W('D','sice',0,0,'Qsw',Qsw,myThid)
ENDIF
#endif
IF (SEAICE_dump_mdsio) THEN
WRITE(suff,'(I10.10)') myIter
_BARRIER
_BEGIN_MASTER( myThid )
CALL WRITE_FLD_XY_RS( 'UWIND.',suff,uwind,myIter,myThid)
CALL WRITE_FLD_XY_RS( 'VWIND.',suff,vwind,myIter,myThid)
CALL WRITE_FLD_XY_RS( 'FU.',suff,fu,myIter,myThid)
CALL WRITE_FLD_XY_RS( 'FV.',suff,fv,myIter,myThid)
CALL WRITE_FLD_XY_RS( 'EmPmR.',suff,EmPmR,myIter,myThid)
CALL WRITE_FLD_XY_RS( 'Qnet.',suff,Qnet,myIter,myThid)
CALL WRITE_FLD_XY_RS( 'Qsw.',suff,Qsw,myIter,myThid)
_END_MASTER( myThid )
_BARRIER
ENDIF
#ifdef SEAICE_DEBUG
CALL PLOT_FIELD_XYRS( uwind , 'Current uwind ', myIter, myThid )
CALL PLOT_FIELD_XYRS( vwind , 'Current vwind ', myIter, myThid )
CALL PLOT_FIELD_XYRS( atemp , 'Current atemp ', myIter, myThid )
CALL PLOT_FIELD_XYRS( aqh , 'Current aqh ', myIter, myThid )
CALL PLOT_FIELD_XYRS( lwdown, 'Current lwdown', myIter, myThid )
CALL PLOT_FIELD_XYRS( swdown, 'Current swdown', myIter, myThid )
CALL PLOT_FIELD_XYRS( precip, 'Current precip', myIter, myThid )
CALL PLOT_FIELD_XYRL( evap , 'Current evap ', myIter, myThid )
CALL PLOT_FIELD_XYRS( runoff, 'Current runoff', myIter, myThid )
CALL PLOT_FIELD_XYRS( SSS , 'Current SSS ', myIter, myThid )
CALL PLOT_FIELD_XYRS( SST , 'Current SST ', myIter, myThid )
CALL PLOT_FIELD_XYRL( fu , 'Current fu ', myIter, myThid )
CALL PLOT_FIELD_XYRL( fv , 'Current fv ', myIter, myThid )
CALL PLOT_FIELD_XYRL( EmPmR , 'Current EmPmR ', myIter, myThid )
CALL PLOT_FIELD_XYRL( Qnet , 'Current Qnet ', myIter, myThid )
CALL PLOT_FIELD_XYRL( Qsw , 'Current Qsw ', myIter, myThid )
#endif
DO bj=myByLo(myThid),myByHi(myThid)
DO bi=myBxLo(myThid),myBxHi(myThid)
DO j=1,sNy
DO i=1,sNx
arr(i,j,bi,bj)=UICE(i,j,1,bi,bj)
ENDDO
ENDDO
ENDDO
ENDDO
_BARRIER
_BEGIN_MASTER( myThid )
CALL WRITE_FLD_XY_RS( 'UICE.',suff,arr,myIter,myThid)
_END_MASTER( myThid )
_BARRIER
#ifdef SEAICE_DEBUG
_EXCH_XY_R4( arr, myThid )
CALL PLOT_FIELD_XYRS( arr , 'Current uice ', myIter, myThid )
#endif
DO bj=myByLo(myThid),myByHi(myThid)
DO bi=myBxLo(myThid),myBxHi(myThid)
DO j=1,sNy
DO i=1,sNx
arr(i,j,bi,bj)=VICE(i,j,1,bi,bj)
ENDDO
ENDDO
ENDDO
ENDDO
_BARRIER
_BEGIN_MASTER( myThid )
CALL WRITE_FLD_XY_RS( 'VICE.',suff,arr,myIter,myThid)
_END_MASTER( myThid )
_BARRIER
#ifdef SEAICE_DEBUG
_EXCH_XY_R4( arr, myThid )
CALL PLOT_FIELD_XYRS( arr , 'Current vice ', myIter, myThid )
#endif
DO bj=myByLo(myThid),myByHi(myThid)
DO bi=myBxLo(myThid),myBxHi(myThid)
DO j=1,sNy
DO i=1,sNx
arr(i,j,bi,bj)=HEFF(i,j,1,bi,bj)
ENDDO
ENDDO
ENDDO
ENDDO
_BARRIER
_BEGIN_MASTER( myThid )
CALL WRITE_FLD_XY_RS( 'HEFF.',suff,arr,myIter,myThid)
_END_MASTER( myThid )
_BARRIER
#ifdef SEAICE_DEBUG
_EXCH_XY_R4( arr, myThid )
CALL PLOT_FIELD_XYRS( arr , 'Current heff ', myIter, myThid )
#endif
DO bj=myByLo(myThid),myByHi(myThid)
DO bi=myBxLo(myThid),myBxHi(myThid)
DO j=1,sNy
DO i=1,sNx
arr(i,j,bi,bj)=AREA(i,j,1,bi,bj)
ENDDO
ENDDO
ENDDO
ENDDO
_BARRIER
_BEGIN_MASTER( myThid )
CALL WRITE_FLD_XY_RS( 'AREA.',suff,arr,myIter,myThid)
_END_MASTER( myThid )
_BARRIER
#ifdef SEAICE_DEBUG
_EXCH_XY_R4( arr, myThid )
CALL PLOT_FIELD_XYRS( arr , 'Current area ', myIter, myThid )
#endif
ENDIF
ENDIF
C----------------------------------------------------------------
C Do SEAICE time averaging.
C----------------------------------------------------------------
#ifdef ALLOW_TIMEAVE
C-- Time-cumulations
DO bj = myByLo(myThid), myByHi(myThid)
DO bi = myBxLo(myThid), myBxHi(myThid)
DO j=1,sNy
DO i=1,sNx
FUtave(i,j,1,bi,bj) =
& FUtave(i,j,1,bi,bj) +FU(i,j,bi,bj) *deltaTclock
FVtave(i,j,1,bi,bj) =
& FVtave(i,j,1,bi,bj) +FV(i,j,bi,bj) *deltaTclock
EmPmRtave(i,j,1,bi,bj)=
& EmPmRtave(i,j,1,bi,bj)+EmPmR(i,j,bi,bj) *deltaTclock
QNETtave(i,j,1,bi,bj) =
& QNETtave(i,j,1,bi,bj) +QNET(i,j,bi,bj) *deltaTclock
QSWtave(i,j,1,bi,bj) =
& QSWtave(i,j,1,bi,bj) +QSW(i,j,bi,bj) *deltaTclock
UICEtave(i,j,1,bi,bj) =
& UICEtave(i,j,1,bi,bj) +UICE(i,j,1,bi,bj)*deltaTclock
VICEtave(i,j,1,bi,bj) =
& VICEtave(i,j,1,bi,bj) +VICE(i,j,1,bi,bj)*deltaTclock
HEFFtave(i,j,1,bi,bj) =
& HEFFtave(i,j,1,bi,bj) +HEFF(i,j,1,bi,bj)*deltaTclock
AREAtave(i,j,1,bi,bj) =
& AREAtave(i,j,1,bi,bj) +AREA(i,j,1,bi,bj)*deltaTclock
ENDDO
ENDDO
DO k=1,Nr
SEAICE_TimeAve(k,bi,bj)=SEAICE_TimeAve(k,bi,bj)+deltaTclock
ENDDO
ENDDO
ENDDO
C Dump files and restart average computation if needed
dumpFiles = .FALSE.
IF ( myIter .NE. nIter0 ) THEN
IF ( DIFFERENT_MULTIPLE(SEAICE_taveFreq,myTime,deltaTClock) )
& dumpFiles = .TRUE.
#ifdef ALLOW_CAL
IF ( calendarDumps .AND. (
& (SEAICE_taveFreq.GE. 2592000.AND.SEAICE_taveFreq.LE. 2678400).OR.
& (SEAICE_taveFreq.GE.31104000.AND.SEAICE_taveFreq.LE.31968000)))
& THEN
C-- Convert approximate months (30-31 days) and years (360-372 days)
C to exact calendar months and years.
C- First determine calendar dates for this and previous time step.
call CAL_GETDATE( myiter ,mytime ,thisdate,mythid )
call CAL_GETDATE( myiter-1,mytime-deltaTClock,prevdate,mythid )
dumpFiles = .FALSE.
C- Monthly SEAICE_taveFreq:
IF( SEAICE_taveFreq.GE. 2592000 .AND. SEAICE_taveFreq.LE. 2678400
& .AND. (thisdate(1)-prevdate(1)).GT.50 ) dumpFiles = .TRUE.
C- Yearly SEAICE_taveFreq:
IF( SEAICE_taveFreq.GE.31104000 .AND. SEAICE_taveFreq.LE.31968000
& .AND. (thisdate(1)-prevdate(1)).GT.5000 ) dumpFiles = .TRUE.
ENDIF
#endif
ENDIF
IF (dumpFiles) THEN
C Normalize by integrated time
DO bj = myByLo(myThid), myByHi(myThid)
DO bi = myBxLo(myThid), myBxHi(myThid)
CALL TIMEAVE_NORMALIZ(FUtave ,SEAICE_timeave, 1,
& bi,bj,myThid)
CALL TIMEAVE_NORMALIZ(FVtave ,SEAICE_timeave, 1,
& bi,bj,myThid)
CALL TIMEAVE_NORMALIZ(EmPmRtave,SEAICE_timeave, 1,
& bi,bj,myThid)
CALL TIMEAVE_NORMALIZ(QNETtave ,SEAICE_timeave, 1,
& bi,bj,myThid)
CALL TIMEAVE_NORMALIZ(QSWtave ,SEAICE_timeave, 1,
& bi,bj,myThid)
CALL TIMEAVE_NORMALIZ(UICEtave ,SEAICE_timeave, 1,
& bi,bj,myThid)
CALL TIMEAVE_NORMALIZ(VICEtave ,SEAICE_timeave, 1,
& bi,bj,myThid)
CALL TIMEAVE_NORMALIZ(HEFFtave ,SEAICE_timeave, 1,
& bi,bj,myThid)
CALL TIMEAVE_NORMALIZ(AREAtave ,SEAICE_timeave, 1,
& bi,bj,myThid)
ENDDO
ENDDO
#ifdef ALLOW_MNC
IF (useMNC .AND. SEAICE_tave_mnc) THEN
CALL MNC_CW_SET_UDIM('sice_tave', -1, myThid)
CALL MNC_CW_RL_W_S('D','sice_tave',0,0,'T', myTime, myThid)
CALL MNC_CW_SET_UDIM('sice_tave', 0, myThid)
CALL MNC_CW_I_W_S('I','sice_tave',0,0,'iter', myIter, myThid)
CALL MNC_CW_RL_W_S('D','sice_tave',0,0,'model_time',
& myTime,myThid)
CALL MNC_CW_RL_W('R','sice_tave',0,0,
& 'UICEtave',UICEtave,myThid)
CALL MNC_CW_RL_W('R','sice_tave',0,0,
& 'VICEtave',VICEtave,myThid)
CALL MNC_CW_RL_W('R','sice_tave',0,0,
& 'FUtave',FUtave,myThid)
CALL MNC_CW_RL_W('R','sice_tave',0,0,
& 'FVtave',FVtave,myThid)
CALL MNC_CW_RL_W('R','sice_tave',0,0,
& 'EmPmRtave',EmPmRtave,myThid)
CALL MNC_CW_RL_W('R','sice_tave',0,0,
& 'QNETtave',QNETtave,myThid)
CALL MNC_CW_RL_W('R','sice_tave',0,0,
& 'QSWtave',QSWtave,myThid)
CALL MNC_CW_RL_W('R','sice_tave',0,0,
& 'HEFFtave',HEFFtave,myThid)
CALL MNC_CW_RL_W('R','sice_tave',0,0,
& 'AREAtave',AREAtave,myThid)
ENDIF
#endif
IF (SEAICE_tave_mdsio) THEN
WRITE(suff,'(I10.10)') myIter
_BARRIER
_BEGIN_MASTER( myThid )
CALL WRITE_FLD_XY_RL('FUtave.' ,suff,FUtave ,myIter,myThid)
CALL WRITE_FLD_XY_RL('FVtave.' ,suff,FVtave ,myIter,myThid)
CALL WRITE_FLD_XY_RL('EmPmRtave.',suff,EmPmRtave,myIter,myThid)
CALL WRITE_FLD_XY_RL('QNETtave.' ,suff,QNETtave ,myIter,myThid)
CALL WRITE_FLD_XY_RL('QSWtave.' ,suff,QSWtave ,myIter,myThid)
CALL WRITE_FLD_XY_RL('UICEtave.' ,suff,UICEtave ,myIter,myThid)
CALL WRITE_FLD_XY_RL('VICEtave.' ,suff,VICEtave ,myIter,myThid)
CALL WRITE_FLD_XY_RL('HEFFtave.' ,suff,HEFFtave ,myIter,myThid)
CALL WRITE_FLD_XY_RL('AREAtave.' ,suff,AREAtave ,myIter,myThid)
_END_MASTER( myThid )
_BARRIER
ENDIF
C Reset averages to zero
DO bj = myByLo(myThid), myByHi(myThid)
DO bi = myBxLo(myThid), myBxHi(myThid)
CALL TIMEAVE_RESET(FUtave ,1,bi,bj,myThid)
CALL TIMEAVE_RESET(FVtave ,1,bi,bj,myThid)
CALL TIMEAVE_RESET(EmPmRtave,1,bi,bj,myThid)
CALL TIMEAVE_RESET(QNETtave ,1,bi,bj,myThid)
CALL TIMEAVE_RESET(QSWtave ,1,bi,bj,myThid)
CALL TIMEAVE_RESET(UICEtave ,1,bi,bj,myThid)
CALL TIMEAVE_RESET(VICEtave ,1,bi,bj,myThid)
CALL TIMEAVE_RESET(HEFFtave ,1,bi,bj,myThid)
CALL TIMEAVE_RESET(AREAtave ,1,bi,bj,myThid)
DO k=1,Nr
SEAICE_TimeAve(k,bi,bj)=ZERO
ENDDO
ENDDO
ENDDO
ENDIF
#endif /* ALLOW_TIMEAVE */
RETURN
END