/[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.15 - (hide annotations) (download)
Thu May 2 20:01:21 2013 UTC (11 years ago) by jmc
Branch: MAIN
CVS Tags: checkpoint64y, checkpoint64x, checkpoint64z, checkpoint64q, checkpoint64p, checkpoint64s, checkpoint64r, checkpoint64u, checkpoint64t, checkpoint64w, checkpoint64v, checkpoint64i, checkpoint64h, checkpoint64k, checkpoint64j, checkpoint64m, checkpoint64l, checkpoint64o, checkpoint64n, checkpoint65, checkpoint66g, checkpoint66f, checkpoint66e, checkpoint66d, checkpoint66c, checkpoint66b, checkpoint66a, checkpoint66o, checkpoint66n, checkpoint66m, checkpoint66l, checkpoint66k, checkpoint66j, checkpoint66i, checkpoint66h, checkpoint65z, checkpoint65x, checkpoint65y, checkpoint65r, checkpoint65s, checkpoint65p, checkpoint65q, checkpoint65v, checkpoint65w, checkpoint65t, checkpoint65u, checkpoint65j, checkpoint65k, checkpoint65h, checkpoint65i, checkpoint65n, checkpoint65o, checkpoint65l, checkpoint65m, checkpoint65b, checkpoint65c, checkpoint65a, checkpoint65f, checkpoint65g, checkpoint65d, checkpoint65e, HEAD
Changes since 1.14: +3 -5 lines 
remove snow-precip output & time-ave
1 jmc 1.15 C $Header: /u/gcmpack/MITgcm/pkg/thsice/thsice_ave.F,v 1.14 2013/04/04 01:59:57 jmc Exp $
2 jmc 1.1 C $Name: $
3    
4     #include "THSICE_OPTIONS.h"
5    
6 jmc 1.2 CBOP
7 jmc 1.1 C !ROUTINE: THSICE_AVE
8     C !INTERFACE:
9 jmc 1.10 SUBROUTINE THSICE_AVE(
10 jmc 1.14 I bi, bj, myTime, myIter, myThid )
11 jmc 1.2 C !DESCRIPTION: \bv
12 jmc 1.1 C *==========================================================*
13 jmc 1.10 C | S/R THSICE_AVE
14 jmc 1.1 C | o save values for timeaveraging
15     C *==========================================================*
16     C \ev
17    
18     C !USES:
19     IMPLICIT NONE
20 jmc 1.2
21 jmc 1.1 C == Global variables ==
22     #include "SIZE.h"
23 jmc 1.2 #include "EEPARAMS.h"
24     #include "PARAMS.h"
25     #include "FFIELDS.h"
26 jmc 1.1 #include "THSICE_PARAMS.h"
27 jmc 1.2 #include "THSICE_VARS.h"
28     #include "THSICE_TAVE.h"
29 jmc 1.1
30     C !INPUT/OUTPUT PARAMETERS:
31     C == Routine Arguments ==
32 jmc 1.2 C frwAtm :: net fresh-water flux (E-P) to the atmosphere [kg/m2/s]
33     C flxAtm :: net heat flux from the atmosphere (+=down) [W/m2]
34 jmc 1.14 C bi,bj :: tile indices
35 jmc 1.2 C myTime :: time counter for this thread
36     C myIter :: iteration counter for this thread
37     C myThid :: thread number for this instance of the routine.
38 jmc 1.10 c _RL frwAtm(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
39     c _RL flxAtm(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
40 jmc 1.1 INTEGER bi,bj
41 jmc 1.2 _RL myTime
42     INTEGER myIter
43 jmc 1.1 INTEGER myThid
44 jmc 1.2 CEOP
45 jmc 1.1
46     #ifdef ALLOW_THSICE
47 jmc 1.2 C == Local variables ==
48 jmc 1.14 c LOGICAL DIFFERENT_MULTIPLE
49     c EXTERNAL DIFFERENT_MULTIPLE
50 jmc 1.10 c INTEGER wrtIter
51 jmc 1.2
52 jmc 1.1 #ifdef ALLOW_TIMEAVE
53    
54 jmc 1.3 IF ( thSIce_taveFreq .GT. 0. _d 0 ) THEN
55 jmc 1.2
56 jmc 1.1 C-- Time-average
57 jmc 1.14 CALL TIMEAVE_CUMULATE( ice_fract_Ave, iceMask,
58     & 1, thSIce_deltaT, bi, bj, myThid )
59 jmc 1.1 C-- area weighted average (with ice-fraction)
60 jmc 1.11 CALL TIMEAVE_CUMUL_2V( ice_iceH_Ave, iceHeight, iceMask,
61 jmc 1.1 & 1, 0, thSIce_deltaT, bi, bj, myThid )
62 jmc 1.11 CALL TIMEAVE_CUMUL_2V( ice_snowH_Ave, snowHeight, iceMask,
63 jmc 1.1 & 1, 0, thSIce_deltaT, bi, bj, myThid )
64 jmc 1.2 CALL TIMEAVE_CUMUL_2V( ice_Tsrf_Ave, Tsrf, iceMask,
65 jmc 1.1 & 1, 0, thSIce_deltaT, bi, bj, myThid )
66 jmc 1.2 CALL TIMEAVE_CUMUL_2V( ice_Tice1_Ave, Tice1, iceMask,
67 jmc 1.1 & 1, 0, thSIce_deltaT, bi, bj, myThid )
68 jmc 1.2 CALL TIMEAVE_CUMUL_2V( ice_Tice2_Ave, Tice2, iceMask,
69 jmc 1.1 & 1, 0, thSIce_deltaT, bi, bj, myThid )
70 jmc 1.15 c CALL TIMEAVE_CUMUL_2V( ice_snowPr_Ave,snowPrc, iceMask,
71     c & 1, 0, thSIce_deltaT, bi, bj, myThid )
72 jmc 1.6 C-- notes: siceAlb is already multiplied by iceFraction
73 jmc 1.14 CALL TIMEAVE_CUMULATE( ice_albedo_Ave, siceAlb,
74     & 1, thSIce_deltaT, bi, bj, myThid )
75 jmc 1.1 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
76     C-- simple average (without ice-fraction)
77 jmc 1.14 c CALL TIMEAVE_CUMULATE( ice_iceH_Ave,iceHeight,
78     c & 1, thSIce_deltaT, bi, bj, myThid )
79     c CALL TIMEAVE_CUMULATE( ice_snowH_Ave,snowHeight,
80     c & 1, thSIce_deltaT, bi, bj, myThid )
81     c CALL TIMEAVE_CUMULATE( ice_Tsrf_Ave, Tsrf,
82     c & 1, thSIce_deltaT, bi, bj, myThid )
83     c CALL TIMEAVE_CUMULATE( ice_Tice1_Ave,Tice1,
84     c & 1, thSIce_deltaT, bi, bj, myThid )
85     c CALL TIMEAVE_CUMULATE( ice_Tice2_Ave,Tice2,
86     c & 1, thSIce_deltaT, bi, bj, myThid )
87 jmc 1.1
88 jmc 1.10 C-- average of atmospheric fluxes :
89 jmc 1.4 IF ( fluidIsWater ) THEN
90 jmc 1.14 CALL TIMEAVE_CUMULATE( ice_flxAtm_Ave, icFlxAtm,
91     & 1, thSIce_deltaT, bi, bj, myThid )
92     CALL TIMEAVE_CUMULATE( ice_frwAtm_Ave, icFrwAtm,
93     & 1, thSIce_deltaT, bi, bj, myThid )
94 jmc 1.3 ENDIF
95 jmc 1.2
96 jmc 1.10 C-- average of oceanic fluxes :
97 jmc 1.2 #ifndef CHECK_ENERGY_CONSERV
98 jmc 1.14 #ifdef RS_IS_REAL8
99     CALL TIMEAVE_CUMULATE( ice_flx2oc_Ave, Qnet,
100     & 1, thSIce_deltaT, bi, bj, myThid )
101     CALL TIMEAVE_CUMULATE( ice_frw2oc_Ave, EmPmR,
102     & 1, thSIce_deltaT, bi, bj, myThid )
103     CALL TIMEAVE_CUMULATE( ice_salFx_Ave, saltFlux,
104     & 1, thSIce_deltaT, bi, bj, myThid )
105 jmc 1.12 #endif
106 jmc 1.2 #endif /* CHECK_ENERGY_CONSERV */
107    
108 jmc 1.13 ice_timeAve(bi,bj) = ice_timeAve(bi,bj)+thSIce_deltaT
109 jmc 1.1
110 jmc 1.3 ENDIF
111 jmc 1.1 #endif /* ALLOW_TIMEAVE */
112 jmc 1.2
113 jmc 1.10 c IF ( DIFFERENT_MULTIPLE(thSIce_diagFreq,
114     c & myTime+deltaTClock,deltaTClock)
115     c & .AND. fluidIsWater ) THEN
116     c wrtIter = myIter+1
117     c CALL WRITE_LOCAL_RL('ice_flxAtm','I10',1,flxAtm,
118     c & bi,bj,1,wrtIter,myThid)
119     c CALL WRITE_LOCAL_RL('ice_frwAtm','I10',1,frwAtm,
120     c & bi,bj,1,wrtIter,myThid)
121     c ENDIF
122 jmc 1.2
123 jmc 1.6 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
124    
125     #ifdef ALLOW_DIAGNOSTICS
126 jmc 1.5 IF ( useDiagnostics ) THEN
127 jmc 1.6
128     C-- Atmospheric fluxes :
129     IF ( fluidIsWater ) THEN
130 jmc 1.10 CALL DIAGNOSTICS_FILL(icFlxAtm,'SIflxAtm',0,1,1,bi,bj,myThid)
131     CALL DIAGNOSTICS_FILL(icFrwAtm,'SIfrwAtm',0,1,1,bi,bj,myThid)
132 jmc 1.6 ENDIF
133    
134     C-- Oceanic fluxes :
135     c #ifndef CHECK_ENERGY_CONSERV
136 jmc 1.14 CALL DIAGNOSTICS_FILL_RS(Qnet, 'SIflx2oc',0,1,1,bi,bj,myThid)
137     CALL DIAGNOSTICS_FILL_RS(EmPmR, 'SIfrw2oc',0,1,1,bi,bj,myThid)
138     CALL DIAGNOSTICS_FILL_RS(saltFlux,'SIsaltFx',0,1,1,bi,bj,myThid)
139 jmc 1.6 c #endif /* CHECK_ENERGY_CONSERV */
140    
141 jmc 1.5 ENDIF
142 jmc 1.6 #endif /* ALLOW_DIAGNOSTICS */
143 jmc 1.5
144 jmc 1.1 #endif /* ALLOW_THSICE */
145    
146     RETURN
147     END
  ViewVC Help
Powered by ViewVC 1.1.22