/[MITgcm]/MITgcm_contrib/enderton/PeriodicCoupling/MakePcForcing.m
ViewVC logotype

Diff of /MITgcm_contrib/enderton/PeriodicCoupling/MakePcForcing.m

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.1 by enderton, Tue Feb 28 23:16:02 2006 UTC revision 1.2 by enderton, Thu Mar 9 00:23:38 2006 UTC
# Line 1  Line 1 
1  function MakePcForcing(cplFiles,atmFiles,gridFiles,prevFocnFiles,period,...  function MakePcForcing(cplFiles,atmFiles,gridFiles,prevFocnFiles,...
2                         outDir,FocnRoot,focnRoot,iter,...                         period,outDir,FocnRoot,focnRoot,iter,...
3                         dtC,gWeight,FocnOnly,InspectFocn,Inspectfocn)                         dtC,gWeight,FocnOnly,InspectFocn,Inspectfocn)
4    
5  % Make fields for forcing ocean component of periodic coupling.  % Make fields for forcing ocean component of periodic coupling.
# Line 22  mkdir(outDir); Line 22  mkdir(outDir);
22  %                                 Load Data                               %  %                                 Load Data                               %
23  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
24    
25  % datacpl = rdmnc(cplFiles,'iter','HFtave','TXtave','TYtave','FWtave','T');  datacpl = rdmnc(cplFiles,'iter','HFtave','TXtave','TYtave','FWtave','T');
26  % dataatm = rdmnc(atmFiles,'iter','THETA'   ,'ALBVISDF','RADSWG'  ,...  dataatm = rdmnc(atmFiles,'iter','THETA'   ,'ALBVISDF','RADSWG'  ,...
27  %                                 'DWNLWG'  ,'TS'      ,'QS'      ,...                                  'DWNLWG'  ,'TS'      ,'QS'      ,...
28  %                                 'WINDS'   ,'PRECON'  ,'PRECLS','T');                                  'WINDS'   ,'PRECON'  ,'PRECLS','T');
29  % datagrd = rdmnc(gridFiles,'rA','XG','YG','XC','YC','HFacC');  datagrd = rdmnc(gridFiles,'rA','XG','YG','XC','YC','HFacC');
30  % save('DataPerCpl.mat','datacpl','dataatm','datagrd');  % save('DataPerCpl.mat','datacpl','dataatm','datagrd');
31  load('DataPerCpl.mat');  % load('DataPerCpl.mat');
32  nc = size(datagrd.rA,2);  nc = size(datagrd.rA,2);
33  RAC = datagrd.rA;  RAC = datagrd.rA;
34  mask = datagrd.HFacC(:,:,1);  mask = datagrd.HFacC(:,:,1);
# Line 57  end Line 57  end
57  %                      Select Coupled Field Data                          %  %                      Select Coupled Field Data                          %
58  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
59    
60    if ~isequal(datacpl.iter,dataatm.iter./8)
61        disp('datacpl.iter');    datacpl.iter
62        disp('dataatm.iter/8');  dataatm.iter./8
63        error('Iterations for coupled and atmospheric data do not match.');
64    end
65  test  = ismember(iter,datacpl.iter);  test  = ismember(iter,datacpl.iter);
66  index = ismember(datacpl.iter,iter);  index = ismember(datacpl.iter,iter);
67  if ~isempty(find(test==0))  if ~isempty(find(test==0))
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22