117 |
#endif |
#endif |
118 |
|
|
119 |
#if ( defined CPL_OCEANCO2 ) |
#if ( defined CPL_OCEANCO2 ) |
120 |
#include "OCM.COM" |
#include "OCM.h" |
121 |
common /Garyflux/pC_atm(jm0),wind_amp,fluxco2(jm0) |
common /Garyflux/pC_atm(jm0),wind_amp,fluxco2(jm0) |
122 |
# if ( defined ML_2D ) |
# if ( defined ML_2D ) |
123 |
common/Garyclim/tggary(jm0),wsgary(jm0),areaml(jm0),arsrf(jm0) |
common/Garyclim/tggary(jm0),wsgary(jm0),areaml(jm0),arsrf(jm0) |
365 |
NSTEP2=NSTEP 29.6 |
NSTEP2=NSTEP 29.6 |
366 |
MRCHT=0. 29.7 |
MRCHT=0. 29.7 |
367 |
ITAU=(NSTEP+NSTEP0)*IDTHR 30. |
ITAU=(NSTEP+NSTEP0)*IDTHR 30. |
368 |
TAU=FLOAT(ITAU)/XINT 31. |
cjrs changed to dfloat 8/2/07 |
369 |
|
TAU=DFLOAT(ITAU)/XINT 31. |
370 |
IDAY=1+ITAU/I24 32. |
IDAY=1+ITAU/I24 32. |
371 |
TOFDAY=(ITAU-(IDAY-1)*I24)/XINT 33. |
TOFDAY=(ITAU-(IDAY-1)*I24)/XINT 33. |
372 |
! if(ISTART.eq.2.or.ISTRT1.eq.0.and..not.CONTRR)then |
! if(ISTART.eq.2.or.ISTRT1.eq.0.and..not.CONTRR)then |
1257 |
C**** 189. |
C**** 189. |
1258 |
500 NSTEP=NSTEP+NDYN 190. |
500 NSTEP=NSTEP+NDYN 190. |
1259 |
ITAU=(NSTEP+NSTEP0)*IDTHR 191. |
ITAU=(NSTEP+NSTEP0)*IDTHR 191. |
1260 |
TAU=FLOAT(ITAU)/XINT 192. |
cJRS fix to DFLOAT 8/2/07 |
1261 |
|
TAU=DFLOAT(ITAU)/XINT 192. |
1262 |
IDAY=1+ITAU/I24 193. |
IDAY=1+ITAU/I24 193. |
1263 |
TOFDAYPR=TOFDAY+1.00 |
TOFDAYPR=TOFDAY+1.00 |
1264 |
TOFDAY=(ITAU-(IDAY-1)*I24)/XINT 194. |
TOFDAY=(ITAU-(IDAY-1)*I24)/XINT 194. |
1880 |
C DTATM time step of atm model in hours |
C DTATM time step of atm model in hours |
1881 |
C precip and evap in mm/day or kg/m**2/day |
C precip and evap in mm/day or kg/m**2/day |
1882 |
do j=1,jm0 |
do j=1,jm0 |
1883 |
#if ( defined OCEAN_3D && defined CPL_OCEANCO2 ) |
Cjrs #if ( defined OCEAN_3D && defined CPL_OCEANCO2 ) |
1884 |
ncallatm=ncallatm+1 |
#ifdef OCEAN_3D |
1885 |
|
!jrs ncallatm=ncallatm+1 |
1886 |
! 020107 |
! 020107 |
1887 |
! co24ocean(j)=pC_atm(j)*1.e6 |
! co24ocean(j)=pC_atm(j)*1.e6 |
1888 |
|
! jrs give CO2 even if ocn carbon off |
1889 |
co24ocean(j)=atm_co2(j) |
co24ocean(j)=atm_co2(j) |
1890 |
! 020107 |
# ifdef CPL_OCEANCO2 |
1891 |
co24ocnan(j)=co24ocnan(j)+co24ocean(j) |
co24ocnan(j)=co24ocnan(j)+co24ocean(j) |
1892 |
|
# endif |
1893 |
#endif |
#endif |
1894 |
#ifdef ML_2D |
#ifdef ML_2D |
1895 |
cjrs block only MD_2D |
cjrs block only MD_2D |