/[MITgcm]/MITgcm/pkg/dic/dic_atmos.F
ViewVC logotype

Diff of /MITgcm/pkg/dic/dic_atmos.F

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

revision 1.1 by stephd, Mon Nov 5 21:01:37 2007 UTC revision 1.2 by stephd, Mon Nov 5 22:28:50 2007 UTC
# Line 146  c         print*,'weights',ntim, aWght, Line 146  c         print*,'weights',ntim, aWght,
146  c interactive atmosphere  c interactive atmosphere
147        if (gchem_int1.eq.3) then        if (gchem_int1.eq.3) then
148    
149         _BEGIN_MASTER(myThid)  c      _BEGIN_MASTER(myThid)
150    
151         total_atmos_moles= 1.5 _d 20         total_atmos_moles= 1.5 _d 20
152    
# Line 179  c interactive atmosphere Line 179  c interactive atmosphere
179           ENDDO           ENDDO
180          ENDDO          ENDDO
181         ENDDO         ENDDO
182    
183         _GLOBAL_SUM_R8(total_flux,myThid)         _GLOBAL_SUM_R8(total_flux,myThid)
184         _GLOBAL_SUM_R8(total_ocean_carbon,myThid)         _GLOBAL_SUM_R8(total_ocean_carbon,myThid)
185    
   
186         if (istate.eq.0) then         if (istate.eq.0) then
187  c read state from output file  c read state from output file
188            DO i = 1,MAX_LEN_FNAM            DO i = 1,MAX_LEN_FNAM
# Line 211  c write out if time for a new pickup Line 211  c write out if time for a new pickup
211              WRITE(fn,'(A,I10.10)') 'dic_atmos.',myIter              WRITE(fn,'(A,I10.10)') 'dic_atmos.',myIter
212  C     Going to really do some IO. Make everyone except master thread wait.  C     Going to really do some IO. Make everyone except master thread wait.
213              _BARRIER              _BARRIER
   
214              open(26,file=fn,status='new')              open(26,file=fn,status='new')
215              write(26,*) total_atmos_carbon, atpco2              write(26,*) total_atmos_carbon, atpco2
216              close(26)              close(26)
# Line 220  C     Going to really do some IO. Make e Line 219  C     Going to really do some IO. Make e
219         endif         endif
220    
221    
 C--   Everyone else must wait  
       _BARRIER  
   
222        atpco2=total_atmos_carbon/total_atmos_moles        atpco2=total_atmos_carbon/total_atmos_moles
223    
224  c     print*,'QQpCO2', total_atmos_carbon, atpco2, total_ocean_carbon,  c     print*,'QQpCO2', total_atmos_carbon, atpco2, total_ocean_carbon,
# Line 282  c Line 278  c
278          endif          endif
279          endif          endif
280    
281        _END_MASTER(myThid)  c     _END_MASTER(myThid)
   
 C--   Everyone else must wait  
       _BARRIER  
   
282    
283         endif         endif
284    
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22