1 |
C $Header: $ |
2 |
C $Name: $ |
3 |
|
4 |
#include "BULK_FORCE_OPTIONS.h" |
5 |
|
6 |
C !ROUTINE: BULKF_AVE |
7 |
C !INTERFACE: |
8 |
SUBROUTINE BULKF_AVE(bi,bj,mythid) |
9 |
C *==========================================================* |
10 |
C | S/R BULKF_AVE |
11 |
C | o save values for timeaveraging |
12 |
C | for sea/ice surface and atmosphere |
13 |
C *==========================================================* |
14 |
|
15 |
C !USES: |
16 |
IMPLICIT NONE |
17 |
C == Global variables == |
18 |
#include "SIZE.h" |
19 |
#include "EEPARAMS.h" |
20 |
#include "PARAMS.h" |
21 |
#include "FFIELDS.h" |
22 |
#include "BULKF.h" |
23 |
#include "BULKF_INT.h" |
24 |
#include "BULKF_DIAG.h" |
25 |
c == routine arguments == |
26 |
|
27 |
integer mythid |
28 |
integer bi,bj |
29 |
|
30 |
#ifdef ALLOW_BULK_FORCE |
31 |
#ifdef ALLOW_TIMEAVE |
32 |
|
33 |
integer i,j,k |
34 |
|
35 |
C-- Time-average |
36 |
do j = 1,sny |
37 |
do i = 1,snx |
38 |
BULK_Qnet_AVE(i,j,bi,bj)=BULK_Qnet_AVE(i,j,bi,bj) |
39 |
& +Qnet(i,j,bi,bj)*deltaTclock |
40 |
c & +lambdaThetaClimRelax* |
41 |
c & (theta(i,j,1,bi,bj)-SST(i,j,bi,bj))/ |
42 |
c & (recip_Cp*recip_rhoNil*recip_dRf(1))*deltaTclock |
43 |
BULK_EmPmR_AVE(i,j,bi,bj)=BULK_EmPmR_AVE(i,j,bi,bj) |
44 |
& +EmPmR(i,j,bi,bj)*deltaTclock |
45 |
BULK_fu_AVE(i,j,bi,bj)=BULK_fu_AVE(i,j,bi,bj) |
46 |
& +fu(i,j,bi,bj)*deltaTclock |
47 |
BULK_fv_AVE(i,j,bi,bj)=BULK_fv_AVE(i,j,bi,bj) |
48 |
& +fv(i,j,bi,bj)*deltaTclock |
49 |
BULK_latent_AVE(i,j,bi,bj)=BULK_latent_AVE(i,j,bi,bj) |
50 |
& +flh(i,j,bi,bj)*deltaTclock |
51 |
BULK_sensible_AVE(i,j,bi,bj)=BULK_sensible_AVE(i,j,bi,bj) |
52 |
& +fsh(i,j,bi,bj)*deltaTclock |
53 |
BULK_evap_AVE(i,j,bi,bj)=BULK_evap_AVE(i,j,bi,bj) |
54 |
& +evap(i,j,bi,bj)*deltaTclock |
55 |
BULK_flwup_AVE(i,j,bi,bj)=BULK_flwup_AVE(i,j,bi,bj) |
56 |
& +flwup(i,j,bi,bj)*deltaTclock |
57 |
BULK_flwupnet_AVE(i,j,bi,bj)=BULK_flwupnet_AVE(i,j,bi,bj) |
58 |
& +flwupnet(i,j,bi,bj)*deltaTclock |
59 |
BULK_solar_AVE(i,j,bi,bj)=BULK_solar_AVE(i,j,bi,bj) |
60 |
& +fswnet(i,j,bi,bj)*deltaTclock |
61 |
BULK_ssq_AVE(i,j,bi,bj)=BULK_ssq_AVE(i,j,bi,bj) |
62 |
& +savssq(i,j,bi,bj)*deltaTclock |
63 |
enddo |
64 |
enddo |
65 |
BULKF_timeave(1,bi,bj)=BULKF_timeave(1,bi,bj)+deltaTclock |
66 |
|
67 |
|
68 |
#endif /*ALLOW_TIMEAVE*/ |
69 |
#endif /*ALLOW_BULK_FORCE*/ |
70 |
RETURN |
71 |
END |
72 |
|