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 | edhill | 1.5 | 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 | cheisey | 1.1 | |
6 | C !ROUTINE: BULKF_AVE | ||
7 | C !INTERFACE: | ||
8 | cheisey | 1.4 | SUBROUTINE BULKF_AVE(bi,bj,mythid) |
9 | cheisey | 1.1 | 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 | cheisey | 1.4 | #include "BULKF_INT.h" |
27 | cheisey | 1.1 | #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 | cheisey | 1.4 | & +flh(i,j,bi,bj)*deltaTclock |
56 | cheisey | 1.1 | BULK_sensible_AVE(i,j,bi,bj)=BULK_sensible_AVE(i,j,bi,bj) |
57 | cheisey | 1.4 | & +fsh(i,j,bi,bj)*deltaTclock |
58 | cheisey | 1.1 | BULK_evap_AVE(i,j,bi,bj)=BULK_evap_AVE(i,j,bi,bj) |
59 | cheisey | 1.4 | & +evap(i,j,bi,bj)*deltaTclock |
60 | cheisey | 1.3 | BULK_flwup_AVE(i,j,bi,bj)=BULK_flwup_AVE(i,j,bi,bj) |
61 | cheisey | 1.4 | & +flwup(i,j,bi,bj)*deltaTclock |
62 | cheisey | 1.3 | BULK_flwupnet_AVE(i,j,bi,bj)=BULK_flwupnet_AVE(i,j,bi,bj) |
63 | cheisey | 1.4 | & +flwupnet(i,j,bi,bj)*deltaTclock |
64 | cheisey | 1.1 | BULK_solar_AVE(i,j,bi,bj)=BULK_solar_AVE(i,j,bi,bj) |
65 | cheisey | 1.4 | & +fswnet(i,j,bi,bj)*deltaTclock |
66 | cheisey | 1.1 | BULK_ssq_AVE(i,j,bi,bj)=BULK_ssq_AVE(i,j,bi,bj) |
67 | cheisey | 1.4 | & +savssq(i,j,bi,bj)*deltaTclock |
68 | cheisey | 1.1 | enddo |
69 | enddo | ||
70 | cheisey | 1.2 | BULKF_timeave(1,bi,bj)=BULKF_timeave(1,bi,bj)+deltaTclock |
71 | cheisey | 1.1 | |
72 | |||
73 | #endif /*ALLOW_TIMEAVE*/ | ||
74 | #endif /*ALLOW_BULK_FORCE*/ | ||
75 | return | ||
76 | end | ||
77 |
ViewVC Help | |
Powered by ViewVC 1.1.22 |