计算字符串运算的结果 如 (10.5*((1+4/2)-5+(9+1)))/2.0

某公司的笔试题,做出来还是花了点时间 用比较笨的方法

/*************************************************************************
    > File Name: numstr_calculate.cpp
    > Author: chunquanL
    > Created Time: 2017-03-22
 ************************************************************************/

#include<iostream>
#include<vector>
#include<cstring>
using namespace std;


//计算字符串运算的结果
//(10.5*((1+4/2)-5+(9+1)))/2.0

int isnum(char ch)
{
    if(ch=='.')
        return 1;
    return ('0'<=ch && '9'>=ch)?1:0;
}
double myatof(const char *str,int *len)
{
    double result=0;
    double base =0.1;
    int length=0;
    int i=0;
    for(int j=0;isnum(str[j]);j++)
        length++;
    *len = length;
    while(i<length && str[i]!='.')
    {
        result*=10;
        result += str[i]-'0';
        i++;
    }
    if(i<length && str[i]=='.')
    {
        i++;
    }
    while(i<length)
    {
        result += (str[i]-'0')*base;
        i++;
        base*=0.1;
    }

    for(int i=0;i<length;i++)
    {

    }
    return result;
}



double getResult_normal(vector<double> &num,vector<char> &sign,int start,int end)
{
    //首先处理乘除法

    int num_of_handle1=0; //乘除运算次数
    int num_of_handle2=0; //加减运算次数

    double result;
    for(int i=start;i<=end;i++)
    {
        if(sign[i]=='*' || sign[i]=='/')
            num_of_handle1++;
        if(sign[i]=='+' || sign[i]=='-')
            num_of_handle2++;
    }
    while(num_of_handle1)
    {
        int i;
        for(i=start;i<=end;i++)
        {
            if(sign[i]=='*' || sign[i]=='/')
                break;
        }
        int before = i;
        int next = i;
        double res=0;
        while(sign[--before]!='n');
        while(sign[++next]!='n');
        if(sign[i] == '*')
            res = num[before]*num[next];
        else if(sign[i]=='/')
            res = num[before]/num[next];
        sign[before]=' ';
        sign[next] = ' ';
        sign[i] = 'n';
        num[i] = res;
        num_of_handle1--;

    }
    while(num_of_handle2)
    {
        int i;
        for(i=start;i<=end;i++)
        {
            if(sign[i]=='+' || sign[i]=='-')
                break;
        }
        int before = i;
        int next = i;
        double res=0;
        while(sign[--before]!='n');
        while(sign[++next]!='n');
        if(sign[i] == '+')
            res = num[before]+num[next];
        else if(sign[i]=='-')
            res = num[before]-num[next];
        sign[before]=' ';
        sign[next] = ' ';
        sign[i] = 'n';
        num[i] = res;
        result = res;
        num_of_handle2--;
    }
    return result;
}

double getResult(char *expr,int length)
{
    if(!expr)
        return 0;

    vector<double> num;
    vector<char> sign;

    for(int i=0;i<length;)
    {
        if(isnum(expr[i]))
        {
            int len=0;
            double res = myatof(&expr[i],&len);
            i+=len;
            num.push_back(res);
            sign.push_back('n');
        }
        else
        {
            sign.push_back(expr[i]);
            i++;
            num.push_back(0.0);
        }
    }
/*
    for(int i=0;i<num.size();i++)
        cout<<num[i]<<"  ";
    cout<<endl;
    for(int i=0;i<sign.size();i++)
        cout<<sign[i]<<" ";
    cout<<endl;
*/
    //先解决括号
    int num_of_kuohao = 0;
    for(int i=0;i<length;i++)
    {
        if(sign[i]=='(')
            num_of_kuohao++;
    }
    if(num_of_kuohao)
    {
        while(num_of_kuohao)
        {
            //找到最后面那个括号的位置
            int before=0;
            int cnt=0;
            for(before=0;before<length;before++)
            {
                if(sign[before]=='(')
                    ++cnt;
                if(cnt==num_of_kuohao)
                    break;
            }
            //找到后面的第一个括号
            int next = before;
            for(;next<length;next++)
            {
                if(sign[next]==')')
                    break;
            }
            getResult_normal(num,sign,before+1,next-1);
            sign[before] = ' ';
            sign[next] =' ';
            num_of_kuohao--;
        }
    }
    getResult_normal(num,sign,0,length-1);
//  cout<<getResult_normal(num,sign,0,num.size()-1);    
    double res=0;
    for(int i=0;i<length;i++)
    {
        if(sign[i]=='n')
            res=num[i];
    }
    return res;
}

int main(void)
{
/*  char a[]="0.12121";
    int len;
    cout<<myatof(a,&len)<<endl;
    cout<<len<<endl;
    cout<<isnum('8');
    cout<<1*1+2*3/4;
*/

    char str[] = "((1.0+2.3)*2)+3/2*(6-5.9)";
    for(int i=0;i<strlen(str);i++)
        cout<<str[i];
    cout<<endl;
    cout<<getResult(str,strlen(str))<<endl;


    return 0;
}


因各个项目中需要使用根据字符串计算数值,这里写出一个算法,专门计算字符串。配有大量常用公式。只有一个人方法,直接调用即可。 类名:CustomMath 函数名:Calculations(string value) 说明:求解算式表达式字符串的值 表达式中包含的符号或函数: truncate, ceiling,floor,round,log10, sign,sinh,sqrt, asin,atan,cosh, tanh, sin,cos,tan ,abs,acos, exp,log,max,min,pow,mod,+,-,*,/,&#39;,&#39;,(,) 函数说明:(不区分大小写) truncate(num) 计算指定数的整数部分 truncate(1.23)=1 ceiling (num) 返回大于或等于指定的双精度浮点数的最小整数值 ceiling(1.23)=2 floor(num) 返回小于或等于指定双精度浮点数的最大整数 floor(1.23)=1 round(num) 将双精度浮点值舍入为最接近的整数值 round(1.23)=1 round(num,num1) 将小数值按指定的小数位数舍入 round(1.23,1)=1.2 log10(num) 返回指定数字以 10 为底的对数 log10(10)=1 sign(num) 返回表示数字符号的值 sign(1.23)=1 sinh(num) 返回指定角度的双曲正弦值 sinh(1.23)=1.5644 sqrt(num) 返回指定数字的平方根 sqrt(9)=3 sqrt(num,num1) 返回指定数字的num1根 sqrt(27,3)=3 asin(num) 返回正弦值为指定数字的角度 asin(0.5)=PI/6 atan(num) 返回正切值为指定数字的角度 atan(1)=45 cosh(num) 返回指定角度的双曲余弦值 cosh(1.23)=1.8567 tanh(num) 返回指定角度的双曲正切值 tanh(1.23)=0.8425 sin(num) 返回指定角度的正弦值 sin(PI/6)=0.5 cos(num) 返回指定角度的余弦值 sin(PI/3)=0.5 tan(num) 返回指定角度的余切值 sin(PI/4)=1 abs(num) 返回数字的绝对值 abs(-12)=12 acos(num) 返回余弦值为指定数字的角度 acos(0.5)=PI/3 exp(num) 返回 e 的指定次幂 exp(1)=2.718 log(num) 返回指定数字的自然对数(底为 e) log(e)=1 log(num,num1) 返回指定数字在使用指定底时的对数 log(e,e)=1 max(num,um1) 返回最大值 max(1,2)=2 min(num,num1) 返回最小值 min(1,2)=1 pow(num,num1) 返回指定数字的指定次幂 pow(2,2)=4 mod(num,num1) 返回余数 mod(3,2)=1 常量: PI 值:3.14159265358979323846 E 值:2.7182818284590452354 YEAR 值:当前年份 MONTH 值:当前月份 DAY 值: 当前日 HOUR 值:当前时 MINUTE 值:当前分 SECOND 值:当前秒 RANDOM 值:一个随机数(0-1 之间) 实例 系统计算1+2*3/4-0.5=2 函数计算1+2*3/4-0.5=2 调用方式:CustomMath.Calculations("1+2*3/4-0.5") 系统计算(1+2)*3/4-0.5=1.75 函数计算(1+2)*3/4-0.5=1.75 调用方式:CustomMath.Calculations("(1+2)*3/4-0.5") 系统计算(sin(pi)+sqrt(3+5*7+(2+8/4*5+2)))/6=1.20185042515466 公式计算(sin(pi)+sqrt(3+5*7+(2+8/4*5+2)))/6=1.20185042515466 调用方式:CustomMath.Calculations("(sin(pi)+sqrt(3+5*7+(2+8/4*5+2)))/6") 系统计算:sin(pow(3,2)/4)+3.5-9*sqrt(81)=-76.7219268031121 函数计算:sin(pow(3,2)/4)+3.5-9*sqrt(81)=-76.7219268031121 调用方式:CustomMath.Calculations("sin(pow(3,2)/4)+3.5-9*sqrt(81)")
! m=20: 纬向格点数 ! n=16: 经向格点数 ! d: 网格距 ! rm: 地图放大系数 ! f: 地转参数 ! w: 工作数组 ! cla,clo: 区域中心纬度和经度 ! dt: 时间步长 ! s: 平滑系数 ! ua,ub,uc:n-1,n,n+1时间层的纬向风速 ! va,vb,vc:n-1,n,n+1时间层的经向风速 ! za,zb,zc:n-1,n,n+1时间层的位势高度 ! na: 控制12小时预报的参数 ! nb: 记录时间积分步数的参数 ! nt2=72: 判别是否积分12小时,是否该做内点平滑; ! nt4=6: 判定是否该做边界平滑; ! nt5=36: 判定是否该做时间平滑。 ! zo: 为了减小重力惯性外波的波速,增加差分格式的稳定性而引入的位势高度。 ! ni: 是否进行初始风场的静力初始化。 ! ni=0为不进行初始化,使用读入的风场和高度场; ! ni=1为进行初始化,需要位势高度场做初值即可。 program Barotropic_Model parameter(m=20,n=16,d=300000.0,cla=51.0,clo=118.0,dt=600.0) dimension ua(m,n),va(m,n),za(m,n),ub(m,n),vb(m,n),zb(m,n),uc(m,n),vc(m,n),zc(m,n),rm(m,n),f(m,n),w(m,n) zo=2500.0 s=0.5 nt2=72 nt4=6 nt5=36 c1=dt/2.0 c2=dt*2.0 write(*,*) write(*,*)&#39;!!!!欢迎使用正压原始方程模式!!!!&#39; write(*,*) write(*,*)&#39;打开文件、读入初始场、准备输出场文件...........&#39; open(11,file=&#39;practice/Input/za.dat&#39;) ! 输入的高度场(文本) open(12,file=&#39;practice/Input/ua.dat&#39;) ! 输入的u风场(文本) open(13,file=&#39;practice/Input/va.dat&#39;) ! 输入的v风场(文本) open(17,file=&#39;practice/Input/za.grd&#39;,form=&#39;binary&#39;) ! 输出的高度场(二进制) open(18,file=&#39;practice/Input/ua.grd&#39;,form=&#39;binary&#39;) ! 输出的u风场(二进制) open(19,file=&#39;practice/Input/va.grd&#39;,form=&#39;binary&#39;) ! 输出的v风场(二进制) open(21,file=&#39;practice/Output/rm.dat&#39;) ! 地图放大系数 open(22,file=&#39;practice/Output/f.dat&#39;) ! 地转参数 open(23,file=&#39;practice/Output/ub.dat&#39;) ! 静力初始化得到的u风场(文本) open(24,file=&#39;practice/Output/ub.grd&#39;,form=&#39;binary&#39;) ! 静力初始化得到的u风场(二进制) open(25,file=&#39;practice/Output/vb.dat&#39;) ! 静力初始化得到的v风场(文本) open(26,file=&#39;practice/Output/vb.grd&#39;,form=&#39;binary&#39;) ! 静力初始化得到的v风场(二进制) open(27,file=&#39;practice/Output/zc.dat&#39;) ! 预报的高度场(文本) open(28,file=&#39;practice/Output/zc.grd&#39;,form=&#39;binary&#39;) ! 预报的高度场(二进制) open(29,file=&#39;practice/Output/uc.dat&#39;) ! 预报的u风场(文本) open(30,file=&#39;practice/Output/uc.grd&#39;,form=&#39;binary&#39;) ! 预报的u风场(二进制) open(31,file=&#39;practice/Output/vc.dat&#39;) ! 预报的v风场(文本) open(32,file=&#39;practice/Output/vc.grd&#39;,form=&#39;binary&#39;) ! 预报的v风场(二进制) ! 读入初始资料场 read(11,&#39;(20f6.0)&#39;)za read(12,&#39;(20f10.5)&#39;)ua read(13,&#39;(20f10.5)&#39;)va write(*,*) write(*,*)&#39;将初始高度场和风场写成二进制文件,便于Grads绘图......&#39; write(17)((za(i,j),i=1,m),j=1,n) write(18)((ua(i,j),i=1,m),j=1,n) write(19)((va(i,j),i=1,m),j=1,n) ! 计算放大系数和地转参数,并写入数据文件中 write(*,*) write(*,*)&#39;计算每个格点上的地图放大系数和地转参数,并写入对应输出文件......&#39; call cmf(rm,f,d,cla,m,n) write(21,&#39;(20f10.5)&#39;)rm write(22,&#39;(20e15.5)&#39;)f write(*,*) write(*,*)&#39;请输入数字01后按回车运行,0表示不进行静力初始化;1表示进行静力初始化。&#39; write(*,*) write(*,*)&#39;注意:如果求地转风的子程序(风场初始化)未完成,则只能输入数字0。&#39; read(*,*)ni if(ni==1)then ! 计算地转风初值 write(*,*)&#39;进行静力初始化,由风场求出高度场.....&#39; call cgw(ua,va,za,rm,f,d,m,n) write(23,&#39;(20f10.5)&#39;)ua write(24)((ua(i,j),i=1,m),j=1,n) write(25,&#39;(20f10.5)&#39;)va write(26)((va(i,j),i=1,m),j=1,n) elseif(ni==0)then write(*,*)&#39;不进行静力初始化使用给出的位势高度场和风场......&#39; else write(*,*)&#39;啊!!!输入了错误字符,请重新运行程序,输入01!&#39; goto 90 endif ! 边值传送子程序 write(*,*) write(*,*)&#39;固定边界条件赋值........&#39; call tbv(ub,vb,zb,ua,va,za,m,n) call tbv(uc,vc,zc,ua,va,za,m,n) ! 开始预报 write(*,*) write(*,*)&#39;开始12小时预报.......&#39; do na=1,2 nb=0 ! 欧拉后差积分1小时 do nn=1,6 call ti(ua,va,za,ua,va,za,ub,vb,zb,rm,f,d,dt,zo,m,n) call ti(ua,va,za,ub,vb,zb,ua,va,za,rm,f,d,dt,zo,m,n) nb=nb+1 enddo ! 边界平滑子程序 call ssbp(za,w,s,m,n) call ssbp(ua,w,s,m,n) call ssbp(va,w,s,m,n) ! 前差积分半步 call ti(ua,va,za,ua,va,za,ub,vb,zb,rm,f,d,c1,zo,m,n) ! 中央差积分半步 call ti(ua,va,za,ub,vb,zb,uc,vc,zc,rm,f,d,dt,zo,m,n) nb=nb+1 ! 数组传送子程序 call ta(ub,vb,zb,uc,vc,zc,m,n) ! 中央差积分一步,共积分11小时 do nn=1,66 call ti(ua,va,za,ub,vb,zb,uc,vc,zc,rm,f,d,c2,zo,m,n) nb=nb+1 ! 打印积分步数,na大循环步,nb小循环步 call pv(na,nb) ! 判断是否积分12小时 if(nb.eq.nt2) go to 80 ! 判断是否做边界平滑 if(nb/nt4*nt4.eq.nb)then call ssbp(zc,w,s,m,n) call ssbp(uc,w,s,m,n) call ssbp(vc,w,s,m,n) else ! 判断是否做时间平滑 if(nb.eq.nt5 .or. nb.eq.nt5+1)then ! 时间平滑子程序 call ts(ua,ub,uc,va,vb,vc,za,zb,zc,s,m,n) else ! 数组传送,为下一轮积分做准备 call ta(ua,va,za,ub,vb,zb,m,n) call ta(ub,vb,zb,uc,vc,zc,m,n) endif endif enddo ! 区域内点平滑 80 call ssip(zc,w,s,m,n,k,2) call ssip(uc,w,s,m,n,k,2) call ssip(vc,w,s,m,n,k,2) ! 打印积分步数 call pv(na,nb) ! 数组传送,为后12小时的积分做准备 call ta(ua,va,za,uc,vc,zc,m,n) enddo ! 存放预报结果 write(*,*) write(*,*)&#39;输出预报结果.......&#39; write(27,&#39;(20f6.0)&#39;) zc write(28) ((zc(i,j),i=1,m),j=1,n) write(29,&#39;(20f10.5)&#39;)uc write(30) ((uc(i,j),i=1,m),j=1,n) write(31,&#39;(20f10.5)&#39;)vc write(32) ((vc(i,j),i=1,m),j=1,n) write(*,*) write(*,*)&#39;!!!预报结束!!!&#39; 90 stop end ! computing map factors and coriolis parameter ! rk: 圆锥常数 ! rlq: 兰勃特投影映像平面上赤道到北极点的距离 ! a: 地球半径 ! sita:标准余纬 ! psx: 区域中心余纬 ! r: 模式中心到北极的距离 subroutine cmf(rm,f,d,cla,m,n) dimension rm(m,n),f(m,n) rk=0.7156 rlq=11423370.0 a=6371000.0 conv=57.29578 w1=2.0/rk sita=30.0/conv psx=(90.0-cla)/conv ! 计算模式中心到北极的距离r cel0=a*sin(sita)/rk cel=(tan(psx/2.0))/(tan(sita/2.0)) r=cel0*cel**rk ! 确定网格坐标原点在地图坐标系中的位置 xi0=-(m-1)/2.0 yj0=r/d+(n-1)/2.0 ! 求各格点至北极点的距离rl,(xj,yi)为模式各格点在地图坐标系中的位置 do i=1,m do j=1,n xi=xi0+(i-1) yj=yj0-(j-1) rl=sqrt(xi2+yj2)*d ! 求放大系数rm和柯氏参数f w2=(rl/rlq)**w1 sinl=(1.0-w2)/(1.0+w2) rm(i,j)=rk*rl/(a*sqrt(1.0-sinl**2)) f(i,j)=1.4584e-4*sinl enddo enddo return end !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! ! computing geostrophic winds ! 请同学编写地转风初值的子程序!!!应用书中(4.134)式 subroutine cgw(ua,va,za,rm,f,d,m,n) dimension ua(m,n),va(m,n),za(m,n),f(m,n),rm(m,n) end !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! ! time integrations subroutine ti(ua,va,za,ub,vb,zb,uc,vc,zc,rm,f,d,dt,zo,m,n) dimension ua(m,n),va(m,n),za(m,n),ub(m,n),vb(m,n),zb(m,n),uc(m,n),vc(m,n),zc(m,n),rm(m,n),f(m,n) c=0.25/d m1=m-1 n1=n-1 do i=2,m1 do j=2,n1 e=-c*rm(i,j)*((ub(i+1,j)+ub(i,j))*(ub(i+1,j)-ub(i,j))+(ub(i,j)+ub(i-1,j))*(ub(i,j)-ub(i-1,j)) & +(vb(I,j-1)+vb(i,j))*(ub(i,j)-ub(i,j-1))+(vb(I,j)+vb(i,j+1))*(ub(i,j+1)-ub(i,j)) & +19.6*(zb(i+1,j)-zb(i-1,j)))+f(i,j)*vb(i,j) uc(i,j)=ua(i,j)+e*dt g=-c*rm(i,j)*((ub(I+1,j)+ub(i,j))*(vb(i+1,j)-vb(i,j))+(ub(I,j)+ub(i-1,j))*(vb(i,j)-vb(i-1,j)) & +(vb(I,j-1)+vb(i,j))*(vb(i,j)-vb(i,j-1))+(vb(I,j)+vb(i,j+1))*(vb(i,j+1)-vb(i,j)) & +19.6*(zb(i,j+1)-zb(i,j-1)))-f(i,j)*ub(i,j) vc(i,j)=va(i,j)+g*dt enddo enddo do i=2,m1 do j=2,n1 h=-c*rm(i,j)*((ub(I+1,j)+ub(i,j))*(zb(i+1,j)-zb(i,j))+(ub(I,j)+ub(i-1,j))*(zb(i,j)-zb(i-1,j)) & +(vb(I,j-1)+vb(i,j))*(zb(i,j)-zb(i,j-1))+(vb(I,j)+vb(i,j+1))*(zb(i,j+1)-zb(i,j)) & +2.0*(zb(i,j)-zo)*(ub(i+1,j)-ub(i-1,j)+vb(i,j+1)-vb(i,j-1))) zc(i,j)=za(i,j)+h*dt enddo enddo return end ! time smoothimg subroutine ts(ua,ub,uc,va,vb,vc,za,zb,zc,s,m,n) dimension ua(m,n),va(m,n),za(m,n),ub(m,n),vb(m,n),zb(m,n),uc(m,n),vc(m,n),zc(m,n) m1=m-1 n1=n-1 do i=2,m1 do j=2,n1 ub(i,j)=ub(i,j)+s*(ua(i,j)+uc(i,j)-2.0*ub(i,j))/2.0 vb(i,j)=vb(i,j)+s*(va(i,j)+vc(i,j)-2.0*vb(i,j))/2.0 zb(i,j)=zb(i,j)+s*(za(i,j)+zc(i,j)-2.0*zb(i,j))/2.0 enddo enddo return end !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! ! space smoothing for internal points 区域内5点平滑(正逆平滑) ! 请同学编写区域内5点平滑(正逆平滑)的子程序!!!应用书中(4.126)式 ! l=1为只执行正平滑,l=2为执行正逆平滑 subroutine ssip(a,w,s,m,n,k,l) dimension a(m,n),w(m,n) end !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! ! space smoothing for boundary points 边界九点平滑 subroutine ssbp(a,w,s,m,n) dimension a(m,n),w(m,n) m1=m-1 m3=m-3 n1=n-1 n2=n-2 n3=n-3 do i=2,m1 do j=2,n1,n3 w(i,j)=a(i,j)+0.5*s*(1.0-s)*(a(i-1,j)+a(i+1,j)+a(i,j-1)+a(i,j+1)-4.0*a(i,j)) & +0.25*s*s*(a(i-1,j-1)+a(i-1,j+1)+a(i+1,j-1)+a(i+1,j+1)-4.0*a(i,j)) enddo enddo do i=2,m1,m3 do j=3,n2 w(i,j)=a(i,j)+0.5*s*(1.0-s)*(a(i-1,j)+a(i+1,j)+a(i,j-1)+a(i,j+1)-4.0*a(i,j)) & +0.25*s*s*(a(i-1,j-1)+a(i-1,j+1)+a(i+1,j-1)+a(i+1,j+1)-4.0*a(i,j)) enddo enddo do i=2,m1 do j=2,n1,n3 a(i,j)=w(i,j) enddo enddo do i=2,m1,m3 do j=3,n2 a(i,j)=w(i,j) enddo enddo return end ! transmiting arrays 数组传送 subroutine ta(ua,va,za,ub,vb,zb,m,n) dimension ua(m,n),va(m,n),za(m,n),ub(m,n),vb(m,n),zb(m,n) do i=1,m do j=1,n ua(i,j)=ub(i,j) va(i,j)=vb(i,j) za(i,j)=zb(i,j) enddo enddo return end ! transmiting boundary valaus 赋固定边界值 subroutine tbv(ua,va,za,ub,vb,zb,m,n) dimension ua(m,n),va(m,n),za(m,n),ub(m,n),vb(m,n),zb(m,n) m1=m-1 n1=n-1 do i=1,m do j=1,n,n1 ua(i,j)=ub(i,j) va(i,j)=vb(i,j) za(i,j)=zb(i,j) enddo enddo do i=1,m,m1 do j=1,n ua(i,j)=ub(i,j) va(i,j)=vb(i,j) za(i,j)=zb(i,j) enddo enddo return end ! printing variables 打印积分步数 subroutine pv(na,nb) write(*,&#39;(5x,3hna=,i3,5x,3hnb=,i2/)&#39;)na,nb return end 将以上Fortran代码改写为对应python代码,同时编写五点平滑子程序、地转风初值子程序, 其中dat文件存放于D:\HuaweiMoveData\Users\80704\Desktop\实习内容中
最新发布
05-17
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值