/[MITgcm]/MITgcm/pkg/aim/aim_calc_diags.F
ViewVC logotype

Contents of /MITgcm/pkg/aim/aim_calc_diags.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph

Revision 1.6 - (show annotations) (download)
Tue Sep 25 19:53:57 2001 UTC (22 years, 8 months ago) by jmc
Branch: MAIN
CVS Tags: icebear2, checkpoint44h_pre, release1_p12, release1_p10, release1_p11, release1_p16, release1_p15, ecco_c44_e17, ecco_c44_e16, ecco_c44_e21, ecco_c44_e26, ecco_c44_e27, ecco_c44_e24, icebear5, icebear4, checkpoint44f_pre, icebear3, checkpoint46f_post, release1_p13_pre, checkpoint46d_pre, checkpoint46e_post, release1-branch_tutorials, release1_p14, checkpoint44g_post, checkpoint46h_pre, checkpoint44h_post, release1_p12_pre, checkpoint44e_post, checkpoint46e_pre, ecco-branch-mod4, checkpoint43a-release1mods, checkpoint45d_post, checkpoint45b_post, checkpoint46b_pre, chkpt44a_pre, release1-branch-end, release1_final_v1, ecco_c44_e19, checkpoint46, ecco_c44_e20, checkpoint44, release1_p13, ecco_c44_e18, checkpoint44f_post, release1_p17, release1_b1, checkpoint44b_post, chkpt44c_post, chkpt44d_post, checkpoint42, release1_p9, release1_p8, checkpoint43, checkpoint46g_pre, release1_p2, release1_p3, release1_p4, release1_p6, checkpoint46a_post, chkpt44a_post, checkpoint44b_pre, release1_p1, checkpoint46a_pre, ecco-branch-mod1, checkpoint45c_post, release1_p5, checkpoint44e_pre, chkpt44c_pre, release1_p7, ecco_ice2, ecco_ice1, checkpoint46b_post, checkpoint46d_post, ecco-branch-mod2, checkpoint46g_post, checkpoint45a_post, checkpoint46c_pre, ecco-branch-mod3, ecco-branch-mod5, ecco_c44_e22, release1_beta1, ecco_c44_e23, release1-branch_branchpoint, checkpoint46c_post, checkpoint45, checkpoint46h_post, release1_chkpt44d_post, ecco_c44_e25
Branch point for: c24_e25_ice, ecco-branch, release1_coupled, icebear, release1_final, release1-branch, release1, release1_50yr
Changes since 1.5: +3 -3 lines 
add a CPP option to turn back to old AIM Interface.
diagnostic of surface stress consistent with dynamical effects.
1 C $Header: /u/gcmpack/models/MITgcmUV/pkg/aim/aim_calc_diags.F,v 1.5 2001/05/29 19:28:53 cnh Exp $
2 C $Name: $
3
4 #include "AIM_OPTIONS.h"
5
6 CStartOfInterFace
7 SUBROUTINE AIM_CALC_DIAGS( bi, bj, myTime, myThid )
8 C /==========================================================\
9 C | SUBROUTINE AIM_CALC_DIAGS |
10 C | o Calculate AIM diagnostics |
11 C \==========================================================/
12 IMPLICIT rEAL*8 (A-H,O-Z)
13
14 C === Global variables ===
15 #include "SIZE.h"
16 #include "EEPARAMS.h"
17 #include "PARAMS.h"
18 #include "AIM_DIAGS.h"
19
20 #ifdef ALLOW_AIM
21 #include "atparam0.h"
22 #include "atparam1.h"
23 INTEGER NGP
24 INTEGER NLON
25 INTEGER NLAT
26 INTEGER NLEV
27 PARAMETER ( NLON=IX, NLAT=IL, NLEV=KX, NGP=NLON*NLAT )
28 #include "com_physvar.h"
29 #endif /* ALLOW_AIM */
30
31 LOGICAL DIFFERENT_MULTIPLE
32 EXTERNAL DIFFERENT_MULTIPLE
33
34 C == Routine arguments ==
35 C myTime - Current time of simulation ( s )
36 C myThid - Number of this instance of the routine
37 C bi,bj - Tile index
38 _RL myTime
39 INTEGER myThid
40 INTEGER bi, bj
41 CEndOfInterface
42
43 #ifdef ALLOW_AIM
44
45 C == Local variables ==
46 INTEGER I,J,I2,K
47 _RL DDTT
48
49 #ifdef ALLOW_TIMEAVE
50
51 DDTT = deltaTclock
52
53 DO J=1,sNy
54 DO I=1,sNx
55 I2 = sNx*(J-1)+I
56 c USTRtave(i,j,bi,bj) = USTRtave(i,j,bi,bj) + USTR(I2,3,myThid)*DDTT
57 c VSTRtave(i,j,bi,bj) = VSTRtave(i,j,bi,bj) + VSTR(I2,3,myThid)*DDTT
58 TSRtave(i,j,bi,bj) = TSRtave(i,j,bi,bj) + TSR(I2,myThid)*DDTT
59 OLRtave(i,j,bi,bj) = OLRtave(i,j,bi,bj) + OLR(I2,myThid)*DDTT
60 SSRtave(i,j,bi,bj) = SSRtave(i,j,bi,bj) + SSR(I2,myThid)*DDTT
61 SLRtave(i,j,bi,bj) = SLRtave(i,j,bi,bj) + SLR(I2,myThid)*DDTT
62 SHFtave(i,j,bi,bj) = SHFtave(i,j,bi,bj) + SHF(I2,3,myThid)*DDTT
63 EVAPtave(i,j,bi,bj) = EVAPtave(i,j,bi,bj) + EVAP(I2,3,myThid)*DDTT
64 PRECNVtave(i,j,bi,bj) = PRECNVtave(i,j,bi,bj) + PRECNV(I2,myThid)*DDTT
65 PRECLStave(i,j,bi,bj) = PRECLStave(i,j,bi,bj) + PRECLS(I2,myThid)*DDTT
66 CLOUDCtave(i,j,bi,bj) = CLOUDCtave(i,j,bi,bj) + CLOUDC(I2,myThid)*DDTT
67 ENDDO
68 ENDDO
69
70 C- Keep record of how much time has been integrated over
71 DO K=1,Nr
72 AIM_TimeAve(k,bi,bj)=AIM_TimeAve(k,bi,bj)+DDTT
73 ENDDO
74
75
76 #endif /* ALLOW_TIMEAVE */
77
78 #endif /* ALLOW_AIM */
79
80 RETURN
81 END
  ViewVC Help
Powered by ViewVC 1.1.22