| 1 |
% MITgcm_contrib Matlab routines of Guillaume MAZE |
% ECCO2: potential vorticity toolbox |
| 2 |
|
% |
| 3 |
|
% This package tries to provide some useful and simple routines to compute, visualize and |
| 4 |
|
% analyze Potential Vorticity from the global high resolution (1/8deg) simulation of the |
| 5 |
|
% MITgcm. |
| 6 |
|
% Routines are as general as possible for extended applications, but note that they were |
| 7 |
|
% developped to focus on the Western Atlantic region for the CLIMODE project. |
| 8 |
|
% Enjoy ! |
| 9 |
% |
% |
| 10 |
% gmaze@mit.edu |
% gmaze@mit.edu |
| 11 |
|
% Last update: Feb1/2007 |
| 12 |
% |
% |
| 13 |
|
% --------------------------------------------------------------------------------------------- |
| 14 |
|
% PROGRAMS LIST (NOT A FUNCTIONS): |
| 15 |
% |
% |
| 16 |
% FUNCTIONS LIST: |
% eg_main_getPV |
| 17 |
|
% This program is an example of how to define global setup and |
| 18 |
|
% to launch the PV computing. |
| 19 |
|
% eg_write_bin2cdf_latlongrid_subdomain |
| 20 |
|
% This program is an example of how to extract a subdomain from |
| 21 |
|
% a lat/lon grid (1/8) binary file and write it into netcdf. A |
| 22 |
|
% directory is created for each time step. |
| 23 |
|
% eg_write_bin2cdf_csgrid_subdomain |
| 24 |
|
% This program is an example of how to extract a subdomain from |
| 25 |
|
% a cube sphere grid (CS510) binary file and write it into netcdf |
| 26 |
|
% and lat/lon grid (1/4). A directory is created for each time step. |
| 27 |
|
% eg_write_UVbin2cdf_csgrid_subdomain |
| 28 |
|
% Idem, except adapted to U and V fields. |
| 29 |
|
% |
| 30 |
|
% --------------------------------------------------------------------------------------------- |
| 31 |
|
% FUNCTIONS LIST 1: NETCDF FILES DIAGNOSTICS |
| 32 |
|
% From netcdf files contained into SNAPSHOT sub-directory of the |
| 33 |
|
% ./netcdf-files/ home folder, these functions ... |
| 34 |
|
% |
| 35 |
|
% A_compute_potential_density(SNAPSHOT) |
| 36 |
|
% Computes potential density SIGMATHETA from potential |
| 37 |
|
% temperature THETA and anomalous salinity SALTanom. |
| 38 |
|
% B_compute_relative_vorticity(SNAPSHOT) |
| 39 |
|
% Computes the 3 components of the relative vorticity from the |
| 40 |
|
% horizontal flow. Take care to the (U,V) grid ! |
| 41 |
|
% C_compute_potential_vorticity(SNAPSHOT,[WANT_SPL_PV]) |
| 42 |
|
% Computes the potential vorticity field from the relative |
| 43 |
|
% vorticity components and the potential density. Option |
| 44 |
|
% WANT_SPL_PV turned 1 (0 by default) makes the function only |
| 45 |
|
% computing the PV based on the planetary vorticity. |
| 46 |
|
% D_compute_potential_vorticity(SNAPSHOT,[WANT_SPL_PV]) |
| 47 |
|
% Multiplies the potential vorticity computed with |
| 48 |
|
% C_COMPUTE_POTENTIAL_VORTICITY by the coefficient: -1/RHO |
| 49 |
|
% Optional flag WANTSPLPV is turned to 0 by default. Turn it to 1 |
| 50 |
|
% if the PV computed was the simple one (f.dSIGMATHETA/dz). It's |
| 51 |
|
% needed for the output netcdf file informations. |
| 52 |
|
% compute_JBz(SNAPSHOT) |
| 53 |
|
% Computes the surface PV flux due to diabatic processes. |
| 54 |
|
% compute_JFz(SNAPSHOT) |
| 55 |
|
% Computes the surface PV flux due to frictionnal forces. |
| 56 |
|
% compute_density(SNAPSHOT) |
| 57 |
|
% Computes density RHO from potential temperature THETA |
| 58 |
|
% and anomalous salinity SALTanom. |
| 59 |
|
% compute_alpha(SNAPSHOT) |
| 60 |
|
% Computes the thermal expansion coefficient ALPHA from potential |
| 61 |
|
% temperature THETA and salinity anomaly SALTanom. |
| 62 |
|
% compute_QEk(SNAPSHOT) |
| 63 |
|
% Computes QEk, the lateral heat flux induced by Ekman currents |
| 64 |
|
% from JFz, the PV flux induced by frictional forces. |
| 65 |
|
% compute_EKL(SNAPSHOT) |
| 66 |
|
% Compute the Ekman Layer Depth from the wind stress and the density |
| 67 |
|
% fields. |
| 68 |
|
% compute_MLD(SNAPSHOT) |
| 69 |
|
% Compute the Mixed Layer Depth from the SST, SSS and potential |
| 70 |
|
% density fields. |
| 71 |
|
% |
| 72 |
|
% --------------------------------------------------------------------------------------------- |
| 73 |
|
% FUNCTIONS LIST 2: ANALYSIS FUNCTIONS |
| 74 |
|
% |
| 75 |
|
% volbet2iso(TRACER,LIMITS,DEPTH,LAT,LONG) |
| 76 |
|
% This function computes the volume embedded between two |
| 77 |
|
% iso-TRACER values and limited eastward, westward and southward |
| 78 |
|
% by fixed limits. |
| 79 |
|
% surfbet2outcrops(TRACER,LIMITS,LAT,LONG) |
| 80 |
|
% This function computes the horizontal surface limited |
| 81 |
|
% by two outcrops of a tracer. |
| 82 |
|
% intbet2outcrops(TRACER,LIMITS,LAT,LONG) |
| 83 |
|
% This function computes the horizontal surface integral |
| 84 |
|
% of the field TRACER on the area limited by two outcrops. |
| 85 |
|
% subfct_getisoS(TRACER,ISO) |
| 86 |
|
% This function determines the iso-surface ISO of the |
| 87 |
|
% 3D field TRACER(Z,Y,X). |
| 88 |
|
% |
| 89 |
|
% --------------------------------------------------------------------------------------------- |
| 90 |
|
% LOWER LEVEL AND SUB-FUNCTIONS LIST: |
| 91 |
|
% |
| 92 |
|
% pv_checkpath |
| 93 |
|
% This function, systematicaly called by the others, ensures that |
| 94 |
|
% all needed sub-directories of the package are in the path. |
| 95 |
% |
% |
| 96 |
% VOLBET2ISO(TRACER,LIMITS,DEPTH,LAT,LONG) |
% --------------------------------------------------------------------------------------------- |
| 97 |
% This function computes the volume limited northward by |
% PS: |
|
% two tracer values and eastward, westward and southward |
|
|
% by fixed limits. |
|
|
% Very simple test available with: TEST_VOLBET2ISO |
|
|
% |
|
|
% SURFBET2OUTCROPS(TRACER,LIMITS,LAT,LONG) |
|
|
% This function computes the horizontal surface limited |
|
|
% by two outcrops of a tracer. |
|
|
% Very simple test available with: TEST_SURFBET2OUTCROPS |
|
|
% |
|
|
% INTBET2OUTCROPS(TRACER,LIMITS,LAT,LONG) |
|
|
% This function computes the horizontal surface integral |
|
|
% of the field TRACER on the area limited by two outcrops. |
|
|
% Very simple test available with: TEST_INTBET2OUTCROPS |
|
| 98 |
% |
% |
| 99 |
|
% > Functions name are case sensitive. |
| 100 |
|
% > See sub-directory "subfct" for further functions. |
| 101 |
|
% > Following packages are required: |
| 102 |
|
% M_MAP: http://www.eos.ubc.ca/~rich/map.html |
| 103 |
|
% SEAWATER: http://www.marine.csiro.au/datacentre/processing.htm |
| 104 |
% |
% |
| 105 |
% PS: |
% --------------------------------------------------------------------------------------------- |
|
% Several functions use the m_map package. |
|
|
% Be sure to get it in the Matlab path. |
|
|
% (M_MAP is available at: http://www.eos.ubc.ca/~rich/map.html) |
|
| 106 |
% |
% |
Parent Directory
| 
Revision Log
| 
Revision Graph
| 
Patch