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/trcnam_medusa.F90

    r9258 r10149  
    8787      &  xsdiss,                                              & 
    8888      &  sedlam,sedlostpoc,jpkb,jdms,jdms_input,jdms_model,   & 
    89       &  scl_chl, chl_out, dmsmin, dmscut, dmsslp 
     89      &  scl_chl, chl_out, dmsmin, dmscut, dmsslp,            & 
     90      &  ln_foam_medusa 
    9091#if defined key_roam 
    9192      NAMELIST/natroam/ xthetaphy,xthetazoo,xthetanit,        & 
     93      &    xobs_xco2a,                                        & 
    9294      &    xthetarem,xo2min  
    9395#endif 
     
    255257      dmscut      = 1.72 !! Anderson DMS default 
    256258      dmsslp      = 8.24 !! Anderson DMS default 
     259!! 
     260      ln_foam_medusa = .FALSE. 
    257261             
    258262      !REWIND(numnatm) 
     
    415419!!       chl_out     : select the chl field to send at the UM: 
    416420!!                     1- Surf Chl ; 2- MLD Chl  
     421!! 
     422!! FOAM - observation operator and data assimilation 
     423!!       ln_foam_medusa : calculate required diagnostics 
    417424 
    418425      IF(lwp) THEN 
     
    975982            ENDIF 
    976983         ENDIF ! IF lk_oasis=true 
     984!! FOAM 
     985         WRITE(numout,*) '=== FOAM-related parameters' 
     986         WRITE(numout,*)     & 
     987         &   ' calculate diagnostics for data assimilation,            ln_foam_medusa = ', ln_foam_medusa 
    977988!! 
    978989      ENDIF 
     
    10301041      xthetarem = 0. 
    10311042      xo2min    = 0. 
     1043      xobs_xco2a = 0. 
    10321044 
    10331045      !READ(numnatm,natroam) 
     
    10491061!!       xthetarem :  oxygen consumption by carbon remineralisation 
    10501062!!       xo2min    :  oxygen minimum concentration 
     1063!!       xobs_xco2a : observed atmospheric xCO2 (not used if <= 0) 
    10511064 
    10521065      IF(lwp) THEN 
     
    10661079          WRITE(numout,*)     & 
    10671080          &   ' oxygen minimum concentration                               xo2min      = ', xo2min 
     1081          WRITE(numout,*)     & 
     1082          &   ' observed atmospheric xCO2 (not used if <= 0)               xobs_xco2a  = ', xobs_xco2a 
    10681083       ENDIF 
    10691084 
Note: See TracChangeset for help on using the changeset viewer.