Changeset 11500
- Timestamp:
- 2019-09-04T15:49:41+02:00 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/UKMO/NEMO_4.0_GC_couple_pkg_cons_checks/src/OCE/SBC/sbccpl.F90
r11436 r11500 1666 1666 ! 1667 1667 INTEGER :: ji, jj, jl ! dummy loop index 1668 INTEGER :: itx ! index of flux over ice 1668 1669 REAL(wp) :: ztri ! local scalar 1669 1670 REAL(wp), DIMENSION(jpi,jpj) :: zcptn, zcptrain, zcptsnw, ziceld, zmsk, zsnw … … 1673 1674 !!---------------------------------------------------------------------- 1674 1675 ! 1676 IF( TRIM( sn_rcv_emp%cldes ) == 'conservative' ) THEN ; itx=jpr_snow ; 1677 ELSE ; itx=jpr_semp ; ENDIF 1678 1679 IF( nrcvinfo(itx) == OASIS_Rcv ) THEN 1680 1675 1681 IF( ln_mixcpl ) zmsk(:,:) = 1. - xcplmask(:,:,0) 1676 1682 ziceld(:,:) = 1._wp - picefr(:,:) … … 1790 1796 ! 1791 1797 #endif 1798 ENDIF 1792 1799 1793 1800 ! outputs … … 1806 1813 ! note: runoff output is done in sbcrnf (which includes icebergs too) and iceshelf output is done in sbcisf 1807 1814 ! 1815 !TG comment - the variable used here should depend on the type of coupling 1816 !used for qns 1817 IF( TRIM( sn_rcv_qns%cldes ) == 'oce only' ) THEN ; itx=jpr_qsnoce ; 1818 ELSE ; itx=jpr_qnsice ; ENDIF 1819 1820 IF( nrcvinfo(itx) == OASIS_Rcv ) THEN 1808 1821 ! ! ========================= ! 1809 1822 SELECT CASE( TRIM( sn_rcv_qns%cldes ) ) ! non solar heat fluxes ! (qns) … … 1923 1936 1924 1937 #endif 1925 ! outputs 1938 ENDIF !End of code only run when qns received 1939 1940 ! outputs - these must be called every time step even if coupling hasn't 1941 ! happened. 1926 1942 IF ( srcv(jpr_cal)%laction ) CALL iom_put('hflx_cal_cea' , - frcv(jpr_cal)%z3(:,:,1) * rLfus ) ! latent heat from calving 1927 1943 IF ( srcv(jpr_icb)%laction ) CALL iom_put('hflx_icb_cea' , - frcv(jpr_icb)%z3(:,:,1) * rLfus ) ! latent heat from icebergs melting … … 1936 1952 ! note: hflx for runoff and iceshelf are done in sbcrnf and sbcisf resp. 1937 1953 ! 1954 1955 IF( TRIM( sn_rcv_qsr%cldes ) == 'oce only') THEN ; itx=jpr_qsroce 1956 ELSE ; itx = jpr_qsrice ; ENDIF 1957 IF( nrcvinfo(itx) == OASIS_Rcv ) THEN 1938 1958 ! ! ========================= ! 1939 1959 SELECT CASE( TRIM( sn_rcv_qsr%cldes ) ) ! solar heat fluxes ! (qsr) … … 2059 2079 ! 2060 2080 #endif 2081 ENDIF !End of test for coupled solar flux 2061 2082 ! 2062 2083 END SUBROUTINE sbc_cpl_ice_flx
Note: See TracChangeset
for help on using the changeset viewer.