Fortran求兰伯特投影地图放大系数

？？/

13top

lz可以参考一下数值天气预报 这本书或者相关课件，弄懂了公式就可以用Fortran编程了

yumi-zihan

可以参考一下这个，思路是先求出坐标点到北极点的距离，再利用兰勃特投影地图放大计算系数公式计算
program test1
implicit none

!输入坐标值
real xq,yq,xw,yw,d,lq,lw,k,le,m1,m2,f1,f2,a,s1,s2,omega
integer xo,yo
write(*,*)'输入Q的横纵坐标'
read*,xq,yq
write(*,*)'输入W的横纵坐标'
read*,xw,yw
write(*,*)'输入网格距'
read*,d

!北极点投影 O 的坐标
xo=-10
yo=4

lq=sqrt((xq-xo)**2+(yq-yo)**2)*d
lw=sqrt((xw-xo)**2+(yw-yo)**2)*d

print*,"Q点到北极点的距离为:",lq,"千米"
print*,"W点到北极点的距离为:",lw,"千米"

k=0.7156
le=11423.37
a=6371.0
omega=7.292E-05
print*,"Lambert投影的圆锥系数为",k,"赤道与北极点的距离为",le,"千米"
print*,"地球半径为",a,"千米"

!Q点地图放大系数和科氏力参数
s1=(le**(2/k)-lq**(2/k))/(le**(2/k)+lq**(2/k))
m1=k*lq/(a*sqrt(1-s1**2))
f1=2*omega*s1

!W点地图放大系数和科氏力参数
s2=(le**(2/k)-lw**(2/k))/(le**(2/k)+lw**(2/k))
m2=k*lw/(a*sqrt(1-s2**2))
f2=2*omega*s2

print*,"Q点地图放大系数为",m1,"科氏力参数为",f1
print*,"W点地图放大系数为",m2,"科氏力参数为",f2

end