Changeset 5084 for LMDZ6/trunk/libf/phylmd/inlandsis/surf_inlandsis_mod.F90

Timestamp:

Jul 19, 2024, 6:40:44 PM (2 months ago)

Author:

Laurent Fairhead

Message:

Reverting to r4065. Updating fortran standard broke too much stuff. Will do it by smaller chunks
AB, LF

File:

• LMDZ6/trunk/libf/phylmd/inlandsis/surf_inlandsis_mod.F90 (modified) (20 diffs)

LMDZ6/trunk/libf/phylmd/inlandsis/surf_inlandsis_mod.F90

@@ -638 +638 @@
             END DO
 
-            IF (sissnow(ikl) <= sn_low) THEN  !add snow
-                IF (isnoSV(ikl)>=1) THEN
+            IF (sissnow(ikl) .LE. sn_low) THEN  !add snow
+                IF (isnoSV(ikl).GE.1) THEN
                     dzsnSV(ikl, 1) = dzsnSV(ikl, 1) + sn_add / max(ro__SV(ikl, 1), epsi)
                     toicSV(ikl) = toicSV(ikl) - sn_add
@@ -657 +657 @@
             END IF
 
-            IF (sissnow(ikl) >= sn_upp) THEN  !thinnen snow layer below
+            IF (sissnow(ikl) .ge. sn_upp) THEN  !thinnen snow layer below
                 dzsnSV(ikl, 1) = dzsnSV(ikl, 1) / sn_div
                 toicSV(ikl) = toicSV(ikl) + dzsnSV(ikl, 1) * ro__SV(ikl, 1) / sn_div
@@ -1049 +1049 @@
         ! Objet: Lecture du fichier de conditions initiales pour SISVAT
         !======================================================================
+        include "netcdf.inc"
         !    include "indicesol.h"
 
@@ -1117 +1118 @@
 
         DO isn = 1, nsno
-            IF (isn<=99) THEN
+            IF (isn.LE.99) THEN
                 WRITE(str2, '(i2.2)') isn
                 CALL get_field("AGESNOW" // str2, &
@@ -1127 +1128 @@
         ENDDO
         DO isn = 1, nsno
-            IF (isn<=99) THEN
+            IF (isn.LE.99) THEN
                 WRITE(str2, '(i2.2)') isn
                 CALL get_field("DZSNOW" // str2, &
@@ -1137 +1138 @@
         ENDDO
         DO isn = 1, nsno
-            IF (isn<=99) THEN
+            IF (isn.LE.99) THEN
                 WRITE(str2, '(i2.2)') isn
                 CALL get_field("G2SNOW" // str2, &
@@ -1147 +1148 @@
         ENDDO
         DO isn = 1, nsno
-            IF (isn<=99) THEN
+            IF (isn.LE.99) THEN
                 WRITE(str2, '(i2.2)') isn
                 CALL get_field("G1SNOW" // str2, &
@@ -1157 +1158 @@
         ENDDO
         DO isn = 1, nsismx
-            IF (isn<=99) THEN
+            IF (isn.LE.99) THEN
                 WRITE(str2, '(i2.2)') isn
                 CALL get_field("ETA" // str2, &
@@ -1167 +1168 @@
         ENDDO
         DO isn = 1, nsismx
-            IF (isn<=99) THEN
+            IF (isn.LE.99) THEN
                 WRITE(str2, '(i2.2)') isn
                 CALL get_field("RO" // str2, &
@@ -1177 +1178 @@
         ENDDO
         DO isn = 1, nsismx
-            IF (isn<=99) THEN
+            IF (isn.LE.99) THEN
                 WRITE(str2, '(i2.2)') isn
                 CALL get_field("TSS" // str2, &
@@ -1187 +1188 @@
         ENDDO
         DO isn = 1, nsno
-            IF (isn<=99) THEN
+            IF (isn.LE.99) THEN
                 WRITE(str2, '(i2.2)') isn
                 CALL get_field("HISTORY" // str2, &
@@ -1286 +1287 @@
         IMPLICIT none
 
+        include "netcdf.inc"
         !    include "indicesol.h"
         !    include "dimsoil.h"
@@ -1401 +1403 @@
 
             DO isn = 1, nsno
-                IF (isn<=99) THEN
+                IF (isn.LE.99) THEN
                     WRITE(str2, '(i2.2)') isn
                     CALL put_field(pass, "AGESNOW" // str2, &
@@ -1412 +1414 @@
             ENDDO
             DO isn = 1, nsno
-                IF (isn<=99) THEN
+                IF (isn.LE.99) THEN
                     WRITE(str2, '(i2.2)') isn
                     CALL put_field(pass, "DZSNOW" // str2, &
@@ -1423 +1425 @@
             ENDDO
             DO isn = 1, nsno
-                IF (isn<=99) THEN
+                IF (isn.LE.99) THEN
                     WRITE(str2, '(i2.2)') isn
                     CALL put_field(pass, "G2SNOW" // str2, &
@@ -1434 +1436 @@
             ENDDO
             DO isn = 1, nsno
-                IF (isn<=99) THEN
+                IF (isn.LE.99) THEN
                     WRITE(str2, '(i2.2)') isn
                     CALL put_field(pass, "G1SNOW" // str2, &
@@ -1445 +1447 @@
             ENDDO
             DO isn = 1, nsismx
-                IF (isn<=99) THEN
+                IF (isn.LE.99) THEN
                     WRITE(str2, '(i2.2)') isn
                     CALL put_field(pass, "ETA" // str2, &
@@ -1456 +1458 @@
             ENDDO
             DO isn = 1, nsismx   !nsno
-                IF (isn<=99) THEN
+                IF (isn.LE.99) THEN
                     WRITE(str2, '(i2.2)') isn
                     CALL put_field(pass, "RO" // str2, &
@@ -1467 +1469 @@
             ENDDO
             DO isn = 1, nsismx
-                IF (isn<=99) THEN
+                IF (isn.LE.99) THEN
                     WRITE(str2, '(i2.2)') isn
                     CALL put_field(pass, "TSS" // str2, &
@@ -1478 +1480 @@
             ENDDO
             DO isn = 1, nsno
-                IF (isn<=99) THEN
+                IF (isn.LE.99) THEN
                     WRITE(str2, '(i2.2)') isn
                     CALL put_field(pass, "HISTORY" // str2, &