Skip to content

Commit

Permalink
Merge branch 'main' into feature/borescan_update_legacy_bcs_path
Browse files Browse the repository at this point in the history
  • Loading branch information
gmao-rreichle committed May 3, 2023
2 parents 07bba17 + 3a4c3f0 commit 5861717
Show file tree
Hide file tree
Showing 13 changed files with 527 additions and 413 deletions.
15 changes: 15 additions & 0 deletions CHANGELOG.md
Original file line number Diff line number Diff line change
Expand Up @@ -17,6 +17,21 @@ and this project adheres to [Semantic Versioning](https://semver.org/spec/v2.0.0

### Deprecated

## [1.1.1] - 2023-03-29

### Fixed

- Fixed bad indentation in remap_restarts.py

## [1.1.0] - 2023-03-10

- Plot updates for MOIST variables.
- Add stretched grid to regrid and remap_restarts

### Fixed

- Plot updates to fix QBO labeling.

## [1.0.2] - 2023-03-07

### Added
Expand Down
22 changes: 11 additions & 11 deletions plots/grads_util/corcmp_plot.gs
Original file line number Diff line number Diff line change
Expand Up @@ -513,7 +513,7 @@ while( m<=mexps )
endwhile
'define zvarave = zvarave / 'mexps
'define se = sqrt( (zvar0 + zvarave)/'numfiles' )'
'define dx = se*'critval90
'define dx = se*'critval95
'define rUave = (exp( 2*(zave0+dx))-1)/(exp( 2*(zave0+dx))+1)'
'define rLave = (exp( 2*(zave0-dx))-1)/(exp( 2*(zave0-dx))+1)'

Expand All @@ -530,7 +530,7 @@ while( m<=mexps )
'define rUp68'm' = 2*(exp( 2*dx)-1)/(exp( 2*dx)+1)'
'define rLp68'm' = 2*(exp(-2*dx)-1)/(exp(-2*dx)+1)'

'define dx = se*'critval90
'define dx = se*'critval95
'define rUp90'm' = 2*(exp( 2*dx)-1)/(exp( 2*dx)+1)'
'define rLp90'm' = 2*(exp(-2*dx)-1)/(exp(-2*dx)+1)'

Expand Down Expand Up @@ -643,7 +643,7 @@ endwhile
'draw ylab Anomaly Correlation'


* Plot 90% Confidence Intervals for Synoptic Variance from Average of All Experiment
* Plot 95% Confidence Intervals for Synoptic Variance from Average of All Experiment
* ----------------------------------------------------------------------------------
'set cmark 0'
'set cthick 2'
Expand Down Expand Up @@ -715,7 +715,7 @@ endif
m = m + 1
endwhile
* Compute Upper and Lower Bounds for 90% Confidence Interval Values for Synoptic Variability
* Compute Upper and Lower Bounds for 95% Confidence Interval Values for Synoptic Variability
* ------------------------------------------------------------------------------------------
'set dfile 'dfile
'set t 'tmin
Expand Down Expand Up @@ -767,16 +767,16 @@ endwhile
* Plot Difference plus Significance
* ---------------------------------
axfac = 1.2
axmax = valrUp90.1*axfac
axmin = valrLp90.1*axfac
axmax = valrUp95.1*axfac
axmin = valrLp95.1*axfac
* Compute Axis Limits based on Error Bars
* ---------------------------------------
if( mexps>1 )
m = 2
while( m<=mexps )
if( valrUp90.m*axfac > axmax ) ; axmax = valrUp90.m*axfac ; endif
if( valrLp90.m*axfac < axmin ) ; axmin = valrLp90.m*axfac ; endif
if( valrUp95.m*axfac > axmax ) ; axmax = valrUp95.m*axfac ; endif
if( valrLp95.m*axfac < axmin ) ; axmin = valrLp95.m*axfac ; endif
m = m + 1
endwhile
else
Expand Down Expand Up @@ -836,7 +836,7 @@ if( mexps=1 )
'd rUp99'm'*1000;rLp99'm'*1000'
else
'set ccolor 'expcol.m
'd rUp90'm'*1000;rLp90'm'*1000'
'd rUp95'm'*1000;rLp95'm'*1000'
endif

m = m + 1
Expand Down Expand Up @@ -1014,9 +1014,9 @@ endwhile
'!/bin/mkdir -p 'SOURCE'/corcmp'
if( nday = ndaymax )
'myprint -name 'SOURCE'/corcmp/stats_'label'_corcmp_'reg'_'level'_'months' -rotate 90 -density 100x100'
'myprint2 -name 'SOURCE'/corcmp/stats_'label'_corcmp_'reg'_'level'_'months' -rotate 90 -density 100x100'
else
'myprint -name 'SOURCE'/corcmp/stats_'label'_corcmp_'reg'_'level'_'months'_'nday'DAY -rotate 90 -density 100x100'
'myprint2 -name 'SOURCE'/corcmp/stats_'label'_corcmp_'reg'_'level'_'months'_'nday'DAY -rotate 90 -density 100x100'
endif
if( debug = "TRUE" )
Expand Down
Loading

0 comments on commit 5861717

Please sign in to comment.