zhanghouxing的个人博客分享 http://blog.sciencenet.cn/u/zhanghouxing

博文

monthly water balance model code

已有 4084 次阅读 2018-7-14 18:08 |系统分类:科研笔记

源代码python编写的,https://github.com/marcae777/Monthly-water-balance-model/blob/master/Monthly_water_balance_model.ipynb

发现这个代码错了,我纠正了,用matlab编写,分享!

data=xlsread('Example_Data.csv');

EP=data(:,3);

P=data(:,2);

Area=862.69;

SC=5000*Area*3.858e-4;

c=8.71259;


 EP=EP*Area*3.858e-4;

 P =Area*P*3.858e-4;

 E=zeros(length(P),1);

 Q=zeros(length(P),1);

 for i=1:length(P)

     if i==1

         S_last=1650*Area*3.858e-4;

     end

     E(i)=c*EP(i)*tanh(P(i)/EP(i));

     Q(i) =(S_last+P(i)-E(i))*tanh((S_last+P(i)-E(i))/SC);

     S_last=S_last+P(i)-E(i)-Q(i);

 end

% xlswrite('output1.xlsx',Q);

     

plot(Q(66:150,:));

hold on;

plot(data(66:150,1)+10,'r');




https://blog.sciencenet.cn/blog-922140-1124050.html

上一篇:DHSVM模型计算流程
下一篇:GSFLOW模型的Codeblocks编译
收藏 IP: 42.91.140.*| 热度|

1 陈昌春

该博文允许注册用户评论 请点击登录 评论 (1 个评论)

数据加载中...
扫一扫,分享此博文

Archiver|手机版|科学网 ( 京ICP备07017567号-12 )

GMT+8, 2024-7-16 05:46

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部