Skip to content
Snippets Groups Projects
Commit 78d81a8c authored by Juan Escobar's avatar Juan Escobar
Browse files

Juan 5/12/2017:num_diff.f90 , Pb SegFault , correct IF(ONUMDIFTH/OZDIFFU) nesting

parent f746cef0
No related branches found
No related tags found
No related merge requests found
...@@ -214,6 +214,7 @@ END MODULE MODI_NUM_DIFF ...@@ -214,6 +214,7 @@ END MODULE MODI_NUM_DIFF
!! 05/07 (C.Lac) Separation between variables !! 05/07 (C.Lac) Separation between variables
!! 07/09 (C.Lac) Correction on budget calls !! 07/09 (C.Lac) Correction on budget calls
!! J.Escobar : 15/09/2015 : WENO5 & JPHEXT <> 1 !! J.Escobar : 15/09/2015 : WENO5 & JPHEXT <> 1
!! J.Escobar : 05/12/2017 : Pb SegFault , correct IF(ONUMDIFTH/OZDIFFU) nesting
!! !!
!------------------------------------------------------------------------------- !-------------------------------------------------------------------------------
! !
...@@ -347,20 +348,20 @@ IF (ONUMDIFTH) THEN ...@@ -347,20 +348,20 @@ IF (ONUMDIFTH) THEN
CALL NUM_DIFF_ALGO_Z(PRTHS, PTHM, IGRID, PRHODJ, & CALL NUM_DIFF_ALGO_Z(PRTHS, PTHM, IGRID, PRHODJ, &
PDK2TH, PDK4TH, PLSTHM, & PDK2TH, PDK4TH, PLSTHM, &
TZHALO2LIST%HALO2, TZHALO2LSLIST%HALO2) TZHALO2LIST%HALO2, TZHALO2LSLIST%HALO2)
!!$ ELSE
!!$ CALL NUM_DIFF_ALGO(PRTHS, PTHM, IGRID, PRHODJ, &
!!$ PDK2TH, PDK4TH, PLSTHM, &
!!$ TZHALO2LIST%HALO2, TZHALO2LSLIST%HALO2)
!!$ ENDIF
ELSE
IF (OZDIFFU) THEN ! call z-diffusion for potential temperature
CALL NUM_DIFF_ALGO_Z(PRTHS, PTHM, IGRID, PRHODJ, &
PDK2TH, PDK4TH, PLSTHM )
ELSE ELSE
CALL NUM_DIFF_ALGO(PRTHS, PTHM, IGRID, PRHODJ, & CALL NUM_DIFF_ALGO(PRTHS, PTHM, IGRID, PRHODJ, &
PDK2TH, PDK4TH, PLSTHM ) PDK2TH, PDK4TH, PLSTHM, &
TZHALO2LIST%HALO2, TZHALO2LSLIST%HALO2)
ENDIF ENDIF
ENDIF !!$ ELSE
!!$ IF (OZDIFFU) THEN ! call z-diffusion for potential temperature
!!$ CALL NUM_DIFF_ALGO_Z(PRTHS, PTHM, IGRID, PRHODJ, &
!!$ PDK2TH, PDK4TH, PLSTHM )
!!$ ELSE
!!$ CALL NUM_DIFF_ALGO(PRTHS, PTHM, IGRID, PRHODJ, &
!!$ PDK2TH, PDK4TH, PLSTHM )
!!$ ENDIF
!!$ ENDIF
! !
IF ( GTKEALLOC ) THEN IF ( GTKEALLOC ) THEN
!!$ IF(NHALO == 1) THEN !!$ IF(NHALO == 1) THEN
...@@ -383,20 +384,20 @@ IF (ONUMDIFTH) THEN ...@@ -383,20 +384,20 @@ IF (ONUMDIFTH) THEN
CALL NUM_DIFF_ALGO_Z(PRRS(:,:,:,1), PRM(:,:,:,1), IGRID, PRHODJ, & CALL NUM_DIFF_ALGO_Z(PRRS(:,:,:,1), PRM(:,:,:,1), IGRID, PRHODJ, &
PDK2TH, PDK4TH, & PDK2TH, PDK4TH, &
PLSRVM, TZHALO2LIST%HALO2, TZHALO2LSLIST%HALO2) PLSRVM, TZHALO2LIST%HALO2, TZHALO2LSLIST%HALO2)
!!$ ELSE
!!$ CALL NUM_DIFF_ALGO(PRRS(:,:,:,1), PRM(:,:,:,1), IGRID, PRHODJ, &
!!$ PDK2TH, PDK4TH, PLSRVM, &
!!$ TZHALO2LIST%HALO2, TZHALO2LSLIST%HALO2)
!!$ ENDIF
ELSE
IF (OZDIFFU) THEN ! call z-diffusion for wv mixing ratio
CALL NUM_DIFF_ALGO_Z(PRRS(:,:,:,1), PRM(:,:,:,1), IGRID, PRHODJ, &
PDK2TH, PDK4TH, PLSRVM)
ELSE ELSE
CALL NUM_DIFF_ALGO(PRRS(:,:,:,1), PRM(:,:,:,1), IGRID, PRHODJ, & CALL NUM_DIFF_ALGO(PRRS(:,:,:,1), PRM(:,:,:,1), IGRID, PRHODJ, &
PDK2TH, PDK4TH, PLSRVM) PDK2TH, PDK4TH, PLSRVM, &
TZHALO2LIST%HALO2, TZHALO2LSLIST%HALO2)
ENDIF ENDIF
ENDIF !!$ ELSE
!!$ IF (OZDIFFU) THEN ! call z-diffusion for wv mixing ratio
!!$ CALL NUM_DIFF_ALGO_Z(PRRS(:,:,:,1), PRM(:,:,:,1), IGRID, PRHODJ, &
!!$ PDK2TH, PDK4TH, PLSRVM)
!!$ ELSE
!!$ CALL NUM_DIFF_ALGO(PRRS(:,:,:,1), PRM(:,:,:,1), IGRID, PRHODJ, &
!!$ PDK2TH, PDK4TH, PLSRVM)
!!$ ENDIF
!!$ ENDIF
ENDIF ENDIF
! !
! In some situations, it makes sense to use the z-diffusion also for cloud water, but it is ! In some situations, it makes sense to use the z-diffusion also for cloud water, but it is
......
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