Changeset 11046
- Timestamp:
- 2019-05-23T18:02:59+02:00 (5 years ago)
- Location:
- branches/UKMO/dev_r5518_GO6_fix_zemp_ice_10681/NEMOGCM/NEMO/OPA_SRC/SBC
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/dev_r5518_GO6_fix_zemp_ice_10681/NEMOGCM/NEMO/OPA_SRC/SBC/sbccpl.F90
r11029 r11046 186 186 LOGICAL :: ln_usecplmask ! use a coupling mask file to merge data received from several models 187 187 ! -> file cplmask.nc with the float variable called cplmask (jpi,jpj,nn_cplmodel) 188 LOGICAL, PUBLIC :: ln_fix_sea_ice_fluxes ! Apply sea ice flux bug fixes (GMED#449)189 188 190 189 TYPE :: DYNARR … … 263 262 & ln_usecplmask, nn_coupled_iceshelf_fluxes, ln_iceshelf_init_atmos, & 264 263 & rn_greenland_total_fw_flux, rn_greenland_calving_fraction, & 265 & rn_antarctica_total_fw_flux, rn_antarctica_calving_fraction, rn_iceshelf_fluxes_tolerance, & 266 & ln_fix_sea_ice_fluxes 264 & rn_antarctica_total_fw_flux, rn_antarctica_calving_fraction, rn_iceshelf_fluxes_tolerance 267 265 !!--------------------------------------------------------------------- 268 266 … … 335 333 WRITE(numout,*)' nn_cplmodel = ', nn_cplmodel 336 334 WRITE(numout,*)' ln_usecplmask = ', ln_usecplmask 337 WRITE(numout,*)' ln_fix_sea_ice_fluxes = ', ln_fix_sea_ice_fluxes338 335 WRITE(numout,*)' nn_coupled_iceshelf_fluxes = ', nn_coupled_iceshelf_fluxes 339 336 WRITE(numout,*)' ln_iceshelf_init_atmos = ', ln_iceshelf_init_atmos … … 1845 1842 IF ( TRIM(sn_rcv_emp%clcat) == 'yes' ) THEN 1846 1843 ! zemp_ice is the sum of frcv(jpr_ievp)%z3(:,:,1) over all layers - snow 1847 IF ( ln_fix_sea_ice_fluxes ) THEN 1848 zemp_ice(:,:) = - frcv(jpr_snow)%z3(:,:,1) * zicefr(:,:) 1849 ELSE 1850 zemp_ice(:,:) = - frcv(jpr_snow)%z3(:,:,1) 1851 ENDIF 1844 zemp_ice(:,:) = - frcv(jpr_snow)%z3(:,:,1) * zicefr(:,:) 1852 1845 DO jl=1,jpl 1853 IF ( ln_fix_sea_ice_fluxes ) THEN 1854 zemp_ice(:,: ) = zemp_ice(:,:) + frcv(jpr_ievp)%z3(:,:,jl) * a_i_last_couple(:,:,jl) 1855 ELSE 1856 zemp_ice(:,: ) = zemp_ice(:,:) + frcv(jpr_ievp)%z3(:,:,jl) 1857 ENDIF 1846 zemp_ice(:,: ) = zemp_ice(:,:) + frcv(jpr_ievp)%z3(:,:,jl) * a_i_last_couple(:,:,jl) 1858 1847 ENDDO 1859 1848 ! latent heat coupled for each category in CICE -
branches/UKMO/dev_r5518_GO6_fix_zemp_ice_10681/NEMOGCM/NEMO/OPA_SRC/SBC/sbcice_cice.F90
r11029 r11046 395 395 ELSE IF (ksbc == jp_purecpl) THEN 396 396 DO jl=1,ncat 397 IF ( ln_fix_sea_ice_fluxes ) THEN 398 ztmpn(:,:,jl)=qla_ice(:,:,jl)*a_i_last_couple(:,:,jl) 399 ELSE 400 ztmpn(:,:,jl)=qla_ice(:,:,jl)*a_i(:,:,jl) 401 ENDIF 397 ztmpn(:,:,jl)=qla_ice(:,:,jl)*a_i_last_couple(:,:,jl) 402 398 ENDDO 403 399 ELSE … … 412 408 ! Convert to GBM 413 409 IF (ksbc == jp_flx .OR. ksbc == jp_purecpl) THEN 414 IF ( ln_fix_sea_ice_fluxes ) THEN 415 ztmp(:,:) = botmelt(:,:,jl)*a_i_last_couple(:,:,jl) 416 ELSE 417 ztmp(:,:) = botmelt(:,:,jl)*a_i(:,:,jl) 418 ENDIF 410 ztmp(:,:) = botmelt(:,:,jl)*a_i_last_couple(:,:,jl) 419 411 ELSE 420 412 ztmp(:,:) = botmelt(:,:,jl) … … 425 417 ! Convert to GBM 426 418 IF (ksbc == jp_flx .OR. ksbc == jp_purecpl) THEN 427 IF ( ln_fix_sea_ice_fluxes ) THEN 428 ztmp(:,:) = (topmelt(:,:,jl)+botmelt(:,:,jl))*a_i_last_couple(:,:,jl) 429 ELSE 430 ztmp(:,:) = (topmelt(:,:,jl)+botmelt(:,:,jl))*a_i(:,:,jl) 431 ENDIF 419 ztmp(:,:) = (topmelt(:,:,jl)+botmelt(:,:,jl))*a_i_last_couple(:,:,jl) 432 420 ELSE 433 421 ztmp(:,:) = (topmelt(:,:,jl)+botmelt(:,:,jl))
Note: See TracChangeset
for help on using the changeset viewer.