Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Adding a new example analysis script. #26

Merged
merged 2 commits into from
Jul 22, 2014

Conversation

nikizadehgfdl
Copy link
Contributor

- EddyKineticEnergy.py calculates the surface eddy kinetic energy annual mean
- Required adding a logarithmic color plot option to 6plot.xyplot
- Serves as an example analysis script that can be called inside the refineDiag script to be done annually
- Motivated by the plot in http://www.gfdl.noaa.gov/cms-filesystem-action/user_files/kd/pdf/c37_dixon_th85b.pdf
- Unaware of the scientific value, just an example.
@nikizadehgfdl
Copy link
Contributor Author

If you're curious what the plot looks like see:
/ptmp/Niki.Zadeh/archive/Niki.Zadeh/tikal_201403_mom6_20140523/OM4_SIS_baseline/gfdl.ncrc2-intel-prod/history/19000101.nc/ptemp/EKE_mean.png

eke_mean

ssv = rootGroup.variables['ssv']
ssv_mean = ssv[:].mean(axis=0)

eke = ((ssu[:] - ssu_mean)**2 + (ssv[:] - ssv_mean)**2)/2
Copy link
Collaborator

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

There should be an i-average of u^2 from u-points to t-points and a j-average of v^2 from v-pointsto t-points. Use numpy.roll().

@adcroft
Copy link
Collaborator

adcroft commented Jun 23, 2014

Nice addition! I noted (in a line comment) a discretization error in the eke calculation which I'll let you fix.

@nikizadehgfdl
Copy link
Contributor Author

I fixed the U to T issue. The new figure is :

eke_mean

@nikizadehgfdl nikizadehgfdl merged commit 70d47e6 into dev/master Jul 22, 2014
@nikizadehgfdl nikizadehgfdl deleted the user/nnz/EddyKineticEnergyAnalysisScript branch July 27, 2016 17:51
gustavo-marques pushed a commit to gustavo-marques/MOM6 that referenced this pull request Sep 7, 2017
Remove dependency on ocean_model_mod; write/read restart files following CESM convention
nikizadehgfdl pushed a commit to nikizadehgfdl/MOM6 that referenced this pull request Oct 9, 2017
- "coupler" is needed for ice-ocean and fully coupled models.
- Using version a61771c "Merge for Ulm release (21.0.0)".
  - Later versions need patches for gnu/pgi compatibility.
- To build ice-ocean models the components land_null and atmos_null
  are still needed which will be added to GitHub sometime soon...
- Closes mom-ocean#26.
- No answer changes.
MJHarrison-GFDL pushed a commit to MJHarrison-GFDL/MOM6 that referenced this pull request Aug 16, 2022
Add optional parameter to register_diag_field "realm" to pass in a string to add a new attribute "modeling_realm".

This new attribute can be used in post-processing to correct the "area"/"volume" field names.

Fixes mom-ocean#26
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

2 participants