#!/usr/bin/env python import os import sys import getpass # just to get user login import shutil import datetime import glob #=============================================================================== #@ Choices #=============================================================================== debug=False #var1='thetao' ; project_list=['woa2009','woa23','glorys12v1','soda3','cmip5','cmip6','trop'] #var1='so' ; project_list=['woa2009','woa23','glorys12v1','soda3','cmip5','cmip6','trop'] #var1='uo'; project_list=['glorys12v1','soda3','cmip5','cmip6','trop'] #var1='vo'; project_list=['glorys12v1','soda3','cmip5','cmip6','trop'] #var1='wo' ; project_list=['cmip5','cmip6'] var1='wod' ; project_list=['cmip5','cmip6'] #var1='hfdsO' ; project_list=['soda3','cmip5','cmip6','trop'] #var1='o2' ; project_list=['woa2009','woa23','cmip5','cmip6','trop'] #var1='no3' ; project_list=['woa2009','woa23','cmip5','cmip6','trop'] html_filename='atlas' title1='PERU '+var1+' atlas' img_pattern_name='img_pattern' F_from_exp_list=False exp_list_name='cmip6_HR' if F_from_exp_list : title1 = title1 + ' (exp_list : '+ exp_list_name+')' else : title1 = title1 + ' (all exp)' img_root_dir='../images/' ### where to generate the exp_list: search_dir='../images/EUC_Upwelling_characterize/PERU.'+var1+'.SAV' if var1 == 'hfdsO' : suffix='_PERU_'+var1+'.pcf025.zip.png' else : suffix='_PERU_'+var1+'_sst.pcf025.zip.png' separateur="_PERU_" num_col=3 # nombre de colonnes pour répéter le nom de l'expérience Klev1="20-26"; lev1='69m-133m' Klev2="24"; lev2="97m" Klev3="18"; lev3="47m" ### Obs if (var1 == 'uo') | (var1 == 'vo') : #obs='GLORYS' # thetao, so, uo, vo obs='glorys12v1' # thetao, so, uo, vo elif (var1 == 'thetao') | (var1 == 'so') | (var1 == 'o2') | (var1 == 'no3'): # obs='woa2009' # thetao, so, o2, no3 obs='woa23' # thetao, so, o2, no3 else : obs='none' #=============================================================================== #@ Definitions #=============================================================================== local_time1=datetime.datetime.strftime(datetime.datetime.now(), '%Y-%m-%d') # '2018-09-19' local_time3=datetime.datetime.strftime(datetime.datetime.now(), 'the %d %b %Y at %H:%M:%S') # '20-Jun-2019 at 16:27:43' local_time4=datetime.datetime.strftime(datetime.datetime.now(), '%Y_%m_%d-%H_%M_%S') # '2019_08_12-18_33_13' script_name=os.path.basename(sys.argv[0]) script_name_root=script_name.split('_3D')[0] script_dir=os.getcwd()+'/' #=============================================================================== #@ Main #=============================================================================== html_filename=html_filename+'.'+var1 if F_from_exp_list : ### read selected/rejected lists infos : with open('exp_list-'+exp_list_name+'.py') as f0: exec(f0.read()) ### we bypass the “for each project” loop because just the selected_list counts project_list=['cmip5'] html_filename=html_filename+'.list_'+exp_list_name else : html_filename=html_filename+'.all_exp' ### to open/create a new html file in the write mode #f1 = open(script_name_root+'_'+var1+'.html', 'w') f1 = open(script_dir+html_filename+'.html', 'w') ### to get id for each exp : with open('id_exp_list.py') as f3: exec(f3.read()) #------------------------------------------------------------------------------- #@@ HTML & javascript header - begin - #------------------------------------------------------------------------------- html_template_begin1 = """ """ f1.write(html_template_begin1) f1.write(' '+title1+'\n') html_template_begin2 = """ """ f1.write(html_template_begin2) f1.write('

'+title1+'

\n') html_template_begin3 = """ """ f1.write(html_template_begin3) #------------------------------------------------------------------------------- #@@ HTML & javascript header - end - #------------------------------------------------------------------------------- os.chdir(search_dir) mm=1 iproject=1 nn=0 #------------------------------------------------------------------------------- #@@ for each project (cmip5, cmip6, woa23)... #------------------------------------------------------------------------------- for project in project_list : #------------------------------------------------------------------------------- #@@ for each exp_file... #------------------------------------------------------------------------------- if F_from_exp_list : exp_file_list=selected_list else : exp_file_list=sorted(glob.glob(project+'*'+suffix)) ii=1 if debug : print('\nfor ',project, ' : exp_file_list : \n',exp_file_list) for exp_file in exp_file_list : exp_name=exp_file.split(separateur)[0] # print(exp_file,':') if debug : print(exp_name) id='none' for jj in range(len(exp_list)) : if exp_list[jj] == exp_name : id=id_list[jj] print(mm,exp_name,':' , id) ### The color definition must be here, because if from_exp_list=true the “for each project” loop is bypassed. color1='#f8e3ff' if 'cmip5' in exp_name : color1='#e2ffd7' if 'cmip6' in exp_name : color1='#d7fffb' if 'trop' in exp_name : color1='#fffed7' ### read : ### Bank of images ### List of image to display for each exp with open(script_dir+img_pattern_name+'.py') as f2: exec(f2.read()) num_col_total=len(img_list) ### create header str1='\n' str2='\n' for num_img in range (len(img_list)) : for jj in range (len(pattern_list)) : header01=pattern_list[jj]['header1'] header02=pattern_list[jj]['header2'] if pattern_list[jj]['name'] == img_list[num_img] : str1=str1+'\n' str2=str2+'\n' str1=str1+'\n' str2=str2+'\n' ### print header (every 4 exp) if nn%4 == 0 : f1.write(str1); f1.write(str2) ### compute the number of the exp in the project : ### ex : n°33/46 if F_from_exp_list : ii_exp=str(ii) else : num_exp_total=len (glob.glob(project+'*'+suffix)) ii_exp=str(ii)+'/'+str(num_exp_total) ### print the number & the name of the exp (every num_col columns) ### ex : n°33/46: cmip5_IPSL-CM5A-MR_historical_19890101-20051231_Omon_1m_r1i1p1_thetao for jj in range(num_col_total) : # print(jj, num_col, jj%num_col) if (jj%num_col == 00) : if (num_col_total-jj) >= num_col : f1.write('\n') f1.write('\n') f1.write('\n') ### display all images for the exp for num_img in range (len(img_list)) : for jj in range (len(pattern_list)) : exp_file00=pattern_list[jj]['exp'] if pattern_list[jj]['name'] == img_list[num_img] : f1.write('\n') f1.write('\n') ii=ii+1 nn=nn+1 mm=mm+1 if ('cmip' in exp_name) | ('trop' in exp_name) : iproject=iproject+1 os.chdir(script_dir) html_template_end = """

'+header01+'	'+header02+'
n°'+ii_exp+': '+'('+id+') '+exp_name+' ('+str(iproject)+')
\n') f1.write('

Webmaster: christophe.hourdin @ locean-ipsl.upmc.fr

Croco-Agrif-Pisces Projects - Home

""" f1.write(html_template_end) # close the file f1.close() print('\n\npattern_list available (generated by img_pattern.py) : \n') for ii in range(len(pattern_list)) : print(pattern_list[ii]['name']) print('\n\nimg_list to be displayed (generated by img_pattern.py) : \n') print(*img_list,sep = '\n')