1 | |
---|
2 | /*************************************************************************** |
---|
3 | module classe sbcgyre_utau_vtau.h - description |
---|
4 | ***************************************************************************/ |
---|
5 | // Mohamed Berrada |
---|
6 | // locean-ipsl.upmc, Paris, Mars 18, 2009 |
---|
7 | //=========================================================================== |
---|
8 | // methode forward |
---|
9 | forward () |
---|
10 | { /* verification norm(xY-xF),norm(xY-xF)/norm(xF) |
---|
11 | utau 6.368455529132161e-16, 3.809412027961928e-16 |
---|
12 | vtau 6.451334920495171e-16, 3.858987902188287e-16 |
---|
13 | */ |
---|
14 | if(Yt==TU){ |
---|
15 | YS1=0.; |
---|
16 | YS2=0.; |
---|
17 | } |
---|
18 | else{ |
---|
19 | // !accumulates days of previous months of this year |
---|
20 | // ! day (in hours) since january the 1st |
---|
21 | /* double ztime = float(Yt-TU)*rdt/(rmmss*rhhmm)// ! incrementation in hour |
---|
22 | - (nyear - 1) * rjjhh * raajj;// ! - nber of hours the precedent years */ |
---|
23 | double ztime = kt*rdt/(rmmss*rhhmm)// ! incrementation in hour |
---|
24 | - (nyear - 1) * rjjhh * raajj;// ! - nber of hours the precedent years |
---|
25 | double ztimemax = ((5.*30.)+21.)* 24.;// ! 21th june in hours |
---|
26 | double ztimemin = ztimemax+rjjhh*raajj/2.;// ! 21th december in hours |
---|
27 | // ! ! NB: rjjhh * raajj / 4 = 1 seasonal cycle in hours |
---|
28 | |
---|
29 | // ! mean intensity at 0.105 ; srqt(2) because projected with 45deg angle |
---|
30 | double ztau = 0.105/sqrt(2.); |
---|
31 | // ! seasonal oscillation intensity |
---|
32 | double ztau_sais = 0.015; |
---|
33 | |
---|
34 | double ztaun = ztau-ztau_sais*cos((ztime-ztimemax)/(ztimemin-ztimemax)*rpi); |
---|
35 | |
---|
36 | |
---|
37 | // ! domain from 15deg to 50deg and 1/2 period along 14deg |
---|
38 | // ! so 5/4 of half period with seasonal cycle |
---|
39 | YS1 = - ztaun*sin(rpi*(gphiu(Yi,Yj)-15.)/(29.-15.)); |
---|
40 | YS2 = ztaun*sin(rpi*(gphiv(Yi,Yj)-15.)/(29.-15.)); |
---|
41 | |
---|
42 | //JUL A MODIFIER APRES TEST |
---|
43 | // YS1= 2; |
---|
44 | //YS2=0; |
---|
45 | //FIN MODIF JUL |
---|
46 | } |
---|
47 | // |
---|
48 | } |
---|
49 | |
---|
50 | //=========================================================================== |
---|
51 | // methode backward |
---|
52 | |
---|
53 | backward () |
---|
54 | { |
---|
55 | |
---|
56 | } |
---|
57 | |
---|
58 | //=========================================================================== |
---|
59 | //********************* FIN DU MODULE sbcgyre_utau_vtau ********************** |
---|
60 | |
---|