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

    r10020 r10149  
    7474                                   jriver_alk, jriver_c,                  & 
    7575                                   jriver_n, jriver_si,                   & 
     76                                   ln_foam_medusa,                        & 
    7677                                   riv_alk, riv_c, riv_n, riv_si,         & 
    7778                                   zn_dms_chd, zn_dms_chn, zn_dms_din,    & 
    7879                                   zn_dms_mld, zn_dms_qsr,                & 
     80                                   f2_pco2w, f2_fco2w,                    & 
    7981                                   xnln, xnld  
    8082      USE trc,               ONLY: med_diag 
     
    8789#  else 
    8890      USE trcco2_medusa,     ONLY: trc_co2_medusa 
     91      USE mocsy_mainmod,     ONLY: p2fCO2 
    8992#  endif 
    9093      USE trcdms_medusa,     ONLY: trc_dms_medusa 
     
    327330                     iters, ' AT (', ji, ', ', jj, ', 1) AT ', kt 
    328331               endif 
     332               IF ( ln_foam_medusa ) THEN 
     333                  !! DAF (Aug 2017): calculate fCO2 for observation operator 
     334                  CALL p2fCO2( f_pco2w, ztmp, f_pp0, 0.0, 1, f_fco2w ) 
     335               ENDIF 
    329336            ENDIF 
    330337         ENDDO 
     
    506513                              CO2flux_conv 
    507514               !! ENDIF 
     515               IF ( ln_foam_medusa ) THEN 
     516                  !! DAF (Aug 2017): Save pCO2 and fCO2 for observation operator 
     517                  f2_pco2w(ji,jj) = f_pco2w(ji,jj) 
     518                  f2_fco2w(ji,jj) = f_pco2w(ji,jj) 
     519               ENDIF 
    508520               IF ( lk_iomput ) THEN 
    509521                  IF( med_diag%ATM_PCO2%dgsave ) THEN 
Note: See TracChangeset for help on using the changeset viewer.