Ignore:
Timestamp:
2020-04-23T10:11:38+02:00 (5 months ago)
Author:
gsamson
Message:

move rhoa computation from sbcabl to sbcblk (to avoid restartability issue over sea-ice); correct prt_ctl calls in ablmod and comment a ctl_stop in iom to allow using 3D atm forcing with bulk mode

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2020/dev_r12563_ASINTER-06_ABL_improvement/src/OCE/SBC/sbcblk.F90

    r12588 r12799  
    658658 
    659659      IF( ln_abl ) THEN         !==  ABL formulation  ==!   multiplication by rho_air and turbulent fluxes computation done in ablstp 
    660          !! FL do we need this multiplication by tmask ... ??? 
    661660         DO_2D_11_11 
    662             zztmp = zU_zu(ji,jj) !* tmask(ji,jj,1) 
     661            zztmp = zU_zu(ji,jj) 
    663662            wndm(ji,jj)   = zztmp                   ! Store zU_zu in wndm to compute ustar2 in ablmod 
    664663            pcd_du(ji,jj) = zztmp * zcd_oce(ji,jj) 
    665664            psen(ji,jj)   = zztmp * zch_oce(ji,jj) 
    666665            pevp(ji,jj)   = zztmp * zce_oce(ji,jj) 
     666            rhoa(ji,jj)   = rho_air( ptair(ji,jj), phumi(ji,jj), pslp(ji,jj) ) 
    667667         END_2D 
    668668      ELSE                      !==  BLK formulation  ==!   turbulent fluxes computation 
Note: See TracChangeset for help on using the changeset viewer.