Parent Directory
| 
Revision Log
| 
Revision Graph
>> clean up and speed up transport computations. - added routines gcmfaces_edge_mask infer edge masks associated with an interior subdomain mask gcmfaces_subset extract subset of unmasked points from field (e.g. for transports) - modified routines gcmfaces_lines_zonal use gcmfaces_edge_mask and make field names in LATS_MASKS more explicit (mskCint, mskCedge, etc.) gcmfaces_lines_transp same thing for section masks; also remove the longer arc masks (we only use the shorter arcs) calc_zonmean_T use new names (see gcmfaces_lines_zonal) calc_zonmedian_T use new names (see gcmfaces_lines_zonal) calc_overturn use new names (see gcmfaces_lines_zonal) and gcmfaces_subset calc_MeridionalTransport use new names (see gcmfaces_lines_zonal) and gcmfaces_subset calc_transports use new names (see gcmfaces_lines_transp) and gcmfaces_subset disp_transport fix case of a single record
| 1 | gforget | 1.1 | function []=gcmfaces_lines_zonal(varargin); |
| 2 | %object: define the set of quasi longitudinal lines | ||
| 3 | % along which transports will integrated; | ||
| 4 | % LATS_MASKS that will be added to mygrid. | ||
| 5 | %optional input: | ||
| 6 | % LATS_VAL is the latitudes vector ([-89:89]' by default) | ||
| 7 | |||
| 8 | global mygrid; | ||
| 9 | |||
| 10 | if nargin>0; LATS_VAL=varargin{1}; else; LATS_VAL=[-89:89]'; end; | ||
| 11 | |||
| 12 | for iy=1:length(LATS_VAL); | ||
| 13 | gforget | 1.2 | |
| 14 | mskCint=1*(mygrid.YC>=LATS_VAL(iy)); | ||
| 15 | [mskCedge,mskWedge,mskSedge]=gcmfaces_edge_mask(mskCint); | ||
| 16 | |||
| 17 | gforget | 1.1 | %store: |
| 18 | if iy==1; | ||
| 19 | gforget | 1.2 | LATS_MASKS=struct('mskCint',mskCint,'mskCedge',mskCedge,'mskWedge',mskWedge,'mskSedge',mskSedge,'lat',LATS_VAL(iy)); |
| 20 | gforget | 1.1 | else; |
| 21 | gforget | 1.2 | LATS_MASKS(iy)=struct('mskCint',mskCint,'mskCedge',mskCedge,'mskWedge',mskWedge,'mskSedge',mskSedge,'lat',LATS_VAL(iy)); |
| 22 | gforget | 1.1 | end; |
| 23 | |||
| 24 | end; | ||
| 25 | |||
| 26 | mygrid.LATS_MASKS=LATS_MASKS; | ||
| 27 | |||
| 28 |
| ViewVC Help | |
| Powered by ViewVC 1.1.22 |