From c05c8cdbfd3e93aad8abfc37701e5bab36a7f5f3 Mon Sep 17 00:00:00 2001 From: Quentin Rodier <quentin.rodier@meteo.fr> Date: Thu, 18 Jul 2024 15:11:23 +0200 Subject: [PATCH] Jean W. 18/07/2024: SURFEX, init of PSOLD to XUNDEF to identify init phase --- src/SURFEX/prep_hor_teb_field.F90 | 2 +- src/SURFEX/prep_teb_buffer.F90 | 2 +- src/SURFEX/prep_teb_extern.F90 | 4 ++-- src/SURFEX/prep_teb_grib.F90 | 2 +- src/SURFEX/prep_teb_unif.F90 | 2 +- 5 files changed, 6 insertions(+), 6 deletions(-) diff --git a/src/SURFEX/prep_hor_teb_field.F90 b/src/SURFEX/prep_hor_teb_field.F90 index b6e625203..d9c1a63e8 100644 --- a/src/SURFEX/prep_hor_teb_field.F90 +++ b/src/SURFEX/prep_hor_teb_field.F90 @@ -329,7 +329,7 @@ IF (YDCTL%LPART5) THEN CASE('PSOLD ') ALLOCATE(B%XPSOLD (SIZE(ZFIELDOUT,1))) ALLOCATE(ZPS(SIZE(ZFIELDOUT,1))) - ZPS = XP00 - ZRHOA * XG * XZS_LS + ZPS = XUNDEF ZFIELDOUT(:,1) = ZPS DEALLOCATE(ZPS) B%XPSOLD(:) = ZFIELDOUT(:,1) diff --git a/src/SURFEX/prep_teb_buffer.F90 b/src/SURFEX/prep_teb_buffer.F90 index 6b7dfa256..5565b0f7c 100644 --- a/src/SURFEX/prep_teb_buffer.F90 +++ b/src/SURFEX/prep_teb_buffer.F90 @@ -223,7 +223,7 @@ SELECT CASE(HSURF) ! CASE('PSOLD ') ALLOCATE(PFIELD(NNI,1)) - PFIELD = 101325.0 + PFIELD = XUNDEF ! CASE('VENTNIG') ALLOCATE(PFIELD(NNI,1)) diff --git a/src/SURFEX/prep_teb_extern.F90 b/src/SURFEX/prep_teb_extern.F90 index 25482e528..8c464625f 100644 --- a/src/SURFEX/prep_teb_extern.F90 +++ b/src/SURFEX/prep_teb_extern.F90 @@ -516,7 +516,7 @@ ELSE YRECFM=ADJUSTL(YRECFM) IF (HSURF=='PSOLD ') THEN IF (GOLD_NAME2 .OR. YBEM == 'DEF') THEN - PFIELD(:,1) = 101325.0 + PFIELD(:,1) = XUNDEF ELSE CALL READ_SURF(HFILETYPE,YRECFM,PFIELD(:,1),IRESP,HDIR='E') ENDIF @@ -608,7 +608,7 @@ ELSE ! CASE('PSOLD ') ALLOCATE(PFIELD(INI,1)) - PFIELD = 101325.0 + PFIELD = XUNDEF ! CASE('VENTNIG') ALLOCATE(PFIELD(INI,1)) diff --git a/src/SURFEX/prep_teb_grib.F90 b/src/SURFEX/prep_teb_grib.F90 index 638e63932..e704a91c7 100644 --- a/src/SURFEX/prep_teb_grib.F90 +++ b/src/SURFEX/prep_teb_grib.F90 @@ -227,7 +227,7 @@ SELECT CASE(HSURF) ! CASE('PSOLD') ALLOCATE(PFIELD(NNI,1)) - PFIELD = 101325.0 + PFIELD = XUNDEF ! CASE('VENTNIG') ALLOCATE(PFIELD(NNI,1)) diff --git a/src/SURFEX/prep_teb_unif.F90 b/src/SURFEX/prep_teb_unif.F90 index 2688c8c03..c81e0e404 100644 --- a/src/SURFEX/prep_teb_unif.F90 +++ b/src/SURFEX/prep_teb_unif.F90 @@ -153,7 +153,7 @@ SELECT CASE(HSURF) ! CASE('PSOLD') ALLOCATE(PFIELD(1,1)) - PFIELD = 101325.0 + PFIELD = XUNDEF CASE('VENTNIG') ALLOCATE(PFIELD(1,1)) -- GitLab