/[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.8 - (show annotations) (download)
Tue May 9 18:59:48 2006 UTC (18 years, 4 months ago) by jmc
Branch: MAIN
CVS Tags: checkpoint59, checkpoint58f_post, checkpoint58y_post, checkpoint58t_post, checkpoint58m_post, checkpoint60, checkpoint61, checkpoint62, checkpoint58w_post, checkpoint58o_post, checkpoint58p_post, checkpoint58q_post, checkpoint58e_post, mitgcm_mapl_00, checkpoint58r_post, checkpoint58n_post, checkpoint59q, checkpoint59p, checkpoint59r, checkpoint59e, checkpoint59d, checkpoint59g, checkpoint59f, checkpoint59a, checkpoint59c, checkpoint59b, checkpoint59m, checkpoint59l, checkpoint59o, checkpoint59n, checkpoint59i, checkpoint59h, checkpoint59k, checkpoint59j, checkpoint58k_post, checkpoint58v_post, checkpoint58l_post, checkpoint61f, checkpoint58g_post, checkpoint58x_post, checkpoint61n, checkpoint58h_post, checkpoint58j_post, checkpoint61q, checkpoint61z, checkpoint61e, checkpoint58i_post, checkpoint58u_post, checkpoint58s_post, checkpoint61g, checkpoint61d, checkpoint61b, checkpoint61c, checkpoint61a, checkpoint61o, checkpoint61l, checkpoint61m, checkpoint61j, checkpoint61k, checkpoint61h, checkpoint61i, checkpoint61v, checkpoint61w, checkpoint61t, checkpoint61u, checkpoint61r, checkpoint61s, checkpoint61p, checkpoint61x, checkpoint61y
Changes since 1.7: +18 -12 lines 
rename BULKF_DIAG.h -> BULKF_TAVE.h (better name for time-average fields)
1 C $Header: /u/gcmpack/MITgcm/pkg/bulk_force/bulkf_ave.F,v 1.7 2005/11/04 01:28:38 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 BULKF_timeave(1,bi,bj)=BULKF_timeave(1,bi,bj)+deltaTclock
71
72
73 #endif /*ALLOW_TIMEAVE*/
74 #endif /*ALLOW_BULK_FORCE*/
75
76 RETURN
77 END
78
  ViewVC Help
Powered by ViewVC 1.1.22