Changeset 5116 for LMDZ6/branches/Amaury_dev/libf/phylmd/inlandsis/sisvat_bsn.f90

Timestamp:

Jul 24, 2024, 2:54:37 PM (2 months ago)

Author:

abarral

Message:

rename modules properly lmdz_*
move ismin, ismax, minmax into new lmdz_libmath.f90
(lint) uppercase fortran keywords

File:

• LMDZ6/branches/Amaury_dev/libf/phylmd/inlandsis/sisvat_bsn.f90 (modified) (4 diffs)

LMDZ6/branches/Amaury_dev/libf/phylmd/inlandsis/sisvat_bsn.f90

@@ -33 +33 @@
 
 
-  integer :: ikl   ,isn
-  real :: h_mmWE                        ! Eroded Snow Layer Min Thickness
-  real :: dbsaux(knonv)                 ! Drift Amount   (Dummy Variable)
-  real :: dzweqo,dzweqn,bsno_x          ! Conversion variables for erosion
-  real :: dz_new,rho_new
-  real :: snofOK                        ! Threshd Snow Fall
-  real :: Fac                           ! Correction factor for erosion
-  real :: densif                        ! Densification rate if erosion
+  INTEGER :: ikl   ,isn
+  REAL :: h_mmWE                        ! Eroded Snow Layer Min Thickness
+  REAL :: dbsaux(knonv)                 ! Drift Amount   (Dummy Variable)
+  REAL :: dzweqo,dzweqn,bsno_x          ! Conversion variables for erosion
+  REAL :: dz_new,rho_new
+  REAL :: snofOK                        ! Threshd Snow Fall
+  REAL :: Fac                           ! Correction factor for erosion
+  REAL :: densif                        ! Densification rate if erosion
 
   ! +--DATA
@@ -73 +73 @@
     if((dzweqo-dzweqn)>0                    .and. &
           dzsnSV(ikl,isn)>0                    .and. &
-          ro__SV(ikl,max(1,isnoSV(ikl)))<roBdSV) then
-
+          ro__SV(ikl,max(1,isnoSV(ikl)))<roBdSV) THEN
     !characteristic time scale for drifting snow compaction set to 24h
     !linear densification rate [kg/m3/s] over 24h
@@ -84 +83 @@
     Fac         = max(0.,min(1.,Fac))
 
-    if (ro__SV(ikl,max(1,isnoSV(ikl)))>roBdSV) then
+    if (ro__SV(ikl,max(1,isnoSV(ikl)))>roBdSV) THEN
       densif=densif*Fac
     endif
@@ -94 +93 @@
     endif
 
-    if(dzsnSV(ikl,isn)>0 .and.dzsnSV(ikl,isn)<0.0001)then
+    IF(dzsnSV(ikl,isn)>0 .and.dzsnSV(ikl,isn)<0.0001)THEN
     dbs_SV(ikl) = dbs_SV(ikl)+ dzsnSV(ikl,isn)*ro__SV(ikl,isn)
     dbs_Er(ikl) = dbs_Er(ikl)+ dzsnSV(ikl,isn)*ro__SV(ikl,isn)