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