C $Header: /home/ubuntu/mnt/e9_copy/MITgcm/pkg/seaice/Attic/seaice_do_diags.F,v 1.20 2006/03/20 21:36:11 mlosch 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 DIFFERENT_MULTIPLE INTEGER i, j, k, bi, bj _RS arr(1-oLx:sNx+oLx,1-oLy:sNy+oLy,nSx,nSy) LOGICAL dumpFiles CHARACTER*(1) pf IF ( writeBinaryPrec .EQ. precFloat64 ) THEN pf(1:1) = 'D' ELSE pf(1:1) = 'R' ENDIF 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) #if ( defined(SEAICE_CGRID) && defined(SEAICE_ALLOW_EVP) ) IF ( SEAICEuseEVP ) THEN CALL MNC_CW_RS_W(pf,'sice',0,0,'si_sigma1', & seaice_sigma1,myThid) CALL MNC_CW_RS_W(pf,'sice',0,0,'si_sigma2', & seaice_sigma2,myThid) CALL MNC_CW_RS_W(pf,'sice',0,0,'si_sigma12', & seaice_sigma12,myThid) ENDIF #endif /* SEAICE_CGRID and SEAICE_ALLOW_EVP */ CALL MNC_CW_RS_W(pf,'sice',0,0,'si_UWIND',uwind,myThid) CALL MNC_CW_RS_W(pf,'sice',0,0,'si_VWIND',vwind,myThid) CALL MNC_CW_RS_W(pf,'sice',0,0,'fu',fu,myThid) CALL MNC_CW_RS_W(pf,'sice',0,0,'fv',fv,myThid) CALL MNC_CW_RS_W(pf,'sice',0,0,'EmPmR',EmPmR,myThid) CALL MNC_CW_RS_W(pf,'sice',0,0,'Qnet',Qnet,myThid) CALL MNC_CW_RS_W(pf,'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) #if ( defined(SEAICE_CGRID) && defined(SEAICE_ALLOW_EVP) ) IF ( SEAICEuseEVP ) THEN CALL WRITE_FLD_XY_Rs('SIGMA1.',suff,seaice_sigma1, & myIter,myThid) CALL WRITE_FLD_XY_Rs('SIGMA2.',suff,seaice_sigma2, & myIter,myThid) CALL WRITE_FLD_XY_Rs('SIGMA12.',suff,seaice_sigma12, & myIter,myThid) ENDIF #endif /* SEAICE_CGRID and SEAICE_ALLOW_EVP */ _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-Oly,sNy+Oly DO i=1-Olx,sNx+Olx arr(i,j,bi,bj)=UICE(i,j,1,bi,bj) ENDDO ENDDO ENDDO ENDDO IF (SEAICE_dump_mdsio) THEN _BARRIER _BEGIN_MASTER( myThid ) CALL WRITE_FLD_XY_RS( 'UICE.',suff,arr,myIter,myThid) _END_MASTER( myThid ) _BARRIER ENDIF #ifdef ALLOW_MNC IF ( useMNC .AND. SEAICE_dump_mnc ) THEN CALL MNC_CW_RS_W(pf,'sice',0,0,'si_UICE',arr,myThid) ENDIF #endif #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-Oly,sNy+Oly DO i=1-Olx,sNx+Olx arr(i,j,bi,bj)=VICE(i,j,1,bi,bj) ENDDO ENDDO ENDDO ENDDO IF (SEAICE_dump_mdsio) THEN _BARRIER _BEGIN_MASTER( myThid ) CALL WRITE_FLD_XY_RS( 'VICE.',suff,arr,myIter,myThid) _END_MASTER( myThid ) _BARRIER ENDIF #ifdef ALLOW_MNC IF ( useMNC .AND. SEAICE_dump_mnc ) THEN CALL MNC_CW_RS_W(pf,'sice',0,0,'si_VICE',arr,myThid) ENDIF #endif #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-Oly,sNy+Oly DO i=1-Olx,sNx+Olx arr(i,j,bi,bj)=HEFF(i,j,1,bi,bj) ENDDO ENDDO ENDDO ENDDO IF (SEAICE_dump_mdsio) THEN _BARRIER _BEGIN_MASTER( myThid ) CALL WRITE_FLD_XY_RS( 'HEFF.',suff,arr,myIter,myThid) _END_MASTER( myThid ) _BARRIER ENDIF #ifdef ALLOW_MNC IF ( useMNC .AND. SEAICE_dump_mnc ) THEN CALL MNC_CW_RS_W(pf,'sice',0,0,'si_HEFF',arr,myThid) ENDIF #endif #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-Oly,sNy+Oly DO i=1-Olx,sNx+Olx arr(i,j,bi,bj)=AREA(i,j,1,bi,bj) ENDDO ENDDO ENDDO ENDDO IF (SEAICE_dump_mdsio) THEN _BARRIER _BEGIN_MASTER( myThid ) CALL WRITE_FLD_XY_RS( 'AREA.',suff,arr,myIter,myThid) _END_MASTER( myThid ) _BARRIER ENDIF #ifdef ALLOW_MNC IF ( useMNC .AND. SEAICE_dump_mnc ) THEN CALL MNC_CW_RS_W(pf,'sice',0,0,'si_AREA',arr,myThid) ENDIF #endif #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 dumpFiles = & DIFFERENT_MULTIPLE(SEAICE_taveFreq,myTime,deltaTClock) #ifdef ALLOW_CAL IF ( useCAL ) THEN CALL CAL_TIME2DUMP( SEAICE_taveFreq, deltaTClock, U dumpFiles, I myTime, myIter, myThid ) 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) C CALL MNC_CW_RL_W_S('D','sice_tave',0,0,'model_time', C & myTime,myThid) CALL MNC_CW_RL_W(pf,'sice_tave',0,0, & 'si_UICEtave',UICEtave,myThid) CALL MNC_CW_RL_W(pf,'sice_tave',0,0, & 'si_VICEtave',VICEtave,myThid) CALL MNC_CW_RL_W(pf,'sice_tave',0,0, & 'si_FUtave',FUtave,myThid) CALL MNC_CW_RL_W(pf,'sice_tave',0,0, & 'si_FVtave',FVtave,myThid) CALL MNC_CW_RL_W(pf,'sice_tave',0,0, & 'si_EmPmRtave',EmPmRtave,myThid) CALL MNC_CW_RL_W(pf,'sice_tave',0,0, & 'si_QNETtave',QNETtave,myThid) CALL MNC_CW_RL_W(pf,'sice_tave',0,0, & 'si_QSWtave',QSWtave,myThid) CALL MNC_CW_RL_W(pf,'sice_tave',0,0, & 'si_HEFFtave',HEFFtave,myThid) CALL MNC_CW_RL_W(pf,'sice_tave',0,0, & 'si_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