Changeset 379 for trunk/UTIL

Timestamp:

Nov 14, 2011, 12:25:49 PM (13 years ago)

Author:

acolaitis

Message:

Fixed a rare but annoying conflict occuring when calling make_gcm_netcdf with a time dimension which length is equal to longitude/or/latitude/or/altitude

File:

• trunk/UTIL/PYTHON/make_netcdf.py (modified) (4 diffs)

trunk/UTIL/PYTHON/make_netcdf.py

@@ -55 +55 @@
    
    zdims={}
-   zdims[nx]='longitude'
-   zdims[ny]='latitude'
-   zdims[nz]='altitude'
-   zdims[nt]='Time'
+   zdims['longitude']=nx
+   zdims['latitude']=ny
+   zdims['altitude']=nz
+   zdims['Time']=nt
    print zdims 
    # Find which variable uses which dimensions
@@ -70 +70 @@
        j=0
        for dim in zvardim:
-           if dim not in zdims.keys():
+           if dim not in zdims.values():
               print "WARNING -----------------------------"
               print "Dimensions given to subroutine do not match variables dimensions :"
@@ -78 +78 @@
               exit()
            else:
-              zzvardim[j]=zdims[dim]
+              a=get_key(zdims,dim)
+              if len(a) is not 1:
+                 if j is 0:                ##this should solve most conflicts with Time
+                    zzvardim[j]=a[1]
+                 else:
+                    zzvardim[j]=a[0]
+              else:
+                 zzvardim[j]=a[0]
               j=j+1
        zvarcarac[zznames[i]]=zzvardim
        i=i+1
-   
+   print zvarcarac  
+ 
    #########################
    ## Clean previous file ##
@@ -130 +138 @@
    ## close file
    file.close()
+
+def find_key(dic, val):
+    """return the key of dictionary dic given the value"""
+    return [k for k, v in dic.iteritems() if v == val][0]
+
+def get_key(self, value):
+    """find the key(s) as a list given a value"""   
+    return [item[0] for item in self.items() if item[1] == value]