148 |
|
|
149 |
c _BEGIN_MASTER(myThid) |
c _BEGIN_MASTER(myThid) |
150 |
|
|
151 |
|
cMass dry atmosphere = (5.1352+/-0.0003)d18 kg (Trenberth & Smith, |
152 |
|
cJournal of Climate 2005) |
153 |
|
cand Mean molecular mass air = 28.97 g/mol (NASA earth fact sheet) |
154 |
|
total_atmos_moles= 1.77 _d 20 |
155 |
|
c for 278ppmv we need total_atmos_carbon=4.9206e+16 |
156 |
|
|
157 |
total_atmos_moles= 1.5 _d 20 |
total_atmos_moles= 1.5 _d 20 |
158 |
|
|
159 |
if (istate.gt.0) then |
if (istate.gt.0) then |
217 |
WRITE(fn,'(A,I10.10)') 'dic_atmos.',myIter |
WRITE(fn,'(A,I10.10)') 'dic_atmos.',myIter |
218 |
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. |
219 |
_BARRIER |
_BARRIER |
220 |
|
c write values to new pickup |
221 |
open(26,file=fn,status='new') |
open(26,file=fn,status='new') |
222 |
write(26,*) total_atmos_carbon, atpco2 |
write(26,*) total_atmos_carbon, atpco2 |
223 |
close(26) |
close(26) |