/[MITgcm]/MITgcm/pkg/therm_seaice/ice_ave.F
ViewVC logotype

Contents of /MITgcm/pkg/therm_seaice/ice_ave.F

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


Revision 1.2 - (show annotations) (download)
Wed Dec 4 14:53:39 2002 UTC (21 years, 5 months ago) by cheisey
Branch: MAIN
CVS Tags: checkpoint47b_post, checkpoint51f_post, checkpoint48i_post, checkpoint51l_post, checkpoint51k_post, checkpoint48d_post, checkpoint50b_post, checkpoint47i_post, checkpoint51o_post, checkpoint48g_post, branchpoint-genmake2, checkpoint50c_pre, checkpoint50, checkpoint51j_post, branch-exfmods-tag, checkpoint51q_post, checkpoint47e_post, checkpoint50f_post, checkpoint50a_post, checkpoint48e_post, checkpoint47c_post, checkpoint50f_pre, checkpoint52a_pre, checkpoint48b_post, checkpoint47j_post, checkpoint47d_pre, checkpoint50d_pre, checkpoint47h_post, checkpoint51d_post, checkpoint51, checkpoint51r_post, checkpoint48c_pre, checkpoint52, checkpoint52b_pre, checkpoint48c_post, checkpoint50d_post, checkpoint51o_pre, checkpoint47f_post, checkpoint51t_post, checkpoint51b_pre, checkpoint52a_post, checkpoint51i_post, checkpoint50e_post, checkpoint47g_post, checkpoint50h_post, checkpoint50c_post, checkpoint51a_post, checkpoint51n_pre, checkpoint47d_post, checkpoint50e_pre, checkpoint50b_pre, checkpoint48d_pre, checkpoint50i_post, checkpoint51p_post, checkpoint51n_post, checkpoint51e_post, checkpoint51b_post, checkpoint48a_post, checkpoint51h_pre, checkpoint48f_post, checkpoint51i_pre, checkpoint51l_pre, checkpoint50g_post, checkpoint51u_post, checkpoint51c_post, checkpoint51g_post, checkpoint51m_post, ecco_c52_e35, checkpoint48, checkpoint49, checkpoint51f_pre, checkpoint51s_post, checkpoint48h_post
Branch point for: checkpoint51n_branch, branch-nonh, tg2-branch, branch-genmake2, branch-exfmods-curt
Changes since 1.1: +1 -1 lines
Fixing bug in which diagnostic quantities are
not properly computed when using multiple tiles
on a single processor (Time cumulation array needs
to be dimensioned over Nr, but only k=1 is actually
used.  This makes routine compatible with routines
in pkg/timeave).

1 #include "CPP_OPTIONS.h"
2
3 CStartofinterface
4 C !ROUTINE: ICE_AVE
5 C !INTERFACE:
6 SUBROUTINE ICE_AVE( bi, bj, iMin, iMax, jMin, jMax,
7 I myThid )
8 C *==========================================================*
9 C | S/R ICE_AVE
10 C | o save values for timeaveraging
11 C *==========================================================*
12 C \ev
13
14 C !USES:
15 IMPLICIT NONE
16 C == Global variables ==
17 #include "SIZE.h"
18 #include "DYNVARS.h"
19 #include "EEPARAMS.h"
20 #include "PARAMS.h"
21 #include "GRID.h"
22 cswdice -- add ---
23 #ifdef ALLOW_THERM_SEAICE
24 #include "ICE.h"
25 #include "ICE_DIAGS.h"
26 #endif
27 cswdice -- end add ---
28
29 C !INPUT/OUTPUT PARAMETERS:
30 C == Routine Arguments ==
31 INTEGER bi,bj,iMin,iMax,jMin,jMax
32 INTEGER myThid
33 CEndofinterface
34
35 #ifdef ALLOW_THERM_SEAICE
36
37 #ifdef ALLOW_TIMEAVE
38
39 c
40 integer i,j
41
42 C-- Time-average
43 DO j=jMin,jMax
44 DO i=iMin,iMax
45 ICE_icemask_AVE(i,j,bi,bj)=ICE_icemask_AVE(i,j,bi,bj)
46 & +iceMask(i,j,bi,bj)*deltaTclock
47 ICE_iceheight_AVE(i,j,bi,bj)=ICE_iceheight_AVE(i,j,bi,bj)
48 & +iceheight(i,j,bi,bj)*deltaTclock
49 ICE_snowheight_AVE(i,j,bi,bj)=ICE_snowheight_AVE(i,j,bi,bj)
50 & +snowHeight(i,j,bi,bj)*deltaTclock
51 ICE_Tsrf_AVE(i,j,bi,bj)=ICE_Tsrf_AVE(i,j,bi,bj)
52 & +Tsrf(i,j,bi,bj)*deltaTclock
53 ICE_Tice1_AVE(i,j,bi,bj)=ICE_Tice1_AVE(i,j,bi,bj)
54 & +Tice1(i,j,bi,bj)*deltaTclock
55 ICE_Tice2_AVE(i,j,bi,bj)=ICE_Tice2_AVE(i,j,bi,bj)
56 & +Tice2(i,j,bi,bj)*deltaTclock
57 ICE_snow_AVE(i,j,bi,bj)=ICE_snow_AVE(i,j,bi,bj)
58 & +snow(i,j,bi,bj)*deltaTclock
59 ENDDO
60 ENDDO
61
62 ICE_timeave(1,bi,bj)=ICE_timeave(1,bi,bj)+deltaTclock
63 c print*,'QQ end ice ave '
64 #endif /* ALLOW_TIMEAVE */
65
66 #endif
67
68 RETURN
69 END

  ViewVC Help
Powered by ViewVC 1.1.22