Changeset 12814 for NEMO/trunk/src/OCE/SBC
- Timestamp:
- 2020-04-24T19:51:20+02:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/trunk/src/OCE/SBC/sbcblk.F90
r12629 r12814 643 643 644 644 IF( ln_abl ) THEN !== ABL formulation ==! multiplication by rho_air and turbulent fluxes computation done in ablstp 645 !! FL do we need this multiplication by tmask ... ???646 645 DO_2D_11_11 647 zztmp = zU_zu(ji,jj) !* tmask(ji,jj,1)646 zztmp = zU_zu(ji,jj) 648 647 wndm(ji,jj) = zztmp ! Store zU_zu in wndm to compute ustar2 in ablmod 649 648 pcd_du(ji,jj) = zztmp * zcd_oce(ji,jj) 650 649 psen(ji,jj) = zztmp * zch_oce(ji,jj) 651 650 pevp(ji,jj) = zztmp * zce_oce(ji,jj) 651 rhoa(ji,jj) = rho_air( ptair(ji,jj), phumi(ji,jj), pslp(ji,jj) ) 652 652 END_2D 653 653 ELSE !== BLK formulation ==! turbulent fluxes computation … … 883 883 884 884 ! local scalars ( place there for vector optimisation purposes) 885 !IF (ln_abl) rhoa (:,:) = rho_air( ptair(:,:), phumi(:,:), pslp(:,:) ) !!GS: rhoa must be (re)computed here with ABL to avoid division by zero after (TBI)886 885 zcd_dui(:,:) = wndm_ice(:,:) * Cd_ice(:,:) 887 886
Note: See TracChangeset
for help on using the changeset viewer.