/[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.25 by molod, Tue May 31 21:33:41 2005 UTC revision 1.28 by molod, Thu Jun 30 23:09:08 2005 UTC
# Line 62  c--------------------------------------- Line 62  c---------------------------------------
62         endif         endif
63  #endif  #endif
64    
65    #ifdef ALLOW_DIAGNOSTICS
66          if ( useDiagnostics ) then
67            call fizhi_diagnostics_init( myThid )
68          endif
69    #endif
70    
71         call fizhi_alarms(nymdb,nhmsb,deltaTClock)         call fizhi_alarms(nymdb,nhmsb,deltaTClock)
72    
73         do bj = myByLo(myThid), myByHi(myThid)         do bj = myByLo(myThid), myByHi(myThid)
# Line 114  C     to send to fizhi init chem Line 120  C     to send to fizhi init chem
120       .      Nrphys,pressure,n2o,methane,co2,cfc11,cfc12,cfc22)       .      Nrphys,pressure,n2o,methane,co2,cfc11,cfc12,cfc22)
121    
122  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)
123          _BEGIN_MASTER( myThid )         CALL READ_REC_XY_RL('topvar19232.data',tempoverlap,1,0,mythid )
124          CALL READ_REC_XY_RL('topvar19232.data',tempoverlap,1,0,mythid )         _BARRIER
         _END_MASTER(myThid)  
         _BARRIER  
125    
126         do bj = myByLo(myThid), myByHi(myThid)         do bj = myByLo(myThid), myByHi(myThid)
127         do bi = myBxLo(myThid), myBxHi(myThid)         do bi = myBxLo(myThid), myBxHi(myThid)
# Line 133  C Read dataset that contains topography Line 137  C Read dataset that contains topography
137         enddo         enddo
138         enddo         enddo
139    
140    C Finally, grab unit numbers for reading sst and sea ice
141    C   (held in common block fizhi ocean coms)
142    
143           call mdsfindunit( kice, myThid )
144           open(kice)
145           call mdsfindunit( ksst, myThid )
146           open(ksst)
147    
148         return         return
149         end         end
Legend:
Removed from v.1.25  
changed lines
  Added in v.1.28
  ViewVC Help
Powered by ViewVC 1.1.22