/[MITgcm]/MITgcm/pkg/frazil/frazil_tendency_apply.F
ViewVC logotype

Annotation of /MITgcm/pkg/frazil/frazil_tendency_apply.F

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

Revision 1.5 - (hide annotations) (download)
Wed Jul 9 17:00:49 2014 UTC (9 years, 11 months ago) by jmc
Branch: MAIN
CVS Tags: checkpoint64z, checkpoint66g, checkpoint66f, checkpoint66e, checkpoint66d, checkpoint66c, checkpoint66b, checkpoint66a, checkpoint66o, checkpoint66n, checkpoint66m, checkpoint66l, checkpoint66k, checkpoint66j, checkpoint66i, checkpoint66h, checkpoint65z, checkpoint65x, checkpoint65y, checkpoint65r, checkpoint65s, checkpoint65p, checkpoint65q, checkpoint65v, checkpoint65w, checkpoint65t, checkpoint65u, checkpoint65j, checkpoint65k, checkpoint65h, checkpoint65i, checkpoint65n, checkpoint65o, checkpoint65l, checkpoint65m, checkpoint65b, checkpoint65c, checkpoint65a, checkpoint65f, checkpoint65g, checkpoint65d, checkpoint65e, checkpoint65, HEAD
Changes since 1.4: +19 -9 lines 
- change argument list of all S/R called from external_forcing.F
  to pass, as argument, the current level dendency array to update
  (instead of a direct update of the common bloc array).
1 jmc 1.5 C $Header: /u/gcmpack/MITgcm/pkg/frazil/frazil_tendency_apply.F,v 1.4 2014/05/06 15:50:14 jmc Exp $
2 dimitri 1.1 C $Name: $
3    
4     #include "FRAZIL_OPTIONS.h"
5     C-- File frazil_tendency_apply.F: Routines to apply FRAZIL tendencies
6     C-- Contents
7     C-- o FRAZIL_TENDENCY_APPLY_T
8    
9     C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
10     CBOP 0
11    
12     C !INTERFACE:
13     SUBROUTINE FRAZIL_TENDENCY_APPLY_T(
14 jmc 1.5 U gT_arr,
15     I iMin,iMax,jMin,jMax, k, bi, bj,
16     I myTime, myIter, myThid )
17 dimitri 1.1
18     C !DESCRIPTION:
19 dimitri 1.3 C Add frazil tendency terms to T (theta) tendency.
20     C Routine works for one level at a time.
21 dimitri 1.1
22     C !USES:
23     IMPLICIT NONE
24 jmc 1.5 #include "SIZE.h"
25 dimitri 1.3 #include "EEPARAMS.h"
26     #include "PARAMS.h"
27 dimitri 1.1 #include "GRID.h"
28 jmc 1.5 c#include "DYNVARS.h"
29 dimitri 1.3 #ifdef ALLOW_FRAZIL
30     # include "FRAZIL.h"
31     #endif
32 dimitri 1.1
33 dimitri 1.3 C !INPUT/OUTPUT PARAMETERS:
34 jmc 1.5 C gT_arr :: the tendency array
35     C iMin,iMax :: Working range of x-index for applying forcing.
36     C jMin,jMax :: Working range of y-index for applying forcing.
37     C k :: Current vertical level index
38     C bi,bj :: Current tile indices
39     C myTime :: Current time in simulation
40     C myIter :: Current iteration number
41     C myThid :: my Thread Id number
42     _RL gT_arr(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
43     INTEGER iMin, iMax, jMin, jMax
44     INTEGER k, bi, bj
45 dimitri 1.3 _RL myTime
46 jmc 1.5 INTEGER myIter
47 dimitri 1.3 INTEGER myThid
48 dimitri 1.1 CEOP
49    
50     #ifdef ALLOW_FRAZIL
51 dimitri 1.3
52 dimitri 1.1 C !LOCAL VARIABLES:
53     INTEGER i, j
54 jmc 1.4 _RL recip_Cp
55 dimitri 1.1
56 jmc 1.4 recip_Cp = 1. _d 0 / HeatCapacity_Cp
57 dimitri 1.2 DO j=jMin,jMax
58     DO i=iMin,iMax
59 jmc 1.5 gT_arr(i,j) = gT_arr(i,j)
60 dimitri 1.3 & + FrazilForcingT(i,j,k,bi,bj)
61     & * recip_Cp * mass2rUnit
62     & * recip_drF(k) * _recip_hFacC(i,j,k,bi,bj)
63 dimitri 1.2 ENDDO
64     ENDDO
65 dimitri 1.3
66 dimitri 1.1 #endif /* ALLOW_FRAZIL */
67    
68     RETURN
69     END
  ViewVC Help
Powered by ViewVC 1.1.22