Changeset 5309 for LMDZ6/trunk/libf/phylmd/inlandsis

Timestamp:

Nov 1, 2024, 12:39:44 PM (6 weeks ago)

Author:

abarral

Message:

Turn YOEGWD.h sisvat_weq into module

File:

• LMDZ6/trunk/libf/phylmd/inlandsis/sisvat_weq_mod.f90 (moved) (moved from LMDZ6/trunk/libf/phylmd/inlandsis/sisvat_weq.f90) (1 diff)

LMDZ6/trunk/libf/phylmd/inlandsis/sisvat_weq_mod.f90

@@ +1 @@
+MODULE sisvat_weq_mod
+  IMPLICIT NONE; PRIVATE
+  PUBLIC sisvat_weq
+
+  LOGICAL :: logWEq
+
+CONTAINS
+
+  SUBROUTINE SISVAT_wEq(labWEq, istart)
+
+    ! +------------------------------------------------------------------------+
+    ! | MAR          SISVAT_wEq                                22-09-2001  MAR |
+    ! |   SubRoutine SISVAT_wEq computes the Snow/Ice  Water  Equivalent       |
+    ! |                                                                        |
+    ! |                                                                        |
+    ! |   Preprocessing  Option: SISVAT IO (not always a standard preprocess.) |
+    ! |   ^^^^^^^^^^^^^^^^^^^^^  ^^^^^^^^^                                     |
+    ! |   FILE                 |      CONTENT                                  |
+    ! |   ~~~~~~~~~~~~~~~~~~~~~+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ |
+    ! | # SISVAT_wEq.ve        | #ve: OUTPUT/Verification: Snow/Ice Water Eqv. |
+    ! |                        |      unit 45, SubRoutine  SISVAT_wEq **ONLY** |
+    ! +------------------------------------------------------------------------+
+
+    ! +--Global Variables
+    ! +  ================
+    USE VARphy
+    USE VAR_SV
+    USE VARxSV
+
+    IMPLICIT NONE
+
+    CHARACTER(LEN = 6) :: labWEq
+    INTEGER :: istart
+
+    ! +--Local  Variables
+    ! +  ================
+
+    INTEGER :: ikl, isn
+    REAL :: SnoWEQ, IceWEQ
 
 
-subroutine SISVAT_wEq( labWEq ,istart)
+    ! +--Switch Initialization
+    ! +  =====================
 
-  ! +------------------------------------------------------------------------+
-  ! | MAR          SISVAT_wEq                                22-09-2001  MAR |
-  ! |   SubRoutine SISVAT_wEq computes the Snow/Ice  Water  Equivalent       |
-  ! |                                                                        |
-  ! |                                                                        |
-  ! |   Preprocessing  Option: SISVAT IO (not always a standard preprocess.) |
-  ! |   ^^^^^^^^^^^^^^^^^^^^^  ^^^^^^^^^                                     |
-  ! |   FILE                 |      CONTENT                                  |
-  ! |   ~~~~~~~~~~~~~~~~~~~~~+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ |
-  ! | # SISVAT_wEq.ve        | #ve: OUTPUT/Verification: Snow/Ice Water Eqv. |
-  ! |                        |      unit 45, SubRoutine  SISVAT_wEq **ONLY** |
-  ! +------------------------------------------------------------------------+
+    IF (.NOT.logWEq)                                              THEN
+      logWEq = .TRUE.
+      open(unit = 45, status = 'unknown', file = 'SISVAT_wEq.ve')
+      rewind    45
+    END IF
 
 
+    ! +--Snow Water Equivalent
+    ! +  =====================
+
+    ikl = 1
+    IF          (isnoSV(ikl)>iiceSV(ikl))                      THEN
+      SnoWEQ = 0.
+      DO isn = iiceSV(ikl) + 1, isnoSV(ikl)
+        SnoWEQ = SnoWEQ + ro__SV(ikl, isn) * dzsnSV(ikl, isn)
+      END DO
+
+    END IF
 
 
-  ! +--Global Variables
-  ! +  ================
+    ! +--Ice  Water Equivalent
+    ! +  =====================
+    IF        (iiceSV(1)>0)                                    THEN
+      IceWEQ = 0.
+      DO isn = 1, iiceSV(ikl)
+        IceWEQ = IceWEQ + ro__SV(ikl, isn) * dzsnSV(ikl, isn)
+      END DO
 
-  use VARphy
-  use VAR_SV
-  use VARxSV
+    END IF
 
-  IMPLICIT NONE
+  END SUBROUTINE sisvat_weq
 
-
-
-  character(len=6) :: labWEq
-  integer :: istart
-
-  logical :: logWEq
-  common/SISVAT_wEq_L/logWEq
-
-
-  ! +--Local  Variables
-  ! +  ================
-
-  integer :: ikl   ,isn
-  real :: SnoWEQ,IceWEQ
-
-
-  ! +--Switch Initialization
-  ! +  =====================
-
-  IF (.NOT.logWEq)                                              THEN
-           logWEq = .true.
-           open(unit=45,status='unknown',file='SISVAT_wEq.ve')
-           rewind    45
-  END IF
-
-
-  ! +--Snow Water Equivalent
-  ! +  =====================
-
-       ikl   = 1
-  IF          (isnoSV(ikl).gt.iiceSV(ikl))                      THEN
-
-      SnoWEQ = 0.
-    DO isn   = iiceSV(ikl)+1 ,isnoSV(ikl)
-      SnoWEQ = SnoWEQ       + ro__SV(ikl,isn) * dzsnSV(ikl,isn)
-    END DO
-
-  END IF
-
-
-  ! +--Ice  Water Equivalent
-  ! +  =====================
-
-  IF        (iiceSV(1).gt.0)                                    THEN
-
-      IceWEQ = 0.
-    DO isn   =             1 ,iiceSV(ikl)
-      IceWEQ = IceWEQ       + ro__SV(ikl,isn) * dzsnSV(ikl,isn)
-    END DO
-
-  END IF
-
-
-  ! +--OUTPUT
-  ! +  ======
-
-  !!      IF (istart.eq.1)                                              THEN
-  !!        write(45,45)dahost,i___SV(lwriSV(1)),j___SV(lwriSV(1)),
-  !!     .              n___SV(lwriSV(1))
-  !! 45     format(a18,10('-'),'Pt.',3i4,60('-'))
-  !!      END IF
-
-  !!      write(45,450) labWEq,IceWEQ,iiceSV(ikl),SnoWEQ
-  !!     .                    ,IceWEQ+SnoWEQ,isnoSV(ikl)
-  !!     .                                  ,drr_SV(ikl)*dt__SV
-  !!     .                                  ,dsn_SV(ikl)*dt__SV
-  !!     .                                  ,BufsSV(ikl)
-  !! 450  format(a6,3x,'  I+S =',f11.4,'(',i2,') +',f11.4,' =',
-  !!     .                       f11.4,'(',i2,')',
-  !!     .             '  drr =', f7.4,
-  !!     .             '  dsn =', f7.4,
-  !!     .             '  Buf =', f7.4)
-
-  return
-end subroutine sisvat_weq
+END MODULE sisvat_weq_mod