New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 12814 for NEMO/trunk/src/OCE – NEMO

Changeset 12814 for NEMO/trunk/src/OCE


Ignore:
Timestamp:
2020-04-24T19:51:20+02:00 (4 years ago)
Author:
gsamson
Message:

ABL restartability bugfixes (ticket #2451); no SETTE change compared to r12794

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/trunk/src/OCE/SBC/sbcblk.F90

    r12629 r12814  
    643643 
    644644      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 ... ??? 
    646645         DO_2D_11_11 
    647             zztmp = zU_zu(ji,jj) !* tmask(ji,jj,1) 
     646            zztmp = zU_zu(ji,jj) 
    648647            wndm(ji,jj)   = zztmp                   ! Store zU_zu in wndm to compute ustar2 in ablmod 
    649648            pcd_du(ji,jj) = zztmp * zcd_oce(ji,jj) 
    650649            psen(ji,jj)   = zztmp * zch_oce(ji,jj) 
    651650            pevp(ji,jj)   = zztmp * zce_oce(ji,jj) 
     651            rhoa(ji,jj)   = rho_air( ptair(ji,jj), phumi(ji,jj), pslp(ji,jj) ) 
    652652         END_2D 
    653653      ELSE                      !==  BLK formulation  ==!   turbulent fluxes computation 
     
    883883 
    884884      ! 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) 
    886885      zcd_dui(:,:) = wndm_ice(:,:) * Cd_ice(:,:) 
    887886 
Note: See TracChangeset for help on using the changeset viewer.