Changeset 1060 in lmdz_wrf for trunk

Timestamp:

Aug 27, 2016, 12:43:55 PM (9 years ago)

Author:

lfita

Message:

Starting to add differences among models

File:

• trunk/tools/model_graphics.py (modified) (3 diffs)

trunk/tools/model_graphics.py

@@ -2896 +2896 @@
 # Models loop
 ##
+
+# dictionary with the experiments of each model
+modexps = {}
+
 for mod in mods:
     print mod
@@ -2910 +2914 @@
     exps, fheaders = exp_headers(mod,cnf)
 
+    modexps[mod] = exps
     # Characteristics of the model
     Modinf = ncvar.model_characteristics(mod,'None','False')
@@ -3162 +3167 @@
           dbg)
 
-#    quit()
 # end of mods loop
 
+### ## #
+# Model differences
+## # ## #
+
+Nmods = len(mods)
+for imod1 in range(0,Nmods-1):
+    mod1 = mods[imod1]
+    exps1 = modexps[mod1]
+    for imod2 in range(imod1+1,Nmods):
+        mod2 = mods[imod2]
+        exps2 = modexps[mod2]
+
+        Smods = mod2 + '-' + mod1
+        print '  ******* ****** ***** **** *** ** * ** *** **** ***** ****** *******'
+        print "  ** '" + Smods + "': Inter models differences  "
+        print '  ******* ****** ***** **** *** ** * ** *** **** ***** ****** *******'
+        print '    experiments mod1:', exps1
+        print '    experiments mod2:', exps2
+
+# Experiments loop
+##
+        difffiles = ['diffop.inf', 'diffvar.inf']
+
+        Nexps1 = len(exps1)
+        Nexps2 = len(exps2)
+        for iexp1 in range(0,Nexps1):
+            exp1 = exps1[iexp1]
+            for iexp2 in range(0,Nexps2):
+                exp2 = exps2[iexp2]
+                Sexps = exp2 + '-' + exp1
+                print '    ' + Sexps + '...'
+                owdir = cnf['ofold'] + '/' + Smods + '/' + exp2 + '-' + exp1
+                sout = sub.call('mkdir -p ' + owdir, shell=True)
+
+                # Removing files with list of differences if should be started from 
+                #   scratch or add differences
+                for fdiff in difffiles:
+                    difff = owdir + '/' + fdiff
+                    if diffscratch:
+                        sub.call('rm -rf' + owdir +' >& /dev/null', shell=True)
+                        sub.call('rm ' + difff +' >& /dev/null', shell=True)
+                        objf = open(owdir+'/all_' + fdiff,'w')
+                        objf.write("## Computation and drawing of differences " +    \
+                          "between '" + mod2+'@'+exp2 + "'-'" + mod1+'@'+exp1 + "'\n")
+                        objf.close()
+                        difff = owdir + '/all_vardiffstatistics.inf'
+                        sub.call('rm ' + difff +' >& /dev/null', shell=True)
+                        objf = open(owdir+'/all_' + fdiff,'w')
+                        objf.write("## Computation of all differences statistics " + \
+                          "between '" + mod2+'@'+exp2 + "'-'" + mod1+'@'+exp1 + "'\n")
+                        objf.close()
+
+                    if adddiffs:
+                        sub.call('rm ' + difff +' >& /dev/null', shell=True)
+
+                # op and var differences
+                diffvarcompf = owdir + '/' + difffiles[0]
+                if not os.path.isfile(diffvarcompf):
+                    alldiffop, alldiffvar, Nopdiffs, Nvardiffs =                     \
+                      diffvarop_listconstruct(cnf, [mod1, mod2], [exp1, exp2], owdir,\
+                      dbg)
+                else: 
+                    print warnmsg
+                    print '  ' + main + ": getting 'op' and 'var'differences to " +  \
+                      "compute already from file !!"
+                    alldiffop, Nopdiffops = read_diff_file(diffvarcompf)
+                    alldiffvar, Nopdiffvars = read_diff_file(owdir+'/'+difffiles[1])
+
+                # End of avoiding to repeat all the experiment search
+
+                print "    For experiments '" + mod2+'@'+exp2+"'-'"+mod+'@'+exp1+    \ 
+                  "' is required to compute:", Nvar, "differences"
+
+                if dbg:
+                    print 'Differences to compute _______'
+                    gen.printing_dictionary(alldiffop)
+                    gen.printing_dictionary(alldiffvar)
+
+# Computing differences
+##
+                print "    Computing operation differences ..."
+                compute_op_diffs(cnf, alldiffop, owdir, diffscratch, dbg)
+                print "    Computing variable differences ..."
+                compute_var_diffs(cnf, alldiffvar, owdir, diffscratch, dbg)
+
+# Plotting operation differences
+##
+                print "  " + main + ": Plotting operation differences' figures ..."
+                dirfigf = owdir + '/diffopplotsdraw.inf'
+                if figscratch:
+                    sout = sub.call('rm ' + dirfigf + ' >& /dev/null', shell=True)
+
+                objf = open(owdir+'/all_diffopfigures.inf','w')
+                objf.write("## Drawing of all operation difference figures \n")
+                objf.close()
+
+                if adddifffigures:
+                    sout = sub.call('rm ' + dirfigf + ' >& /dev/null', shell=True)
+
+                if not os.path.isfile(dirfigf):
+                    listplots, Nplt = plots_listconstruct(cnf, 'PLOTDIFFOP', dirfigf,\
+                      owdir, dbg)
+                else:
+                    print warnmsg
+                    print '  ' + main + ": getting plots to draw already from file !!"
+                        listplots, Nplt = read_plot_file(dirfigf)
+
+                # End of avoiding to repeat all the plots search
+
+                print "  For experiment 'operation' differences '" + Sexps + "' is "+\
+                  "required to plot:" , Nplt, "plots"
+
+                if dbg:
+                    print 'Plots to draw _______'
+                    gen.printing_dictionary(listplots)
+
+                draw_diff_plots(cnf, listplots, owdir, alldiffop, 'diffop',          \
+                  figdiffscratch, dbg)
+                draw_diff_plots(cnf, listplots, owdir, alldiffvar, 'diffvar',        \
+                  figdiffscratch, dbg)
+
+            # end of exp2 loop
+        # end of exp1 loop
+    # end of mod2 loop
+# end of mod1 loop
+