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

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

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

revision 1.23 by molod, Tue May 24 14:57:35 2005 UTC revision 1.27 by molod, Tue Jun 14 18:14:21 2005 UTC
# Line 45  c--------------------------------------- Line 45  c---------------------------------------
45         _RL pressure(Nrphys)         _RL pressure(Nrphys)
46         _RL lats(sNx,sNy,Nsx,Nsy), lons(sNx,sNy,Nsx,Nsy)         _RL lats(sNx,sNy,Nsx,Nsy), lons(sNx,sNy,Nsx,Nsy)
47         _RL fracland(sNx,sNy,Nsx,Nsy)         _RL fracland(sNx,sNy,Nsx,Nsy)
48           _RL tempoverlap(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nsx,Nsy)
49    
50         idim2 = sNx+OLx         idim2 = sNx+OLx
51         jdim2 = sNy+OLy         jdim2 = sNy+OLy
# Line 113  C     to send to fizhi init chem Line 114  C     to send to fizhi init chem
114       .      Nrphys,pressure,n2o,methane,co2,cfc11,cfc12,cfc22)       .      Nrphys,pressure,n2o,methane,co2,cfc11,cfc12,cfc22)
115    
116  C Read dataset that contains topography variance (10m by 10m in a 2x2.5 box)  C Read dataset that contains topography variance (10m by 10m in a 2x2.5 box)
117          _BEGIN_MASTER( myThid )         CALL READ_REC_XY_RL('topvar19232.data',tempoverlap,1,0,mythid )
118          CALL READ_REC_XY_RL('topvar19232.data',phis_var,1,0,mythid )         _BARRIER
119          _END_MASTER(myThid)  
120          _BARRIER         do bj = myByLo(myThid), myByHi(myThid)
121           do bi = myBxLo(myThid), myBxHi(myThid)
122            do j=jm1,jm2
123            do i=im1,im2
124             if(fracland(i,j,bi,bj).gt.0.3) then
125              phis_var(i,j,bi,bj) = tempoverlap(i,j,bi,bj)
126             else
127              phis_var(i,j,bi,bj) = 0.
128             endif
129            enddo
130            enddo
131           enddo
132           enddo
133    
134    C Finally, grab unit numbers for reading sst and sea ice
135    C   (held in common block fizhi ocean coms)
136    
137           call mdsfindunit( kice, myThid )
138           open(kice)
139           call mdsfindunit( ksst, myThid )
140           open(ksst)
141    
142         return         return
143         end         end
Legend:
Removed from v.1.23  
changed lines
  Added in v.1.27
  ViewVC Help
Powered by ViewVC 1.1.22