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

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

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

revision 1.6 by molod, Tue Aug 23 18:25:33 2005 UTC revision 1.7 by molod, Wed Apr 19 21:45:46 2006 UTC
# Line 96  C Only do I/O if I am the master thread Line 96  C Only do I/O if I am the master thread
96          tempgrid(i,j) = globalarr(i+bislot,j+bjslot,1)          tempgrid(i,j) = globalarr(i+bislot,j+bjslot,1)
97         enddo         enddo
98         enddo         enddo
99         call grd2msc(tempgrid,sNx,sNy,igrd,temptile,nchp,nchptot(bi,bj))         call grd2msc(tempgrid,sNx,sNy,igrd(1,bi,bj),
100         .                 temptile,nchp,nchptot(bi,bj))
101         do n = 1,nchp         do n = 1,nchp
102          tcanopy(n,bi,bj) = temptile(n)          tcanopy(n,bi,bj) = temptile(n)
103         enddo         enddo
# Line 108  C Only do I/O if I am the master thread Line 109  C Only do I/O if I am the master thread
109       .    tempgrid(i,j) = globalarr(i+bislot,j+bjslot,1) - 0.5       .    tempgrid(i,j) = globalarr(i+bislot,j+bjslot,1) - 0.5
110         enddo         enddo
111         enddo         enddo
112         call grd2msc(tempgrid,sNx,sNy,igrd,temptile,nchp,nchptot(bi,bj))         call grd2msc(tempgrid,sNx,sNy,igrd(1,bi,bj),
113         .                    temptile,nchp,nchptot(bi,bj))
114         do n = 1,nchp         do n = 1,nchp
115          tdeep(n,bi,bj) = temptile(n)          tdeep(n,bi,bj) = temptile(n)
116         enddo         enddo
# Line 120  C Only do I/O if I am the master thread Line 122  C Only do I/O if I am the master thread
122       .    tempgrid(i,j) = 0.01       .    tempgrid(i,j) = 0.01
123         enddo         enddo
124         enddo         enddo
125         call grd2msc(tempgrid,sNx,sNy,igrd,temptile,nchp,nchptot(bi,bj))         call grd2msc(tempgrid,sNx,sNy,igrd(1,bi,bj),
126         .                    temptile,nchp,nchptot(bi,bj))
127         do n = 1,nchp         do n = 1,nchp
128          ecanopy(n,bi,bj) = temptile(n)          ecanopy(n,bi,bj) = temptile(n)
129         enddo         enddo
# Line 132  C Only do I/O if I am the master thread Line 135  C Only do I/O if I am the master thread
135       .    tempgrid(i,j) = 0.7       .    tempgrid(i,j) = 0.7
136         enddo         enddo
137         enddo         enddo
138         call grd2msc(tempgrid,sNx,sNy,igrd,temptile,nchp,nchptot(bi,bj))         call grd2msc(tempgrid,sNx,sNy,igrd(1,bi,bj),
139         .                    temptile,nchp,nchptot(bi,bj))
140         do n = 1,nchp         do n = 1,nchp
141          swetshal(n,bi,bj) = temptile(n)          swetshal(n,bi,bj) = temptile(n)
142         enddo         enddo
# Line 144  C Only do I/O if I am the master thread Line 148  C Only do I/O if I am the master thread
148       .    tempgrid(i,j) = 0.5       .    tempgrid(i,j) = 0.5
149         enddo         enddo
150         enddo         enddo
151         call grd2msc(tempgrid,sNx,sNy,igrd,temptile,nchp,nchptot(bi,bj))         call grd2msc(tempgrid,sNx,sNy,igrd(1,bi,bj),
152         .                    temptile,nchp,nchptot(bi,bj))
153         do n = 1,nchp         do n = 1,nchp
154          swetroot(n,bi,bj) = temptile(n)          swetroot(n,bi,bj) = temptile(n)
155         enddo         enddo
# Line 156  C Only do I/O if I am the master thread Line 161  C Only do I/O if I am the master thread
161       .    tempgrid(i,j) = 0.3       .    tempgrid(i,j) = 0.3
162         enddo         enddo
163         enddo         enddo
164         call grd2msc(tempgrid,sNx,sNy,igrd,temptile,nchp,nchptot(bi,bj))         call grd2msc(tempgrid,sNx,sNy,igrd(1,bi,bj),
165         .                    temptile,nchp,nchptot(bi,bj))
166         do n = 1,nchp         do n = 1,nchp
167          swetdeep(n,bi,bj) = temptile(n)          swetdeep(n,bi,bj) = temptile(n)
168         enddo         enddo
# Line 168  C Only do I/O if I am the master thread Line 174  C Only do I/O if I am the master thread
174       .    tempgrid(i,j) = 0.       .    tempgrid(i,j) = 0.
175         enddo         enddo
176         enddo         enddo
177         call grd2msc(tempgrid,sNx,sNy,igrd,temptile,nchp,nchptot(bi,bj))         call grd2msc(tempgrid,sNx,sNy,igrd(1,bi,bj),
178         .                    temptile,nchp,nchptot(bi,bj))
179         do n = 1,nchp         do n = 1,nchp
180          snodep(n,bi,bj) = temptile(n)          snodep(n,bi,bj) = temptile(n)
181         enddo         enddo
# Line 180  C Only do I/O if I am the master thread Line 187  C Only do I/O if I am the master thread
187       .    tempgrid(i,j) = 0.       .    tempgrid(i,j) = 0.
188         enddo         enddo
189         enddo         enddo
190         call grd2msc(tempgrid,sNx,sNy,igrd,temptile,nchp,nchptot(bi,bj))         call grd2msc(tempgrid,sNx,sNy,igrd(1,bi,bj),
191         .                    temptile,nchp,nchptot(bi,bj))
192         do n = 1,nchp         do n = 1,nchp
193          capac(n,bi,bj) = temptile(n)          capac(n,bi,bj) = temptile(n)
194         enddo         enddo
Legend:
Removed from v.1.6  
changed lines
  Added in v.1.7
  ViewVC Help
Powered by ViewVC 1.1.22