爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

搜索
查看: 10276|回复: 1

FVCOM4.3 遇到的问题

[复制链接]
发表于 2020-3-13 20:51:30 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册

x
I encountered some problems when running FVCOM4.3:

1952724-20200226131936834-355417528.png

I experimented with two examples (River Plume and channel) .

River plume works, and Estuary reports an error.

River Plume does not have OBC_ELEVATION_FILE, and Channel have a obc file (idealised_channel_02_spectide_M2.nc).

Is there a problem with open boundary conditions?

I experimented with my example, and I have this hint:

1952724-20200226132012065-283782510.png

Should I add the 'components' to the OBC file?

My commonly used version is FVCOM3.x.

Here is some information about my computer:
system:  centos 6.6
compiler: intel 2013
mpi:    openmpi-1.6.5

我认为的两个值得尝试的地方:
(1)查看源码,看看格式,把components加入到开边界文件里试一试
(2)我的开边界条件用的是各个分潮的振幅迟角,换成水位的时间序列能不能用?


密码修改失败请联系微信:mofangbao
 楼主| 发表于 2020-4-24 12:23:34 | 显示全部楼层
来填坑了

在FVCOM_source/mod_force.F里
把line 1795行: ! tide_elv_p%curr_stkcnt = 0; tide_elv_n%curr_stkcnt = 0 !!yw 0311
移动到
line 15164行:

CASE(TIDE_FORCING_TIMESERIES)

       FTM => TIDE_FILE%FTIME

       ! PRECIPITATION
       CALL UPDATE_VAR_BRACKET(TIDE_FILE,TIDE_ELV_P,TIDE_ELV_N,NOW,STATUS)
       IF (STATUS /= 0) THEN
          CALL FATAL_ERROR("COULD NOT UPATE TIDE ELVATION BRACKET: BOUNDS EXCEEDED?")
       end if

       CALL NC_POINT_VAR(TIDE_ELV_N,VNP)
       CALL NC_POINT_VAR(TIDE_ELV_P,VPP)   
       BND_ELV = FTM%NEXT_WGHT * VNP + FTM%PREV_WGHT * VPP
       tide_elv_p%curr_stkcnt = 0; tide_elv_n%curr_stkcnt = 0 !yw 0311

因为两个指针在这里赋值
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表