Changeset 1106 for LMDZ4/branches/LMDZ4-dev/libf/phylmd/radiation_AR4.F
- Timestamp:
- Feb 24, 2009, 3:46:23 PM (16 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ4/branches/LMDZ4-dev/libf/phylmd/radiation_AR4.F
r998 r1106 2123 2123 ZALPHA = 3.*ZWCP*ZRM2*(1.+ZGP*ZWM)/ZX2 2124 2124 ZBETA = 3.*ZWCP* PRMUZ(JL) *(1.+3.*ZGP*ZRM2*ZWM)/ZX2 2125 CMAF ZARG=MIN(ZTOP/PRMUZ(JL),200.) 2126 ZARG=MIN(ZTOP/PRMUZ(JL),2.0d+2) 2125 ZARG=MIN(ZTOP/PRMUZ(JL),200._8) 2127 2126 ZEXMU0=EXP(-ZARG) 2128 CMAF ZARG2=MIN(ZRK*ZTOP,200.) 2129 ZARG2=MIN(ZRK*ZTOP,2.0d+2) 2127 ZARG2=MIN(ZRK*ZTOP,200._8) 2130 2128 ZEXKP=EXP(ZARG2) 2131 2129 ZEXKM = 1./ZEXKP … … 2898 2896 ZTX2=ZTX*ZTX 2899 2897 ZZABLY = ZABLY(JL,6,JAE1)+ZABLY(JL,6,JAE2)+ZABLY(JL,6,JAE3) 2900 CMAF ZUP=MIN( MAX( 0.5*R10E*LOG( ZZABLY ) + 5., 0.), 6.0) 2901 ZUP=MIN( MAX( 0.5*R10E*LOG( ZZABLY ) + 5., 0.d+0), 6.d+0) 2898 ZUP=MIN( MAX( 0.5*R10E*LOG( ZZABLY ) + 5., 0._8), 6._8) 2902 2899 ZCAH1=AT(1,1)+ZUP*(AT(1,2)+ZUP*(AT(1,3))) 2903 2900 ZCBH1=BT(1,1)+ZUP*(BT(1,2)+ZUP*(BT(1,3))) … … 2930 2927 ZZABLY = ZABLY(JL,9,JAE1)+ZABLY(JL,9,JAE2)+ZABLY(JL,9,JAE3) 2931 2928 ZALUP = R10E * LOG ( ZZABLY ) 2932 CMAF ZUP = MAX( 0.0 , 5.0 + 0.5 * ZALUP ) 2933 ZUP = MAX( 0.d+0 , 5.0 + 0.5 * ZALUP ) 2929 ZUP = MAX( 0._8, 5.0 + 0.5 * ZALUP ) 2934 2930 ZPSC2(JL) = (ZTAVI(JL)/TREF) ** ZUP 2935 2931 ZCAC8=AT(8,1)+ZUP*(AT(8,2)+ZUP*(AT(8,3)))
Note: See TracChangeset
for help on using the changeset viewer.