Changeset 10196 for branches/NERC/dev_r5518_GO6_under_ice_relax/NEMOGCM/NEMO/TOP_SRC/MEDUSA/air_sea.F90
- Timestamp:
- 2018-10-16T12:15:14+02:00 (5 years ago)
- 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 7 7 !! - ! 2017-04 (M. Stringer) Code taken from trcbio_medusa.F90 8 8 !! - ! 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 10 11 !!---------------------------------------------------------------------- 11 12 #if defined key_medusa … … 43 44 dms_surf2d, dms_andm, dms_andm2d, & 44 45 dms_nlim, dms_wtkn, & 46 dms_flux, dms_flux2d, & 45 47 f_co2flux, f_co2flux2d, & 46 48 f_co2starair_2d, f_co3, & … … 105 107 USE trcco2_medusa, ONLY: trc_co2_medusa 106 108 # endif 107 USE trcdms_medusa, ONLY: trc_dms_medusa 109 USE trcdms_medusa, ONLY: trc_dms_medusa, dms_flux_ocn 108 110 USE trcoxy_medusa, ONLY: trc_oxy_medusa 109 111 # endif … … 540 542 endif 541 543 !! 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 !! 542 551 !! 2D diag through iom_use 543 552 IF( med_diag%DMS_SURF%dgsave ) THEN … … 558 567 IF( med_diag%DMS_ANDM%dgsave ) THEN 559 568 dms_andm2d(ji,jj) = dms_andm 569 ENDIF 570 IF( med_diag%DMS_FLUX%dgsave ) THEN 571 dms_flux2d(ji,jj) = dms_flux 560 572 ENDIF 561 573 ENDIF
Note: See TracChangeset
for help on using the changeset viewer.