- Timestamp:
- 2020-07-27T17:51:59+02:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/UKMO/NEMO_4.0.1_penetrating_solar/src/OCE/SBC/sbccpl.F90
r11715 r13347 116 116 INTEGER, PARAMETER :: jpr_tauwy = 56 ! y component of the ocean stress from waves 117 117 INTEGER, PARAMETER :: jpr_ts_ice = 57 ! Sea ice surface temp 118 119 INTEGER, PARAMETER :: jprcv = 57 ! total number of fields received 118 INTEGER, PARAMETER :: jpr_qtr = 58 ! Transmitted solar 119 120 INTEGER, PARAMETER :: jprcv = 58 ! total number of fields received 120 121 121 122 INTEGER, PARAMETER :: jps_fice = 1 ! ice fraction sent to the atmosphere … … 174 175 TYPE(FLD_C) :: sn_rcv_w10m, sn_rcv_taumod, sn_rcv_tau, sn_rcv_tauw, sn_rcv_dqnsdt, sn_rcv_qsr, & 175 176 & sn_rcv_qns , sn_rcv_emp , sn_rcv_rnf, sn_rcv_ts_ice 176 TYPE(FLD_C) :: sn_rcv_cal, sn_rcv_iceflx, sn_rcv_co2, sn_rcv_mslp, sn_rcv_icb, sn_rcv_isf 177 TYPE(FLD_C) :: sn_rcv_cal, sn_rcv_iceflx, sn_rcv_co2, sn_rcv_mslp, sn_rcv_icb, sn_rcv_isf, sn_rcv_qtr 177 178 ! Send to waves 178 179 TYPE(FLD_C) :: sn_snd_ifrac, sn_snd_crtw, sn_snd_wlev … … 256 257 & sn_rcv_iceflx, sn_rcv_co2 , nn_cplmodel , ln_usecplmask, sn_rcv_mslp , & 257 258 & sn_rcv_icb , sn_rcv_isf , sn_rcv_wfreq , sn_rcv_tauw, nn_cats_cpl , & 258 & sn_rcv_ts_ice 259 & sn_rcv_ts_ice, sn_rcv_qtr 259 260 260 261 !!--------------------------------------------------------------------- … … 295 296 WRITE(numout,*)' ice shelf = ', TRIM(sn_rcv_isf%cldes ), ' (', TRIM(sn_rcv_isf%clcat ), ')' 296 297 WRITE(numout,*)' sea ice heat fluxes = ', TRIM(sn_rcv_iceflx%cldes), ' (', TRIM(sn_rcv_iceflx%clcat), ')' 298 WRITE(numout,*)' transmitted solar = ', TRIM(sn_rcv_qtr%cldes ), ' (', TRIM(sn_rcv_qtr%clcat ), ')' 297 299 WRITE(numout,*)' atm co2 = ', TRIM(sn_rcv_co2%cldes ), ' (', TRIM(sn_rcv_co2%clcat ), ')' 298 300 WRITE(numout,*)' significant wave heigth = ', TRIM(sn_rcv_hsig%cldes ), ' (', TRIM(sn_rcv_hsig%clcat ), ')' … … 566 568 ENDIF 567 569 ! ! ------------------------- ! 570 ! ! transmitted solar ! 571 ! ! ------------------------- ! 572 srcv(jpr_topm )%clname = 'OQtr' 573 IF( TRIM(sn_rcv_qtr%cldes) == 'coupled' ) THEN 574 IF ( TRIM( sn_rcv_qtr%clcat ) == 'yes' ) THEN 575 srcv(jpr_qtr)%nct = nn_cats_cpl 576 ELSE 577 CALL ctl_stop( 'sbc_cpl_init: sn_rcv_qtr%clcat should always be set to yes currently' ) 578 ENDIF 579 srcv(jpr_qtr)%laction = .TRUE. 580 ENDIF 581 582 ! ! ------------------------- ! 568 583 ! ! ice skin temperature ! 569 584 ! ! ------------------------- ! … … 2025 2040 ELSEIF( ln_cndflx .AND. .NOT.ln_cndemulate ) THEN !== conduction flux as surface forcing ==! 2026 2041 ! 2027 ! ! ===> here we must receive the qtr_ice_top array from the coupler 2028 ! for now just assume zero (fully opaque ice) 2029 qtr_ice_top(:,:,:) = 0._wp 2042 SELECT CASE( TRIM( sn_rcv_qtr%cldes ) ) 2043 ! 2044 ! ! ===> here we receive the qtr_ice_top array from the coupler 2045 CASE ('coupled') 2046 qtr_ice_top(:,:,:) = frcv(jpr_qtr)%z3(:,:,:) 2047 2048 ! if we are not getting this data from the coupler then assume zero (fully opaque ice) 2049 CASE ('none') 2050 qtr_ice_top(:,:,:) = 0._wp 2051 END SELECT 2030 2052 ! 2031 2053 ENDIF
Note: See TracChangeset
for help on using the changeset viewer.