/[MITgcm]/MITgcm/model/src/cycle_tracer.F
ViewVC logotype

Annotation of /MITgcm/model/src/cycle_tracer.F

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

Revision 1.4 - (hide annotations) (download)
Wed Sep 19 02:43:27 2001 UTC (23 years ago) by adcroft
Branch: MAIN
CVS Tags: checkpoint40
Changes since 1.3: +3 -4 lines 
Re-arranged sequence of operations for Adams-Bashforth
 o this does not change numbers
 o this makes it very easy to extract forcing/diffusion out of ABII
   by changing calling sequence in calc_gt, calc_gs,...

Key modifications:
 o new s/r: ADAMS_BASHFORTH2  gT=3/2*gT-1/2*gTnm1
 o changed TIMESTEP_TRACER from gTnm1=t+dt*(3/2*gT-1/2*gTnm1)
   to  gT=T+dt*gT
 o changed CALC_GT,CALC_GS & CALC_GTR1 to calcuate "gT" defined
   by new timestep_tracer (ie. including forcing, ABII, N-L F-S, etc...)
   now calls ADAMS_BASHFORTH2 and FREESURF_RESCALE_G
 o changed CYCLE_TRACER appropriately  T=gT only

Other associated mods:
 o new s/r: FREESURF_RESCALE_G applies non-linear free-surface term
   this used to be in TIMESTEP_TRACER
 o added myIter as argument to CALC_GS,CALC_GT,CALC_GTR1
1 adcroft 1.4 C $Header: /u/gcmpack/models/MITgcmUV/model/src/cycle_tracer.F,v 1.3 2001/02/04 14:38:46 cnh Exp $
2     C $Name: $
3 adcroft 1.2
4     #include "CPP_OPTIONS.h"
5    
6     C /==========================================================\
7     C | S/R CYCLE_TRACER |
8     C |==========================================================|
9     C | o Cycles the time-stepping arrays for a tracer field |
10     C \==========================================================/
11     SUBROUTINE CYCLE_TRACER(
12     I bi, bj, iMin, iMax, jMin, jMax, K,
13     U tracer, gTracer, gTrNm1,
14     I myTime, myThid )
15     IMPLICIT NONE
16     C Common
17     #include "SIZE.h"
18     #include "EEPARAMS.h"
19     #include "PARAMS.h"
20     C == Routine Arguments ==
21     C bi,bj,iMin,iMax,jMin,jMax, K - Loop counters
22     C tracer, gTracer, gTrNm1 - Tracer field, with tendencies
23     C myThid - Instance number for
24     C this call to S/R CORRECTION_STEP
25     C myTime - Current simulation time for this instance.
26     INTEGER bi,bj,iMin,iMax,jMin,jMax,K
27     _RL tracer (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
28     _RL gTracer(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
29     _RL gTrNm1 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
30     INTEGER myThid
31     _RL myTime
32    
33     C == Local variables ==
34     INTEGER i,j
35    
36     C Rotate tracer/gTracer/gTrNm1
37     DO j=jMin,jMax
38     DO i=iMin,iMax
39 adcroft 1.4 tracer(i,j,k,bi,bj)=gTracer(i,j,k,bi,bj)
40 adcroft 1.2 ENDDO
41     ENDDO
42    
43     RETURN
44     END
  ViewVC Help
Powered by ViewVC 1.1.22