/[MITgcm]/MITgcm/pkg/fizhi/fizhi_wrapper.F
ViewVC logotype

Diff of /MITgcm/pkg/fizhi/fizhi_wrapper.F

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

revision 1.28 by molod, Thu Jun 30 23:09:08 2005 UTC revision 1.29 by molod, Thu Feb 16 23:03:53 2006 UTC
# Line 30  c--------------------------------------- Line 30  c---------------------------------------
30  #include "fizhi_earth_coms.h"  #include "fizhi_earth_coms.h"
31  #include "fizhi_ocean_coms.h"  #include "fizhi_ocean_coms.h"
32  #include "fizhi_chemistry_coms.h"  #include "fizhi_chemistry_coms.h"
33    #ifdef ALLOW_DIAGNOSTICS
34    #include "fizhi_SHP.h"
35    #endif
36    
37         integer myIter, myThid         integer myIter, myThid
38         _RL myTime         _RL myTime
# Line 56  c pe on dynamics and physics grid refers Line 59  c pe on dynamics and physics grid refers
59         jm1 = 1         jm1 = 1
60         jm2 = sNy         jm2 = sNy
61    
62    #ifdef ALLOW_DIAGNOSTICS
63          if ( useDiagnostics ) then
64           if(diagnostics_is_on('TENDUFIZ',myThid) .or.
65         .       diagnostics_is_on('CORRDU  ',myThid) ) then
66            do bj = myByLo(myThid), myByHi(myThid)
67            do bi = myBxLo(myThid), myBxHi(myThid)
68            do L = 1,Nrphys
69            do j = 1,sNy
70            do i = 1,sNx
71             ubef(i,j,L,bi,bj) = uphy(i,j,L,bi,bj)
72            enddo
73            enddo
74            enddo
75            do L = 1,Nr
76            do j = 1,sNy
77            do i = 1,sNx+1
78             udynbef(i,j,L,bi,bj) = uvel(i,j,L,bi,bj)
79            enddo
80            enddo
81            enddo
82            enddo
83            enddo
84           endif
85           if(diagnostics_is_on('TENDVFIZ',myThid) .or.
86         .       diagnostics_is_on('CORRDV  ',myThid) ) then
87            do bj = myByLo(myThid), myByHi(myThid)
88            do bi = myBxLo(myThid), myBxHi(myThid)
89            do L = 1,Nrphys
90            do j = 1,sNy
91            do i = 1,sNx
92             vbef(i,j,L,bi,bj) = vphy(i,j,L,bi,bj)
93            enddo
94            enddo
95            enddo
96            do L = 1,Nr
97            do j = 1,sNy+1
98            do i = 1,sNx
99             vdynbef(i,j,L,bi,bj) = vvel(i,j,L,bi,bj)
100            enddo
101            enddo
102            enddo
103            enddo
104            enddo
105           endif
106           if(diagnostics_is_on('TENDTFIZ',myThid) .or.
107         .       diagnostics_is_on('CORRDT  ',myThid) ) then
108            do bj = myByLo(myThid), myByHi(myThid)
109            do bi = myBxLo(myThid), myBxHi(myThid)
110            do L = 1,Nrphys
111            do j = 1,sNy
112            do i = 1,sNx
113             thbef(i,j,L,bi,bj) = thphy(i,j,L,bi,bj)
114            enddo
115            enddo
116            enddo
117            do L = 1,Nr
118            do j = 1,sNy
119            do i = 1,sNx
120             thdynbef(i,j,L,bi,bj) = theta(i,j,L,bi,bj)
121            enddo
122            enddo
123            enddo
124            enddo
125            enddo
126           endif
127           if(diagnostics_is_on('TENDQFIZ',myThid) .or.
128         .       diagnostics_is_on('CORRDQ  ',myThid) ) then
129            do bj = myByLo(myThid), myByHi(myThid)
130            do bi = myBxLo(myThid), myBxHi(myThid)
131            do L = 1,Nrphys
132            do j = 1,sNy
133            do i = 1,sNx
134             sbef(i,j,L,bi,bj) = sphy(i,j,L,bi,bj)
135            enddo
136            enddo
137            enddo
138            do L = 1,Nr
139            do j = 1,sNy
140            do i = 1,sNx
141             sdynbef(i,j,L,bi,bj) = salt(i,j,L,bi,bj)
142            enddo
143            enddo
144            enddo
145            enddo
146            enddo
147           endif
148          endif
149    #endif
150    
151         do bj = myByLo(myThid), myByHi(myThid)         do bj = myByLo(myThid), myByHi(myThid)
152         do bi = myBxLo(myThid), myBxHi(myThid)         do bi = myBxLo(myThid), myBxHi(myThid)
153    
Legend:
Removed from v.1.28  
changed lines
  Added in v.1.29
  ViewVC Help
Powered by ViewVC 1.1.22