| 307 |
Pa2Atm = 1.01325 _d 5 |
Pa2Atm = 1.01325 _d 5 |
| 308 |
#endif |
#endif |
| 309 |
|
|
| 310 |
|
#ifdef ALLOW_CDOM |
| 311 |
|
fraccdom=2. _d 0 / 100. _d 0 |
| 312 |
|
cdomdegrd= 1. _d 0 / (200 _d 0 *pday) |
| 313 |
|
cdombleach = 1. _d 0 / (15 _d 0 *pday) |
| 314 |
|
PARcdom = 20. _d 0 |
| 315 |
|
rnp_cdom = 16. _d 0 |
| 316 |
|
rfep_cdom = 1. _d -3 |
| 317 |
|
rcp_cdom = 120. _d 0 |
| 318 |
|
cdomcoeff = .1 _d -1 / 1.d -4 |
| 319 |
|
#endif |
| 320 |
|
|
| 321 |
C make sure we have reserved enough space in Ptracers |
C make sure we have reserved enough space in Ptracers |
| 322 |
IF ( nCompZooMax .LT. 4 ) THEN |
IF ( nCompZooMax .LT. 4 ) THEN |
| 323 |
WRITE(msgBuf,'(A,A,I3)') |
WRITE(msgBuf,'(A,A,I3)') |
| 328 |
ENDIF |
ENDIF |
| 329 |
DO nz = 1,nzmax |
DO nz = 1,nzmax |
| 330 |
iZooP (nz) = iZoo + (nz-1)*strideTypeZoo |
iZooP (nz) = iZoo + (nz-1)*strideTypeZoo |
| 331 |
|
#ifdef ONLY_P_CYCLE |
| 332 |
|
iZooN (nz) = nptot |
| 333 |
|
iZooFe(nz) = nptot |
| 334 |
|
iZooSi(nz) = nptot |
| 335 |
|
#else |
| 336 |
iZooN (nz) = iZoo + 1*strideCompZoo + (nz-1)*strideTypeZoo |
iZooN (nz) = iZoo + 1*strideCompZoo + (nz-1)*strideTypeZoo |
| 337 |
iZooFe(nz) = iZoo + 2*strideCompZoo + (nz-1)*strideTypeZoo |
iZooFe(nz) = iZoo + 2*strideCompZoo + (nz-1)*strideTypeZoo |
| 338 |
iZooSi(nz) = iZoo + 3*strideCompZoo + (nz-1)*strideTypeZoo |
iZooSi(nz) = iZoo + 3*strideCompZoo + (nz-1)*strideTypeZoo |
| 339 |
|
#endif |
| 340 |
ENDDO |
ENDDO |
| 341 |
#ifdef ALLOW_CARBON |
#ifdef ALLOW_CARBON |
| 342 |
DO nz = 1,nzmax |
DO nz = 1,nzmax |