Changeset 10049


Ignore:
Timestamp:
2018-08-10T17:45:44+02:00 (2 years ago)
Author:
jpalmier
Message:

correct PI CO2 flux bug

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/NERC/dev_r5518_GO6_under_ice_relax/NEMOGCM/NEMO/TOP_SRC/MEDUSA/air_sea.F90

    r10047 r10049  
    455455               !! mmol/m2/s -> mmol/m3/d; correct for sea-ice; divide  
    456456               !! through by layer thickness 
    457                f_co2flux(ji,jj) = (1. - fr_i(ji,jj)) * f_co2flux(ji,jj) *    & 
     457               f_co2flux(ji,jj) = (1. - fr_i(ji,jj)) * f_co2flux(ji,jj) *       & 
    458458                                  86400. / fse3t(ji,jj,1) 
    459                !! 
     459#   if defined key_omip_dic 
     460               !! PI DIC - for omip runs 
     461               f_pi_co2flux(ji,jj) = (1. - fr_i(ji,jj)) * f_pi_co2flux(ji,jj) * & 
     462                                      86400. / fse3t(ji,jj,1) 
     463#   endif 
    460464               !! oxygen (O2); OCMIP-2 code 
    461465               !! AXY (23/06/15): amend input list for oxygen to account  
    462466               !!                 for common gas transfer velocity 
    463                CALL trc_oxy_medusa(ztmp(ji,jj),zsal(ji,jj),f_kw660(ji,jj),   & 
    464                                    f_pp0(ji,jj),zoxy(ji,jj),                 & 
    465                                    f_kwo2_dum,f_o2flux(ji,jj),               & 
     467               CALL trc_oxy_medusa(ztmp(ji,jj),zsal(ji,jj),f_kw660(ji,jj),     & 
     468                                   f_pp0(ji,jj),zoxy(ji,jj),                   & 
     469                                   f_kwo2_dum,f_o2flux(ji,jj),                 & 
    466470                                   f_o2sat(ji,jj)) 
    467471               !! 
    468472               !! mmol/m2/s -> mol/m3/d; correct for sea-ice; divide  
    469473               !! through by layer thickness 
    470                f_o2flux(ji,jj)  = (1. - fr_i(ji,jj)) * f_o2flux(ji,jj) *     & 
     474               f_o2flux(ji,jj)  = (1. - fr_i(ji,jj)) * f_o2flux(ji,jj) *       & 
    471475                                  86400. / fse3t(ji,jj,1) 
    472476            ENDIF 
Note: See TracChangeset for help on using the changeset viewer.