Changeset 993 for codes/icosagcm/devel/src/output/set_bounds.f90

Timestamp:

11/15/19 11:03:17 (5 years ago)

Author:

rpennel

Message:

devel : temporary hack to use both cellset and read_metric

File:

• codes/icosagcm/devel/src/output/set_bounds.f90 (modified) (3 diffs)

codes/icosagcm/devel/src/output/set_bounds.f90

@@ -32 +32 @@
     cells%ncell = n
     
-    ! now set bounds
-    ALLOCATE(cells%ij(n), cells%lon(n), cells%lat(n), cells%ind_glo(n))
-    ALLOCATE(cells%bnds_lon(0:5,n), cells%bnds_lat(0:5,n))
+    IF ( .NOT. ALLOCATED(cells%ij) ) THEN 
+        ! now set bounds
+        ALLOCATE(cells%ij(n), cells%lon(n), cells%lat(n), cells%ind_glo(n))
+        ALLOCATE(cells%bnds_lon(0:5,n), cells%bnds_lat(0:5,n))
+    END IF
 
     n=0
@@ -74 +76 @@
     cells%ncell = n
 
-    ! now set bounds
-    ALLOCATE(cells%ind_glo(n)) ! not set but must be allocated
-    ALLOCATE(cells%ij(n), cells%lon(n), cells%lat(n))
-    ALLOCATE(cells%bnds_lon(0:2,n), cells%bnds_lat(0:2,n))
+    IF ( .NOT. ALLOCATED( cells%ind_glo ) ) THEN 
+        ! now set bounds
+        ALLOCATE(cells%ind_glo(n)) ! not set but must be allocated
+        ALLOCATE(cells%ij(n), cells%lon(n), cells%lat(n))
+        ALLOCATE(cells%bnds_lon(0:2,n), cells%bnds_lat(0:2,n))
+    END IF
 
     n=0
@@ -138 +142 @@
     cells%ncell = n
 
-    ! now set bounds
-    ALLOCATE(cells%ij(n), cells%lon(n), cells%lat(n), cells%ind_glo(n))
-    ALLOCATE(cells%sgn(n)) ! flip sign when reading/writing
-    ALLOCATE(cells%bnds_lon(2,n), cells%bnds_lat(2,n))
+    IF ( .NOT. ALLOCATED( cells%ij ) ) THEN 
+        ! now set bounds
+        ALLOCATE(cells%ij(n), cells%lon(n), cells%lat(n), cells%ind_glo(n))
+        ALLOCATE(cells%sgn(n)) ! flip sign when reading/writing
+        ALLOCATE(cells%bnds_lon(2,n), cells%bnds_lat(2,n))
+    END IF 
 
     CALL swap_dimensions(ind)