38 |
integer myiter |
integer myiter |
39 |
_RL mytime |
_RL mytime |
40 |
|
|
|
#ifdef ALLOW_BULKFORMULAE |
|
|
|
|
41 |
c == local variables == |
c == local variables == |
42 |
|
|
43 |
integer bi,bj |
integer bi,bj |
51 |
integer ilnblnk |
integer ilnblnk |
52 |
external ilnblnk |
external ilnblnk |
53 |
|
|
54 |
#ifndef ALLOW_ATM_WIND |
_RL tmp1 |
55 |
_RL TMP1 |
_RL tmp2 |
56 |
_RL TMP2 |
_RL tmp3 |
57 |
_RL TMP3 |
_RL tmp4 |
58 |
_RL TMP4 |
_RL tmp5 |
|
_RL TMP5 |
|
|
#endif |
|
59 |
|
|
60 |
c == end of interface == |
c == end of interface == |
61 |
|
|
89 |
endif |
endif |
90 |
#else /* ifndef ALLOW_ATM_WIND */ |
#else /* ifndef ALLOW_ATM_WIND */ |
91 |
c |
c |
|
#ifdef ALLOW_ATM_TEMP |
|
|
c |
|
92 |
c The variables us, sh and rdn have to be computed from |
c The variables us, sh and rdn have to be computed from |
93 |
c given wind stresses inverting relationship for neutral |
c given wind stresses inverting relationship for neutral |
94 |
c drag coeff. cdn. |
c drag coeff. cdn. |
121 |
us(i,j,bi,bj) = (tmp4 + tmp5)**(1/3) + |
us(i,j,bi,bj) = (tmp4 + tmp5)**(1/3) + |
122 |
& tmp3**2 * (tmp4 + tmp5)**(-1/3) - tmp3 |
& tmp3**2 * (tmp4 + tmp5)**(-1/3) - tmp3 |
123 |
endif |
endif |
124 |
c |
uwind(i,j,bi,bj) = us(i,j,bi,bj)*cw(i,j,bi,bj) |
125 |
#endif /* ALLOW_ATM_TEMP */ |
vwind(i,j,bi,bj) = us(i,j,bi,bj)*sw(i,j,bi,bj) |
126 |
c |
c |
127 |
#endif /* ifndef ALLOW_ATM_WIND */ |
#endif /* ifndef ALLOW_ATM_WIND */ |
128 |
|
|
143 |
enddo |
enddo |
144 |
enddo |
enddo |
145 |
|
|
146 |
#endif /* ALLOW_BULKFORMULAE */ |
return |
|
|
|
147 |
end |
end |