Changeset 501
- Timestamp:
- Jan 18, 2012, 4:38:50 PM (13 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/LMDZ.MARS/libf/phymars/vdifc.F
r499 r501 320 320 IF (callrichsl) THEN 321 321 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 323 324 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) 325 327 326 328 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) 328 331 tst(:)=(ph(:,1)-ptsrf(:))*zcdh(:)/ust(:) 329 332 ! ust(:)=sqrt(zcdv_true(:))*sqrt(zu2(:)) … … 332 335 zcdv(:)=zcdv_true(:)*sqrt(zu2(:)) ! 1 / bulk aerodynamic momentum conductance 333 336 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(:)) 334 339 ENDIF 335 340 336 341 337 342 ! Some usefull diagnostics for the new surface layer parametrization : 338 343 339 344 ! call WRITEDIAGFI(ngridmx,'pcdv', 340 345 ! & 'momentum drag','no units',
Note: See TracChangeset
for help on using the changeset viewer.