/[MITgcm]/MITgcm/pkg/layers/layers_calc.F
ViewVC logotype

Diff of /MITgcm/pkg/layers/layers_calc.F

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

revision 1.5 by dfer, Sat Dec 4 23:50:32 2010 UTC revision 1.6 by jmc, Sun Dec 5 14:17:39 2010 UTC
# Line 2  C $Header$ Line 2  C $Header$
2  C $Name$  C $Name$
3    
4  #include "LAYERS_OPTIONS.h"  #include "LAYERS_OPTIONS.h"
5    #ifdef ALLOW_GMREDI
6  #include "GMREDI_OPTIONS.h"  #include "GMREDI_OPTIONS.h"
7    #endif
8    
9  CBOP 0  CBOP 0
10  C !ROUTINE: LAYERS_CALC  C !ROUTINE: LAYERS_CALC
# Line 26  C !USES: Line 28  C !USES:
28  #include "DYNVARS.h"  #include "DYNVARS.h"
29  #include "LAYERS_SIZE.h"  #include "LAYERS_SIZE.h"
30  #include "LAYERS.h"  #include "LAYERS.h"
31  #if (defined ALLOW_GMREDI) && (defined GM_BOLUS_ADVEC)  #ifdef ALLOW_GMREDI
32  #include "GMREDI.h"  # include "GMREDI.h"
33  #endif  #endif
34    
35  C !INPUT PARAMETERS:  C !INPUT PARAMETERS:
# Line 53  C     TatV     :: temperature at U point Line 55  C     TatV     :: temperature at U point
55  C     TatV     :: temperature at V point  C     TatV     :: temperature at V point
56    
57        INTEGER bi, bj        INTEGER bi, bj
58        INTEGER i,j,k,kk,kg,kci,kcip1        INTEGER i,j,k,kk,kg,kci
59        INTEGER kgu(sNx+1,sNy+1), kgv(sNx+1,sNy+1)        INTEGER kgu(sNx+1,sNy+1), kgv(sNx+1,sNy+1)
60        _RL TatU, TatV, delPsi, maskp1        _RL TatU, TatV
61        CHARACTER*(MAX_LEN_MBUF) msgBuf        CHARACTER*(MAX_LEN_MBUF) msgBuf
62    #if (defined ALLOW_GMREDI) && (defined GM_BOLUS_ADVEC)
63          INTEGER kcip1
64          _RL delPsi, maskp1
65    #endif
66    
67  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
68    
# Line 164  C ------ Augment the bin values Line 170  C ------ Augment the bin values
170             delPsi = GM_PsiX(i,j,kcip1,bi,bj)*maskp1             delPsi = GM_PsiX(i,j,kcip1,bi,bj)*maskp1
171       &            - GM_PsiX(i,j, kci, bi,bj)       &            - GM_PsiX(i,j, kci, bi,bj)
172             layers_UFlux(i,j,kgu(i,j),bi,bj) =             layers_UFlux(i,j,kgu(i,j),bi,bj) =
173       &      layers_UFlux(i,j,kgu(i,j),bi,bj)       &      layers_UFlux(i,j,kgu(i,j),bi,bj)
174       &      + delPsi*recip_drF(kci)*_recip_hFacW(i,j,kci,bi,bj)       &      + delPsi*recip_drF(kci)*_recip_hFacW(i,j,kci,bi,bj)
175       &      * dZZf(kk)*hFacW(i,j,kci,bi,bj)       &      * dZZf(kk)*hFacW(i,j,kci,bi,bj)
176           ENDIF           ENDIF
Legend:
Removed from v.1.5  
changed lines
  Added in v.1.6
  ViewVC Help
Powered by ViewVC 1.1.22