Changeset 13376
- Timestamp:
- 2020-08-05T10:43:58+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
r13353 r13376 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 INTEGER, PARAMETER :: jpr_qtr = 58 ! Transmitted solar 119 120 INTEGER, PARAMETER :: jprcv = 58 ! total number of fields received 118 INTEGER, PARAMETER :: jpr_grnm = 58 ! Greenland ice mass 119 INTEGER, PARAMETER :: jpr_antm = 59 ! Antarctic ice mass 120 INTEGER, PARAMETER :: jpr_rnf_1d = 60 ! 1D river runoff 121 INTEGER, PARAMETER :: jpr_qtr = 61 ! Transmitted solar 122 123 INTEGER, PARAMETER :: jprcv = 61 ! total number of fields received 121 124 122 125 INTEGER, PARAMETER :: jps_fice = 1 ! ice fraction sent to the atmosphere … … 175 178 TYPE(FLD_C) :: sn_rcv_w10m, sn_rcv_taumod, sn_rcv_tau, sn_rcv_tauw, sn_rcv_dqnsdt, sn_rcv_qsr, & 176 179 & sn_rcv_qns , sn_rcv_emp , sn_rcv_rnf, sn_rcv_ts_ice 177 TYPE(FLD_C) :: sn_rcv_cal, sn_rcv_iceflx, sn_rcv_co2, sn_rcv_mslp, sn_rcv_icb, sn_rcv_isf , sn_rcv_qtr180 TYPE(FLD_C) :: sn_rcv_cal, sn_rcv_iceflx, sn_rcv_co2, sn_rcv_mslp, sn_rcv_icb, sn_rcv_isf 178 181 ! Send to waves 179 182 TYPE(FLD_C) :: sn_snd_ifrac, sn_snd_crtw, sn_snd_wlev … … 181 184 TYPE(FLD_C) :: sn_rcv_hsig, sn_rcv_phioc, sn_rcv_sdrfx, sn_rcv_sdrfy, sn_rcv_wper, sn_rcv_wnum, sn_rcv_tauwoc, & 182 185 sn_rcv_wdrag, sn_rcv_wfreq 186 ! Transmitted solar 187 TYPE(FLD_C) :: sn_rcv_qtr 183 188 ! ! Other namelist parameters 184 189 INTEGER :: nn_cplmodel ! Maximum number of models to/from which NEMO is potentialy sending/receiving data … … 254 259 & sn_rcv_w10m , sn_rcv_taumod, sn_rcv_tau , sn_rcv_dqnsdt, sn_rcv_qsr , & 255 260 & sn_rcv_sdrfx , sn_rcv_sdrfy , sn_rcv_wper , sn_rcv_wnum , sn_rcv_tauwoc, & 261 & sn_rcv_qtr, & 256 262 & sn_rcv_wdrag , sn_rcv_qns , sn_rcv_emp , sn_rcv_rnf , sn_rcv_cal , & 257 263 & sn_rcv_iceflx, sn_rcv_co2 , nn_cplmodel , ln_usecplmask, sn_rcv_mslp , & 258 264 & sn_rcv_icb , sn_rcv_isf , sn_rcv_wfreq , sn_rcv_tauw, nn_cats_cpl , & 259 & sn_rcv_ts_ice , sn_rcv_qtr265 & sn_rcv_ts_ice 260 266 261 267 !!--------------------------------------------------------------------- … … 2045 2051 CASE ('coupled') 2046 2052 IF (ln_scale_ice_fluxes) THEN 2047 WHERE( a_i(:,:,:) > 0.0_wp ) qtr_ice_top(:,:,:) = frcv(jpr_qtr)%z3(:,:,:) * a_i_last_couple(:,:,:) / a_i(:,:,:)2048 WHERE( a_i(:,:,:) <= 0.0_wp ) qtr_ice_top(:,:,:) = 0.0_wp2053 WHERE( a_i(:,:,:) > 0.0_wp ) zqtr_ice_top(:,:,:) = frcv(jpr_qtr)%z3(:,:,:) * a_i_last_couple(:,:,:) / a_i(:,:,:) 2054 WHERE( a_i(:,:,:) <= 0.0_wp ) zqtr_ice_top(:,:,:) = 0.0_wp 2049 2055 ELSE 2050 qtr_ice_top(:,:,:) = frcv(jpr_qtr)%z3(:,:,:)2056 zqtr_ice_top(:,:,:) = frcv(jpr_qtr)%z3(:,:,:) 2051 2057 ENDIF 2052 2058 2053 2059 ! if we are not getting this data from the coupler then assume zero (fully opaque ice) 2054 2060 CASE ('none') 2055 qtr_ice_top(:,:,:) = 0._wp2061 zqtr_ice_top(:,:,:) = 0._wp 2056 2062 END SELECT 2063 ! 2064 ENDIF 2065 ! 2066 IF( ln_mixcpl ) THEN 2067 DO jl=1,jpl 2068 qtr_ice_top(:,:,jl) = qtr_ice_top(:,:,jl) * xcplmask(:,:,0) + zqtr_ice_top(:,:,jl) * zmsk(:,:) 2069 ENDDO 2070 ELSE 2071 qtr_ice_top(:,:,:) = zqtr_ice_top(:,:,:) 2072 ENDIF 2073 ! ! ================== ! 2074 ! ! ice skin temp. ! 2075 ! ! ================== ! 2076 ! needed by Met Office 2077 IF( srcv(jpr_ts_ice)%laction ) THEN 2078 WHERE ( frcv(jpr_ts_ice)%z3(:,:,:) > 0.0 ) ; ztsu(:,:,:) = 0.0 + rt0 2079 ELSEWHERE( frcv(jpr_ts_ice)%z3(:,:,:) < -60. ) ; ztsu(:,:,:) = -60. + rt0 2080 ELSEWHERE ; ztsu(:,:,:) = frcv(jpr_ts_ice)%z3(:,:,:) + rt0 2081 END WHERE 2082 ! 2083 IF( ln_mixcpl ) THEN 2084 DO jl=1,jpl 2085 pist(:,:,jl) = pist(:,:,jl) * xcplmask(:,:,0) + ztsu(:,:,jl) * zmsk(:,:) 2086 ENDDO 2087 ELSE 2088 pist(:,:,:) = ztsu(:,:,:) 2089 ENDIF 2057 2090 ! 2058 2091 ENDIF
Note: See TracChangeset
for help on using the changeset viewer.