- Timestamp:
- Nov 10, 2022, 3:10:32 PM (19 months ago)
- Location:
- BOL/Multi_atlas
- Files:
-
- 2 edited
-
Utils/atlas0.sh (modified) (3 diffs)
-
atlas/atlas_none.py (modified) (3 diffs)
Legend:
- Unmodified
- Added
- Removed
-
BOL/Multi_atlas/Utils/atlas0.sh
r4340 r4345 21 21 #myCC=$DODSDIR/$login/IMclimafCache 22 22 #if [ 1 == 0 ]; then 23 myCC=/scratchu/${login}/climafCache041122run 24 #\rm -rf ${myCC} 23 #myCC=/scratchu/${login}/climafCache041122run 24 myCC=/scratchu/${login}/climafCache$$ 25 \rm -rf ${myCC} 25 26 mkdir -p ${myCC} 26 27 export CLIMAF_CACHE=${myCC} … … 159 160 160 161 echo Fabrication des cartes $var pour $ref $sim $sim-$ref 161 162 ### 163 echo 4 call to climaf instead of 1 to get rid of cache issue 164 ### 162 165 if [ "$var" == "ua" ] || [ "$var" == "va" ] || [ "$var" == "ta" ] || [ "$var" == "hur" ] || [ "$var" == "hus" ]; then 163 echo $python ${local}/$pyatlas -v $var -s $sim -t $ref -p $seas -r $proj --dirpng ${dirpng} 164 $python ${local}/$pyatlas -v $var -s $sim -t $ref -p $seas -r $proj --dirpng ${dirpng} 165 if [ ! -f ${mapok} ]; then 166 echo $python ${local}/$pyatlas -v $var -s $sim -t $ref -p $seas -r $proj --dirpng ${dirpng} 167 $python ${local}/$pyatlas -v $var -s $sim -t $ref -p $seas -r $proj --dirpng ${dirpng} 168 echo $python ${local}/$pyatlas -v $var -s $sim -t $ref -p $seas -r $proj --dirpng ${dirpng} 169 $python ${local}/$pyatlas -v $var -s $sim -t $ref -p $seas -r $proj --dirpng ${dirpng} 170 elif [ ! -f ${refok} ]; then 171 echo $python ${local}/$pyatlas -v $var -s OBS -t $ref -p $seas -r $proj --dirpng ${dirpng} 172 $python ${local}/$pyatlas -v $var -s OBS -t $ref -p $seas -r $proj --dirpng ${dirpng} 173 echo $python ${local}/$pyatlas -v $var -s OBS -t $ref -p $seas -r $proj --dirpng ${dirpng} 174 $python ${local}/$pyatlas -v $var -s OBS -t $ref -p $seas -r $proj --dirpng ${dirpng} 175 elif [ ! -f ${dmapok} ]; then 176 echo $python ${local}/$pyatlas -v $var -s $sim -t $ref -p $seas -r $proj --dirpng ${dirpng} 177 $python ${local}/$pyatlas -v $var -s $sim -t $ref -p $seas -r $proj --dirpng ${dirpng} 178 echo $python ${local}/$pyatlas -v $var -s $sim -t $ref -p $seas -r $proj --dirpng ${dirpng} 179 $python ${local}/$pyatlas -v $var -s $sim -t $ref -p $seas -r $proj --dirpng ${dirpng} 180 fi 166 ### 167 echo $python ${local}/$pyatlas -v $var -s $sim -t $ref -p $seas -r $proj --dirpng ${dirpng} 168 $python ${local}/$pyatlas -v $var -s $sim -t $ref -p $seas -r $proj --dirpng ${dirpng} 169 ### 170 echo $python ${local}/$pyatlas -v $var -s $sim -t $ref -p $seas -r $proj --dirpng ${dirpng} 171 $python ${local}/$pyatlas -v $var -s $sim -t $ref -p $seas -r $proj --dirpng ${dirpng} 172 # 173 echo $python ${local}/$pyatlas -v $var -s $sim -t $ref -p $seas -r $proj --dirpng ${dirpng} 174 $python ${local}/$pyatlas -v $var -s $sim -t $ref -p $seas -r $proj --dirpng ${dirpng} 175 ### 176 echo $python ${local}/$pyatlas -v $var -s $sim -t $ref -p $seas -r $proj --dirpng ${dirpng} 177 $python ${local}/$pyatlas -v $var -s $sim -t $ref -p $seas -r $proj --dirpng ${dirpng} 178 181 179 else 182 echo $python ${local}/$pyatlas -v $var -s $sim -t $ref -p $seas -r $proj --dirpng ${dirpng} --lonmin ${lonmin} --lonmax ${lonmax} --latmin ${latmin} --latmax ${latmax} 183 $python ${local}/$pyatlas -v $var -s $sim -t $ref -p $seas -r $proj --dirpng ${dirpng} --lonmin ${lonmin} --lonmax ${lonmax} --latmin ${latmin} --latmax ${latmax} 184 if [ ! -f ${mapok} ]; then 185 echo $python ${local}/$pyatlas -v $var -s $sim -t $ref -p $seas -r $proj --dirpng ${dirpng} 186 $python ${local}/$pyatlas -v $var -s $sim -t $ref -p $seas -r $proj --dirpng ${dirpng} 187 echo $python ${local}/$pyatlas -v $var -s $sim -t $ref -p $seas -r $proj --dirpng ${dirpng} 188 $python ${local}/$pyatlas -v $var -s $sim -t $ref -p $seas -r $proj --dirpng ${dirpng} 189 elif [ ! -f ${refok} ]; then 190 echo $python ${local}/$pyatlas -v $var -s OBS -t $ref -p $seas -r $proj --dirpng ${dirpng} 191 $python ${local}/$pyatlas -v $var -s OBS -t $ref -p $seas -r $proj --dirpng ${dirpng} 192 echo $python ${local}/$pyatlas -v $var -s OBS -t $ref -p $seas -r $proj --dirpng ${dirpng} 193 $python ${local}/$pyatlas -v $var -s OBS -t $ref -p $seas -r $proj --dirpng ${dirpng} 194 elif [ ! -f ${dmapok} ]; then 195 echo $python ${local}/$pyatlas -v $var -s $sim -t $ref -p $seas -r $proj --dirpng ${dirpng} 196 $python ${local}/$pyatlas -v $var -s $sim -t $ref -p $seas -r $proj --dirpng ${dirpng} 197 echo $python ${local}/$pyatlas -v $var -s $sim -t $ref -p $seas -r $proj --dirpng ${dirpng} 198 $python ${local}/$pyatlas -v $var -s $sim -t $ref -p $seas -r $proj --dirpng ${dirpng} 199 fi 180 echo $python ${local}/$pyatlas -v $var -s $sim -t $ref -p $seas -r $proj --dirpng ${dirpng} --lonmin ${lonmin} --lonmax ${lonmax} --latmin ${latmin} --latmax ${latmax} --pmin ${pmin} 181 $python ${local}/$pyatlas -v $var -s $sim -t $ref -p $seas -r $proj --dirpng ${dirpng} --lonmin ${lonmin} --lonmax ${lonmax} --latmin ${latmin} --latmax ${latmax} --pmin ${pmin} 182 # 183 echo $python ${local}/$pyatlas -v $var -s $sim -t $ref -p $seas -r $proj --dirpng ${dirpng} --lonmin ${lonmin} --lonmax ${lonmax} --latmin ${latmin} --latmax ${latmax} --pmin ${pmin} 184 $python ${local}/$pyatlas -v $var -s $sim -t $ref -p $seas -r $proj --dirpng ${dirpng} --lonmin ${lonmin} --lonmax ${lonmax} --latmin ${latmin} --latmax ${latmax} --pmin ${pmin} 185 # 186 echo $python ${local}/$pyatlas -v $var -s $sim -t $ref -p $seas -r $proj --dirpng ${dirpng} --lonmin ${lonmin} --lonmax ${lonmax} --latmin ${latmin} --latmax ${latmax} --pmin ${pmin} 187 $python ${local}/$pyatlas -v $var -s $sim -t $ref -p $seas -r $proj --dirpng ${dirpng} --lonmin ${lonmin} --lonmax ${lonmax} --latmin ${latmin} --latmax ${latmax} --pmin ${pmin} 188 # 189 echo $python ${local}/$pyatlas -v $var -s $sim -t $ref -p $seas -r $proj --dirpng ${dirpng} --lonmin ${lonmin} --lonmax ${lonmax} --latmin ${latmin} --latmax ${latmax} --pmin ${pmin} 190 $python ${local}/$pyatlas -v $var -s $sim -t $ref -p $seas -r $proj --dirpng ${dirpng} --lonmin ${lonmin} --lonmax ${lonmax} --latmin ${latmin} --latmax ${latmax} --pmin ${pmin} 191 # 200 192 fi 201 193 # … … 227 219 228 220 #IM041122 \rm -rf $DODSDIR/$login/climafCache$$ 221 \rm -rf ${myCC} -
BOL/Multi_atlas/atlas/atlas_none.py
r4330 r4345 272 272 if not is3d(variable): 273 273 # -- apply a mask corresponding to the reference 274 mask = div(reff,reff)275 msimu = mul(wsimu,mask)274 #IMorig mask = div(reff,reff) 275 #IMorig msimu = mul(wsimu,mask) 276 276 # -- Compute the zonal mean 277 zmean=ccdo(msimu, operator='zonmean')278 ref_zmean=ccdo(wreff, operator='zonmean')277 #IMorig zmean=ccdo(msimu, operator='zonmean') 278 #IMorig ref_zmean=ccdo(wreff, operator='zonmean') 279 279 # 280 sim=opts.simulation280 #IMorig sim=opts.simulation 281 281 #if variable in ['zg500']: 282 282 # ref_zmean = ccdo(ref_zmean,operator='-b F32 mulc,1') 283 283 # zmean = ccdo(zmean,operator='-b F32 mulc,1') 284 zmean_fig=curves(cens([sim,ref],zmean,ref_zmean),285 title="",286 lgcols=3,287 options=#'tiYAxisString=""|'+\288 #'+\'+\289 'tmYROn=True|'+\290 'tmYRBorderOn=True|'+\291 'tmYLOn=False|'+\292 'tmYUseRight=True|'+\293 'vpXF=0|'+\294 'vpWidthF=0.66|'+\295 'vpHeightF=0.33|'+\296 'tmYRLabelsOn=True|'+\297 'tmXBLabelFontHeightF=0.018|'+\298 'tmYLLabelFontHeightF=0.016|'+\299 'lgLabelFontHeightF=0.018|'+\300 #'pmLegendSide=Bottom|'+\301 'pmLegendOrthogonalPosF=-0.32|'+\302 'pmLegendParallelPosF=1.0|'+\303 'tmXMajorGrid=True|'+\304 'tmYMajorGrid=True|'+\305 'tmXMajorGridLineDashPattern=2|'+\306 'tmYMajorGridLineDashPattern=2|'+\307 'xyLineThicknessF=8|'+\308 'gsnLeftString='+variable+'|'+\309 'gsnCenterString='+opts.simulation+' vs '+ref+'|'+\310 'gsnRightString='+opts.season+'|'+\311 'gsnStringFontHeightF='+str(stringFontHeight))284 #IMorig zmean_fig=curves(cens([sim,ref],zmean,ref_zmean), 285 #IMorig title="", 286 #IMorig lgcols=3, 287 #IMorig options=#'tiYAxisString=""|'+\ 288 #IMorig #'+\'+\ 289 #IMorig 'tmYROn=True|'+\ 290 #IMorig 'tmYRBorderOn=True|'+\ 291 #IMorig 'tmYLOn=False|'+\ 292 #IMorig 'tmYUseRight=True|'+\ 293 #IMorig 'vpXF=0|'+\ 294 #IMorig 'vpWidthF=0.66|'+\ 295 #IMorig 'vpHeightF=0.33|'+\ 296 #IMorig 'tmYRLabelsOn=True|'+\ 297 #IMorig 'tmXBLabelFontHeightF=0.018|'+\ 298 #IMorig 'tmYLLabelFontHeightF=0.016|'+\ 299 #IMorig 'lgLabelFontHeightF=0.018|'+\ 300 #IMorig #'pmLegendSide=Bottom|'+\ 301 #IMorig 'pmLegendOrthogonalPosF=-0.32|'+\ 302 #IMorig 'pmLegendParallelPosF=1.0|'+\ 303 #IMorig 'tmXMajorGrid=True|'+\ 304 #IMorig 'tmYMajorGrid=True|'+\ 305 #IMorig 'tmXMajorGridLineDashPattern=2|'+\ 306 #IMorig 'tmYMajorGridLineDashPattern=2|'+\ 307 #IMorig 'xyLineThicknessF=8|'+\ 308 #IMorig 'gsnLeftString='+variable+'|'+\ 309 #IMorig 'gsnCenterString='+opts.simulation+' vs '+ref+'|'+\ 310 #IMorig 'gsnRightString='+opts.season+'|'+\ 311 #IMorig 'gsnStringFontHeightF='+str(stringFontHeight)) 312 312 # # Single figures 313 313 page=cpage([[simu_fig]],orientation='landscape', page_trim=True, fig_trim=True) … … 333 333 page_trim=True, fig_trim=True, format='pdf', title=variable+" "+opts.simulation) 334 334 else: 335 page=cpage([[simu_fig,ref_fig],[diff_fig,zmean_fig]],orientation='landscape', page_trim=True, fig_trim=True) 336 pdf_page=cpage([[simu_fig,ref_fig],[diff_fig,zmean_fig]],orientation='landscape', 335 #IMorig page=cpage([[simu_fig,ref_fig],[diff_fig,zmean_fig]],orientation='landscape', page_trim=True, fig_trim=True) 336 page=cpage([[simu_fig,ref_fig],[diff_fig,diff_fig]],orientation='landscape', page_trim=True, fig_trim=True) 337 #IMorig pdf_page=cpage([[simu_fig,ref_fig],[diff_fig,zmean_fig]],orientation='landscape', 338 pdf_page=cpage([[simu_fig,ref_fig],[diff_fig,diff_fig]],orientation='landscape', 337 339 page_trim=True, fig_trim=True, format='pdf', title=variable+" "+opts.simulation) 338 340 pdfargs.append(cfile(pdf_page)) … … 357 359 cell(ref,cfile(ref_fig),thumbnail=thumbnail_size,hover=False)+\ 358 360 cell(text_diff,cfile(diff_fig),thumbnail=thumbnail_size,hover=False)+\ 359 cell('zonal mean',cfile( zmean_fig),thumbnail=thumbnail_size,hover=False)+\361 cell('zonal mean',cfile(diff_fig),thumbnail=thumbnail_size,hover=False)+\ 360 362 cell('page',cfile(page),thumbnail=thumbnail_size,hover=False)+\ 361 363 cell('Pdf',cfile(pdf_page)) 362 364 close_line() 363 # 365 #IMorig: au-dessus cell('zonal mean',cfile(zmean_fig),thumbnail=thumbnail_size,hover=False)+\ 366 364 367 # Finalisons l'index html 365 368 index += close_table()
Note: See TracChangeset
for help on using the changeset viewer.
