计算大气热源Q1遇到的问题

良辰

dx=cos((-87.5+2.5*(j-1))*RAD)*RR*RAD*2.5*2没太理解。弧长等于弧度乘以半径吧。dx=cos((-87.5+2.5*(j-1))*RAD)*RR，就可以了吧，不太懂，求指教

LiuJiawei

良辰 发表于 2018-3-8 21:56
dx=cos((-87.5+2.5*(j-1))*RAD)*RR*RAD*2.5*2没太理解。弧长等于弧度乘以半径吧。dx=cos((-87.5+2.5*(j-1)) ...

RAD*2.5*2是弧度，cos((-87.5+2.5*(j-1))*RAD)*RR是不同维度的半径。

良辰

LiuJiawei 发表于 2018-3-10 16:25
RAD*2.5*2是弧度，cos((-87.5+2.5*(j-1))*RAD)*RR是不同维度的半径。

谢谢，理解了！您的加热率的单位是多少呀？

LiuJiawei

良辰 发表于 2018-3-11 09:23
谢谢，理解了！您的加热率的单位是多少呀？

整层的是W/m^2

xgm324

正需要该程序，学习

Miko_

您好，我也用这个程序进行了计算，感觉最后算出来的值相比于文献，是其三倍左右我看您也提到了这个问题，请问后来解决了吗？

半生铁马半生花

Miko_ 发表于 2018-4-17 16:05
您好，我也用这个程序进行了计算，感觉最后算出来的值相比于文献，是其三倍左右我看您也提到了这 ...

请问您后来解决了吗

雪寒月明

厉害厉害哦，楼主不错。

爆浆炸鸡

楼主用的欧洲中心的资料，w是m/s的单位，但是原fortran程序是omega,应该是pa/s的单位，不知道楼主后来怎么解决的？

清宵逝水

1. ug(i,j,l)=(u(i,j,l)*(t(i+1,j,l)-t(i-1,j,l))/dx +(v(i,j,l)*(t(i,j+1,l)-t(i,j-1,l))))/dy
   
2. 
   

ug(i,j,l)=(u(i,j,l)*(t(i+1,j,l)-t(i-1,j,l))/dx +(v(i,j,l)*(t(i,j+1,l)-t(i,j-1,l))))/dy<br /> <br />

我下载的程序中，平流项都是这么写的。。。
dy难道不是写在里面？？

1. (u(i,j,l)*(t(i+1,j,l)-t(i-1,j,l))/dx + (v(i,j,l)*(t(i,j+1,l)-t(i,j-1,l))/dy))

(u(i,j,l)*(t(i+1,j,l)-t(i-1,j,l))/dx + (v(i,j,l)*(t(i,j+1,l)-t(i,j-1,l))/dy))