Jean W. 18/07/2024: SURFEX, init of PSOLD to XUNDEF to identify init phase

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)

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))

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))

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))

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))