/[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.3 by molod, Sat May 7 14:15:09 2005 UTC revision 1.4 by molod, Wed May 11 19:04:01 2005 UTC
# Line 47  C     !LOCAL VARIABLES: Line 47  C     !LOCAL VARIABLES:
47        Real*8 globalarr(Nx,Ny,8)        Real*8 globalarr(Nx,Ny,8)
48        _RL tempgrid(sNx,sNy)        _RL tempgrid(sNx,sNy)
49        _RL temptile(nchp)        _RL temptile(nchp)
50          _RL fracland(sNx,sNy,Nsx,Nsy)
51    
52        ihour = nhms/10000        ihour = nhms/10000
53        WRITE(fn,'(a,I8,a,I2.2,a)')'vegtiles_cs32.d',nymd,'z',ihour,'.bin'        WRITE(fn,'(a,I8,a,I2.2,a)')'vegtiles_cs32.d',nymd,'z',ihour,'.bin'
# Line 79  C Only do I/O if I am the master thread Line 80  C Only do I/O if I am the master thread
80         bjslot = myYGlobalLo-1+(bj-1)*sNy         bjslot = myYGlobalLo-1+(bj-1)*sNy
81  #endif /* ALLOW_EXCH2 */  #endif /* ALLOW_EXCH2 */
82    
83          call get_landfrac(sNx,sNy,Nsx,Nsy,bi,bj,maxtyp,
84         .        surftype,tilefrac,fracland)
85    
86         do j = 1,sNy         do j = 1,sNy
87         do i = 1,sNx         do i = 1,sNx
88          tempgrid(i,j) = globalarr(i+bislot,j+bjslot,1)          tempgrid(i,j) = globalarr(i+bislot,j+bjslot,1)
# Line 92  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         do j = 1,sNy         do j = 1,sNy
97         do i = 1,sNx         do i = 1,sNx
98          tempgrid(i,j) = globalarr(i+bislot,j+bjslot,2)          tempgrid(i,j) = globalarr(i+bislot,j+bjslot,2)
99            if (tempgrid(i,j).gt.1.e14 .and. fracland(i,j,bi,bj).gt.0.0001)
100         .    tempgrid(i,j) = globalarr(i+bislot,j+bjslot,1) - 0.5
101         enddo         enddo
102         enddo         enddo
103         call grd2msc(tempgrid,sNx,sNy,igrd,temptile,nchp,nchptot(bi,bj))         call grd2msc(tempgrid,sNx,sNy,igrd,temptile,nchp,nchptot(bi,bj))
# Line 102  C Only do I/O if I am the master thread Line 108  C Only do I/O if I am the master thread
108         do j = 1,sNy         do j = 1,sNy
109         do i = 1,sNx         do i = 1,sNx
110          tempgrid(i,j) = globalarr(i+bislot,j+bjslot,3)          tempgrid(i,j) = globalarr(i+bislot,j+bjslot,3)
111            if (tempgrid(i,j).gt.1.e14 .and. fracland(i,j,bi,bj).gt.0.0001)
112         .    tempgrid(i,j) = 0.01
113         enddo         enddo
114         enddo         enddo
115         call grd2msc(tempgrid,sNx,sNy,igrd,temptile,nchp,nchptot(bi,bj))         call grd2msc(tempgrid,sNx,sNy,igrd,temptile,nchp,nchptot(bi,bj))
# Line 112  C Only do I/O if I am the master thread Line 120  C Only do I/O if I am the master thread
120         do j = 1,sNy         do j = 1,sNy
121         do i = 1,sNx         do i = 1,sNx
122          tempgrid(i,j) = globalarr(i+bislot,j+bjslot,4)          tempgrid(i,j) = globalarr(i+bislot,j+bjslot,4)
123            if (tempgrid(i,j).gt.1.e14 .and. fracland(i,j,bi,bj).gt.0.0001)
124         .    tempgrid(i,j) = 0.7
125         enddo         enddo
126         enddo         enddo
127         call grd2msc(tempgrid,sNx,sNy,igrd,temptile,nchp,nchptot(bi,bj))         call grd2msc(tempgrid,sNx,sNy,igrd,temptile,nchp,nchptot(bi,bj))
# Line 122  C Only do I/O if I am the master thread Line 132  C Only do I/O if I am the master thread
132         do j = 1,sNy         do j = 1,sNy
133         do i = 1,sNx         do i = 1,sNx
134          tempgrid(i,j) = globalarr(i+bislot,j+bjslot,5)          tempgrid(i,j) = globalarr(i+bislot,j+bjslot,5)
135            if (tempgrid(i,j).gt.1.e14 .and. fracland(i,j,bi,bj).gt.0.0001)
136         .    tempgrid(i,j) = 0.5
137         enddo         enddo
138         enddo         enddo
139         call grd2msc(tempgrid,sNx,sNy,igrd,temptile,nchp,nchptot(bi,bj))         call grd2msc(tempgrid,sNx,sNy,igrd,temptile,nchp,nchptot(bi,bj))
# Line 132  C Only do I/O if I am the master thread Line 144  C Only do I/O if I am the master thread
144         do j = 1,sNy         do j = 1,sNy
145         do i = 1,sNx         do i = 1,sNx
146          tempgrid(i,j) = globalarr(i+bislot,j+bjslot,6)          tempgrid(i,j) = globalarr(i+bislot,j+bjslot,6)
147            if (tempgrid(i,j).gt.1.e14 .and. fracland(i,j,bi,bj).gt.0.0001)
148         .    tempgrid(i,j) = 0.3
149         enddo         enddo
150         enddo         enddo
151         call grd2msc(tempgrid,sNx,sNy,igrd,temptile,nchp,nchptot(bi,bj))         call grd2msc(tempgrid,sNx,sNy,igrd,temptile,nchp,nchptot(bi,bj))
# Line 142  C Only do I/O if I am the master thread Line 156  C Only do I/O if I am the master thread
156         do j = 1,sNy         do j = 1,sNy
157         do i = 1,sNx         do i = 1,sNx
158          tempgrid(i,j) = globalarr(i+bislot,j+bjslot,7)          tempgrid(i,j) = globalarr(i+bislot,j+bjslot,7)
159            if (tempgrid(i,j).gt.1.e14 .and. fracland(i,j,bi,bj).gt.0.0001)
160         .    tempgrid(i,j) = 0.
161         enddo         enddo
162         enddo         enddo
163         call grd2msc(tempgrid,sNx,sNy,igrd,temptile,nchp,nchptot(bi,bj))         call grd2msc(tempgrid,sNx,sNy,igrd,temptile,nchp,nchptot(bi,bj))
# Line 152  C Only do I/O if I am the master thread Line 168  C Only do I/O if I am the master thread
168         do j = 1,sNy         do j = 1,sNy
169         do i = 1,sNx         do i = 1,sNx
170          tempgrid(i,j) = globalarr(i+bislot,j+bjslot,8)          tempgrid(i,j) = globalarr(i+bislot,j+bjslot,8)
171            if (tempgrid(i,j).gt.1.e14 .and. fracland(i,j,bi,bj).gt.0.0001)
172         .    tempgrid(i,j) = 0.
173         enddo         enddo
174         enddo         enddo
175         call grd2msc(tempgrid,sNx,sNy,igrd,temptile,nchp,nchptot(bi,bj))         call grd2msc(tempgrid,sNx,sNy,igrd,temptile,nchp,nchptot(bi,bj))
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22