ubuntu下运行POM2K模型中GRID.f生成正交曲线网格

jiao_yu

平台：ubuntu.14.10
功能：为POM模型生产IC、grid文件
主程序：GRID.f  
头文件：gridcom.txt  gridborder.txt  params.txt
输入数据文件：TS.dat（温度、盐度） TOPO.dat（海底高程） WIND.dat（风）
运行文件：rungrid.txt
用户修改部分：输入数据，rungrid.txt等

操作：
1.将主程序、头文件、输入数据放在同一目录下（如：/home/jiao/文档/debug/test_GRID）。头文件可放在该目录新建的include（如：/home/jiao/文档/debug/test_GRID/include）文件夹下。
2.修改rungrid.txt（可与主程序放在同一目录）。
gfortran -O3 GRID.f -I /home/jiao/文档/debug/test_GRID/include
（将编译时需要用到的头文件路径更改）
3.打开终端（Ctrl+Atl+T）进入目录:cd /home/jiao/文档/debug/test_GRID
4.修改rungrid.txt的权限：chmod +x rungrid.txt （增加可执行权限）
5.运行：./rungrid.txt

生成结果：
IC.dat           POM2K的输入初始文件,formatted IC,grid & constant wind
ijk.dat          IM, JM, KB, Z
plt.dat          formatted grid
bnd.dat          boundary points (gridborder)
wnd.dat          wind velocity
dxy.dat          dx,dx (km)

ijk.dat、plt.dat、bnd.dat、wnd.dat、dxy.dat用于matlab绘图（MATgrid.m）。

GRID.f的说明：
gridcom文件是自定义的边界文件，若在rungrid文件中设置echo '       IGRID=0' >  params  # 0= specified grid, 1= gridborder
则将不用到gridborder文件，自己生成均匀矩形网格，而GRID.f文件中的      
    ASOUTH=32.
        ANORTH=39.
        AWEST=-77.
        AEAST=-71.  是指矩形网格的边界经纬度



C*********************************************************
C                           说明            'gridborder'
C
C---- Selected boundary points. The values I=100 or J= 10     ----
C---- are scaled to I=IM or J=JM respectively.                ----
C---- LINE1: set model resolution                             ----
C---- LINE2: X (longitude) of boundary points                 ----
C---- LINE3: Y (latitude)  of boundary points                 ----
C
C--------------- points on north boundary ------------------------
C----NT、NB、NL、NR中N为number，T为TOP（指绘图时图的顶部），B为Bottom，L为Left，R为Right。
C----ITOP用来指出点的位置（第几个网格），TOPX、TOPY对应为该店的具体经纬度，其余的边界类似
      NT=2
      DATA ITOP(1),ITOP(2)/  1,         100/
      DATA TOPX(1),TOPX(2)/-75.5,      -70.5/
      DATA TOPY(1),TOPY(2)/ 39.5,       39.5/
C
C--------------- points on south boundary ------------------------
      NB=2
      DATA IBOT(1),IBOT(2)/  1,         100/
      DATA BOTX(1),BOTX(2)/-77.5,      -71.5/
      DATA BOTY(1),BOTY(2)/ 31.5,       31.5/
C
C--------------- points on west boundary -------------------------
      NL=4
      DATA JLFT(1),JLFT(2),JLFT(3),JLFT(4)/  1,   30,   50,   100/
      DATA LFTX(1),LFTX(2),LFTX(3),LFTX(4)/-77.5,-77.5,-76.5,-75.5/
      DATA LFTY(1),LFTY(2),LFTY(3),LFTY(4)/ 31.5, 34.0, 36.0, 39.5/   
C
C--------------- points on east boundary -------------------------
      NR=2
      DATA JRHT(1),JRHT(2)/  1,         100/
      DATA RHTX(1),RHTX(2)/-71.5,      -70.5/
      DATA RHTY(1),RHTY(2)/ 31.5,       39.5/
C
      RHTX(NR)=TOPX(NT)
      RHTY(NR)=TOPY(NT)
C*********************************************************

生成的bnd.dat说明： 由GIRD.f中CALL BORDER(IM,JM,X,Y)生成，输出格式'2F10.3'
TOPX   TOPY  NT个
BOTX   BOTY  NB个
LEFX   TEFY  NL个
RHTX   RHTY  NR个

cibaouc

很不错，谢谢楼主

ningm3

楼主你好~因为最近刚开始学习POM，我运行完了./rungrid.txt和./runpom2k.txt之后，得到的是fort.71，不是nc格式的文件，过程也没用出现错误，请问最后的输出文件就是fort.71吗
谢谢！

Lionhart

学习了。。感谢楼主

lovewwt

ningm3 发表于 2016-3-30 13:53
楼主你好~因为最近刚开始学习POM，我运行完了./rungrid.txt和./runpom2k.txt之后，得到的是fort.71，不是nc ...

楼主最后怎么解决的呀，菜鸟一只也遇到这个问题了

fourteen

ningm3 发表于 2016-3-30 13:53
楼主你好~因为最近刚开始学习POM，我运行完了./rungrid.txt和./runpom2k.txt之后，得到的是fort.71，不是nc ...

解决了吗？还需要把pom.n里面做一些更改，将fort.71转为.nc文件