/[MITgcm]/MITgcm/pkg/bulk_force/bulkf_ave.F
ViewVC logotype

Contents of /MITgcm/pkg/bulk_force/bulkf_ave.F

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

Revision 1.9 - (show annotations) (download)
Sat Jan 2 22:46:29 2010 UTC (14 years, 4 months ago) by jmc
Branch: MAIN
CVS Tags: checkpoint65z, checkpoint65x, checkpoint65y, checkpoint65r, checkpoint65s, checkpoint65p, checkpoint65q, checkpoint65v, checkpoint65w, checkpoint65t, checkpoint65u, checkpoint65j, checkpoint65k, checkpoint65h, checkpoint65i, checkpoint65n, checkpoint65l, checkpoint65m, checkpoint65b, checkpoint65c, checkpoint65a, checkpoint65f, checkpoint65g, checkpoint65d, checkpoint65e, checkpoint62c, checkpoint62a, checkpoint62g, checkpoint62f, checkpoint62e, checkpoint62d, checkpoint62k, checkpoint62j, checkpoint62i, checkpoint62h, checkpoint62o, checkpoint62n, checkpoint62m, checkpoint62l, checkpoint62s, checkpoint62r, checkpoint62q, checkpoint62p, checkpoint62w, checkpoint62v, checkpoint62u, checkpoint62t, checkpoint62z, checkpoint62y, checkpoint62x, checkpoint64, checkpoint65, checkpoint63, checkpoint66g, checkpoint66f, checkpoint66e, checkpoint66d, checkpoint66c, checkpoint66b, checkpoint66a, checkpoint66o, checkpoint66n, checkpoint66m, checkpoint66l, checkpoint66k, checkpoint66j, checkpoint66i, checkpoint66h, checkpoint63p, checkpoint63q, checkpoint63r, checkpoint63s, checkpoint63l, checkpoint63m, checkpoint63n, checkpoint63o, checkpoint63h, checkpoint63i, checkpoint63j, checkpoint63k, checkpoint63d, checkpoint63e, checkpoint63f, checkpoint63g, checkpoint63a, checkpoint63b, checkpoint63c, checkpoint65o, checkpoint62b, checkpoint64y, checkpoint64x, checkpoint64z, checkpoint64q, checkpoint64p, checkpoint64s, checkpoint64r, checkpoint64u, checkpoint64t, checkpoint64w, checkpoint64v, checkpoint64i, checkpoint64h, checkpoint64k, checkpoint64j, checkpoint64m, checkpoint64l, checkpoint64o, checkpoint64n, checkpoint64a, checkpoint64c, checkpoint64b, checkpoint64e, checkpoint64d, checkpoint64g, checkpoint64f, HEAD
Changes since 1.8: +24 -25 lines 
time-ave: use simpler (no level index) cumulative-time counter
1 C $Header: /u/gcmpack/MITgcm/pkg/bulk_force/bulkf_ave.F,v 1.8 2006/05/09 18:59:48 jmc Exp $
2 C $Name: $
3
4 #include "BULK_FORCE_OPTIONS.h"
5
6 CBOP
7 C !ROUTINE: BULKF_AVE
8 C !INTERFACE:
9 SUBROUTINE BULKF_AVE(bi,bj,myThid)
10
11 C !DESCRIPTION: \bv
12 C *==========================================================*
13 C | S/R BULKF_AVE
14 C | o save values for timeaveraging
15 C | for sea/ice surface and atmosphere
16 C *==========================================================*
17 C \ev
18
19 C !USES:
20 IMPLICIT NONE
21 C == Global variables ==
22 #include "SIZE.h"
23 #include "EEPARAMS.h"
24 #include "PARAMS.h"
25 #include "FFIELDS.h"
26 #include "BULKF.h"
27 #include "BULKF_INT.h"
28 #include "BULKF_TAVE.h"
29
30 C !INPUT/OUTPUT PARAMETERS:
31 INTEGER bi,bj
32 INTEGER myThid
33 CEOP
34
35 #ifdef ALLOW_BULK_FORCE
36 #ifdef ALLOW_TIMEAVE
37
38 INTEGER i,j
39
40 C-- Time-average
41 DO j = 1,sNy
42 DO i = 1,sNx
43 BULK_Qnet_Ave(i,j,bi,bj) = BULK_Qnet_Ave(i,j,bi,bj)
44 & +Qnet(i,j,bi,bj)*deltaTclock
45 c & +lambdaThetaClimRelax*
46 c & (theta(i,j,1,bi,bj)-SST(i,j,bi,bj))/
47 c & (recip_Cp*recip_rhoNil*recip_dRf(1))*deltaTclock
48 BULK_EmPmR_Ave(i,j,bi,bj) = BULK_EmPmR_Ave(i,j,bi,bj)
49 & +EmPmR(i,j,bi,bj)*deltaTclock
50 BULK_fu_Ave(i,j,bi,bj) = BULK_fu_Ave(i,j,bi,bj)
51 & +fu(i,j,bi,bj)*deltaTclock
52 BULK_fv_Ave(i,j,bi,bj) = BULK_fv_Ave(i,j,bi,bj)
53 & +fv(i,j,bi,bj)*deltaTclock
54 BULK_latent_Ave(i,j,bi,bj) = BULK_latent_Ave(i,j,bi,bj)
55 & +flh(i,j,bi,bj)*deltaTclock
56 BULK_sensible_Ave(i,j,bi,bj) = BULK_sensible_Ave(i,j,bi,bj)
57 & +fsh(i,j,bi,bj)*deltaTclock
58 BULK_evap_Ave(i,j,bi,bj) = BULK_evap_Ave(i,j,bi,bj)
59 & +evap(i,j,bi,bj)*deltaTclock
60 BULK_flwup_Ave(i,j,bi,bj) = BULK_flwup_Ave(i,j,bi,bj)
61 & +flwup(i,j,bi,bj)*deltaTclock
62 BULK_flwupnet_Ave(i,j,bi,bj) = BULK_flwupnet_Ave(i,j,bi,bj)
63 & +flwupnet(i,j,bi,bj)*deltaTclock
64 BULK_solar_Ave(i,j,bi,bj) = BULK_solar_Ave(i,j,bi,bj)
65 & +fswnet(i,j,bi,bj)*deltaTclock
66 BULK_ssq_Ave(i,j,bi,bj) = BULK_ssq_Ave(i,j,bi,bj)
67 & +savssq(i,j,bi,bj)*deltaTclock
68 ENDDO
69 ENDDO
70 BULK_timeAve(bi,bj) = BULK_timeAve(bi,bj)+deltaTclock
71
72 #endif /*ALLOW_TIMEAVE*/
73 #endif /*ALLOW_BULK_FORCE*/
74
75 RETURN
76 END
77
  ViewVC Help
Powered by ViewVC 1.1.22