合成分析，异常年份的提取，以及grsds处理nc格式的环流资料

冬日的阳光

    !----合成分析年份提取子程序-------!
        !input：气象序列x(n);
        !n:样本个数或者年份
        !temp:控制变量
        !yr:起始年份
        !x(n):异常气象要素
        !output:k1,k2为偏多偏少年的年份总数
        !by 冬日的阳光  2014.5.15
         subroutine com(x,n,temp,k1,k2,yr)
                  implicit none
                 real temp
             integer::i,n,k1,k2,yr
                 real,dimension(n)::x
         k1=0
                 k2=0
                 open (1,file='composite_unnormal.txt')
                 write(1,*)"气候异常年份如下："
         do i=1,n
                   if(x(i).GE.temp) then
                   k1=k1+1
                   !print *,"正异常年份",i+yr
                           write(1,*) "正异常年份",i+yr,  "      负异常年份","         无"

                else if(x(i).LE.-temp) then
                   k2=k2+1
                   !print *,"负异常年份",i+yr
                           write(1,*) "正异常年份","        无","         负异常年份",i+yr
                endif
                 enddo
                 !print *,"正异常年份共有",k1,"    负异常年份共有",k2
                 write(1,*) "正异常年份共有",k1,"    负异常年份共有",k2
                 close(1)

        end subroutine com

来自群组: 广东海洋大学大气科学

yysunnyboy

赞一个！！！   

hzgloria

真的很好呢！请问你强年和弱年是如何判定的？

/xin村儿/

write(1,*) "正异常年份",i+yr, 应改为
write(1,*) "正异常年份",i+yr-1,

冬日的阳光

hzgloria 发表于 2015-9-15 08:43
真的很好呢！请问你强年和弱年是如何判定的？

temp控制变量，就是可以设置的强年和弱年分解值

冬日的阳光

/xin村儿/ 发表于 2015-9-15 12:10
write(1,*) "正异常年份",i+yr, 应改为
write(1,*) "正异常年份",i+yr-1,

很抱歉，我可能写错了，等我有时间一下。

明明~~

楼主棒棒哒~~

卡卡的暴风眼

学习中……

孟泽华

有脚本好感动

yhsmith

感谢分享 雪中送炭！！