Changeset 8159
- Timestamp:
- 2017-06-09T15:38:39+02:00 (7 years ago)
- Location:
- branches/UKMO/dev_r5518_GO6_MEDUSA_dummyrun/NEMOGCM
- Files:
-
- 4 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/dev_r5518_GO6_MEDUSA_dummyrun/NEMOGCM/CONFIG/SHARED/namelist_ref
r8046 r8159 363 363 sn_snd_crt = 'none' , 'no' , 'spherical' , 'eastward-northward' , 'T' 364 364 sn_snd_co2 = 'coupled' , 'no' , '' , '' , '' 365 sn_snd_bio_co2 = 'medusa' , 'no' , '' , '' , '' 366 sn_snd_bio_dms = 'medusa' , 'no' , '' , '' , '' 367 sn_snd_bio_chloro = 'medusa' , 'no' , '' , '' , '' 368 sn_snd_cond = 'weighted ice' , 'no' , '' , '' , '' 369 sn_snd_mpnd = 'ice only' , 'no' , '' , '' , '' 370 sn_snd_sstfrz = 'coupled' , 'no' , '' , '' , '' 371 sn_snd_thick1 = 'ice and snow' , 'no' , '' , '' , '' 365 372 ! receive 366 373 sn_rcv_w10m = 'none' , 'no' , '' , '' , '' … … 374 381 sn_rcv_cal = 'coupled' , 'no' , '' , '' , '' 375 382 sn_rcv_co2 = 'coupled' , 'no' , '' , '' , '' 383 sn_rcv_antm = 'coupled' , 'no' , '' , '' , '' 384 sn_rcv_grnm = 'coupled' , 'no' , '' , '' , '' 385 sn_rcv_iceflx = 'coupled' , 'no' , '' , '' , '' 386 sn_rcv_ts_ice = 'ice' , 'no' , '' , '' , '' 387 sn_rcv_atm_dust = 'medusa' , 'no' , '' , '' , '' 388 sn_rcv_atm_pco2 = 'medusa' , 'no' , '' , '' , '' 389 376 390 ! 377 391 nn_cplmodel = 1 ! Maximum number of models to/from which NEMO is potentialy sending/receiving data -
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 -
branches/UKMO/dev_r5518_GO6_MEDUSA_dummyrun/NEMOGCM/NEMO/OPA_SRC/nemogcm.F90
r8155 r8159 610 610 IF( numdct_heat /= -1 ) CLOSE( numdct_heat ) ! heat transports 611 611 IF( numdct_salt /= -1 ) CLOSE( numdct_salt ) ! salt transports 612 #if defined key_top 612 613 IF( numstr /= -1 ) CLOSE( numstr ) ! tracer statistics 613 614 #endif 614 615 ! 615 616 numout = 6 ! redefine numout in case it is used after this point... -
branches/UKMO/dev_r5518_GO6_MEDUSA_dummyrun/NEMOGCM/NEMO/OPA_SRC/oce.F90
r7770 r8159 77 77 REAL(wp), PUBLIC, ALLOCATABLE, SAVE :: CO2Flux_out_cpl(:,:) ! Output coupling CO2 flux 78 78 REAL(wp), PUBLIC, ALLOCATABLE, SAVE :: DMS_out_cpl(:,:) ! Output coupling DMS 79 REAL(wp), PUBLIC, ALLOCATABLE, SAVE :: chloro_out_cpl(:,:) ! Output coupling chlorophyll 80 ! (expected in Kg/M3) 79 81 80 82 REAL(wp), PUBLIC, ALLOCATABLE, SAVE :: PCO2a_in_cpl(:,:) ! Input coupling CO2 partial pressure … … 138 140 ! are enabled 139 141 ALLOCATE( CO2Flux_out_cpl(jpi,jpj), DMS_out_cpl(jpi,jpj), & 142 chloro_out_cpl(jpi,jpj), & 140 143 PCO2a_in_cpl(jpi,jpj), Dust_in_cpl(jpi,jpj), STAT=ierr(5) ) 141 144
Note: See TracChangeset
for help on using the changeset viewer.