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

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

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

Revision 1.4 - (show annotations) (download)
Wed Sep 19 02:43:27 2001 UTC (22 years, 9 months 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 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
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 tracer(i,j,k,bi,bj)=gTracer(i,j,k,bi,bj)
40 ENDDO
41 ENDDO
42
43 RETURN
44 END
  ViewVC Help
Powered by ViewVC 1.1.22