/[lmdze]/trunk/Sources/phylmd/Radlwsw/lwb.f
ViewVC logotype

Diff of /trunk/Sources/phylmd/Radlwsw/lwb.f

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 154 by guez, Wed Apr 29 15:47:56 2015 UTC revision 155 by guez, Wed Jul 8 17:03:45 2015 UTC
# Line 1260  Line 1260 
1260  !*         1.0     PLANCK FUNCTIONS AND GRADIENTS  !*         1.0     PLANCK FUNCTIONS AND GRADIENTS
1261  !                  ------------------------------  !                  ------------------------------
1262  !  !
  100  CONTINUE  
 !  
1263        DO 102 JK = 1 , KFLEV+1        DO 102 JK = 1 , KFLEV+1
1264        DO 101 JL = 1, KDLON        DO 101 JL = 1, KDLON
1265        PBINT(JL,JK) = 0.        PBINT(JL,JK) = 0.
# Line 1277  Line 1275 
1275  !*         1.1   LEVELS FROM SURFACE TO KFLEV  !*         1.1   LEVELS FROM SURFACE TO KFLEV
1276  !                ----------------------------  !                ----------------------------
1277  !  !
  110  CONTINUE  
 !  
1278        DO 112 JK = 1 , KFLEV        DO 112 JK = 1 , KFLEV
1279        DO 111 JL = 1, KDLON        DO 111 JL = 1, KDLON
1280        ZTI(JL)=(PTL(JL,JK)-TSTAND)/TSTAND        ZTI(JL)=(PTL(JL,JK)-TSTAND)/TSTAND
# Line 1300  Line 1296 
1296  !*         1.2   TOP OF THE ATMOSPHERE AND SURFACE  !*         1.2   TOP OF THE ATMOSPHERE AND SURFACE
1297  !                ---------------------------------  !                ---------------------------------
1298  !  !
  120  CONTINUE  
 !  
1299        DO 121 JL = 1, KDLON        DO 121 JL = 1, KDLON
1300        ZTI(JL)=(PTL(JL,KFLEV+1)-TSTAND)/TSTAND        ZTI(JL)=(PTL(JL,KFLEV+1)-TSTAND)/TSTAND
1301        ZTI2(JL) = (PTL(JL,1) + PDT0(JL) - TSTAND) / TSTAND        ZTI2(JL) = (PTL(JL,1) + PDT0(JL) - TSTAND) / TSTAND
# Line 1323  Line 1317 
1317  !*         1.3   GRADIENTS IN SUB-LAYERS  !*         1.3   GRADIENTS IN SUB-LAYERS
1318  !                -----------------------  !                -----------------------
1319  !  !
  130  CONTINUE  
 !  
1320        DO 132 JK = 1 , KFLEV        DO 132 JK = 1 , KFLEV
1321        JK2 = 2 * JK        JK2 = 2 * JK
1322        JK1 = JK2 - 1        JK1 = JK2 - 1
# Line 1339  Line 1331 
1331  !*         2.0   CHOOSE THE RELEVANT SETS OF PADE APPROXIMANTS  !*         2.0   CHOOSE THE RELEVANT SETS OF PADE APPROXIMANTS
1332  !                ---------------------------------------------  !                ---------------------------------------------
1333  !  !
  200  CONTINUE  
 !  
 !  
  210  CONTINUE  
 !  
1334        DO 211 JL=1, KDLON        DO 211 JL=1, KDLON
1335        ZDSTO1 = (PTL(JL,KFLEV+1)-TINTP(1)) / TSTP        ZDSTO1 = (PTL(JL,KFLEV+1)-TINTP(1)) / TSTP
1336        IXTOX = MAX( 1, MIN( MXIXT, INT( ZDSTO1 + 1. ) ) )        IXTOX = MAX( 1, MIN( MXIXT, INT( ZDSTO1 + 1. ) ) )
# Line 1378  Line 1365 
1365   213  CONTINUE   213  CONTINUE
1366   214  CONTINUE   214  CONTINUE
1367  !  !
  220  CONTINUE  
 !  
1368        DO 225 JK=1,KFLEV        DO 225 JK=1,KFLEV
1369        DO 221 JL=1, KDLON        DO 221 JL=1, KDLON
1370        ZDST1 = (PTAVE(JL,JK)-TINTP(1)) / TSTP        ZDST1 = (PTAVE(JL,JK)-TINTP(1)) / TSTP
Legend:
Removed from v.154  
changed lines
  Added in v.155
  ViewVC Help
Powered by ViewVC 1.1.21