Changeset 10149 for branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/TOP_SRC/MEDUSA/trcbio_medusa.F90
- Timestamp:
- 2018-09-24T11:47:08+02:00 (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/TOP_SRC/MEDUSA/trcbio_medusa.F90
r10020 r10149 105 105 USE sbc_oce, ONLY: lk_oasis 106 106 USE sms_medusa, ONLY: hist_pco2, co2_yinit, co2_yend, & 107 lk_pi_co2 107 # if defined key_roam 108 xobs_xco2a, & 109 # endif 110 pgrow_avg, & 111 ploss_avg, phyt_avg, mld_max, & 112 lk_pi_co2, ln_foam_medusa 108 113 USE trc, ONLY: ln_rsttr, nittrc000, trn 109 114 USE bio_medusa_init_mod, ONLY: bio_medusa_init … … 330 335 !! f_xco2a(:,:) = 284.725 !! CMIP5 pre-industrial pCO2 331 336 f_xco2a(:,:) = 284.317 !! CMIP6 pre-industrial pCO2 337 ELSEIF ( xobs_xco2a > 0.0 ) THEN 338 IF(lwp) WRITE(numout,*) ' using observed atm pCO2 = ', xobs_xco2a 339 f_xco2a(:,:) = xobs_xco2a 332 340 ELSE 333 341 !! xCO2 from file … … 418 426 !! now use the NEMO calendar tool : nsec_month to be sure to call 419 427 !! at the beginning of a new month . 428 !! DAF: For FOAM we want to run daily 420 429 IF ( (kt == nittrc000 .AND. .NOT.ln_rsttr) .OR. & 421 ( nsec_month .LE. INT(rdt) ) ) THEN 430 ( nsec_month .LE. INT(rdt) ) .OR. & 431 ( nsec_day .LE. INT(rdt) .AND. ln_foam_medusa ) ) THEN 422 432 IF ( lwp ) WRITE(numout,*) & 423 433 ' *** 3D carb chem call *** -- kt:', kt, &
Note: See TracChangeset
for help on using the changeset viewer.