Ignore:
Timestamp:
2019-06-14T11:50:40+02:00 (18 months ago)
Author:
frrh
Message:

Commit changes from Dan Copsey's sea ice heat coupling
flux fixes in branch:
branches/UKMO/dev_r5518_GO6_fix_zemp_ice_10681
revisions 11028:11088.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/SBC/sbcice_cice.F90

    r11101 r11107  
    401401         ELSE IF (ksbc == jp_purecpl) THEN 
    402402            DO jl=1,ncat 
    403                ztmpn(:,:,jl)=qla_ice(:,:,jl)*a_i(:,:,jl) 
     403               IF ( ln_fix_sea_ice_fluxes ) THEN 
     404                  ztmpn(:,:,jl)=qla_ice(:,:,jl)*a_i_last_couple(:,:,jl) 
     405               ELSE 
     406                  ztmpn(:,:,jl)=qla_ice(:,:,jl)*a_i(:,:,jl) 
     407               ENDIF 
    404408            ENDDO 
    405409    ELSE 
     
    414418!  Convert to GBM 
    415419            IF (ksbc == jp_flx .OR. ksbc == jp_purecpl) THEN 
    416                ztmp(:,:) = botmelt(:,:,jl)*a_i(:,:,jl) 
     420               IF ( ln_fix_sea_ice_fluxes ) THEN 
     421                  ztmp(:,:) = botmelt(:,:,jl)*a_i_last_couple(:,:,jl) 
     422               ELSE 
     423                  ztmp(:,:) = botmelt(:,:,jl)*a_i(:,:,jl) 
     424               ENDIF 
    417425            ELSE 
    418426               ztmp(:,:) = botmelt(:,:,jl) 
     
    423431!  Convert to GBM 
    424432            IF (ksbc == jp_flx .OR. ksbc == jp_purecpl) THEN 
    425                ztmp(:,:) = (topmelt(:,:,jl)+botmelt(:,:,jl))*a_i(:,:,jl)  
     433               IF ( ln_fix_sea_ice_fluxes ) THEN 
     434                  ztmp(:,:) = (topmelt(:,:,jl)+botmelt(:,:,jl))*a_i_last_couple(:,:,jl)  
     435               ELSE 
     436                  ztmp(:,:) = (topmelt(:,:,jl)+botmelt(:,:,jl))*a_i(:,:,jl) 
     437               ENDIF 
    426438            ELSE 
    427439               ztmp(:,:) = (topmelt(:,:,jl)+botmelt(:,:,jl)) 
Note: See TracChangeset for help on using the changeset viewer.