Ignore:
Timestamp:
Feb 24, 2009, 3:46:23 PM (16 years ago)
Author:
lguez
Message:

"comconst.h", and "comgeom2.h" are now both fixed and free form.
Removed calls to procedure "flush".
Corrected kinds of constants which appeared as arguments to "min" or
"max" (all arguments are now of the same type and kind).

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ4/branches/LMDZ4-dev/libf/phylmd/radiation_AR4.F

    r998 r1106  
    21232123      ZALPHA = 3.*ZWCP*ZRM2*(1.+ZGP*ZWM)/ZX2
    21242124      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)
    21272126      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)
    21302128      ZEXKP=EXP(ZARG2)
    21312129      ZEXKM = 1./ZEXKP
     
    28982896      ZTX2=ZTX*ZTX
    28992897      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)
    29022899      ZCAH1=AT(1,1)+ZUP*(AT(1,2)+ZUP*(AT(1,3)))
    29032900      ZCBH1=BT(1,1)+ZUP*(BT(1,2)+ZUP*(BT(1,3)))
     
    29302927      ZZABLY = ZABLY(JL,9,JAE1)+ZABLY(JL,9,JAE2)+ZABLY(JL,9,JAE3)
    29312928      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 )
    29342930      ZPSC2(JL) = (ZTAVI(JL)/TREF) ** ZUP
    29352931      ZCAC8=AT(8,1)+ZUP*(AT(8,2)+ZUP*(AT(8,3)))
Note: See TracChangeset for help on using the changeset viewer.