Changeset 8039 for branches/UKMO/dev_r5518_medusa_chg_trc_bio_medusa/NEMOGCM/NEMO/TOP_SRC/MEDUSA/air_sea.F90
- Timestamp:
- 2017-05-18T11:14:31+02:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/dev_r5518_medusa_chg_trc_bio_medusa/NEMOGCM/NEMO/TOP_SRC/MEDUSA/air_sea.F90
r8023 r8039 39 39 dms_aran2d, dms_hall, dms_hall2d, & 40 40 dms_simo, dms_simo2d, dms_surf, & 41 dms_surf2d, f_BetaD,&42 f_co2flux, f_co2flux2d, f_co2starair,&41 dms_surf2d, & 42 f_co2flux, f_co2flux2d, & 43 43 f_co2starair_2d, f_co3, & 44 f_dcf, f_dpco2, f_fco2a_2d, f_fco2atm, & 45 f_fco2w, f_fco2w_2d, f_h2co3, & 46 f_hco3, f_henry, f_insitut, f_K0, & 47 f_kw660, f_kw6602d, f_kwco2, f_kwo2, & 48 f_o2flux, f_o2flux2d, & 49 f_o2sat, f_o2sat2d, f_ocndpco2_2d, & 44 f_dcf, f_fco2a_2d, f_fco2w_2d, & 45 f_h2co3, f_hco3, f_henry, & 46 f_kw660, f_kw6602d, & 47 f_o2flux, f_o2flux2d, f_o2sat, & 48 f_o2sat2d, f_ocndpco2_2d, & 50 49 f_ocnk0_2d, f_ocnkwco2_2d, & 51 50 f_ocnrhosw_2d, f_ocnschco2_2d, & 52 f_omarg, f_omcal, f_opres, & 53 f_pco2a2d, f_pco2atm, & 54 f_pco2w, f_pco2w2d, f_ph, f_pp0, & 55 f_pp02d, f_rhosw, & 56 f_schmidtco2, f_TALK, f_TALK2d, & 57 f_TDIC, f_TDIC2d, f_xco2a, f_xco2a_2d, & 58 iters, & 51 f_omarg, f_omcal, & 52 f_pco2a2d, f_pco2atm, f_pco2w, & 53 f_pco2w2d, f_ph, f_pp0, f_pp02d, & 54 f_TALK, f_TALK2d, f_TDIC, f_TDIC2d, & 55 f_xco2a, f_xco2a_2d, & 59 56 zalk, zdic, zoxy, zsal, ztmp, & 60 57 # endif … … 69 66 USE in_out_manager, ONLY: lwp, numout 70 67 USE oce, ONLY: PCO2a_in_cpl 71 USE par_oce, ONLY: jpim1, jpjm1 68 USE par_kind, ONLY: wp 69 USE par_oce, ONLY: jpi, jpim1, jpj, jpjm1 72 70 USE sbc_oce, ONLY: fr_i, lk_oasis, qsr, wndm 73 71 USE sms_medusa, ONLY: jdms, jdms_input, jdms_model, & … … 97 95 INTEGER, INTENT( in ) :: kt 98 96 97 !! Loop variables 98 INTEGER :: ji, jj 99 99 100 # if defined key_roam 100 101 !! jpalm 14-07-2016: convert CO2flux diag from mmol/m2/d to kg/m2/s … … 102 103 REAL, PARAMETER :: secs_in_day = 86400.0 !! s / d 103 104 REAL, PARAMETER :: CO2flux_conv = (1.e-6 * weight_CO2_mol) / secs_in_day 105 106 INTEGER :: iters 107 108 !! AXY (23/06/15): additional diagnostics for MOCSY and oxygen 109 REAL(wp), DIMENSION(jpi,jpj) :: f_fco2w, f_rhosw 110 REAL(wp), DIMENSION(jpi,jpj) :: f_fco2atm 111 REAL(wp), DIMENSION(jpi,jpj) :: f_schmidtco2, f_kwco2, f_K0 112 REAL(wp), DIMENSION(jpi,jpj) :: f_co2starair, f_dpco2 113 !! Output arguments from mocsy_interface, which aren't used 114 REAL(wp) :: f_BetaD_dum, f_opres_dum 115 REAL(wp) :: f_insitut_dum 116 REAL(wp) :: f_kwo2_dum 104 117 # endif 105 118 106 INTEGER :: ji, jj107 119 108 120 # if defined key_roam … … 171 183 f_h2co3(ji,jj),f_hco3(ji,jj), & 172 184 f_co3(ji,jj),f_omarg(ji,jj), & 173 f_omcal(ji,jj),f_BetaD (ji,jj),&174 f_rhosw(ji,jj),f_opres (ji,jj),&175 f_insitut (ji,jj),f_pco2atm(ji,jj),&185 f_omcal(ji,jj),f_BetaD_dum, & 186 f_rhosw(ji,jj),f_opres_dum, & 187 f_insitut_dum,f_pco2atm(ji,jj), & 176 188 f_fco2atm(ji,jj),f_schmidtco2(ji,jj), & 177 189 f_kwco2(ji,jj),f_K0(ji,jj), & … … 192 204 DO ji = 2,jpim1 193 205 if (tmask(ji,jj,1) == 1) then 194 iters (ji,jj)= 0206 iters = 0 195 207 !! 196 208 !! carbon dioxide (CO2); Jerry Blackford code (ostensibly … … 205 217 f_co2flux(ji,jj),f_TDIC(ji,jj), & 206 218 f_TALK(ji,jj),f_dcf(ji,jj), & 207 f_henry(ji,jj),iters (ji,jj))219 f_henry(ji,jj),iters) 208 220 !! 209 221 !! AXY (09/01/14): removed iteration and NaN checks; these have … … 213 225 !! output warnings are retained here so that 214 226 !! failure position can be determined 215 if (iters (ji,jj).eq. 25) then227 if (iters .eq. 25) then 216 228 IF(lwp) WRITE(numout,*) ' trc_bio_medusa: ITERS WARNING, ', & 217 iters (ji,jj), ' AT (', ji, ', ', jj, ', 1) AT ', kt229 iters, ' AT (', ji, ', ', jj, ', 1) AT ', kt 218 230 endif 219 231 ENDIF … … 244 256 !! AXY (23/06/15): add in some extra MOCSY diagnostics 245 257 f_fco2w(ji,jj) = f_xco2a(ji,jj) 246 f_BetaD(ji,jj) = 1. 258 ! This doesn't seem to be used - marc 16/5/17 259 ! f_BetaD(ji,jj) = 1. 247 260 f_rhosw(ji,jj) = 1.026 248 f_opres(ji,jj) = 0. 249 f_insitut(ji,jj) = ztmp(ji,jj) 261 ! This doesn't seem to be used - marc 16/5/17 262 ! f_opres(ji,jj) = 0. 263 ! f_insitut(ji,jj) = ztmp(ji,jj) 250 264 f_pco2atm(ji,jj) = f_xco2a(ji,jj) 251 265 f_fco2atm(ji,jj) = f_xco2a(ji,jj) … … 271 285 !! AXY (23/06/15): amend input list for oxygen to account 272 286 !! for common gas transfer velocity 273 !! Note that f_kwo2 is an about from the subroutine below,274 !! which doesn't seem to be used - marc 10/4/17275 287 CALL trc_oxy_medusa(ztmp(ji,jj),zsal(ji,jj),f_kw660(ji,jj), & 276 288 f_pp0(ji,jj),zoxy(ji,jj), & 277 f_kwo2 (ji,jj),f_o2flux(ji,jj),&289 f_kwo2_dum,f_o2flux(ji,jj), & 278 290 f_o2sat(ji,jj)) 279 291 !!
Note: See TracChangeset
for help on using the changeset viewer.