Changeset 501


Ignore:
Timestamp:
Jan 18, 2012, 4:38:50 PM (13 years ago)
Author:
acolaitis
Message:

More robust gustiness speed computation. Sorry for changing it so often, this time it should be pretty final... This formulation make the 1D model match LES result very well on high vertical resolution. The offset in mixed layer temperature observed on normal resolution is definitely due to vertical resolution (mixing does not occur high enough in the PBL) and not to surface layer underestimating sensible heat flux.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/LMDZ.MARS/libf/phymars/vdifc.F

    r499 r501  
    320320        IF (callrichsl) THEN
    321321       zcdv(:)=zcdv_true(:)*sqrt(zu2(:)+
    322      &  (14.75*sqrt(hfmax(:)) - 10.*hfmax(:) + 1.7*hfmax(:)**2)**2)     ! subgrid gustiness added by quadratic interpolation with a coeff beta
     322     &  (log(1.+0.7*wstar(:) + 2.3*wstar(:)**2))**2)
     323!     &  (14.75*sqrt(hfmax(:)) - 10.*hfmax(:) + 1.7*hfmax(:)**2)**2)     ! subgrid gustiness added by quadratic interpolation with a coeff beta
    323324       zcdh(:)=zcdh_true(:)*sqrt(zu2(:)+
    324      &  (20.6*hfmax(:) - 22.8*hfmax(:)**2 + 9.*hfmax(:)**4)**2)     ! LES comparisons. This parameter is linked to the thermals settings)
     325     &  (log(1.+0.7*wstar(:) + 2.3*wstar(:)**2))**2)
     326!     &  (8.*hfmax(:) + 0.*hfmax(:)**2 + 0.*hfmax(:)**4)**2)     ! LES comparisons. This parameter is linked to the thermals settings)
    325327
    326328        ust(:)=sqrt(zcdv_true(:))*sqrt(zu2(:)+
    327      &  (14.75*sqrt(hfmax(:)) - 10.*hfmax(:) + 1.7*hfmax(:)**2)**2)
     329     &  (log(1.+0.7*wstar(:) + 2.3*wstar(:)**2))**2)
     330!     &  (14.75*sqrt(hfmax(:)) - 10.*hfmax(:) + 1.7*hfmax(:)**2)**2)
    328331        tst(:)=(ph(:,1)-ptsrf(:))*zcdh(:)/ust(:)
    329332!        ust(:)=sqrt(zcdv_true(:))*sqrt(zu2(:))
     
    332335        zcdv(:)=zcdv_true(:)*sqrt(zu2(:))     ! 1 / bulk aerodynamic momentum conductance
    333336        zcdh(:)=zcdh_true(:)*sqrt(zu2(:))     ! 1 / bulk aerodynamic heat conductance
     337        ust(:)=sqrt(zcdv_true(:))*sqrt(zu2(:))
     338        tst(:)=(ph(:,1)-ptsrf(:))*zcdh_true(:)/sqrt(zcdv_true(:))
    334339        ENDIF
    335340
    336341
    337342! Some usefull diagnostics for the new surface layer parametrization :
    338 
     343 
    339344!         call WRITEDIAGFI(ngridmx,'pcdv',
    340345!     &              'momentum drag','no units',
Note: See TracChangeset for help on using the changeset viewer.