/[MITgcm]/MITgcm/pkg/matrix/matrix_init.F
ViewVC logotype

Contents of /MITgcm/pkg/matrix/matrix_init.F

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

Revision 1.3 - (show annotations) (download)
Mon Nov 5 18:58:00 2007 UTC (16 years, 6 months ago) by jmc
Branch: MAIN
CVS Tags: checkpoint63a, checkpoint60, checkpoint61, checkpoint62, checkpoint63, checkpoint61l, checkpoint59q, checkpoint59p, checkpoint59r, checkpoint59m, checkpoint59l, checkpoint59o, checkpoint59n, checkpoint59k, checkpoint59j, checkpoint62c, checkpoint62b, checkpoint62a, checkpoint62g, checkpoint62f, checkpoint62e, checkpoint62d, checkpoint62k, checkpoint62j, checkpoint62i, checkpoint62h, checkpoint62o, checkpoint62n, checkpoint62m, checkpoint62l, checkpoint62s, checkpoint62r, checkpoint62q, checkpoint62p, checkpoint62w, checkpoint62v, checkpoint62u, checkpoint62t, checkpoint62z, checkpoint62y, checkpoint62x, checkpoint61f, checkpoint61g, checkpoint61d, checkpoint61e, checkpoint61b, checkpoint61c, checkpoint61a, checkpoint61n, checkpoint61o, checkpoint61m, checkpoint61j, checkpoint61k, checkpoint61h, checkpoint61i, checkpoint61v, checkpoint61w, checkpoint61t, checkpoint61u, checkpoint61r, checkpoint61s, checkpoint61p, checkpoint61q, checkpoint61z, checkpoint61x, checkpoint61y
Changes since 1.2: +27 -26 lines 
split PTRACERS.h in 2 header files: PTRACERS_FIELDS.h & PTRACERS_PARAMS
 declare Namelist in matrix_readparms.F (was previously in MATRIX.h)
1 C $Header: /u/gcmpack/MITgcm/pkg/matrix/matrix_init.F,v 1.2 2005/04/20 18:05:08 spk Exp $
2 C $Name: $
3
4 #include "PACKAGES_CONFIG.h"
5 #include "CPP_OPTIONS.h"
6
7 SUBROUTINE MATRIX_INIT( myThid )
8
9 IMPLICIT NONE
10
11 #include "EEPARAMS.h"
12 #include "SIZE.h"
13 #include "PTRACERS_SIZE.h"
14 #include "PTRACERS_PARAMS.h"
15 #include "PTRACERS_FIELDS.h"
16 #include "MATRIX.h"
17 #include "PARAMS.h"
18
19 INTEGER myThid
20
21 #ifdef ALLOW_MATRIX
22
23 C local variables:
24 INTEGER i, j, k, bi, bj
25 INTEGER iTracer
26
27 C Initialize counters
28 expMatrixCounter=0
29 impMatrixCounter=0
30 DO iTracer=1,PTRACERS_num
31 DO bj=myByLo(myThid), myByHi (myThid)
32 DO bi=myBxLo (myThid), myBxHi (myThid)
33 DO k=1,Nr
34 DO j=1-Oly, sNy+Oly
35 DO i=1-Olx, sNx+Olx
36 PTRACERS_initial(i,j,k,bi,bj,iTracer)=
37 & pTracer(i,j,k,bi,bj,iTracer)
38 MATRIX(i,j,k,bi,bj,iTracer,1)= 0. _d 0
39 MATRIX(i,j,k,bi,bj,iTracer,2)= 0. _d 0
40 ENDDO
41 ENDDO
42 ENDDO
43 ENDDO
44 ENDDO
45 ENDDO
46
47 C DO iTracer=1,PTRACERS_numInUse
48 C DO bj=myByLo(myThid), myByHi (myThid)
49 C DO bi=myBxLo (myThid), myBxHi (myThid)
50 C DO k=1,Nr
51 C DO j=1-Oly, sNy+Oly
52 C DO i=1-Olx, sNx+Olx
53 C MATRIX(i,j,k,bi,bj,iTracer,1)= 0. _d 0
54 C MATRIX(i,j,k,bi,bj,iTracer,2)= 0. _d 0
55 C ENDDO
56 C ENDDO
57 C ENDDO
58 C ENDDO
59 C ENDDO
60 C ENDDO
61
62 #endif /* ALLOW_MATRIX */
63 RETURN
64 END
  ViewVC Help
Powered by ViewVC 1.1.22