Changeset 12733
- Timestamp:
- 2020-04-10T13:07:50+02:00 (5 years ago)
- Location:
- NEMO/branches/2020/r4.0-HEAD_r12713_clem_dan_fixcpl
- Files:
-
- 5 edited
-
cfgs/SHARED/namelist_ref (modified) (1 diff)
-
doc/namelists/namsbc_cpl (modified) (1 diff)
-
src/ICE/ice.F90 (modified) (3 diffs)
-
src/OCE/SBC/sbc_ice.F90 (modified) (1 diff)
-
src/OCE/SBC/sbccpl.F90 (modified) (21 diffs)
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2020/r4.0-HEAD_r12713_clem_dan_fixcpl/cfgs/SHARED/namelist_ref
r12720 r12733 286 286 &namsbc_cpl ! coupled ocean/atmosphere model ("key_oasis3") 287 287 !----------------------------------------------------------------------- 288 nn_cplmodel = 1 ! Maximum number of models to/from which NEMO is potentially sending/receiving data 289 ln_usecplmask = .false. ! use a coupling mask file to merge data received from several models 290 ! ! -> file cplmask.nc with the float variable called cplmask (jpi,jpj,nn_cplmodel) 291 nn_cats_cpl = 5 ! Number of sea ice categories over which coupling is to be carried out (if not 1) 288 nn_cplmodel = 1 ! Maximum number of models to/from which NEMO is potentially sending/receiving data 289 ln_usecplmask = .false. ! use a coupling mask file to merge data received from several models 290 ! ! -> file cplmask.nc with the float variable called cplmask (jpi,jpj,nn_cplmodel) 291 ln_scale_ice_flux = .false. ! use ice fluxes that are already "ice weighted" ( i.e. multiplied ice concentration) 292 nn_cats_cpl = 5 ! Number of sea ice categories over which coupling is to be carried out (if not 1) 292 293 !_____________!__________________________!____________!_____________!______________________!________! 293 294 ! ! description ! multiple ! vector ! vector ! vector ! -
NEMO/branches/2020/r4.0-HEAD_r12713_clem_dan_fixcpl/doc/namelists/namsbc_cpl
r10075 r12733 2 2 &namsbc_cpl ! coupled ocean/atmosphere model ("key_oasis3") 3 3 !----------------------------------------------------------------------- 4 nn_cplmodel = 1 ! Maximum number of models to/from which NEMO is potentially sending/receiving data5 ln_usecplmask = .false. ! use a coupling mask file to merge data received from several models6 ! ! -> file cplmask.nc with the float variable called cplmask (jpi,jpj,nn_cplmodel)7 nn_cats_cpl = 5 ! Number of sea ice categories over which coupling is to be carried out (if not 1)8 4 nn_cplmodel = 1 ! Maximum number of models to/from which NEMO is potentially sending/receiving data 5 ln_usecplmask = .false. ! use a coupling mask file to merge data received from several models 6 ! ! -> file cplmask.nc with the float variable called cplmask (jpi,jpj,nn_cplmodel) 7 ln_scale_ice_flux = .false. ! use ice fluxes that are already "ice weighted" ( i.e. multiplied ice concentration) 8 nn_cats_cpl = 5 ! Number of sea ice categories over which coupling is to be carried out (if not 1) 9 9 !_____________!__________________________!____________!_____________!______________________!________! 10 10 ! ! description ! multiple ! vector ! vector ! vector ! -
NEMO/branches/2020/r4.0-HEAD_r12713_clem_dan_fixcpl/src/ICE/ice.F90
r12726 r12733 398 398 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: qcn_ice_top !: Surface conduction flux (W/m2) 399 399 400 !!---------------------------------------------------------------------- 401 !! * Only for atmospheric coupling 402 !!---------------------------------------------------------------------- 403 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: a_i_last_couple !: Ice fractional area at last coupling time 400 404 ! 401 405 !!---------------------------------------------------------------------- … … 412 416 INTEGER :: ice_alloc 413 417 ! 414 INTEGER :: ierr(1 6), ii418 INTEGER :: ierr(17), ii 415 419 !!----------------------------------------------------------------- 416 420 ierr(:) = 0 … … 494 498 ALLOCATE( t_si(jpi,jpj,jpl) , tm_si(jpi,jpj) , qcn_ice_bot(jpi,jpj,jpl) , qcn_ice_top(jpi,jpj,jpl) , STAT = ierr(ii) ) 495 499 500 ! * For atmospheric coupling 501 ii = ii + 1 502 ALLOCATE( a_i_last_couple(jpi,jpj,jpl) , STAT=ierr(ii) ) 503 496 504 ice_alloc = MAXVAL( ierr(:) ) 497 505 IF( ice_alloc /= 0 ) CALL ctl_stop( 'STOP', 'ice_alloc: failed to allocate arrays.' ) 498 506 ! 507 499 508 END FUNCTION ice_alloc 500 509 -
NEMO/branches/2020/r4.0-HEAD_r12713_clem_dan_fixcpl/src/OCE/SBC/sbc_ice.F90
r12395 r12733 89 89 ! variables used in the coupled interface 90 90 INTEGER , PUBLIC, PARAMETER :: jpl = ncat 91 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: u_ice, v_ice ! jpi, jpj91 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: u_ice, v_ice 92 92 93 93 ! already defined in ice.F90 for SI3 94 94 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: a_i 95 95 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: h_i, h_s 96 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: a_i_last_couple !: Sea ice fraction on categories at the last coupling point 96 97 97 98 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: tatm_ice !: air temperature [K] -
NEMO/branches/2020/r4.0-HEAD_r12713_clem_dan_fixcpl/src/OCE/SBC/sbccpl.F90
r12725 r12733 48 48 USE lib_mpp ! distribued memory computing library 49 49 USE lbclnk ! ocean lateral boundary conditions (or mpp link) 50 51 #if defined key_oasis3 52 USE mod_oasis, ONLY : OASIS_Sent, OASIS_ToRest, OASIS_SentOut, OASIS_ToRestOut 53 #endif 50 54 51 55 IMPLICIT NONE … … 152 156 INTEGER, PARAMETER :: jps_wlev = 32 ! water level 153 157 INTEGER, PARAMETER :: jps_fice1 = 33 ! first-order ice concentration (for semi-implicit coupling of atmos-ice fluxes) 154 INTEGER, PARAMETER :: jps_a_p = 34 ! meltpond area 158 INTEGER, PARAMETER :: jps_a_p = 34 ! meltpond area fraction 155 159 INTEGER, PARAMETER :: jps_ht_p = 35 ! meltpond thickness 156 160 INTEGER, PARAMETER :: jps_kice = 36 ! sea ice effective conductivity … … 159 163 160 164 INTEGER, PARAMETER :: jpsnd = 38 ! total number of fields sent 165 166 #if ! defined key_oasis3 167 ! Dummy variables to enable compilation when oasis3 is not being used 168 INTEGER :: OASIS_Sent = -1 169 INTEGER :: OASIS_SentOut = -1 170 INTEGER :: OASIS_ToRest = -1 171 INTEGER :: OASIS_ToRestOut = -1 172 #endif 161 173 162 174 ! !!** namelist namsbc_cpl ** … … 184 196 LOGICAL :: ln_usecplmask ! use a coupling mask file to merge data received from several models 185 197 ! -> file cplmask.nc with the float variable called cplmask (jpi,jpj,nn_cplmodel) 198 LOGICAL :: ln_scale_ice_flux ! use ice fluxes that are already "ice weighted" ( i.e. multiplied ice concentration) 199 186 200 TYPE :: DYNARR 187 201 REAL(wp), POINTER, DIMENSION(:,:,:) :: z3 … … 248 262 REAL(wp), DIMENSION(jpi,jpj) :: zacs, zaos 249 263 !! 250 NAMELIST/namsbc_cpl/ sn_snd_temp , sn_snd_alb , sn_snd_thick, sn_snd_crt , sn_snd_co2 , & 264 NAMELIST/namsbc_cpl/ nn_cplmodel , ln_usecplmask, nn_cats_cpl , ln_scale_ice_flux, & 265 & sn_snd_temp , sn_snd_alb , sn_snd_thick, sn_snd_crt , sn_snd_co2 , & 251 266 & sn_snd_ttilyr, sn_snd_cond , sn_snd_mpnd , sn_snd_sstfrz, sn_snd_thick1, & 252 & sn_snd_ifrac , sn_snd_crtw , sn_snd_wlev , sn_rcv_hsig , sn_rcv_phioc ,&253 & sn_rcv_w10m , sn_rcv_taumod, sn_rcv_tau , sn_rcv_dqnsdt, sn_rcv_qsr ,&267 & sn_snd_ifrac , sn_snd_crtw , sn_snd_wlev , sn_rcv_hsig , sn_rcv_phioc , & 268 & sn_rcv_w10m , sn_rcv_taumod, sn_rcv_tau , sn_rcv_dqnsdt, sn_rcv_qsr , & 254 269 & sn_rcv_sdrfx , sn_rcv_sdrfy , sn_rcv_wper , sn_rcv_wnum , sn_rcv_tauwoc, & 255 & sn_rcv_wdrag , sn_rcv_qns , sn_rcv_emp , sn_rcv_rnf , sn_rcv_cal ,&256 & sn_rcv_iceflx, sn_rcv_co2 , nn_cplmodel , ln_usecplmask, sn_rcv_mslp ,&257 & sn_rcv_icb , sn_rcv_isf , sn_rcv_wfreq , sn_rcv_tauw, nn_cats_cpl ,&270 & sn_rcv_wdrag , sn_rcv_qns , sn_rcv_emp , sn_rcv_rnf , sn_rcv_cal , & 271 & sn_rcv_iceflx, sn_rcv_co2 , sn_rcv_mslp , & 272 & sn_rcv_icb , sn_rcv_isf , sn_rcv_wfreq, sn_rcv_tauw , & 258 273 & sn_rcv_ts_ice 259 260 274 !!--------------------------------------------------------------------- 261 275 ! … … 279 293 ENDIF 280 294 IF( lwp .AND. ln_cpl ) THEN ! control print 295 WRITE(numout,*)' nn_cplmodel = ', nn_cplmodel 296 WRITE(numout,*)' ln_usecplmask = ', ln_usecplmask 297 WRITE(numout,*)' ln_scale_ice_flux = ', ln_scale_ice_flux 298 WRITE(numout,*)' nn_cats_cpl = ', nn_cats_cpl 281 299 WRITE(numout,*)' received fields (mutiple ice categogies)' 282 300 WRITE(numout,*)' 10m wind module = ', TRIM(sn_rcv_w10m%cldes ), ' (', TRIM(sn_rcv_w10m%clcat ), ')' … … 327 345 WRITE(numout,*)' - orientation = ', sn_snd_crtw%clvor 328 346 WRITE(numout,*)' - mesh = ', sn_snd_crtw%clvgrd 329 WRITE(numout,*)' nn_cplmodel = ', nn_cplmodel330 WRITE(numout,*)' ln_usecplmask = ', ln_usecplmask331 WRITE(numout,*)' nn_cats_cpl = ', nn_cats_cpl332 347 ENDIF 333 348 … … 821 836 END SELECT 822 837 838 ! Initialise ice fractions from last coupling time to zero (needed by Met-Office) 839 a_i_last_couple(:,:,:) = 0._wp 840 823 841 ! ! ------------------------- ! 824 842 ! ! Ice Meltponds ! … … 1639 1657 REAL(wp), DIMENSION(jpi,jpj) :: zemp_tot, zemp_ice, zemp_oce, ztprecip, zsprecip , zevap_oce, zdevap_ice 1640 1658 REAL(wp), DIMENSION(jpi,jpj) :: zqns_tot, zqns_oce, zqsr_tot, zqsr_oce, zqprec_ice, zqemp_oce, zqemp_ice 1659 REAL(wp), DIMENSION(jpi,jpj) :: zevap_ice_total 1641 1660 REAL(wp), DIMENSION(jpi,jpj,jpl) :: zqns_ice, zqsr_ice, zdqns_ice, zqevap_ice, zevap_ice, zqtr_ice_top, ztsu 1642 1661 !!---------------------------------------------------------------------- … … 1659 1678 ztprecip(:,:) = frcv(jpr_rain)%z3(:,:,1) + zsprecip(:,:) ! May need to ensure positive here 1660 1679 zemp_tot(:,:) = frcv(jpr_tevp)%z3(:,:,1) - ztprecip(:,:) 1661 zemp_ice(:,:) = ( frcv(jpr_ievp)%z3(:,:,1) - frcv(jpr_snow)%z3(:,:,1) ) * picefr(:,:)1662 1680 CASE( 'oce and ice' ) ! received fields: jpr_sbpr, jpr_semp, jpr_oemp, jpr_ievp 1663 1681 zemp_tot(:,:) = ziceld(:,:) * frcv(jpr_oemp)%z3(:,:,1) + picefr(:,:) * frcv(jpr_sbpr)%z3(:,:,1) … … 1671 1689 1672 1690 #if defined key_si3 1691 1692 ! --- evaporation over ice (kg/m2/s) --- ! 1693 IF (ln_scale_ice_flux) THEN ! typically met-office requirements 1694 IF (sn_rcv_emp%clcat == 'yes') THEN 1695 WHERE( a_i(:,:,:) > 1.e-10) ; zevap_ice(:,:,:) = frcv(jpr_ievp)%z3(:,:,:) * a_i_last_couple(:,:,:) / a_i(:,:,:) 1696 ELSEWHERE ; zevap_ice(:,:,:) = 0._wp 1697 END WHERE 1698 zevap_ice_total(:,:) = SUM( zevap_ice(:,:,:), dim=3 ) 1699 ELSE 1700 WHERE( picefr(:,:) > 1.e-10) ; zevap_ice(:,:,1) = frcv(jpr_ievp)%z3(:,:,1) * SUM( a_i_last_couple, dim=3 ) / picefr(:,:) 1701 ELSEWHERE ; zevap_ice(:,:,1) = 0._wp 1702 END WHERE 1703 zevap_ice_total(:,:) = zevap_ice(:,:,1) 1704 ENDIF 1705 ELSE 1706 IF (sn_rcv_emp%clcat == 'yes') THEN 1707 zevap_ice(:,:,1:jpl) = frcv(jpr_ievp)%z3(:,:,1:jpl) 1708 zevap_ice_total(:,:) = SUM( zevap_ice(:,:,:), dim=3 ) 1709 ELSE 1710 zevap_ice(:,:,1) = frcv(jpr_ievp)%z3(:,:,1) 1711 zevap_ice_total(:,:) = zevap_ice(:,:,1) 1712 ENDIF 1713 ENDIF 1714 1715 IF ( TRIM( sn_rcv_emp%cldes ) == 'conservative' ) THEN 1716 ! For conservative case zemp_ice has not been defined yet. Do it now. 1717 zemp_ice(:,:) = zevap_ice_total(:,:) * picefr(:,:) - frcv(jpr_snow)%z3(:,:,1) * picefr(:,:) 1718 ENDIF 1719 1673 1720 ! zsnw = snow fraction over ice after wind blowing (=picefr if no blowing) 1674 1721 zsnw(:,:) = 0._wp ; CALL ice_thd_snwblow( ziceld, zsnw ) … … 1679 1726 1680 1727 ! --- evaporation over ocean (used later for qemp) --- ! 1681 zevap_oce(:,:) = frcv(jpr_tevp)%z3(:,:,1) - frcv(jpr_ievp)%z3(:,:,1) * picefr(:,:) 1682 1683 ! --- evaporation over ice (kg/m2/s) --- ! 1684 DO jl=1,jpl 1685 IF (sn_rcv_emp%clcat == 'yes') THEN ; zevap_ice(:,:,jl) = frcv(jpr_ievp)%z3(:,:,jl) 1686 ELSE ; zevap_ice(:,:,jl) = frcv(jpr_ievp)%z3(:,:,1 ) ; ENDIF 1687 ENDDO 1728 zevap_oce(:,:) = frcv(jpr_tevp)%z3(:,:,1) - zevap_ice_total(:,:) * picefr(:,:) 1688 1729 1689 1730 ! since the sensitivity of evap to temperature (devap/dT) is not prescribed by the atmosphere, we set it to 0 … … 1763 1804 !! IF( srcv(jpr_rnf)%laction ) CALL iom_put( 'runoffs' , rnf(:,:) * tmask(:,:,1) ) ! runoff 1764 1805 !! IF( srcv(jpr_isf)%laction ) CALL iom_put( 'iceshelf_cea', -fwfisf(:,:) * tmask(:,:,1) ) ! iceshelf 1765 IF( srcv(jpr_cal)%laction ) CALL iom_put( 'calving_cea' , frcv(jpr_cal)%z3(:,:,1) * tmask(:,:,1) ) ! calving1766 IF( srcv(jpr_icb)%laction ) CALL iom_put( 'iceberg_cea' , frcv(jpr_icb)%z3(:,:,1) * tmask(:,:,1) ) ! icebergs1767 CALL iom_put( 'snowpre' , sprecip(:,:) ) ! Snow1768 CALL iom_put( 'precip' , tprecip(:,:) ) ! total precipitation1769 IF ( iom_use('rain') )CALL iom_put( 'rain' , tprecip(:,:) - sprecip(:,:) ) ! liquid precipitation1770 IF ( iom_use('snow_ao_cea') )CALL iom_put( 'snow_ao_cea' , sprecip(:,:) * ( 1._wp - zsnw(:,:) ) ) ! Snow over ice-free ocean (cell average)1771 IF ( iom_use('snow_ai_cea') )CALL iom_put( 'snow_ai_cea' , sprecip(:,:) * zsnw(:,:) ) ! Snow over sea-ice (cell average)1772 IF ( iom_use('rain_ao_cea') )CALL iom_put( 'rain_ao_cea' , ( tprecip(:,:) - sprecip(:,:) ) * picefr(:,:) ) ! liquid precipitation over ocean (cell average)1773 IF ( iom_use('subl_ai_cea') ) CALL iom_put( 'subl_ai_cea' , frcv(jpr_ievp)%z3(:,:,1) * picefr(:,:) * tmask(:,:,1) )! Sublimation over sea-ice (cell average)1774 IF ( iom_use('evap_ao_cea') )CALL iom_put( 'evap_ao_cea' , ( frcv(jpr_tevp)%z3(:,:,1) &1775 & - frcv(jpr_ievp)%z3(:,:,1) * picefr(:,:) ) * tmask(:,:,1) )! ice-free oce evap (cell average)1806 IF( srcv(jpr_cal)%laction ) CALL iom_put( 'calving_cea' , frcv(jpr_cal)%z3(:,:,1) * tmask(:,:,1) ) ! calving 1807 IF( srcv(jpr_icb)%laction ) CALL iom_put( 'iceberg_cea' , frcv(jpr_icb)%z3(:,:,1) * tmask(:,:,1) ) ! icebergs 1808 IF( iom_use('snowpre') ) CALL iom_put( 'snowpre' , sprecip(:,:) ) ! Snow 1809 IF( iom_use('precip') ) CALL iom_put( 'precip' , tprecip(:,:) ) ! total precipitation 1810 IF( iom_use('rain') ) CALL iom_put( 'rain' , tprecip(:,:) - sprecip(:,:) ) ! liquid precipitation 1811 IF( iom_use('snow_ao_cea') ) CALL iom_put( 'snow_ao_cea' , sprecip(:,:) * ( 1._wp - zsnw(:,:) ) ) ! Snow over ice-free ocean (cell average) 1812 IF( iom_use('snow_ai_cea') ) CALL iom_put( 'snow_ai_cea' , sprecip(:,:) * zsnw(:,:) ) ! Snow over sea-ice (cell average) 1813 IF( iom_use('rain_ao_cea') ) CALL iom_put( 'rain_ao_cea' , ( tprecip(:,:) - sprecip(:,:) ) * picefr(:,:) ) ! liquid precipitation over ocean (cell average) 1814 IF( iom_use('subl_ai_cea') ) CALL iom_put( 'subl_ai_cea' , frcv(jpr_ievp)%z3(:,:,1) * picefr(:,:) * tmask(:,:,1) ) ! Sublimation over sea-ice (cell average) 1815 IF( iom_use('evap_ao_cea') ) CALL iom_put( 'evap_ao_cea' , ( frcv(jpr_tevp)%z3(:,:,1) & 1816 & - frcv(jpr_ievp)%z3(:,:,1) * picefr(:,:) ) * tmask(:,:,1) ) ! ice-free oce evap (cell average) 1776 1817 ! note: runoff output is done in sbcrnf (which includes icebergs too) and iceshelf output is done in sbcisf 1777 1818 ! … … 1781 1822 CASE( 'oce only' ) ! the required field is directly provided 1782 1823 zqns_tot(:,:) = frcv(jpr_qnsoce)%z3(:,:,1) 1824 ! For Met Office sea ice non-solar fluxes are already delt with by JULES so setting to zero 1825 ! here so the only flux is the ocean only one. 1826 zqns_ice(:,:,:) = 0._wp 1783 1827 CASE( 'conservative' ) ! the required fields are directly provided 1784 1828 zqns_tot(:,:) = frcv(jpr_qnsmix)%z3(:,:,1) … … 1811 1855 zqns_ice(:,:,jl) = frcv(jpr_qnsmix)%z3(:,:,jl) & 1812 1856 & + frcv(jpr_dqnsdt)%z3(:,:,jl) * ( pist(:,:,jl) - ( ( rt0 + psst(:,:) ) * ziceld(:,:) & 1813 & 1857 & + pist(:,:,jl) * picefr(:,:) ) ) 1814 1858 END DO 1815 1859 ELSE … … 1817 1861 zqns_ice(:,:,jl) = frcv(jpr_qnsmix)%z3(:,:, 1) & 1818 1862 & + frcv(jpr_dqnsdt)%z3(:,:, 1) * ( pist(:,:,jl) - ( ( rt0 + psst(:,:) ) * ziceld(:,:) & 1819 & 1863 & + pist(:,:,jl) * picefr(:,:) ) ) 1820 1864 END DO 1821 1865 ENDIF … … 1921 1965 CASE( 'oce only' ) 1922 1966 zqsr_tot(:,: ) = MAX( 0._wp , frcv(jpr_qsroce)%z3(:,:,1) ) 1967 ! For Met Office sea ice solar fluxes are already delt with by JULES so setting to zero 1968 ! here so the only flux is the ocean only one. 1969 zqsr_ice(:,:,:) = 0._wp 1923 1970 CASE( 'conservative' ) 1924 1971 zqsr_tot(:,: ) = frcv(jpr_qsrmix)%z3(:,:,1) … … 2007 2054 ENDDO 2008 2055 ENDIF 2056 CASE( 'none' ) 2057 zdqns_ice(:,:,:) = 0._wp 2009 2058 END SELECT 2010 2059 … … 2022 2071 ! ! ========================= ! 2023 2072 CASE ('coupled') 2024 IF( ln_mixcpl ) THEN 2025 DO jl=1,jpl 2026 qml_ice(:,:,jl) = qml_ice(:,:,jl) * xcplmask(:,:,0) + frcv(jpr_topm)%z3(:,:,jl) * zmsk(:,:) 2027 qcn_ice(:,:,jl) = qcn_ice(:,:,jl) * xcplmask(:,:,0) + frcv(jpr_botm)%z3(:,:,jl) * zmsk(:,:) 2028 ENDDO 2073 IF (ln_scale_ice_flux) THEN 2074 WHERE( a_i(:,:,:) > 1.e-10_wp ) 2075 qml_ice(:,:,:) = frcv(jpr_topm)%z3(:,:,:) * a_i_last_couple(:,:,:) / a_i(:,:,:) 2076 qcn_ice(:,:,:) = frcv(jpr_botm)%z3(:,:,:) * a_i_last_couple(:,:,:) / a_i(:,:,:) 2077 ELSEWHERE 2078 qml_ice(:,:,:) = 0.0_wp 2079 qcn_ice(:,:,:) = 0.0_wp 2080 END WHERE 2029 2081 ELSE 2030 2082 qml_ice(:,:,:) = frcv(jpr_topm)%z3(:,:,:) … … 2245 2297 ENDIF 2246 2298 2299 ! If this coupling was successful then save ice fraction for use between coupling points. 2300 ! This is needed for some calculations where the ice fraction at the last coupling point 2301 ! is needed. 2302 IF( info == OASIS_Sent .OR. info == OASIS_ToRest .OR. & 2303 & info == OASIS_SentOut .OR. info == OASIS_ToRestOut ) THEN 2304 IF ( sn_snd_thick%clcat == 'yes' ) THEN 2305 a_i_last_couple(:,:,1:jpl) = a_i(:,:,1:jpl) 2306 ENDIF 2307 ENDIF 2308 2247 2309 IF( ssnd(jps_fice1)%laction ) THEN 2248 2310 SELECT CASE( sn_snd_thick1%clcat ) … … 2314 2376 ztmp4(:,:,:) = 0.0 2315 2377 DO jl=1,jpl 2316 ztmp3(:,:,1) = ztmp3(:,:,1) + a_ip_ eff(:,:,jpl)2317 ztmp4(:,:,1) = ztmp4(:,:,1) + h_ip(:,:,jpl) 2378 ztmp3(:,:,1) = ztmp3(:,:,1) + a_ip_frac(:,:,jpl) 2379 ztmp4(:,:,1) = ztmp4(:,:,1) + h_ip(:,:,jpl) 2318 2380 ENDDO 2319 2381 CASE default ; CALL ctl_stop( 'sbc_cpl_snd: wrong definition of sn_snd_mpnd%clcat' )
Note: See TracChangeset
for help on using the changeset viewer.
