/[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.2 by molod, Thu May 5 21:23:27 2005 UTC revision 1.3 by molod, Sat May 7 14:15:09 2005 UTC
# Line 58  C     !LOCAL VARIABLES: Line 58  C     !LOCAL VARIABLES:
58  C Only do I/O if I am the master thread  C Only do I/O if I am the master thread
59        _BEGIN_MASTER( myThid )        _BEGIN_MASTER( myThid )
60    
       print *,' Opening ',fn  
61        open(iUnit,file=fn,status='old',access='direct',recl=recl)        open(iUnit,file=fn,status='old',access='direct',recl=recl)
62        read(iunit,rec=1) globalarr        read(iunit,rec=1) globalarr
63        close( iunit )        close( iunit )
# Line 80  C Only do I/O if I am the master thread Line 79  C Only do I/O if I am the master thread
79         bjslot = myYGlobalLo-1+(bj-1)*sNy         bjslot = myYGlobalLo-1+(bj-1)*sNy
80  #endif /* ALLOW_EXCH2 */  #endif /* ALLOW_EXCH2 */
81    
82         do j = 1,sNx         do j = 1,sNy
83         do i = 1,sNx         do i = 1,sNx
84          tempgrid(i,j) = globalarr(i+bislot,j+bjslot,1)          tempgrid(i,j) = globalarr(i+bislot,j+bjslot,1)
85         enddo         enddo
# Line 90  C Only do I/O if I am the master thread Line 89  C Only do I/O if I am the master thread
89          tcanopy(n,bi,bj) = temptile(n)          tcanopy(n,bi,bj) = temptile(n)
90         enddo         enddo
91    
92         do j = 1,sNx         do j = 1,sNy
93         do i = 1,sNx         do i = 1,sNx
94          tempgrid(i,j) = globalarr(i+bislot,j+bjslot,2)          tempgrid(i,j) = globalarr(i+bislot,j+bjslot,2)
95         enddo         enddo
# Line 100  C Only do I/O if I am the master thread Line 99  C Only do I/O if I am the master thread
99          tdeep(n,bi,bj) = temptile(n)          tdeep(n,bi,bj) = temptile(n)
100         enddo         enddo
101    
102         do j = 1,sNx         do j = 1,sNy
103         do i = 1,sNx         do i = 1,sNx
104          tempgrid(i,j) = globalarr(i+bislot,j+bjslot,3)          tempgrid(i,j) = globalarr(i+bislot,j+bjslot,3)
105         enddo         enddo
# Line 110  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          ecanopy(n,bi,bj) = temptile(n)          ecanopy(n,bi,bj) = temptile(n)
110         enddo         enddo
111    
112         do j = 1,sNx         do j = 1,sNy
113         do i = 1,sNx         do i = 1,sNx
114          tempgrid(i,j) = globalarr(i+bislot,j+bjslot,4)          tempgrid(i,j) = globalarr(i+bislot,j+bjslot,4)
115         enddo         enddo
# Line 120  C Only do I/O if I am the master thread Line 119  C Only do I/O if I am the master thread
119          swetshal(n,bi,bj) = temptile(n)          swetshal(n,bi,bj) = temptile(n)
120         enddo         enddo
121    
122         do j = 1,sNx         do j = 1,sNy
123         do i = 1,sNx         do i = 1,sNx
124          tempgrid(i,j) = globalarr(i+bislot,j+bjslot,5)          tempgrid(i,j) = globalarr(i+bislot,j+bjslot,5)
125         enddo         enddo
# Line 130  C Only do I/O if I am the master thread Line 129  C Only do I/O if I am the master thread
129          swetroot(n,bi,bj) = temptile(n)          swetroot(n,bi,bj) = temptile(n)
130         enddo         enddo
131    
132         do j = 1,sNx         do j = 1,sNy
133         do i = 1,sNx         do i = 1,sNx
134          tempgrid(i,j) = globalarr(i+bislot,j+bjslot,6)          tempgrid(i,j) = globalarr(i+bislot,j+bjslot,6)
135         enddo         enddo
# Line 140  C Only do I/O if I am the master thread Line 139  C Only do I/O if I am the master thread
139          swetdeep(n,bi,bj) = temptile(n)          swetdeep(n,bi,bj) = temptile(n)
140         enddo         enddo
141    
142         do j = 1,sNx         do j = 1,sNy
143         do i = 1,sNx         do i = 1,sNx
144          tempgrid(i,j) = globalarr(i+bislot,j+bjslot,7)          tempgrid(i,j) = globalarr(i+bislot,j+bjslot,7)
145         enddo         enddo
# Line 150  C Only do I/O if I am the master thread Line 149  C Only do I/O if I am the master thread
149          snodep(n,bi,bj) = temptile(n)          snodep(n,bi,bj) = temptile(n)
150         enddo         enddo
151    
152         do j = 1,sNx         do j = 1,sNy
153         do i = 1,sNx         do i = 1,sNx
154          tempgrid(i,j) = globalarr(i+bislot,j+bjslot,8)          tempgrid(i,j) = globalarr(i+bislot,j+bjslot,8)
155         enddo         enddo
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22