/[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.6 by jmc, Mon Jan 31 19:42:59 2005 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 125  C-- average of Ocean Mixed-Layer temp. & Line 128  C-- average of Ocean Mixed-Layer temp. &
128       &                      bi,bj,1,wrtIter,myThid)       &                      bi,bj,1,wrtIter,myThid)
129        ENDIF        ENDIF
130    
131    C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
132    
133    #ifdef ALLOW_DIAGNOSTICS
134        IF ( useDiagnostics ) THEN        IF ( useDiagnostics ) THEN
135          CALL THSICE_DIAGNOSTICS_FILL(  
136       I                   frwAtm, flxAtm,          CALL DIAGNOSTICS_FILL(siceAlb,'SIalbedo',0,1,1,bi,bj,myThid)
137       I                   myTime, myIter, bi,bj, myThid )  
138    C-- Atmospheric fluxes :
139           IF ( fluidIsWater ) THEN
140            CALL DIAGNOSTICS_FILL(flxAtm, 'SIflxAtm',0,1,2,bi,bj,myThid)
141            CALL DIAGNOSTICS_FILL(frwAtm, 'SIfrwAtm',0,1,2,bi,bj,myThid)
142           ENDIF
143    
144    C-- Oceanic fluxes :
145    c #ifndef CHECK_ENERGY_CONSERV
146            CALL DIAGNOSTICS_FILL(Qnet,    'SIflx2oc',0,1,1,bi,bj,myThid)
147            CALL DIAGNOSTICS_FILL(EmPmR,   'SIfrw2oc',0,1,1,bi,bj,myThid)
148            CALL DIAGNOSTICS_FILL(saltFlux,'SIsaltFx',0,1,1,bi,bj,myThid)
149    c #endif /* CHECK_ENERGY_CONSERV */
150    
151        ENDIF        ENDIF
152    #endif /* ALLOW_DIAGNOSTICS */
153    
154  #endif  /* ALLOW_THSICE */  #endif  /* ALLOW_THSICE */
155    
Legend:
Removed from v.1.5  
changed lines
  Added in v.1.6
  ViewVC Help
Powered by ViewVC 1.1.22