/[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.3 by jmc, Thu Apr 8 18:54:26 2004 UTC revision 1.8 by jmc, Sun May 15 03:04:57 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 86  c      CALL TIMEAVE_CUMULATE(ice_snowPr_ Line 89  c      CALL TIMEAVE_CUMULATE(ice_snowPr_
89  c    &                       1, thSIce_deltaT, bi, bj, myThid )  c    &                       1, thSIce_deltaT, bi, bj, myThid )
90    
91  C-- average of atmospheric fluxes :  C-- average of atmospheric fluxes :
92  #ifdef ALLOW_TIMEAVE        IF ( fluidIsWater ) THEN
       IF ( buoyancyRelation(1:7) .EQ. 'OCEANIC' ) THEN  
93         CALL TIMEAVE_CUMUL_1T(ice_flxAtm_Ave, flxAtm,         CALL TIMEAVE_CUMUL_1T(ice_flxAtm_Ave, flxAtm,
94       &                       1, thSIce_deltaT, bi, bj, myThid )       &                       1, thSIce_deltaT, bi, bj, myThid )
95         CALL TIMEAVE_CUMUL_1T(ice_frwAtm_Ave, frwAtm,         CALL TIMEAVE_CUMUL_1T(ice_frwAtm_Ave, frwAtm,
96       &                       1, thSIce_deltaT, bi, bj, myThid )       &                       1, thSIce_deltaT, bi, bj, myThid )
97        ENDIF        ENDIF
 #endif /*ALLOW_TIMEAVE*/  
98    
99  C-- average of oceanic fluxes :  C-- average of oceanic fluxes :
100  #ifndef CHECK_ENERGY_CONSERV  #ifndef CHECK_ENERGY_CONSERV
# Line 118  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       &     .AND. buoyancyRelation(1:7) .EQ. 'OCEANIC' ) THEN       &                        myTime+deltaTClock,deltaTClock)
124         &     .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,
127       &                      bi,bj,1,wrtIter,myThid)       &                      bi,bj,1,wrtIter,myThid)
# Line 127  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
136    
137            CALL DIAGNOSTICS_FILL(siceAlb,'SIalbedo',0,1,1,bi,bj,myThid)
138    
139    C-- Atmospheric fluxes :
140           IF ( fluidIsWater ) THEN
141            CALL DIAGNOSTICS_FILL(flxAtm, 'SIflxAtm',0,1,2,bi,bj,myThid)
142            CALL DIAGNOSTICS_FILL(frwAtm, 'SIfrwAtm',0,1,2,bi,bj,myThid)
143           ENDIF
144    
145    C-- Oceanic fluxes :
146    c #ifndef CHECK_ENERGY_CONSERV
147            CALL DIAGNOSTICS_FILL(Qnet,    'SIflx2oc',0,1,1,bi,bj,myThid)
148            CALL DIAGNOSTICS_FILL(EmPmR,   'SIfrw2oc',0,1,1,bi,bj,myThid)
149            CALL DIAGNOSTICS_FILL(saltFlux,'SIsaltFx',0,1,1,bi,bj,myThid)
150    c #endif /* CHECK_ENERGY_CONSERV */
151    
152          ENDIF
153    #endif /* ALLOW_DIAGNOSTICS */
154    
155  #endif  /* ALLOW_THSICE */  #endif  /* ALLOW_THSICE */
156    
157        RETURN        RETURN
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.8
  ViewVC Help
Powered by ViewVC 1.1.22