Changeset 1105 in lmdz_wrf for trunk

Timestamp:

Sep 13, 2016, 10:52:42 AM (9 years ago)

Author:

lfita

Message:

Adding `julday_360d': Function to provide the julian day of a date in a 360 days/yr (or 12 30-days months) calendar

File:

• trunk/tools/generic_tools.py (modified) (6 diffs)

trunk/tools/generic_tools.py

@@ -55 +55 @@
 # index_mat: Function to provide the coordinates of a given value inside a matrix
 # index_vec: Function to provide the coordinates of a given value inside a vector
+# julday_360d: Function to provide the julian day of a date in a 360 days/yr (or 12 30-days months) calendar
 # list_combos: Function to construct a new list with all possible N-combinations of the list-values
 # list_coincidences: Function to provide the coincidences between two lists
@@ -280 +281 @@
     Nlvs = len(listv)
     for ic in range(Ncombos):
-        print fname + ' Lluis: ic:', ic, 'Nlvs:', Nlvs
         prevlist = list(newlist)
         for lv in prevlist:
@@ -850 +850 @@
     return newdate
 
+def julday_360d(date):
+    """ Function to provide the julian day of a date in a 360 days/yr (or 12 30-days months) calendar
+      date=  date [Y, m, d, h, mi, s]
+      >>> julday_360d([1976,2,17,8,30,2])
+      47
+    """
+    fname = 'julday_360d'
+
+    julday = (date[1]-1)*30 + date[2]
+
+    return julday
+print julday_360d([1976,2,17,8,30,2])
+
+quit()
 class dtsec360dyr(object):
     """ Class to operate a number of seconds to a date in a 360 days/yr (or 12 30-days months) calendar
@@ -9198 +9212 @@
     matype = mat.dtype
 
-    print fname + '; Lluis: matype:', matype,'shape:',mat.shape,'opdims:',opdims
-
     if oper.find(',') != -1:
         opern = oper.split(',')[0]
@@ -9227 +9239 @@
         idim = idim + 1
 
+    if len(ldimop) == 0:
+        print errormsg
+        print '  ' + fname + ': no dimensions found for the operation !!'
+        print '    dimensions in file:', dimns
+        print '    dimensions for the operation:', opdims
+        quit(-1)
+
     Ndop = len(opdims)
     Ntotdims = np.prod(ldimop)
-    print fname + '; Lluis ldimop:', ldimop, 'Ntotdims:', Ntotdims
 
     # Checking for coincidence in shapes
@@ -9251 +9269 @@
             if slicevals[idn] == -1:
                 ijkop = ijkop + 1
+
                 opslice.append(idcalc[ijkop])
             else: