Skip to content
Snippets Groups Projects
Commit b4f557ab authored by RODIER Quentin's avatar RODIER Quentin
Browse files

Didier R. 24/11/2023: fix computation of diag theta_es

parent 9812bf52
No related branches found
No related tags found
No related merge requests found
Hide whitespace changes
Inline Side-by-side
src/MNH/write_lfifm1_for_diag.f90
+ 1
− 1
View file @ b4f557ab
...@@ -2454,7 +2454,7 @@ IF (LMOIST_ES .AND. (NRR>0)) THEN ...@@ -2454,7 +2454,7 @@ IF (LMOIST_ES .AND. (NRR>0)) THEN
(3.5*ALOG(XTHT(:,:,:)*( XPABST(:,:,:)/XP00 )**(XRD/XCPD) ) & (3.5*ALOG(XTHT(:,:,:)*( XPABST(:,:,:)/XP00 )**(XRD/XCPD) ) &
- ALOG( XPABST(:,:,:)*0.01*ZWORK31(:,:,:) / ( 0.622+ZWORK31(:,:,:) ) ) & - ALOG( XPABST(:,:,:)*0.01*ZWORK31(:,:,:) / ( 0.622+ZWORK31(:,:,:) ) ) &
-4.805 ) ) + 55. -4.805 ) ) + 55.
ZTHETAES(:,:,:)= XTHT(:,:,:) * EXP( (3376. / ZTHETAE(:,:,:) - 2.54) & ZTHETAES(:,:,:)= XTHT(:,:,:) * EXP( (3376. / ZTHETAES(:,:,:) - 2.54) &
*ZWORK31(:,:,:) *(1. +0.81 *ZWORK31(:,:,:)) ) *ZWORK31(:,:,:) *(1. +0.81 *ZWORK31(:,:,:)) )
TZFIELD = TFIELDMETADATA( & TZFIELD = TFIELDMETADATA( &
CMNHNAME = 'THETAES', & CMNHNAME = 'THETAES', &
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment

Documentation | Mentions légales | CGU | Accessibilité : non conforme