[1] | 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 ztimemax = ((5.*30.)+21.)* 24.;// ! 21th june in hours |
---|
| 24 | double ztimemin = ztimemax+rjjhh*raajj/2.;// ! 21th december in hours |
---|
| 25 | // ! ! NB: rjjhh * raajj / 4 = 1 seasonal cycle in hours |
---|
| 26 | |
---|
| 27 | // ! mean intensity at 0.105 ; srqt(2) because projected with 45deg angle |
---|
| 28 | double ztau = 0.105/sqrt(2.); |
---|
| 29 | // ! seasonal oscillation intensity |
---|
| 30 | double ztau_sais = 0.015; |
---|
| 31 | |
---|
| 32 | double ztaun = ztau-ztau_sais*cos((ztime-ztimemax)/(ztimemin-ztimemax)*rpi); |
---|
| 33 | |
---|
| 34 | |
---|
| 35 | // ! domain from 15deg to 50deg and 1/2 period along 14deg |
---|
| 36 | // ! so 5/4 of half period with seasonal cycle |
---|
| 37 | YS1 = - ztaun*sin(rpi*(gphiu(Yi,Yj)-15.)/(29.-15.)); |
---|
| 38 | YS2 = ztaun*sin(rpi*(gphiv(Yi,Yj)-15.)/(29.-15.)); |
---|
| 39 | |
---|
| 40 | //JUL A MODIFIER APRES TEST |
---|
[9] | 41 | // YS1= 2; |
---|
| 42 | //YS2=0; |
---|
[1] | 43 | //FIN MODIF JUL |
---|
| 44 | } |
---|
| 45 | // |
---|
| 46 | } |
---|
| 47 | |
---|
| 48 | //=========================================================================== |
---|
| 49 | // methode backward |
---|
| 50 | |
---|
| 51 | backward () |
---|
| 52 | { |
---|
| 53 | |
---|
| 54 | } |
---|
| 55 | |
---|
| 56 | //=========================================================================== |
---|
| 57 | //********************* FIN DU MODULE sbcgyre_utau_vtau ********************** |
---|
| 58 | |
---|