/[MITgcm]/MITgcm_contrib/darwin2/pkg/darwin/dic_alk_surfforcing.F
ViewVC logotype

Annotation of /MITgcm_contrib/darwin2/pkg/darwin/dic_alk_surfforcing.F

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

Revision 1.2 - (hide annotations) (download)
Tue Sep 20 17:44:42 2011 UTC (14 years, 3 months ago) by jahn
Branch: MAIN
CVS Tags: ctrb_darwin2_ckpt64k_20130723, ctrb_darwin2_ckpt65w_20160512, ctrb_darwin2_ckpt65j_20150225, ctrb_darwin2_ckpt63l_20120405, ctrb_darwin2_ckpt66g_20170424, ctrb_darwin2_ckpt64h_20130528, ctrb_darwin2_ckpt66k_20171025, ctrb_darwin2_ckpt66n_20180118, ctrb_darwin2_ckpt65v_20160409, ctrb_darwin2_ckpt65s_20160114, ctrb_darwin2_ckpt65_20140718, ctrb_darwin2_ckpt64m_20130820, ctrb_darwin2_ckpt66d_20170214, ctrb_darwin2_ckpt64r_20131210, ctrb_darwin2_ckpt65m_20150615, ctrb_darwin2_ckpt65q_20151118, ctrb_darwin2_ckpt65o_20150914, ctrb_darwin2_ckpt64f_20130405, ctrb_darwin2_ckpt63f_20111201, ctrb_darwin2_ckpt64a_20121116, ctrb_darwin2_ckpt65p_20151023, ctrb_darwin2_ckpt64n_20130826, ctrb_darwin2_ckpt65e_20140929, ctrb_darwin2_ckpt64o_20131024, ctrb_darwin2_ckpt64v_20140411, ctrb_darwin2_ckpt64z_20140711, ctrb_darwin2_ckpt65l_20150504, ctrb_darwin2_ckpt65z_20160929, ctrb_darwin2_ckpt65n_20150729, ctrb_darwin2_ckpt64y_20140622, ctrb_darwin2_ckpt65d_20140915, ctrb_darwin2_ckpt64t_20140202, ctrb_darwin2_ckpt66h_20170602, ctrb_darwin2_ckpt64i_20130622, ctrb_darwin2_ckpt64s_20140105, ctrb_darwin2_ckpt64x_20140524, ctrb_darwin2_ckpt63o_20120629, ctrb_darwin2_ckpt64e_20130305, ctrb_darwin2_ckpt65x_20160612, ctrb_darwin2_ckpt66f_20170407, ctrb_darwin2_ckpt63c_20111011, ctrb_darwin2_ckpt63i_20120124, ctrb_darwin2_ckpt65g_20141120, ctrb_darwin2_ckpt63m_20120506, ctrb_darwin2_ckpt63s_20120908, ctrb_darwin2_ckpt65k_20150402, ctrb_darwin2_ckpt63e_20111107, ctrb_darwin2_ckpt64w_20140502, ctrb_darwin2_ckpt63j_20120217, ctrb_darwin2_ckpt66a_20161020, ctrb_darwin2_ckpt63r_20120817, ctrb_darwin2_ckpt64g_20130503, ctrb_darwin2_ckpt64l_20130806, ctrb_darwin2_ckpt63g_20111220, ctrb_darwin2_ckpt65f_20141014, ctrb_darwin2_ckpt64c_20130120, ctrb_darwin2_ckpt66b_20161219, ctrb_darwin2_ckpt64u_20140308, ctrb_darwin2_ckpt64j_20130704, ctrb_darwin2_ckpt65i_20150123, ctrb_darwin2_ckpt66j_20170815, ctrb_darwin2_ckpt65y_20160801, ctrb_darwin2_ckpt63h_20111230, ctrb_darwin2_ckpt63p_20120707, ctrb_darwin2_ckpt66c_20170121, ctrb_darwin2_ckpt65a_20140728, ctrb_darwin2_ckpt65b_20140812, ctrb_darwin2_ckpt65t_20160221, ctrb_darwin2_ckpt64p_20131118, ctrb_darwin2_ckpt63d_20111107, ctrb_darwin2_ckpt63q_20120731, ctrb_darwin2_ckpt64b_20121224, ctrb_darwin2_ckpt64d_20130219, ctrb_darwin2_ckpt66o_20180209, ctrb_darwin2_ckpt66e_20170314, ctrb_darwin2_ckpt64_20121012, ctrb_darwin2_ckpt64q_20131118, ctrb_darwin2_ckpt64p_20131024, ctrb_darwin2_ckpt65u_20160315, ctrb_darwin2_ckpt65r_20151221, ctrb_darwin2_ckpt66i_20170718, ctrb_darwin2_ckpt63n_20120604, ctrb_darwin2_ckpt63k_20120317, ctrb_darwin2_ckpt65c_20140830, ctrb_darwin2_ckpt66l_20171025, ctrb_darwin2_ckpt65h_20141217, ctrb_darwin2_ckpt66m_20171213, HEAD
Changes since 1.1: +6 -0 lines 
avoid uninitialized variable
1 jahn 1.1 #include "CPP_OPTIONS.h"
2     #include "PTRACERS_OPTIONS.h"
3     #include "DARWIN_OPTIONS.h"
4    
5     #ifdef ALLOW_PTRACERS
6     #ifdef ALLOW_DARWIN
7    
8     #ifdef ALLOW_CARBON
9    
10     CBOP
11     C !ROUTINE: ALK_SURFFORCING
12    
13     C !INTERFACE: ==========================================================
14     SUBROUTINE DIC_ALK_SURFFORCING( PTR_ALK , GALK,
15     I bi,bj,imin,imax,jmin,jmax,
16     I myIter,myTime,myThid)
17    
18     C !DESCRIPTION:
19     C Calculate the alkalinity change due to freshwater flux
20     C Update alkalinity tendency term
21    
22     C !USES: ===============================================================
23     IMPLICIT NONE
24     #include "SIZE.h"
25     #include "DYNVARS.h"
26     #include "EEPARAMS.h"
27     #include "PARAMS.h"
28     #include "GRID.h"
29     #include "FFIELDS.h"
30     #include "DARWIN_FLUX.h"
31    
32     C !INPUT PARAMETERS: ===================================================
33     C myThid :: thread number
34     C myIter :: current timestep
35     C myTime :: current time
36     C bi,bj :: tile indices
37     C PTR_ALK :: alkalinity field
38     INTEGER myIter, myThid
39     _RL myTime
40     _RL PTR_ALK(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)
41     INTEGER iMin,iMax,jMin,jMax, bi, bj
42    
43     C !OUTPUT PARAMETERS: ==================================================
44     C GALK :: tendency term of alkalinity
45     _RL GALK(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
46    
47    
48     C !LOCAL VARIABLES: ====================================================
49     C i,j :: loop indices
50     INTEGER i,j
51     CEOP
52    
53     #ifdef ALLOW_OLD_VIRTUALFLUX
54     DO j=jmin,jmax
55     DO i=imin,imax
56     IF (maskC(i,j,1,bi,bj).NE.0. _d 0) THEN
57     c calculate virtual flux
58     c EminusPforV = dS/dt*(1/Sglob)
59     C NOTE: Be very careful with signs here!
60     C Positive EminusPforV => loss of water to atmos and increase
61     C in salinity. Thus, also increase in other surface tracers
62     C (i.e. positive virtual flux into surface layer)
63     GALK(i,j)=gsm_ALK*surfaceForcingS(i,j,bi,bj)*
64     & recip_drF(1)/gsm_S
65     c OR
66     c let virtual flux be zero
67     c GALK(i,j)=0.d0
68     c
69     ELSE
70     GALK(i,j)=0. _d 0
71     ENDIF
72     ENDDO
73     ENDDO
74 jahn 1.2 #else
75     DO j=jmin,jmax
76     DO i=imin,imax
77     GALK(i,j)=0. _d 0
78     ENDDO
79     ENDDO
80 jahn 1.1 #endif /* ALLOW_OLD_VIRTUALFLUX */
81    
82     RETURN
83     END
84     #endif /*ALLOW_CARBON*/
85    
86     #endif /*DARWIN*/
87     #endif /*ALLOW_PTRACERS*/
88     c ==================================================================
  ViewVC Help
Powered by ViewVC 1.1.22