- Timestamp:
- 2017-06-09T15:38:39+02:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/dev_r5518_GO6_MEDUSA_dummyrun/NEMOGCM/NEMO/OPA_SRC/SBC/sbccpl.F90
r8046 r8159 34 34 USE geo2ocean ! 35 35 USE oce , ONLY : tsn, un, vn, sshn, ub, vb, sshb, fraqsr_1lev, & 36 CO2Flux_out_cpl, DMS_out_cpl, PCO2a_in_cpl, Dust_in_cpl, & 36 CO2Flux_out_cpl, DMS_out_cpl, chloro_out_cpl, & 37 PCO2a_in_cpl, Dust_in_cpl, & 37 38 ln_medusa 38 39 USE albedo ! … … 145 146 INTEGER, PARAMETER :: jps_sstfrz = 32 ! sea-surface freezing temperature 146 147 INTEGER, PARAMETER :: jps_fice1 = 33 ! first-order ice concentration (for time-travelling ice coupling) 147 INTEGER, PARAMETER :: jps_bio_co2 = 34 ! MEDUSA air-sea CO2 flux in 148 INTEGER, PARAMETER :: jps_bio_dms = 35 ! MEDUSA DMS surface concentration in 149 INTEGER, PARAMETER :: jpsnd = 35 ! total number of fields sent 148 INTEGER, PARAMETER :: jps_bio_co2 = 34 ! MEDUSA air-sea CO2 flux 149 INTEGER, PARAMETER :: jps_bio_dms = 35 ! MEDUSA DMS surface concentration 150 INTEGER, PARAMETER :: jps_bio_chloro = 36 ! MEDUSA chlorophyll surface concentration 151 INTEGER, PARAMETER :: jpsnd = 36 ! total number of fields sent 150 152 151 153 REAL(wp), PARAMETER :: dms_unit_conv = 1.0e+6 ! Coversion factor to get outgong DMS in standard units for coupling … … 162 164 ! Send to the atmosphere ! 163 165 TYPE(FLD_C) :: sn_snd_temp, sn_snd_alb, sn_snd_thick, sn_snd_crt, sn_snd_co2, sn_snd_cond, sn_snd_mpnd, sn_snd_sstfrz, sn_snd_thick1 164 TYPE(FLD_C) :: sn_snd_bio_co2, sn_snd_bio_dms 166 TYPE(FLD_C) :: sn_snd_bio_co2, sn_snd_bio_dms, sn_snd_bio_chloro 165 167 166 168 ! Received from the atmosphere ! … … 246 248 247 249 ! Add MEDUSA related fields to namelist 248 NAMELIST/namsbc_cpl/ sn_snd_bio_co2, sn_snd_bio_dms, 250 NAMELIST/namsbc_cpl/ sn_snd_bio_co2, sn_snd_bio_dms, sn_snd_bio_chloro, & 249 251 & sn_rcv_atm_pco2, sn_rcv_atm_dust 250 252 … … 304 306 WRITE(numout,*)' bio co2 flux = ', TRIM(sn_snd_bio_co2%cldes), ' (', TRIM(sn_snd_bio_co2%clcat), ')' 305 307 WRITE(numout,*)' bio dms flux = ', TRIM(sn_snd_bio_dms%cldes), ' (', TRIM(sn_snd_bio_dms%clcat), ')' 308 WRITE(numout,*)' bio dms chlorophyll = ', TRIM(sn_snd_bio_chloro%cldes), ' (', TRIM(sn_snd_bio_chloro%clcat), ')' 306 309 WRITE(numout,*)' oce co2 flux = ', TRIM(sn_snd_co2%cldes ), ' (', TRIM(sn_snd_co2%clcat ), ')' 307 310 WRITE(numout,*)' ice effective conductivity = ', TRIM(sn_snd_cond%cldes ), ' (', TRIM(sn_snd_cond%clcat ), ')' … … 826 829 IF( TRIM(sn_snd_bio_co2%cldes) == 'medusa' ) ssnd(jps_bio_co2 )%laction = .TRUE. 827 830 831 ! Surface chlorophyll from Medusa 832 ssnd(jps_bio_chloro)%clname = 'OBioChlo' 833 IF( TRIM(sn_snd_bio_chloro%cldes) == 'medusa' ) ssnd(jps_bio_chloro )%laction = .TRUE. 834 828 835 ! ! ------------------------- ! 829 836 ! ! Sea surface freezing temp ! … … 1045 1052 REAL(wp), POINTER, DIMENSION(:,:) :: ztx, zty, zmsk, zemp, zqns, zqsr 1046 1053 !!---------------------------------------------------------------------- 1047 1048 ! RSRH temporary arrays for testing, just to recieve incoming MEDUSA related fields1049 ! until we know where they need to go.1050 REAL(wp), ALLOCATABLE :: atm_pco2(:,:)1051 REAL(wp), ALLOCATABLE :: atm_dust(:,:)1052 1054 1053 1055 ! … … 2347 2349 2348 2350 IF (ln_medusa) THEN 2349 ! ! --------------------------------- !2350 ! ! CO2 flux and DMSfrom MEDUSA !2351 ! ! --------------------------------- !2351 ! ! ---------------------------------------------- ! 2352 ! ! CO2 flux, DMS and chlorophyll from MEDUSA ! 2353 ! ! ---------------------------------------------- ! 2352 2354 IF ( ssnd(jps_bio_co2)%laction ) THEN 2353 2355 CALL cpl_snd( jps_bio_co2, isec, RESHAPE( CO2Flux_out_cpl, (/jpi,jpj,1/) ), info ) … … 2356 2358 IF ( ssnd(jps_bio_dms)%laction ) THEN 2357 2359 CALL cpl_snd( jps_bio_dms, isec, RESHAPE( DMS_out_cpl, (/jpi,jpj,1/) ), info ) 2360 ENDIF 2361 2362 IF ( ssnd(jps_bio_chloro)%laction ) THEN 2363 CALL cpl_snd( jps_bio_chloro, isec, RESHAPE( chloro_out_cpl, (/jpi,jpj,1/) ), info ) 2358 2364 ENDIF 2359 2365 ENDIF
Note: See TracChangeset
for help on using the changeset viewer.