v4.4.8

Author: SiqiLiOcean

src/bcond_gcn.F

@@ -225,13 +225,21 @@ SUBROUTINE BCOND_GCN(IDX,K_RK)
          K300_CD = KBM1
          K300_ND = KBM1
          DO K=1,KBM1
-            IF (H(CDX)*(-ZZ(CDX,K)) > 300.0_SP) THEN
+!@---> Qichun Xu, 300to200@20250618
+!@     Changed the top layer thickness from 300 m to 200 m (for Scotain Shelf)
+!@            IF (H(CDX)*(-ZZ(CDX,K)) > 300.0_SP) THEN  
+            IF (H(CDX)*(-ZZ(CDX,K)) > 200.0_SP) THEN  
+!@<---
                 K300_CD = K - 1
                 EXIT
             ENDIF
          ENDDO
          DO K=1,KBM1
-            IF (H(NDX)*(-ZZ(NDX,K)) > 300.0_SP) THEN
+!@---> Qichun Xu, 300to200@20250618
+!@     Changed the top layer thickness from 300 m to 200 m (for Scotain Shelf)
+!@            IF (H(NDX)*(-ZZ(NDX,K)) > 300.0_SP) THEN
+            IF (H(NDX)*(-ZZ(NDX,K)) > 200.0_SP) THEN
+!@<---
                 K300_ND = K - 1
                 EXIT
             ENDIF

src/bcond_gcy.F

@@ -182,13 +182,21 @@ SUBROUTINE BCOND_GCY(IDX,K_RK)
          K300_CD = KBM1
          K300_ND = KBM1
          DO K=1,KBM1
-            IF (H(CDX)*(-ZZ(CDX,K)) > 300.0_SP) THEN
+!@---> Qichun Xu, 300to200@20250618
+!@     Changed the top layer thickness from 300 m to 200 m (for Scotain Shelf)
+!@            IF (H(CDX)*(-ZZ(CDX,K)) > 300.0_SP) THEN
+            IF (H(CDX)*(-ZZ(CDX,K)) > 200.0_SP) THEN
+!@<---
                 K300_CD = K - 1
                 EXIT
             ENDIF
          ENDDO
          DO K=1,KBM1
-            IF (H(NDX)*(-ZZ(NDX,K)) > 300.0_SP) THEN
+!@---> Qichun Xu, 300to200@20250618
+!@     Changed the top layer thickness from 300 m to 200 m (for Scotain Shelf)
+!@            IF (H(NDX)*(-ZZ(NDX,K)) > 300.0_SP) THEN
+            IF (H(NDX)*(-ZZ(NDX,K)) > 200.0_SP) THEN
+!@<---
                 K300_ND = K - 1
                 EXIT
             ENDIF

src/fvcom.F

@@ -27,7 +27,7 @@
 !/===========================================================================/
 
 !==============================================================================!
-!  VERSION 4.4.7
+!  VERSION 4.4.8
 !==============================================================================!
 
 PROGRAM FVCOM

src/mod_startup.F

@@ -729,10 +729,19 @@ SUBROUTINE READ_AIR_PRESS_ELE
 
     ! LOAD Air pressure induced elevation
     VAR => FIND_VAR(NC_START,'el_press',FOUND)
-    IF(.not. FOUND) CALL FATAL_ERROR&
-         &("COULD NOT FIND VARIABLE 'el_press' IN THE HOTSTART FILE OBJECT")
-    CALL NC_CONNECT_AVAR(VAR, EL_AIR)
-    CALL NC_READ_VAR(VAR,STKCNT)
+!@---> Lianyuan Zheng, el_press@20250618
+!@      To avoid crash when el_press is missing in the restart file
+!@    IF(.not. FOUND) CALL FATAL_ERROR&
+!@         &("COULD NOT FIND VARIABLE 'el_press' IN THE HOTSTART FILE OBJECT")
+!@    CALL NC_CONNECT_AVAR(VAR, EL_AIR)
+!@    CALL NC_READ_VAR(VAR,STKCNT)
+    IF(.not. FOUND) THEN
+      EL_AIR = ZERO
+    ELSE
+      CALL NC_CONNECT_AVAR(VAR, EL_AIR)
+      CALL NC_READ_VAR(VAR,STKCNT)
+    END IF
+!@<---
 
     IF(DBG_SET(DBG_SBR)) write(ipt,*) "End: READ_AIR_PRESS_ELE" 
 

src/mod_utils.F

@@ -141,7 +141,7 @@ SUBROUTINE INITIALIZE_CONTROL(NAME)
 !  FVCOM VERSION                                                               !
 !==============================================================================!
 
-   FVCOM_VERSION     = 'FVCOM_4.4.7'
+   FVCOM_VERSION     = 'FVCOM_4.4.8'
    FVCOM_WEBSITE     = 'http://fvcom.smast.umassd.edu, https://github.com/FVCOM-GitHub/FVCOM'
    INSTITUTION       = 'School for Marine Science and Technology'
 
@@ -2324,7 +2324,7 @@ SUBROUTINE WRITE_BANNER(PAR,NP,ID)
 WRITE(IPT,*)'         |  ___)  | |   | || |      | |   | || ||_|| |            '
 WRITE(IPT,*)'         | |       \ \ / / | |_____ | |___| || |   | |            '
 WRITE(IPT,*)'         |_|        \___/   \______) \_____/ |_|   |_|            '
-WRITE(IPT,*)' -- Version 4.4.7 Release'
+WRITE(IPT,*)' -- Version 4.4.8 Release'
 WRITE(IPT,*)'!================================================================!'
 WRITE(IPT,*)'!                                                                !'
 WRITE(IPT,*)'!========DOMAIN DECOMPOSITION USING: METIS 4.0.1 ================!'