Ignore:
Timestamp:
2017-09-04T18:50:10+02:00 (3 years ago)
Author:
dford
Message:

Merge in changes from dev_r5518_GO6_package_asm_surf_bgc, and adapt to the updated MEDUSA structure.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/UKMO/dev_r5518_GO6_package_asm_surf_bgc_v2/NEMOGCM/NEMO/TOP_SRC/MEDUSA/trcbio_medusa.F90

    r8442 r8495  
    9898      !! JPALM (27-06-2016): add lk_oasis for CO2 and DMS coupling with atm 
    9999      USE sbc_oce,                    ONLY: lk_oasis 
     100# if defined key_foam_medusa 
     101      USE sms_medusa,                 ONLY: hist_pco2, xobs_xco2a,          & 
     102                                            pgrow_avg, ploss_avg,           & 
     103                                            phyt_avg, mld_max 
     104# else 
    100105      USE sms_medusa,                 ONLY: hist_pco2 
     106# endif 
    101107      USE trc,                        ONLY: ln_rsttr, nittrc000, trn 
    102108      USE bio_medusa_init_mod,        ONLY: bio_medusa_init 
     
    309315         f_xco2a(:,:) = fq4 
    310316      endif 
     317#  if defined key_foam_medusa 
     318      IF ( xobs_xco2a > 0.0 ) THEN 
     319         IF(lwp) WRITE(numout,*) ' using observed atm pCO2 = ', xobs_xco2a 
     320         f_xco2a(:,:) = xobs_xco2a 
     321      ELSE 
     322         IF(lwp) WRITE(numout,*) ' xobs_xco2a <= 0 so using default atm pCO2' 
     323      ENDIF 
     324#  endif 
    311325#  if defined key_axy_pi_co2 
    312326      !! OCMIP pre-industrial pCO2 
     
    348362      !!          x * 30d + 1*rdt(i.e: mod = rdt)    
    349363      !!          ++ need to pass carb-chem output var through restarts 
     364#if defined key_foam_medusa 
     365      !! DAF (Aug 2017): For FOAM we want to run daily 
     366      If ( (kt == nittrc000 .AND. .NOT.ln_rsttr) .OR.                        & 
     367           (mod(kt*rdt,86400.) == rdt) ) THEN 
     368#else 
    350369      If ( (kt == nittrc000 .AND. .NOT.ln_rsttr) .OR.                        & 
    351370           ((86400*mod(nn_date0,100) + mod(kt*rdt,2592000.)) == rdt) ) THEN 
    352          !!--------------------------------------------------------------- 
     371#endif 
     372         !!---------------------------------------------------------------------- 
    353373         !! Calculate the carbonate chemistry for the whole ocean on the first 
    354374         !! simulation timestep and every month subsequently; the resulting 3D 
Note: See TracChangeset for help on using the changeset viewer.