Changeset 1890 in lmdz_wrf

Timestamp:

Apr 4, 2018, 5:24:51 PM (7 years ago)

Author:

lfita

Message:

Adding calculation of conform frequency in `itime'

File:

• trunk/tools/nc_var_tools.py (modified) (1 diff)

trunk/tools/nc_var_tools.py

@@ -24300 +24300 @@
     print 'diffdate_2nd-1st:', diffdate
     print 'diff_seconds_2nd-1st:', totsecs
+    secfreq = totsecs
+
+    if np.mod(secfreq, 100 * 365 * 24 * 3600.) == 0.:
+        Sfreq = str(int(secfreq/(100 * 365 * 24 * 3600.))) + 'C'
+    elif np.mod(secfreq, 365 * 24 * 3600.) == 0.: 
+        Sfreq = str(int(secfreq/(365 * 24 * 3600.))) + 'Y'
+    elif np.mod(secfreq, 7 * 24 * 3600.) == 0.: 
+        Sfreq = str(int(secfreq/(7 * 24 * 3600.))) + 'W'
+    elif np.mod(secfreq, 24 * 3600.) == 0.:
+        Sfreq = str(int(secfreq/(24.*3600.))) + 'D'
+    elif np.mod(secfreq,3600.) == 0.: 
+        Sfreq = str(int(secfreq/3600.)) + 'H'
+    elif np.mod(secfreq,60.) == 0.:
+        Sfreq = str(int(secfreq/60.)) + 'M'
+    else:
+        print '  ' + fname + ": data from file does not conform "+\
+          "any know time-units!!"
+        print '    frequency from file:', secfreq, ' seconds _______'
+        print '      century:', secfreq/(100 * 365 * 24 * 3600.)
+        print '      year (365d):', secfreq/(365 * 24 * 3600.)
+        print '      week:', secfreq/(7 * 24 * 3600.)
+        print '      day:', secfreq/(24 * 3600.)
+        print '      hour:', secfreq/(3600.)
+        print '      minute:', secfreq/(60.)
+        print '      second:', secfreq
+    print 'diff_conform:', Sfreq
 
     for it in range(1,dimt-1):