Changeset 6683 for branches/UKMO/dev_r5518_RH_MEDUSA_Stable
- Timestamp:
- 2016-06-10T16:48:58+02:00 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/dev_r5518_RH_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC/MEDUSA/trcbio_medusa.F90
r6639 r6683 58 58 !! trc_bio_medusa : 59 59 !!---------------------------------------------------------------------- 60 USE oce, ONLY: CO2Flux_out_cpl, DMS_out_cpl, PCO2a_in_cpl 60 61 USE oce_trc 61 62 USE trc … … 96 97 PUBLIC trc_bio_medusa ! called in ??? 97 98 99 REAL(wp), PARAMETER :: dms_unit_conv = 1.0e+6 ! Coversion factor to get 100 ! outgong DMS in standard 101 ! units for coupling, i.e. 102 ! specifically nmol/L 103 !(= umol/m3) 98 104 !!* Substitution 99 105 # include "domzgr_substitute.h90" … … 357 363 REAL(wp), POINTER, DIMENSION(:,: ) :: oben_fe2d, oben_c2d, oben_si2d, oben_ca2d, sfr_ocal2d 358 364 REAL(wp), POINTER, DIMENSION(:,: ) :: sfr_oarg2d, lyso_ca2d 365 359 366 # endif 360 367 !! 2D var for diagnostics. … … 1538 1545 ENDIF 1539 1546 # if defined key_debug_medusa 1540 IF (lwp) write (numout,*) 'trc_bio_medusa: fin nish calculating dms'1547 IF (lwp) write (numout,*) 'trc_bio_medusa: finish calculating dms' 1541 1548 CALL flush(numout) 1542 1549 # endif … … 4278 4285 !!----------------------------------------- 4279 4286 IF (jk.eq.1) THEN 4287 4288 #if defined key_oasis3 4289 ! Capture outgoing CO2 flux coupling field 4290 CO2Flux_out_cpl(:,:) = f_co2flux2d(:,:) 4291 ! Ensure halos are populated in the couping field 4292 CALL lbc_lnk(CO2Flux_out_cpl(:,:),'T',1. ) 4293 4294 ! Capture outgoing DMS coupling field and convert 4295 ! values to standard units 4296 DMS_out_cpl(:,:) = dms_surf2d(:,:) * dms_unit_conv 4297 ! Ensure halos are populated in the couping field 4298 CALL lbc_lnk(DMS_out_cpl(:,:),'T',1. ) 4299 #endif 4300 4280 4301 # if defined key_debug_medusa 4281 4302 IF (lwp) write (numout,*) 'trc_bio_medusa: diag jk = 1' … … 4368 4389 IF( med_diag%DMS_SURF%dgsave ) THEN 4369 4390 CALL iom_put( "DMS_SURF" , dms_surf2d ) 4391 4392 4370 4393 CALL wrk_dealloc( jpi, jpj, dms_surf2d ) 4371 4394 ENDIF
Note: See TracChangeset
for help on using the changeset viewer.