Changeset 2353 in lmdz_wrf

Timestamp:

Feb 19, 2019, 6:47:45 PM (6 years ago)

Author:

lfita

Message:

Fixing subroutines when Nin == 0!

File:

• trunk/tools/module_scientific.f90 (modified) (8 diffs)

trunk/tools/module_scientific.f90

@@ -5888 +5888 @@
             varout(s1,s2,s3,6) = medv
             varout(s1,s2,s3,7) = Nin*1.
-          ELSE
+          ELSE IF (Nin == 1) THEN
             i1 = gridsin(s1,s2,s3,1,idv)
             IF (i1 > 0 .AND. i1 <= dmaxvarin(1)) THEN
@@ -5898 +5898 @@
               varout(s1,s2,s3,6) = varin(i1)
               varout(s1,s2,s3,7) = Nin*1.
+            ELSE
+              varout(s1,s2,s3,1) = fillval64
+              varout(s1,s2,s3,2) = fillval64
+              varout(s1,s2,s3,3) = fillval64
+              varout(s1,s2,s3,4) = fillval64
+              varout(s1,s2,s3,5) = fillval64
+              varout(s1,s2,s3,6) = fillval64
+              varout(s1,s2,s3,7) = zeroRK
             END IF
+          ELSE
+            varout(s1,s2,s3,1) = fillval64
+            varout(s1,s2,s3,2) = fillval64
+            varout(s1,s2,s3,3) = fillval64
+            varout(s1,s2,s3,4) = fillval64
+            varout(s1,s2,s3,5) = fillval64
+            varout(s1,s2,s3,6) = fillval64
+            varout(s1,s2,s3,7) = zeroRK
           END IF
         END DO
@@ -6018 +6034 @@
             varout(s1,s2,s3,6) = medv
             varout(s1,s2,s3,7) = Nin*1.
-          ELSE
+          ELSE IF (Nin == 1) THEN
             i1 = gridsin(s1,s2,s3,1,1)
             i2 = gridsin(s1,s2,s3,1,2)
@@ -6028 +6044 @@
             varout(s1,s2,s3,6) = varin(i1,i2)
             varout(s1,s2,s3,7) = Nin*1.
+          ELSE
+            varout(s1,s2,s3,1) = fillval64
+            varout(s1,s2,s3,2) = fillval64
+            varout(s1,s2,s3,3) = fillval64
+            varout(s1,s2,s3,4) = fillval64
+            varout(s1,s2,s3,5) = fillval64
+            varout(s1,s2,s3,6) = fillval64
+            varout(s1,s2,s3,7) = zeroRK
           END IF
         END DO
@@ -6149 +6173 @@
               varout(s1,s2,s3,i3,7) = Nin*1.
             END DO
-          ELSE
+          ELSE IF (Nin == 1) THEN
             i1 = gridsin(s1,s2,s3,1,1)
             i2 = gridsin(s1,s2,s3,1,2)
@@ -6159 +6183 @@
             varout(s1,s2,s3,:,6) = varin(i1,i2,:)
             varout(s1,s2,s3,:,7) = Nin*1.
+          ELSE
+            varout(s1,s2,s3,:,1) = fillval64
+            varout(s1,s2,s3,:,2) = fillval64
+            varout(s1,s2,s3,:,3) = fillval64
+            varout(s1,s2,s3,:,4) = fillval64
+            varout(s1,s2,s3,:,5) = fillval64
+            varout(s1,s2,s3,:,6) = fillval64
+            varout(s1,s2,s3,:,7) = zeroRK
           END IF
         END DO
@@ -6414 +6446 @@
                 varout(s1,s2,s3,s4,i3,7) = Nin*1.
               END DO
-            ELSE
+            ELSE IF (Nin == 1) THEN
                 i1 = gridsin(s1,s2,s3,s4,1,1)
                 i2 = gridsin(s1,s2,s3,s4,1,2)
@@ -6424 +6456 @@
                 varout(s1,s2,s3,s4,:,6) = varin(i1,i2,:)
                 varout(s1,s2,s3,s4,:,7) = Nin*1.
+            ELSE
+                varout(s1,s2,s3,s4,:,1) = fillval64
+                varout(s1,s2,s3,s4,:,2) = fillval64
+                varout(s1,s2,s3,s4,:,3) = fillval64
+                varout(s1,s2,s3,s4,:,4) = fillval64
+                varout(s1,s2,s3,s4,:,5) = fillval64
+                varout(s1,s2,s3,s4,:,6) = fillval64
+                varout(s1,s2,s3,s4,:,7) = zeroRK
             END IF
           END DO