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

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

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


Revision 1.1 - (hide annotations) (download)
Sun Nov 23 01:20:13 2003 UTC (20 years, 5 months ago) by jmc
Branch: MAIN
CVS Tags: checkpoint52l_pre, checkpoint52e_pre, hrcube4, checkpoint52j_post, checkpoint52e_post, hrcube_1, branch-netcdf, checkpoint52d_pre, checkpoint52l_post, checkpoint52k_post, checkpoint52d_post, checkpoint52b_post, checkpoint52f_post, checkpoint52c_post, hrcube5, checkpoint52i_post, checkpoint52j_pre, checkpoint52i_pre, checkpoint52h_pre, checkpoint52f_pre, hrcube_2, hrcube_3
Branch point for: netcdf-sm0
new pkg "thSIce" (replace therm_seaice).

1 jmc 1.1 C $Header: $
2     C $Name: $
3    
4     #include "THSICE_OPTIONS.h"
5    
6     CStartofinterface
7     C !ROUTINE: THSICE_AVE
8     C !INTERFACE:
9     SUBROUTINE THSICE_AVE( bi, bj, myThid )
10     C *==========================================================*
11     C | S/R THSICE_AVE
12     C | o save values for timeaveraging
13     C *==========================================================*
14     C \ev
15    
16     C !USES:
17     IMPLICIT NONE
18     C == Global variables ==
19     #include "SIZE.h"
20     c #include "EEPARAMS.h"
21     c #include "PARAMS.h"
22     #include "THSICE_PARAMS.h"
23     #include "THSICE.h"
24     #include "THSICE_DIAGS.h"
25    
26     C !INPUT/OUTPUT PARAMETERS:
27     C == Routine Arguments ==
28     INTEGER bi,bj
29     INTEGER myThid
30     CEndofinterface
31    
32     #ifdef ALLOW_THSICE
33     #ifdef ALLOW_TIMEAVE
34    
35     C-- Time-average
36     CALL TIMEAVE_CUMULATE(ICE_iceMask_AVE,iceMask,
37     & 1, thSIce_deltaT, bi, bj, myThid )
38     C-- area weighted average (with ice-fraction)
39     CALL TIMEAVE_CUMUL_2V( ICE_iceH_AVE, iceheight, iceMask,
40     & 1, 0, thSIce_deltaT, bi, bj, myThid )
41     CALL TIMEAVE_CUMUL_2V( ICE_snowH_AVE, snowheight, iceMask,
42     & 1, 0, thSIce_deltaT, bi, bj, myThid )
43     CALL TIMEAVE_CUMUL_2V( ICE_Tsrf_AVE, Tsrf, iceMask,
44     & 1, 0, thSIce_deltaT, bi, bj, myThid )
45     CALL TIMEAVE_CUMUL_2V( ICE_Tice1_AVE, Tice1, iceMask,
46     & 1, 0, thSIce_deltaT, bi, bj, myThid )
47     CALL TIMEAVE_CUMUL_2V( ICE_Tice2_AVE, Tice2, iceMask,
48     & 1, 0, thSIce_deltaT, bi, bj, myThid )
49     CALL TIMEAVE_CUMUL_2V( ICE_snow_AVE, snow, iceMask,
50     & 1, 0, thSIce_deltaT, bi, bj, myThid )
51     C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
52     C-- simple average (without ice-fraction)
53     c CALL TIMEAVE_CUMULATE(ICE_iceH_AVE,iceheight,
54     c & 1, thSIce_deltaT, bi, bj, myThid )
55     c CALL TIMEAVE_CUMULATE(ICE_snowH_AVE,snowHeight,
56     c & 1, thSIce_deltaT, bi, bj, myThid )
57     c CALL TIMEAVE_CUMULATE(ICE_Tsrf_AVE, Tsrf,
58     c & 1, thSIce_deltaT, bi, bj, myThid )
59     c CALL TIMEAVE_CUMULATE(ICE_Tice1_AVE,Tice1,
60     c & 1, thSIce_deltaT, bi, bj, myThid )
61     c CALL TIMEAVE_CUMULATE(ICE_Tice2_AVE,Tice2,
62     c & 1, thSIce_deltaT, bi, bj, myThid )
63     c CALL TIMEAVE_CUMULATE(ICE_snow_AVE, snow,
64     c & 1, thSIce_deltaT, bi, bj, myThid )
65    
66     ICE_timeave(1,bi,bj)=ICE_timeave(1,bi,bj)+thSIce_deltaT
67    
68     #endif /* ALLOW_TIMEAVE */
69     #endif /* ALLOW_THSICE */
70    
71     RETURN
72     END

  ViewVC Help
Powered by ViewVC 1.1.22