Skip to content
Snippets Groups Projects
Commit b320324e authored by ESCOBAR MUNOZ Juan's avatar ESCOBAR MUNOZ Juan Committed by WAUTELET Philippe
Browse files

Juan 28/03/2014: for pb with emissivity/aerosol reset X/MNH/SURF_TINY=1.0e-80 in real8 case

parent 8abd4fd1
No related branches found
No related tags found
No related merge requests found
...@@ -66,6 +66,7 @@ END MODULE MODI_INI_CST ...@@ -66,6 +66,7 @@ END MODULE MODI_INI_CST
!! V. Masson 05/10/98 add XRHOLI !! V. Masson 05/10/98 add XRHOLI
!! C. Mari 31/10/00 add NDAYSEC !! C. Mari 31/10/00 add NDAYSEC
!! V. Masson 01/03/03 add XCONDI !! V. Masson 01/03/03 add XCONDI
!! J. Escobar 28/03/2014 for pb with emissivity/aerosol reset XMNH_TINY=1.0e-80 in real8 case
!! !!
!------------------------------------------------------------------------------- !-------------------------------------------------------------------------------
! !
...@@ -156,22 +157,25 @@ XALPI = LOG(XESTT) + (XBETAI /XTT) + (XGAMI *LOG(XTT)) ...@@ -156,22 +157,25 @@ XALPI = LOG(XESTT) + (XBETAI /XTT) + (XGAMI *LOG(XTT))
! Some machine precision value depending of real4/8 use ! Some machine precision value depending of real4/8 use
! !
XMNH_TINY = TINY (XMNH_TINY )
XMNH_TINY_12 = SQRT (XMNH_TINY )
XMNH_EPSILON = EPSILON (XMNH_EPSILON ) XMNH_EPSILON = EPSILON (XMNH_EPSILON )
XMNH_HUGE = HUGE (XMNH_HUGE ) XMNH_HUGE = HUGE (XMNH_HUGE )
#ifdef MNH_MPI_DOUBLE_PRECISION #ifdef MNH_MPI_DOUBLE_PRECISION
XMNH_TINY = 1.0e-80
XEPS_DT = 1.0e-5 XEPS_DT = 1.0e-5
XRES_FLAT_CART = 1.0e-12 XRES_FLAT_CART = 1.0e-12
XRES_OTHER = 1.0e-9 XRES_OTHER = 1.0e-9
XRES_PREP = 1.0e-8 XRES_PREP = 1.0e-8
#else #else
XMNH_TINY = TINY (XMNH_TINY )
XEPS_DT = 1.0e-4 XEPS_DT = 1.0e-4
XRES_FLAT_CART = 1.0e-12 XRES_FLAT_CART = 1.0e-12
XRES_OTHER = 1.0e-7 XRES_OTHER = 1.0e-7
XRES_PREP = 1.0e-4 XRES_PREP = 1.0e-4
#endif #endif
XMNH_TINY_12 = SQRT (XMNH_TINY )
! !
......
...@@ -46,6 +46,7 @@ ...@@ -46,6 +46,7 @@
!! C. Mari 31/10/00 add NDAYSEC !! C. Mari 31/10/00 add NDAYSEC
!! V. Masson 01/03/03 add XCONDI !! V. Masson 01/03/03 add XCONDI
!! A. Voldoire 01/12/09 add XTTSI, XICEC, XTTS for ESM !! A. Voldoire 01/12/09 add XTTSI, XICEC, XTTS for ESM
!! J. Escobar 28/03/2014 for pb with emissivity/aerosol reset XSURF_TINY=1.0e-80 in real8 case
!! !!
!------------------------------------------------------------------------------- !-------------------------------------------------------------------------------
! !
...@@ -76,7 +77,15 @@ REAL(KIND=JPRB) :: ZHOOK_HANDLE ...@@ -76,7 +77,15 @@ REAL(KIND=JPRB) :: ZHOOK_HANDLE
IF (LHOOK) CALL DR_HOOK('INI_CSTS',0,ZHOOK_HANDLE) IF (LHOOK) CALL DR_HOOK('INI_CSTS',0,ZHOOK_HANDLE)
#ifdef MNH
#ifdef MNH_MPI_DOUBLE_PRECISION
XSURF_TINY = 1.0e-80
#else
XSURF_TINY = TINY (XSURF_TINY ) XSURF_TINY = TINY (XSURF_TINY )
#endif
#else
XSURF_TINY = 1.0e-80
#endif
XSURF_TINY_12 = SQRT (XSURF_TINY ) XSURF_TINY_12 = SQRT (XSURF_TINY )
XSURF_EPSILON = EPSILON (XSURF_EPSILON ) * 10.0 XSURF_EPSILON = EPSILON (XSURF_EPSILON ) * 10.0
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment