- Timestamp:
- 2017-04-13T12:07:16+02:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/r6232_HZG_WAVE-coupling/NEMOGCM/NEMO/OPA_SRC/SBC/sbccpl.F90
r7878 r7905 168 168 TYPE(FLD_C) :: sn_snd_ifrac, sn_snd_crtw, sn_snd_wlev 169 169 ! Received from waves 170 TYPE(FLD_C) :: sn_rcv_hsig,sn_rcv_phioc,sn_rcv_sdrf x,sn_rcv_sdrfy,sn_rcv_wper, &170 TYPE(FLD_C) :: sn_rcv_hsig,sn_rcv_phioc,sn_rcv_sdrft,sn_rcv_wper, & 171 171 sn_rcv_wfreq,sn_rcv_wnum,sn_rcv_tauoc,sn_rcv_wdrag 172 172 ! Other namelist parameters ! … … 244 244 & sn_rcv_w10m, sn_rcv_taumod, sn_rcv_tau , sn_rcv_dqnsdt, sn_rcv_qsr, & 245 245 & sn_snd_ifrac, sn_snd_crtw , sn_snd_wlev , sn_rcv_hsig , sn_rcv_phioc , & 246 & sn_rcv_sdrfx, sn_rcv_sdrfy, sn_rcv_wper , sn_rcv_wnum , sn_rcv_wfreq, & 247 & sn_rcv_tauoc, sn_rcv_wdrag, sn_rcv_qns , sn_rcv_emp , sn_rcv_rnf, & 248 & sn_rcv_cal , sn_rcv_iceflx, sn_rcv_co2 , sn_rcv_mslp , nn_cplmodel, & 249 & ln_usecplmask 246 & sn_rcv_sdrft, sn_rcv_wper , sn_rcv_wnum , sn_rcv_wfreq, sn_rcv_tauoc, & 247 & sn_rcv_wdrag, sn_rcv_qns , sn_rcv_emp , sn_rcv_rnf , sn_rcv_cal , & 248 & sn_rcv_iceflx, sn_rcv_co2 , sn_rcv_mslp , nn_cplmodel, ln_usecplmask 250 249 !!--------------------------------------------------------------------- 251 250 ! … … 290 289 WRITE(numout,*)' significant wave heigth = ', TRIM(sn_rcv_hsig%cldes ), ' (', TRIM(sn_rcv_hsig%clcat ), ')' 291 290 WRITE(numout,*)' wave to oce energy flux = ', TRIM(sn_rcv_phioc%cldes ), ' (', TRIM(sn_rcv_phioc%clcat ), ')' 292 WRITE(numout,*)' Surface Stokes drift grid u = ', TRIM(sn_rcv_sdrfx%cldes ), ' (', TRIM(sn_rcv_sdrfx%clcat ), ')' 293 WRITE(numout,*)' Surface Stokes drift grid v = ', TRIM(sn_rcv_sdrfy%cldes ), ' (', TRIM(sn_rcv_sdrfy%clcat ), ')' 291 WRITE(numout,*)' Surface Stokes drift u,v = ', TRIM(sn_rcv_sdrft%cldes ), ' (', TRIM(sn_rcv_sdrft%clcat ), ')' 294 292 WRITE(numout,*)' Mean wave period = ', TRIM(sn_rcv_wper%cldes ), ' (', TRIM(sn_rcv_wper%clcat ), ')' 295 293 WRITE(numout,*)' Mean wave number = ', TRIM(sn_rcv_wnum%cldes ), ' (', TRIM(sn_rcv_wnum%clcat ), ')' … … 549 547 ENDIF 550 548 srcv(jpr_sdrftx)%clname = 'O_Sdrfx' ! Stokes drift in the u direction 551 IF( TRIM(sn_rcv_sdrfx%cldes ) == 'coupled' ) THEN 549 srcv(jpr_sdrfty)%clname = 'O_Sdrfy' ! Stokes drift in the v direction 550 IF( TRIM(sn_rcv_sdrft%cldes ) == 'coupled' ) THEN 552 551 srcv(jpr_sdrftx)%laction = .TRUE. 553 cpl_sdrftx = .TRUE.554 ENDIF555 srcv(jpr_sdrfty)%clname = 'O_Sdrfy' ! Stokes drift in the v direction556 IF( TRIM(sn_rcv_sdrfy%cldes ) == 'coupled' ) THEN557 552 srcv(jpr_sdrfty)%laction = .TRUE. 558 cpl_sdrft y= .TRUE.553 cpl_sdrft = .TRUE. 559 554 ENDIF 560 555 srcv(jpr_wper)%clname = 'O_WPer' ! mean wave period … … 964 959 !! emp upward mass flux [evap. - precip. (- runoffs) (- calving)] (ocean only case) 965 960 !!---------------------------------------------------------------------- 966 USE zdf_oce, ONLY : ln_zdfqiao967 961 USE sbcflx , ONLY : ln_shelf_flx 968 962 … … 1171 1165 IF( ln_sdw ) THEN ! Stokes Drift correction activated 1172 1166 ! ! ========================= ! 1173 ! ! Stokes drift u!1167 ! ! Stokes drift u,v ! 1174 1168 ! ! ========================= ! 1175 IF( srcv(jpr_sdrftx)%laction ) ut0sd(:,:) = frcv(jpr_sdrftx)%z3(:,:,1) 1176 ! 1177 ! ! ========================= ! 1178 ! ! Stokes drift v ! 1179 ! ! ========================= ! 1180 IF( srcv(jpr_sdrfty)%laction ) vt0sd(:,:) = frcv(jpr_sdrfty)%z3(:,:,1) 1169 IF( srcv(jpr_sdrftx)%laction .AND. srcv(jpr_sdrfty)%laction ) THEN 1170 ut0sd(:,:) = frcv(jpr_sdrftx)%z3(:,:,1) 1171 vt0sd(:,:) = frcv(jpr_sdrfty)%z3(:,:,1) 1172 ENDIF 1181 1173 ! 1182 1174 ! ! ========================= ! … … 1201 1193 1202 1194 ! Calculate the 3D Stokes drift both in coupled and not fully uncoupled mode 1203 IF( srcv(jpr_sdrftx)%laction .OR. srcv(jpr_sdrfty)%laction.OR. srcv(jpr_wper)%laction &1204 .OR. srcv(jpr_hsig)%laction .OR. srcv(jpr_wfreq)%laction) &1195 IF( (srcv(jpr_sdrftx)%laction .AND. srcv(jpr_sdrfty)%laction) .OR. srcv(jpr_wper)%laction & 1196 .OR. srcv(jpr_hsig)%laction .OR. srcv(jpr_wfreq)%laction) & 1205 1197 CALL sbc_stokes() 1206 1198 ENDIF … … 1215 1207 IF( srcv(jpr_phioc)%laction .AND. ln_phioc ) THEN 1216 1208 rn_crban(:,:) = 29.0 * frcv(jpr_phioc)%z3(:,:,1) 1209 WHERE( rn_crban < 0.0 ) rn_crban = 0.0 1210 WHERE( rn_crban > 1000.0 ) rn_crban = 1000.0 1217 1211 ENDIF 1218 1212
Note: See TracChangeset
for help on using the changeset viewer.