New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 10196 for branches/NERC/dev_r5518_GO6_under_ice_relax/NEMOGCM/NEMO/TOP_SRC/MEDUSA/air_sea.F90 – NEMO

Ignore:
Timestamp:
2018-10-16T12:15:14+02:00 (5 years ago)
Author:
jpalmier
Message:

add DMS flux --

File:
1 edited

Legend:

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

    r10049 r10196  
    77   !!   -   ! 2017-04 (M. Stringer)        Code taken from trcbio_medusa.F90 
    88   !!   -   ! 2017-08 (A. Yool)            Add air-sea flux kill switch 
    9    !!   -   ! 2018-08 (A. Yool)            add OMIP preindustrial DIC 
     9   !!   -   ! 2018-08 (A. Yool)            add OMIP preindustrial DIC  
     10   !!   -   ! 2018-10 (A. Yool)            Add air-sea DMS flux 
    1011   !!---------------------------------------------------------------------- 
    1112#if defined key_medusa 
     
    4344                                   dms_surf2d, dms_andm, dms_andm2d,      & 
    4445                                   dms_nlim, dms_wtkn,                    & 
     46               dms_flux, dms_flux2d,                  & 
    4547                                   f_co2flux, f_co2flux2d,                & 
    4648                                   f_co2starair_2d, f_co3,                & 
     
    105107      USE trcco2_medusa,     ONLY: trc_co2_medusa 
    106108#  endif  
    107       USE trcdms_medusa,     ONLY: trc_dms_medusa 
     109      USE trcdms_medusa,     ONLY: trc_dms_medusa, dms_flux_ocn 
    108110      USE trcoxy_medusa,     ONLY: trc_oxy_medusa 
    109111# endif 
     
    540542                  endif 
    541543                  !! 
     544                  !! calculate air-sea DMS flux using UM code 
     545                  CALL dms_flux_ocn(wndm(ji,jj),ztmp(ji,jj),dms_surf,1, & 
     546                                    dms_flux) 
     547                  !! 
     548                  !! Correct DMS flux for sea-ice, but leave units as kg S m-2 s-1 
     549                  dms_flux = (1. - fr_i(ji,jj)) * dms_flux 
     550                  !! 
    542551                  !! 2D diag through iom_use 
    543552                  IF( med_diag%DMS_SURF%dgsave ) THEN 
     
    558567                  IF( med_diag%DMS_ANDM%dgsave ) THEN  
    559568                     dms_andm2d(ji,jj) = dms_andm 
     569                  ENDIF 
     570                  IF( med_diag%DMS_FLUX%dgsave ) THEN  
     571                     dms_flux2d(ji,jj) = dms_flux 
    560572                  ENDIF 
    561573               ENDIF 
Note: See TracChangeset for help on using the changeset viewer.