Parent Directory
|
Revision Log
|
Revision Graph
o first check-in for the "branch-genmake2" merge o verification suite as run on shelley (gcc 3.2.2): Wed Oct 8 23:42:29 EDT 2003 T S U V G D M c m s m s m s m s E p a R g m m e . m m e . m m e . m m e . N n k u 2 i a a d i a a d i a a d i a a d 2 d e n d n x n . n x n . n x n . n x n . OPTFILE=NONE Y Y Y Y 13 16 16 16 0 16 16 16 16 16 16 16 16 13 12 0 0 pass adjustment.128x64x1 Y Y Y Y 16 16 16 16 0 16 16 16 16 16 16 0 0 16 16 0 0 pass adjustment.cs-32x32x1 Y Y Y Y 16 16 16 16 0 16 16 16 16 16 16 22 0 16 16 22 0 pass adjust_nlfs.cs-32x32x1 Y Y Y Y -- 13 13 16 16 13 13 13 13 16 16 16 16 16 16 16 16 N/O advect_cs Y Y Y Y -- 22 16 16 16 16 16 16 13 16 16 16 16 16 16 16 16 N/O advect_xy Y Y Y Y -- 13 16 13 16 16 16 16 16 16 16 22 16 16 16 16 16 N/O advect_xz Y Y Y Y 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 pass aim.5l_cs Y Y Y Y 14 16 16 16 16 16 16 16 16 13 16 16 16 16 16 13 16 pass aim.5l_Equatorial_Channel Y Y Y Y 16 16 16 16 16 16 16 16 16 16 16 13 16 16 13 13 16 pass aim.5l_LatLon Y Y Y Y 13 16 16 16 16 16 16 16 16 16 13 12 13 13 16 13 16 pass exp0 Y Y Y Y 14 16 16 16 16 16 16 16 22 16 16 16 13 16 16 22 16 pass exp1 Y Y Y Y 13 13 16 13 16 16 16 16 16 13 13 16 16 13 13 13 13 pass exp2 Y Y Y Y 16 16 16 16 16 16 16 16 22 16 16 16 16 16 16 16 16 pass exp4 Y Y Y Y 16 16 16 16 16 16 16 16 16 16 16 22 16 16 16 22 16 pass exp5 Y Y Y Y 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 pass front_relax Y Y Y Y 14 16 16 13 13 16 16 13 13 16 13 13 16 12 13 13 16 pass global_ocean.90x40x15 Y Y Y Y 10 16 16 13 13 16 13 16 16 13 13 13 13 16 16 13 16 FAIL global_ocean.cs32x15 Y Y Y Y 6 11 12 13 13 12 13 16 13 9 9 9 9 10 9 9 11 FAIL global_ocean_pressure Y Y Y Y 14 16 16 13 16 16 16 13 13 13 13 13 16 12 16 13 16 pass global_with_exf Y Y Y Y 14 16 16 16 16 16 16 16 16 11 13 22 13 16 16 9 16 pass hs94.128x64x5 Y Y Y Y 13 16 16 16 16 16 16 16 16 11 16 16 16 13 16 22 13 pass hs94.1x64x5 Y Y Y Y 14 16 16 16 16 16 16 16 16 13 16 13 13 16 16 22 13 pass hs94.cs-32x32x5 Y Y Y Y 10 10 16 13 13 16 16 16 22 16 13 13 13 13 13 22 13 FAIL ideal_2D_oce Y Y Y Y 8 16 16 16 16 16 16 16 16 13 13 8 16 16 16 16 16 FAIL internal_wave Y Y Y Y 14 16 16 16 16 16 16 16 16 13 13 22 13 13 13 22 16 pass inverted_barometer Y Y Y Y 12 16 16 16 16 16 16 16 16 16 13 12 13 13 13 13 13 FAIL lab_sea Y Y Y Y 11 16 16 16 16 16 16 16 13 13 13 12 13 16 13 12 13 FAIL natl_box Y Y Y Y 16 16 16 16 16 16 16 16 22 16 16 16 16 16 16 16 16 pass plume_on_slope Y Y Y Y 13 16 16 16 16 13 16 16 16 16 16 16 16 13 16 16 16 pass solid-body.cs-32x32x1
1 | C $Header: /u/u3/gcmpack/MITgcm/pkg/bulk_force/bulkf_ave.F,v 1.4.4.1 2003/10/02 18:18:33 adcroft Exp $ |
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 "DYNVARS.h" |
20 | #include "EEPARAMS.h" |
21 | #include "PARAMS.h" |
22 | #include "GRID.h" |
23 | #ifdef ALLOW_BULK_FORCE |
24 | #include "FFIELDS.h" |
25 | #include "BULKF.h" |
26 | #include "BULKF_INT.h" |
27 | #include "BULKF_DIAG.h" |
28 | #endif |
29 | c == routine arguments == |
30 | |
31 | integer mythid |
32 | integer bi,bj |
33 | integer i,j,k |
34 | |
35 | |
36 | #ifdef ALLOW_BULK_FORCE |
37 | #ifdef ALLOW_TIMEAVE |
38 | |
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 | return |
76 | end |
77 |
ViewVC Help | |
Powered by ViewVC 1.1.22 |