Changeset 1695 for trunk/NEMO/OPA_SRC/SBC/sbccpl.F90
- Timestamp:
- 2009-10-30T17:42:38+01:00 (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/OPA_SRC/SBC/sbccpl.F90
r1694 r1695 545 545 !! 546 546 !! ** Action : update utau, vtau ocean stress at U,V grid 547 !! taum, wndm wind stres and wind speed module at T-point 547 548 !! qns , qsr non solar and solar ocean heat fluxes ('ocean only case) 548 549 !! emp = emps evap. - precip. (- runoffs) (- calving) ('ocean only case) … … 557 558 REAL(wp) :: zcumulneg, zcumulpos ! temporary scalars 558 559 REAL(wp) :: zcoef ! temporary scalar 560 REAL(wp) :: zrhoa = 1.22 ! Air density kg/m3 561 REAL(wp) :: zcdrag = 1.5e-3 ! drag coefficient 562 REAL(wp) :: zzx, zzy ! temporary variables 559 563 REAL(wp), DIMENSION(jpi,jpj) :: ztx, zty ! 2D workspace 560 564 !!---------------------------------------------------------------------- … … 625 629 ! u(v)tau will be modified by ice model -> need to be reset before each call of the ice/fsbc 626 630 IF( MOD( kt-1, k_fsbc ) == 0 ) THEN 631 ! 627 632 utau(:,:) = frcv(:,:,jpr_otx1) 628 633 vtau(:,:) = frcv(:,:,jpr_oty1) 629 IF( .NOT. srcv(jpr_w10m)%laction ) CALL sbc_tau2wnd 634 ! 635 ! wind stress module at T-point 636 !CDIR NOVERRCHK 637 DO jj = 2, jpjm1 638 !CDIR NOVERRCHK 639 DO ji = fs_2, fs_jpim1 ! vect. opt. 640 zzx = utau(ji-1,jj ) + utau(ji,jj) 641 zzy = vtau(ji ,jj-1) + vtau(ji,jj) 642 taum(ji,jj) = 0.5 * SQRT( zzx * zzx + zzy * zzy ) 643 END DO 644 END DO 645 CALL lbc_lnk( taum(:,:), 'T', 1. ) 646 ! 647 ! wind speed module at T-point 648 IF( .NOT. srcv(jpr_w10m)%laction ) THEN 649 zcoef = 1. / ( zrhoa * zcdrag ) 650 !CDIR NOVERRCHK 651 DO jj = 1, jpj 652 !CDIR NOVERRCHK 653 DO ji = 1, jpi 654 wndm(ji,jj) = SQRT( taum(ji,jj) * zcoef ) 655 END DO 656 END DO 657 ENDIF 658 ! 630 659 ENDIF 631 660 ! ! ========================= ! … … 678 707 ! ! 10 m wind speed 679 708 IF( srcv(jpr_w10m)%laction ) wndm(:,:) = frcv(:,:,jpr_w10m) 680 ! it not, we call sbc_tau2wnd in sbc_cpl_rcv (or later, after the ice???)681 709 ! 682 710 #if defined key_cpl_carbon_cycle
Note: See TracChangeset
for help on using the changeset viewer.