Documentation of coup_windstress

Global Index (short | long) | Local contents | Local Index (short | long)

Help text

 use c = 60m/s, or rossby radius ~10degrees ==> 2c=12000.

Listing of script coup_windstress

im = 180;
jm = 102;
nypt = [45 58]; nxpt = [36 147];
dd = 100e5;
global OMEGA RADUS; 
rearth = RADUS*100;
bcor = 2*OMEGA./rearth;
taux = zeros(jm, im);
      for j=1:jm
         for i = nxpt(1):nxpt(2)
            taux(j,i) = cos(dd*180.*(i-im/2.-1.)/(100.*rearth))* ...
                exp(-1*bcor*((dd*(j-jm/2.+0.5))^2)/12000.);
         end
      end