Changeset 8434 for branches/NERC/dev_r5518_GO6_split_trcbiomedusa/NEMOGCM/NEMO/TOP_SRC/MEDUSA/air_sea.F90
- Timestamp:
- 2017-08-11T12:47:19+02:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/NERC/dev_r5518_GO6_split_trcbiomedusa/NEMOGCM/NEMO/TOP_SRC/MEDUSA/air_sea.F90
r8395 r8434 124 124 !! Air-sea gas exchange 125 125 !!----------------------------------------------------------- 126 127 # if defined key_debug_medusa 128 IF (lwp) write (numout,*) & 129 'air-sea: gas_transfer kt = ', kt 130 CALL flush(numout) 131 # endif 126 132 DO jj = 2,jpjm1 127 133 DO ji = 2,jpim1 … … 145 151 !! Wanninkhof (2014), option 7 146 152 !! 153 CALL gas_transfer( wndm(ji,jj), 1, 7, & ! inputs 154 f_kw660(ji,jj) ) ! outputs 155 ENDIF 156 ENDDO 157 ENDDO 158 147 159 # if defined key_debug_medusa 148 IF (lwp) write (numout,*) 'trc_bio_medusa: entering gas_transfer' 160 IF (lwp) write (numout,*) & 161 'air-sea: carb-chem kt = ', kt 149 162 CALL flush(numout) 150 163 # endif 151 CALL gas_transfer( wndm(ji,jj), 1, 7, & ! inputs152 f_kw660(ji,jj) ) ! outputs153 # if defined key_debug_medusa154 IF (lwp) write (numout,*) 'trc_bio_medusa: exiting gas_transfer'155 CALL flush(numout)156 # endif157 ENDIF158 ENDDO159 ENDDO160 161 164 DO jj = 2,jpjm1 162 165 DO ji = 2,jpim1 … … 228 231 !! failure position can be determined 229 232 if (iters .eq. 25) then 230 IF(lwp) WRITE(numout,*) ' trc_bio_medusa: ITERS WARNING, ',&233 IF(lwp) WRITE(numout,*) 'air-sea: ITERS WARNING, ', & 231 234 iters, ' AT (', ji, ', ', jj, ', 1) AT ', kt 232 235 endif … … 335 338 hmld(ji,jj),qsr(ji,jj), & 336 339 zdin(ji,jj), dms_nlim(ji,jj), & 337 dms_andr(ji,jj),dms_simo(ji,jj), & 338 dms_aran(ji,jj),dms_hall(ji,jj), & 339 dms_andm(ji,jj)) 340 dms_andr,dms_simo,dms_aran,dms_hall, & 341 dms_andm) 340 342 else 341 343 !! use diel-average inputs … … 346 348 zn_dms_mld(ji,jj),zn_dms_qsr(ji,jj), & 347 349 zn_dms_din(ji,jj),dms_nlim(ji,jj), & 348 dms_andr(ji,jj),dms_simo(ji,jj), & 349 dms_aran(ji,jj),dms_hall(ji,jj), & 350 dms_andm(ji,jj)) 350 dms_andr,dms_simo,dms_aran,dms_hall, & 351 dms_andm) 351 352 endif 352 353 !! 353 354 !! assign correct output to variable passed to atmosphere 354 355 if (jdms_model .eq. 1) then 355 dms_surf (ji,jj) = dms_andr(ji,jj)356 dms_surf = dms_andr 356 357 elseif (jdms_model .eq. 2) then 357 dms_surf (ji,jj) = dms_simo(ji,jj)358 dms_surf = dms_simo 358 359 elseif (jdms_model .eq. 3) then 359 dms_surf (ji,jj) = dms_aran(ji,jj)360 dms_surf = dms_aran 360 361 elseif (jdms_model .eq. 4) then 361 dms_surf (ji,jj) = dms_hall(ji,jj)362 dms_surf = dms_hall 362 363 elseif (jdms_model .eq. 5) then 363 dms_surf (ji,jj) = dms_andm(ji,jj)364 dms_surf = dms_andm 364 365 endif 365 366 !! 366 367 !! 2D diag through iom_use 367 IF( lk_iomput ) THEN 368 IF( med_diag%DMS_SURF%dgsave ) THEN 369 dms_surf2d(ji,jj) = dms_surf(ji,jj) 370 ENDIF 371 IF( med_diag%DMS_ANDR%dgsave ) THEN 372 dms_andr2d(ji,jj) = dms_andr(ji,jj) 373 ENDIF 374 IF( med_diag%DMS_SIMO%dgsave ) THEN 375 dms_simo2d(ji,jj) = dms_simo(ji,jj) 376 ENDIF 377 IF( med_diag%DMS_ARAN%dgsave ) THEN 378 dms_aran2d(ji,jj) = dms_aran(ji,jj) 379 ENDIF 380 IF( med_diag%DMS_HALL%dgsave ) THEN 381 dms_hall2d(ji,jj) = dms_hall(ji,jj) 382 ENDIF 383 IF( med_diag%DMS_ANDM%dgsave ) THEN 384 dms_andm2d(ji,jj) = dms_andm(ji,jj) 385 ENDIF 386 # if defined key_debug_medusa 387 IF (lwp) write (numout,*) & 388 'trc_bio_medusa: finish calculating dms' 389 CALL flush(numout) 390 # endif 391 ENDIF !! End iom 368 IF( med_diag%DMS_SURF%dgsave ) THEN 369 dms_surf2d(ji,jj) = dms_surf 370 ENDIF 371 IF( med_diag%DMS_ANDR%dgsave ) THEN 372 dms_andr2d(ji,jj) = dms_andr 373 ENDIF 374 IF( med_diag%DMS_SIMO%dgsave ) THEN 375 dms_simo2d(ji,jj) = dms_simo 376 ENDIF 377 IF( med_diag%DMS_ARAN%dgsave ) THEN 378 dms_aran2d(ji,jj) = dms_aran 379 ENDIF 380 IF( med_diag%DMS_HALL%dgsave ) THEN 381 dms_hall2d(ji,jj) = dms_hall 382 ENDIF 383 IF( med_diag%DMS_ANDM%dgsave ) THEN 384 dms_andm2d(ji,jj) = dms_andm 385 ENDIF 392 386 ENDIF 393 387 ENDDO 394 388 ENDDO 389 # if defined key_debug_medusa 390 IF (lwp) write (numout,*) & 391 'air-sea: finish calculating dms kt = ',kt 392 CALL flush(numout) 393 # endif 395 394 ENDIF !! End IF (jdms == 1) 396 397 395 398 396 !!
Note: See TracChangeset
for help on using the changeset viewer.