/[MITgcm]/MITgcm/pkg/thsice/thsice_ave.F
ViewVC logotype

Diff of /MITgcm/pkg/thsice/thsice_ave.F

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

revision 1.5 by jmc, Tue Jan 4 00:25:31 2005 UTC revision 1.9 by jmc, Tue Feb 21 01:15:15 2006 UTC
# Line 70  C-- area weighted average (with ice-frac Line 70  C-- area weighted average (with ice-frac
70       &                    1, 0, thSIce_deltaT, bi, bj, myThid )       &                    1, 0, thSIce_deltaT, bi, bj, myThid )
71         CALL TIMEAVE_CUMUL_2V( ice_snowPr_Ave,snowPrc, iceMask,         CALL TIMEAVE_CUMUL_2V( ice_snowPr_Ave,snowPrc, iceMask,
72       &                    1, 0, thSIce_deltaT, bi, bj, myThid )       &                    1, 0, thSIce_deltaT, bi, bj, myThid )
73    C--    notes: siceAlb is already multiplied by iceFraction
74           CALL TIMEAVE_CUMULATE(ice_albedo_Ave, siceAlb,
75         &                       1, thSIce_deltaT, bi, bj, myThid )
76  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
77  C-- simple average (without ice-fraction)  C-- simple average (without ice-fraction)
78  c      CALL TIMEAVE_CUMULATE(ice_iceH_Ave,iceheight,  c      CALL TIMEAVE_CUMULATE(ice_iceH_Ave,iceheight,
# Line 116  C-- average of Ocean Mixed-Layer temp. & Line 119  C-- average of Ocean Mixed-Layer temp. &
119        ENDIF        ENDIF
120  #endif /* ALLOW_TIMEAVE */  #endif /* ALLOW_TIMEAVE */
121    
122        IF ( DIFFERENT_MULTIPLE(thSIce_diagFreq,myTime+deltaTClock,myTime)        IF ( DIFFERENT_MULTIPLE(thSIce_diagFreq,
123         &                        myTime+deltaTClock,deltaTClock)
124       &     .AND. fluidIsWater ) THEN       &     .AND. fluidIsWater ) THEN
125          wrtIter = myIter+1          wrtIter = myIter+1
126          CALL WRITE_LOCAL_RL('ice_flxAtm','I10',1,flxAtm,          CALL WRITE_LOCAL_RL('ice_flxAtm','I10',1,flxAtm,
# Line 125  C-- average of Ocean Mixed-Layer temp. & Line 129  C-- average of Ocean Mixed-Layer temp. &
129       &                      bi,bj,1,wrtIter,myThid)       &                      bi,bj,1,wrtIter,myThid)
130        ENDIF        ENDIF
131    
132    C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
133    
134    #ifdef ALLOW_DIAGNOSTICS
135        IF ( useDiagnostics ) THEN        IF ( useDiagnostics ) THEN
136          CALL THSICE_DIAGNOSTICS_FILL(  
137       I                   frwAtm, flxAtm,  C-- Atmospheric fluxes :
138       I                   myTime, myIter, bi,bj, myThid )         IF ( fluidIsWater ) THEN
139            CALL DIAGNOSTICS_FILL(flxAtm, 'SIflxAtm',0,1,2,bi,bj,myThid)
140            CALL DIAGNOSTICS_FILL(frwAtm, 'SIfrwAtm',0,1,2,bi,bj,myThid)
141           ENDIF
142    
143    C-- Oceanic fluxes :
144    c #ifndef CHECK_ENERGY_CONSERV
145            CALL DIAGNOSTICS_FILL(Qnet,    'SIflx2oc',0,1,1,bi,bj,myThid)
146            CALL DIAGNOSTICS_FILL(EmPmR,   'SIfrw2oc',0,1,1,bi,bj,myThid)
147            CALL DIAGNOSTICS_FILL(saltFlux,'SIsaltFx',0,1,1,bi,bj,myThid)
148    c #endif /* CHECK_ENERGY_CONSERV */
149    
150        ENDIF        ENDIF
151    #endif /* ALLOW_DIAGNOSTICS */
152    
153  #endif  /* ALLOW_THSICE */  #endif  /* ALLOW_THSICE */
154    
Legend:
Removed from v.1.5  
changed lines
  Added in v.1.9
  ViewVC Help
Powered by ViewVC 1.1.22