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

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

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

Revision 1.6 - (hide 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 jmc 1.6 C $Header: /u/gcmpack/models/MITgcmUV/pkg/aim/aim_calc_diags.F,v 1.5 2001/05/29 19:28:53 cnh Exp $
2 jmc 1.3 C $Name: $
3 adcroft 1.2
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 cnh 1.5 IMPLICIT rEAL*8 (A-H,O-Z)
13 adcroft 1.2
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 cnh 1.5 INTEGER I,J,I2,K
47 adcroft 1.2 _RL DDTT
48    
49 jmc 1.3 #ifdef ALLOW_TIMEAVE
50 adcroft 1.2
51     DDTT = deltaTclock
52    
53     DO J=1,sNy
54     DO I=1,sNx
55     I2 = sNx*(J-1)+I
56 jmc 1.6 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 cnh 1.5 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 adcroft 1.2 ENDDO
68     ENDDO
69 adcroft 1.4
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 adcroft 1.2
75    
76 jmc 1.3 #endif /* ALLOW_TIMEAVE */
77 adcroft 1.2
78     #endif /* ALLOW_AIM */
79    
80     RETURN
81     END
  ViewVC Help
Powered by ViewVC 1.1.22