Ignore:
Timestamp:
2018-09-24T11:47:08+02:00 (23 months ago)
Author:
frrh
Message:

Met Office GMED ticket 379: Merged David Ford's MEDUSA assimilation changes
using command:

svn merge -r 10054:10141 svn+ssh://forge.ipsl.jussieu.fr/ipsl/forge/projets/nemo/svn/branches/UKMO/dev_r5518_GO6_package_asm_3d_bgc_v3

File:
1 edited

Legend:

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

    r10020 r10149  
    105105      USE sbc_oce,                    ONLY: lk_oasis 
    106106      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 
    108113      USE trc,                        ONLY: ln_rsttr, nittrc000, trn 
    109114      USE bio_medusa_init_mod,        ONLY: bio_medusa_init 
     
    330335         !! f_xco2a(:,:) = 284.725       !! CMIP5 pre-industrial pCO2 
    331336         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 
    332340      ELSE       
    333341         !! xCO2 from file 
     
    418426      !! now use the NEMO calendar tool : nsec_month to be sure to call  
    419427      !! at the beginning of a new month . 
     428      !! DAF: For FOAM we want to run daily 
    420429      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 
    422432           IF ( lwp )  WRITE(numout,*)                                       &          
    423433                              ' *** 3D carb chem call *** -- kt:', kt,       & 
Note: See TracChangeset for help on using the changeset viewer.