想不通啊想不通

onenuister

我又来请教各位大佬了
因论文需要，正在画超前滞后系数图，Fotran计算程序如下：
program main
    implicit none
    integer k
    parameter(k=40)
    integer i,j,ii,it
    real ESWJ_u(k,6),ESWJ_v(k,6),PWJ_u(k,17),PWJ_v(k,17)
    real ESWJ(k,6),PWJ(k,17),z_ESWJ(k,6),z_PWJ(k,17)!!!1-6为6月1-6侯，7-12为7月1-6侯，13-17为8月1-5侯
    real ave,s,rr,r(6,12),t(6,12)
    open(10,file='f:\002\04\ESWJ_uwnd1979-2018.grd',form='binary')
    open(11,file='f:\002\04\ESWJ_vwnd1979-2018.grd',form='binary')
    do i=1,6
        read(10)(ESWJ_u(it,i),it=1,k)
        read(11)(ESWJ_v(it,i),it=1,k)
    enddo
    do i=1,6
        do it=1,k
            ESWJ(it,i)=sqrt(ESWJ_u(it,i)**2+ESWJ_v(it,i)**2)
        enddo
    enddo
    open(12,file='f:\002\04\PWJ_uwnd1979-2018.grd',form='binary')
    open(13,file='f:\002\04\PWJ_vwnd1979-2018.grd',form='binary')
    do i=1,17
        read(12)(PWJ_u(it,i),it=1,k)
        read(13)(PWJ_v(it,i),it=1,k)
    enddo
    do i=1,17
        do it=1,k
            PWJ(it,i)=sqrt(PWJ_u(it,i)**2+PWJ_v(it,i)**2)
        enddo
    enddo
    !!!标准化
    do i=1,6
        ave=0
        do it=1,k
            ave=ave+ESWJ(it,i)
        enddo
        ave=ave/real(k)
        s=0
        do it=1,k
            s=s+(ESWJ(it,i)-ave)**2
        enddo
        s=sqrt(s/real(k))
        do it=1,k
            z_ESWJ(it,i)=(ESWJ(it,i)-ave)/s
        enddo
    enddo
    do i=1,17
        ave=0
        do it=1,k
            ave=ave+PWJ(it,i)
        enddo
        ave=ave/real(k)
        s=0
        do it=1,k
            s=s+(PWJ(it,i)-ave)**2
        enddo
        s=sqrt(s/real(k))
        do it=1,k
            z_PWJ(it,i)=(PWJ(it,i)-ave)/s
        enddo
    enddo
    !!!计算相关系数
    do i=1,6
        ii=1
        do j=i,i+11
            rr=0
            do it=1,k
                rr=rr+z_ESWJ(it,i)*z_PWJ(it,j)
            enddo
            r(i,ii)=rr/real(k)
            t(i,ii)=r(i,ii)*sqrt(real(k-2))/sqrt(1-r(i,ii)**2)
            ii=ii+1
        enddo
    enddo
    open(14,file='f:\002\04\r&t_uv.grd',form='binary')
    write(14)((r(i,j),i=1,6),j=1,12)
    write(14)((t(i,j),i=1,6),j=1,12)
    end

想用grads画一个等值线图，阴影是显著性检验，ctl和gs文件如下：
dset f:\002\04\r&t_uv.grd
undef -999
title corrle
xdef 6 linear 0 1
ydef 12 linear 0 1
zdef 1 linear 300 1
tdef 1 linear 00z1jan1979 1yr
vars 2
r 1 99
t 1 99
endvars


'reinit'
'open f:\002\04\r&t.ctl'
'set grads off'
'set timelab off'
'set grid off'
'set xaxis 1 6 1'
'set yaxis -6 5 1'
'set xlevs 1 2 3 4 5 6'
'set ylevs -6 -5 -4 -3 -2 -1 0 1 2 3 4 5'
'set gxout shaded'
'define_colors'
'set ccols 42 0 62'
'set clevs -2.02 2.02'
'd t'
'set gxout contour'
'set cthick 6'
'set cint 0.1'
'd r'
'draw title correl_uv(SJ,PJ)'
'draw xlab SJ'
'draw ylab t'
'printim f:\002\04\corrle_uv.png x1366 y768 white'
;

结果画出来的图上总有一条折线，就算r和t都用set gxout shaded表示，图上还是存在那条折线，想请教大家问题出在哪里呀？

ljchen1989

额，看错了~楼主请忽略吧

river

你在画图之前加上 'set mpdraw off' 试试

onenuister

river 发表于 2020-3-21 19:14
你在画图之前加上 'set mpdraw off' 试试

太有帮助了，感谢

river

onenuister 发表于 2020-3-21 19:36
太有帮助了，感谢

加上就没有了？那我猜对了，应该是默认画了个不知道是哪的地图······