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 |