Skip to content
Snippets Groups Projects
Commit c51fc8fd authored by Gaelle DELAUTIER's avatar Gaelle DELAUTIER
Browse files

Gaelle 18/5/2018 : bug with khko

parent 6d4195d6
No related branches found
No related tags found
No related merge requests found
Hide whitespace changes
Inline Side-by-side
Showing
with 5 additions and 3 deletions
src/MNH/shallow_mf.f90
+ 5
− 3
View file @ c51fc8fd
...@@ -312,9 +312,11 @@ ENDIF ...@@ -312,9 +312,11 @@ ENDIF
! Thermodynamics functions ! Thermodynamics functions
ZFRAC_ICE(:,:) = 0. ZFRAC_ICE(:,:) = 0.
WHERE(PRM(:,:,2)+PRM(:,:,4) > 1.E-20) IF (SIZE(PRM,3).GE.4) THEN
ZFRAC_ICE(:,:) = PRM(:,:,4) / (PRM(:,:,2)+PRM(:,:,4)) WHERE(PRM(:,:,2)+PRM(:,:,4) > 1.E-20)
ENDWHERE ZFRAC_ICE(:,:) = PRM(:,:,4) / (PRM(:,:,2)+PRM(:,:,4))
ENDWHERE
ENDIF
CALL COMPUTE_FRAC_ICE(HFRAC_ICE,ZFRAC_ICE(:,:),PTHM(:,:)*PEXNM(:,:)) CALL COMPUTE_FRAC_ICE(HFRAC_ICE,ZFRAC_ICE(:,:),PTHM(:,:)*PEXNM(:,:))
! Conservative variables at t-dt ! Conservative variables at t-dt
......
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