--- MITgcm_contrib/enderton/PeriodicCoupling/MakePcForcing.m 2006/06/16 16:57:05 1.3 +++ MITgcm_contrib/enderton/PeriodicCoupling/MakePcForcing.m 2006/06/16 16:59:04 1.4 @@ -22,13 +22,13 @@ % Load Data % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% -% datacpl = rdmnc(cplFiles,'iter','HFtave','TXtave','TYtave','FWtave','T'); -% dataatm = rdmnc(atmFiles,'iter','THETA' ,'ALBVISDF','RADSWG' ,... -% 'DWNLWG' ,'TS' ,'QS' ,... -% 'WINDS' ,'PRECON' ,'PRECLS','T'); -% datagrd = rdmnc(gridFiles,'rA','XG','YG','XC','YC','HFacC'); +datacpl = rdmnc(cplFiles,'iter','HFtave','TXtave','TYtave','FWtave','T'); +dataatm = rdmnc(atmFiles,'iter','THETA' ,'ALBVISDF','RADSWG' ,... + 'DWNLWG' ,'TS' ,'QS' ,... + 'WINDS' ,'PRECON' ,'PRECLS','T'); +datagrd = rdmnc(gridFiles,'rA','XG','YG','XC','YC','HFacC'); % save('DataPerCpl.mat','datacpl','dataatm','datagrd'); -load('DataPerCpl.mat'); +% load('DataPerCpl.mat'); nc = size(datagrd.rA,2); RAC = datagrd.rA; mask = datagrd.HFacC(:,:,1);