Thanks! I implemented it and you can find the code as follows. It does
not have proper error checking mechanisms but it could be used as a
base. In this case, i used the code from airsea Matlab toolbox and
reimplemented it in NCL.
undef("cdnlp")
function cdnlp(sp[*]:numeric, z:numeric)
begin
kappa = 0.4
tol = 0.001
a = log(z/10.0)/kappa
cd = sp
cd = 1.15e-3
u10 = sp/(1.0+a*sqrt(cd))
u10o = sp*0.0
ii = (abs(u10-u10o) .gt. tol)
do while(any(ii))
u10o = u10
cd = 4.9e-4+6.5e-5*u10o
cd = where(u10o .lt. 10.15385, 1.15e-3, cd)
u10 = sp/(1.0+a*sqrt(cd))
ii = (abs(u10-u10o) .gt. tol)
end do
opt = True
opt@cd = cd
opt@wspd = u10
return(opt)
end
Regards,
--ufuk
On 10/12/13 02:27, Dennis Shea wrote:
> Hi Ufuk,
>
> To my knowledge,no one has implemented the Large & Pond drag coefficient
> approach within NCL.
>
> Regards
> D
>
>
> On 12/7/13, 3:11 AM, u.utku.turuncoglu@be.itu.edu.tr wrote:
>> Hi,
>>
>> I just wonder that is there any NCL function to calculate the drag
>> coefficient based on Large&Pond or other formulation?
>>
>> Regards,
>>
>> --ufuk
>>
>> _______________________________________________
>> ncl-talk mailing list
>> List instructions, subscriber options, unsubscribe:
>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Tue Dec 10 00:35:53 2013
This archive was generated by hypermail 2.1.8 : Fri Dec 13 2013 - 11:39:30 MST