湿位涡和假相当位温相同？？？》

shuling216

湿位涡和假相当位温相同？？？》？？？？

shuling216

湿位涡和假相当位温相同？？？

shuling216

湿位涡和假相当位温相同？？？

shuling216

湿位涡和假相当位温相同？？？

shuling216

湿位涡和假相当位温相同？？？

言深深

此前有人传的一个fortran程序，有计算的，参考如下：

       SUBROUTINE CAL_THSE(P,T,TD,LEV,THSE) !! 计算假相当位温C     THIS SUBROUTINE CALCULATE PSEUDO-EQUIVALENT POTENTIAL TEMPERATURE
C        AND IT IS USED WITH CAL_TC_PC,CAL_TA,CAL_TE,CAL_ENGY,
C        IN WHICH SUBROUTINE CAL_TA AND CAL_TC_PC ARE USED WITH CAL_THSE
C    INPUT   T ======> TEMPERATURE (CELSIUS DEGREE)
C            TD =====> DEW POINT (CELSIUS DEGREE)
C            P ======> PRESSURE   (HPA)
C            LEV ====> LEVEL
C    OUTPUT  THSE ===> PSEUDO EQUIVALENT POTENTIAL TEMPERATURE   (KELVIN)
C    TEMP ARRAY
      DIMENSION P(LEV),T(LEV),TD(LEV),THSE(LEV)
      REAL A,B,THSE,TCM,PCM,ECM
        DO K=1,LEV
                IF(T(K).LT.80.AND.TD(K).LT.80)THEN
                    TCM=T(K)-(T(K)-TD(K))*0.976/(0.976-0.000833
        1            *(237.3+TD(K))**2/(273.15+TD(K)))
              PCM=P(K)*((273.15+TCM)/(273.15+T(K)))**3.5
              IF(TCM.GE.0)A=7.5*TCM/(237.3+TCM)
              IF(TCM.LT.0)A=9.5*TCM/(265.5+TCM)
              ECM=6.112*10.**A
              B=1000.0/(PCM-ECM)
!        WRITE(*,*)PCM
              THSE(K)=(273.16+TCM)*B**0.286*EXP((2.5*10.**6-2368.*TCM)
     1                *0.622*ECM/(1004.*(273.16+TCM)*(PCM-ECM)))
            ELSE
                THSE(K)=9999.0
            ENDIF
        ENDDO
        RETURN
        END

原帖地址记不得了，现将完整程序粘贴如下：

Cp.f (69.38 KB, 下载次数: 1)

2017-3-6 09:34 上传

点击文件名下载附件

言深深

另外两个帖子重复，不必要发那么多，因此删掉了。