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

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

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

revision 1.1 by enderton, Wed Oct 5 02:27:21 2005 UTC revision 1.2 by enderton, Tue Oct 18 15:51:19 2005 UTC
# Line 14  mkdir(outDir); Line 14  mkdir(outDir);
14  %                                 Load Data                               %  %                                 Load Data                               %
15  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
16    
17  % datacpl = rdmnc(cplFiles,'iter','HFtave','TXtave',...  datacpl = rdmnc(cplFiles,'iter','HFtave','TXtave',...
18  %                          'TYtave','FWtave','SICtave');                           'TYtave','FWtave','SICtave');
19  % datagrd = rdmnc(gridFiles,'rA','XG','YG','XC','YC');  datagrd = rdmnc(gridFiles,'rA','XG','YG','XC','YC');
20  % RAC = datagrd.rA;  RAC = datagrd.rA;
21  % XC = datagrd.XC; XG = datagrd.XG(1:size(RAC,1),1:size(RAC,2));  XC = datagrd.XC; XG = datagrd.XG(1:size(RAC,1),1:size(RAC,2));
22  % YC = datagrd.YC; YG = datagrd.YG(1:size(RAC,1),1:size(RAC,2));  YC = datagrd.YC; YG = datagrd.YG(1:size(RAC,1),1:size(RAC,2));
23  % save('DataPerCpl.mat','datacpl','RAC','XG','YG','XC','YC');  %save('DataPerCpl.mat','datacpl','RAC','XG','YG','XC','YC');
24  load('DataPerCpl.mat');  %load('DataPerCpl.mat');
25    
26    
27  % If a new focn is to be calculated, the old Focn must be loaded.  % If a new focn is to be calculated, the old Focn must be loaded.
# Line 62  end Line 62  end
62    
63  for ifield = 1:length(mncfld)  for ifield = 1:length(mncfld)
64      FocnOld = dataFocnOld.(mncfld{ifield});      FocnOld = dataFocnOld.(mncfld{ifield});
65        if isequal(mncfld{ifield},'FW')
66            FocnOld = DataCorrections(FocnOld,mncfld{ifield},gravity,RAC);
67        end
68      FocnNew = mean(dataCplSelect.([mncfld{ifield},'tave']),3);      FocnNew = mean(dataCplSelect.([mncfld{ifield},'tave']),3);
69      FocnNew = DataCorrections(FocnNew,mncfld{ifield},gravity,RAC);      FocnNew = DataCorrections(FocnNew,mncfld{ifield},gravity,RAC);
70      if FocnOnly      if FocnOnly
# Line 171  function CorrectedData = ... Line 174  function CorrectedData = ...
174      else      else
175          CorrectedData = DataToCorrect;          CorrectedData = DataToCorrect;
176      end      end
           
177            
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22