/[MITgcm]/MITgcm/pkg/ptracers/ptracers_forcing.F
ViewVC logotype

Contents of /MITgcm/pkg/ptracers/ptracers_forcing.F

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

Revision 1.5 - (show annotations) (download)
Sun Oct 26 01:10:34 2003 UTC (20 years, 8 months ago) by heimbach
Branch: MAIN
CVS Tags: checkpoint52l_pre, hrcube4, hrcube5, checkpoint52d_pre, checkpoint52j_pre, checkpoint52k_post, checkpoint54, checkpoint53, checkpoint52, checkpoint52f_post, checkpoint51t_post, checkpoint52i_pre, hrcube_1, hrcube_2, hrcube_3, checkpoint51s_post, checkpoint52e_pre, checkpoint52e_post, checkpoint53d_post, checkpoint52b_pre, checkpoint54b_post, checkpoint52m_post, checkpoint51q_post, checkpoint52b_post, checkpoint52c_post, checkpoint52h_pre, checkpoint52f_pre, checkpoint54a_pre, checkpoint53c_post, checkpoint54a_post, checkpoint51r_post, checkpoint53a_post, checkpoint52d_post, checkpoint53g_post, checkpoint52a_pre, checkpoint52i_post, checkpoint53f_post, checkpoint52j_post, branch-netcdf, checkpoint52l_post, checkpoint52n_post, checkpoint53b_pre, checkpoint51o_post, checkpoint53b_post, checkpoint52a_post, ecco_c52_e35, checkpoint53d_pre, checkpoint51p_post, checkpoint51u_post
Branch point for: branch-nonh, netcdf-sm0
Changes since 1.4: +6 -7 lines 
o initialisation of rFlx extended to full array (required by TAF)
  and shifted to thermodynamics
o removed PTRACERS.h in ptracers routine
o added surfacetendencyPtr to S/R parameter list pracers_forcing
1 C $Header: /u/gcmpack/MITgcm/pkg/ptracers/ptracers_forcing.F,v 1.4 2003/09/24 04:52:39 dimitri Exp $
2 C $Name: $
3
4 #include "PTRACERS_OPTIONS.h"
5
6 CBOP
7 C !ROUTINE: PTRACERS_FORCING
8
9 C !INTERFACE: ==========================================================
10 SUBROUTINE PTRACERS_FORCING(
11 I bi,bj,iMin,iMax,jMin,jMax,k,
12 U gPtracer,surfPtracer,
13 I myIter,myTime,myThid )
14
15 C !DESCRIPTION:
16 C Adds sources and sinks of passive tracers to the tendancy arrays
17
18 C !USES: ===============================================================
19 IMPLICIT NONE
20 #include "SIZE.h"
21 #include "EEPARAMS.h"
22 #include "PARAMS.h"
23 #include "FFIELDS.h"
24 #include "DYNVARS.h"
25 #include "GRID.h"
26
27 C !INPUT PARAMETERS: ===================================================
28 C bi,bj :: tile indices
29 C iMin iMax jMin jMax :: working range of tile for applying forcing
30 C k :: vertical level number
31 C gPtracer :: the tendancy array
32 C myIter :: time-step number
33 C myTime :: model time
34 C myThid :: thread number
35 INTEGER bi,bj,iMin, iMax, jMin, jMax,k
36 _RL gPtracer(1-Olx:sNx+Olx,1-Oly:sNy+Oly,Nr,nSx,nSy)
37 _RL surfPtracer(1-Olx:sNx+Olx,1-Oly:sNy+Oly,nSx,nSy)
38 INTEGER myIter
39 _RL myTime
40 INTEGER myThid
41
42 C !OUTPUT PARAMETERS: ==================================================
43 C gPtracer :: updates tendancy array
44
45 #ifdef ALLOW_PTRACERS
46
47 C !LOCAL VARIABLES: ====================================================
48 C i,j :: loop indices
49 INTEGER i,j
50 C number of surface interface layer
51 INTEGER kSurface
52 CEOP
53
54 if ( buoyancyRelation .eq. 'OCEANICP' ) then
55 kSurface = Nr
56 else
57 kSurface = 1
58 endif
59
60 C-- Surface forcing term is surfaceTendencyPtr, previously computed
61 C by PTRACERS_FORCING because it is needed by KPP_TRANSPORT_PTR.
62
63 IF ( k .EQ. kSurface ) THEN
64 DO j=jMin,jMax
65 DO i=iMin,iMax
66 gPtracer(i,j,k,bi,bj) = gPtracer(i,j,k,bi,bj) +
67 & maskC(i,j,k,bi,bj) *
68 & surfPtracer(i,j,bi,bj)
69 ENDDO
70 ENDDO
71 ENDIF
72
73 #endif /* ALLOW_PTRACERS */
74
75 RETURN
76 END
  ViewVC Help
Powered by ViewVC 1.1.22