【分享贴】把ncl中的gc_inout改写成f90子程序

立方光年的团子

ncl中的gc_inout：判定目标点在闭合曲线内（含边界上）/外

官网：
F77源码包含在ncl安装包的这个文件中：
/ncl_ncarg-6.6.2/ni/src/lib/nfpfort/sg_tools.f


提取与gcinout有关的子函数，并改写成f90子程序

使用示例：

1. !--------------------------------------------------
   
2. !-------------- USAGE -----------------------------
   
3. !--------------------------------------------------
   
4. 
   
5. !----GCINOUT=0   : in
   
6. !----GCINOUT=1   : out
   
7. !----GCINOUT=-1  : error
   
8. 
   
9. 
   
10. include "./gc_inout.f90"
    
11. program ex1
    
12.     implicit none
    
13.     integer,parameter :: NPTS=5
    
14.     real(kind=8) :: blat(NPTS),blon(NPTS)
    
15.     real(kind=8) :: plat,plon
    
16.     real(kind=8) :: WORK(4,NPTS)
    
17.     integer GCINOUT
    
18. 
    
19.     blat=DBLE((/20,20,40,40,20/))         !边界：闭合曲线（头尾同一点）注意双精度
    
20.     blon=DBLE((/100,120,120,110,100/))   
    
21. 
    
22.     plat=DBLE(30)                         !目标点
    
23.     plon=DBLE(110)
    
24. 
    
25.     print*,GCINOUT(plat, plon, blat, blon, 5, WORK)
    
26. 
    
27.     stop
    
28. end
    

!--------------------------------------------------<br /> !-------------- USAGE -----------------------------<br /> !--------------------------------------------------<br /> <br /> !----GCINOUT=0   : in<br /> !----GCINOUT=1   : out<br /> !----GCINOUT=-1  : error<br /> <br /> <br /> include "./gc_inout.f90"<br /> program ex1<br />     implicit none<br />     integer,parameter :: NPTS=5<br />     real(kind=8) :: blat(NPTS),blon(NPTS)<br />     real(kind=8) :: plat,plon<br />     real(kind=8) :: WORK(4,NPTS)<br />     integer GCINOUT<br /> <br />     blat=DBLE((/20,20,40,40,20/))         !边界：闭合曲线（头尾同一点）注意双精度<br />     blon=DBLE((/100,120,120,110,100/))    <br /> <br />     plat=DBLE(30)                         !目标点<br />     plon=DBLE(110)<br /> <br />     print*,GCINOUT(plat, plon, blat, blon, 5, WORK)<br /> <br />     stop <br /> end<br />

注意事项：
1.边界和目标点都需要双精度
2.边界需闭合，若想看在环装区域的内/外，可使用两次gcinout函数来实现


附件包含
1.打包以下两个文件，由于付出了时间和精力，想要1贡献
①.改写后的可以直接使用的f90子程序 gc_inout.f90
②.源码（F77）sg_tools.f
gcinout.zip (10.81 KB, 下载次数: 3, 售价: 1 贡献)

2021-5-20 10:54 上传

点击文件名下载附件


2.或者也可以自己改写，无需贡献
sg_tools.f (25.6 KB, 下载次数: 1)

2021-5-20 10:54 上传

点击文件名下载附件

【补充原理】如有错误，欢迎指正和补充通过计算目标点P，和边界点Q(i)和Q(i+1)组成的圆弧的弧度之和的绝对值
（1）计算弧度：余弦定理


（2）进行判断点在闭合边界外还是内。


【渔】
转载自https://www.cnblogs.com/jiangleads/p/9461143.html
查找某字符串在某文件下某文本的位置
find ./ncl_ncarg-6.6.2 -name '*.f' -exec grep -sl -i 'gcinout' {} \;

stcya

谢谢分享