Hi,
When I was trying to understand the effects of ocean currents velocity on ICEPACK simulations, I found myself in an awkward situation that no matter what values of uocn and vocn I give to icepack as the oceanic forcing, they seem to produce the same results of, for example, sea ice thickness. I then read the ICEPACK codes and try to understand the reason behind it and I found out that fbot = cpchr * deltaT * ustar; fbot = max (fbot, frzmlt), in which, ustar = sqrt (sqrt(strocnxT**2+strocnyT**2)/rhow); ustar = max (ustar,ustar_min); My problem is that from my understanding, different uocn and vocn should lead to different strocnxT and strocnyT. However, it seems to me that in the ICEPACK codes, the strocnxT and strocnyT are not calculated explicitly? Are they prescribed in the codes or did I missing something? Thank you very much for your help!
Best,
Shengzhe
When I was trying to understand the effects of ocean currents velocity on ICEPACK simulations, I found myself in an awkward situation that no matter what values of uocn and vocn I give to icepack as the oceanic forcing, they seem to produce the same results of, for example, sea ice thickness. I then read the ICEPACK codes and try to understand the reason behind it and I found out that fbot = cpchr * deltaT * ustar; fbot = max (fbot, frzmlt), in which, ustar = sqrt (sqrt(strocnxT**2+strocnyT**2)/rhow); ustar = max (ustar,ustar_min); My problem is that from my understanding, different uocn and vocn should lead to different strocnxT and strocnyT. However, it seems to me that in the ICEPACK codes, the strocnxT and strocnyT are not calculated explicitly? Are they prescribed in the codes or did I missing something? Thank you very much for your help!
Best,
Shengzhe