Ignore:
Timestamp:
Nov 1, 2016, 11:19:45 AM (8 years ago)
Author:
oboucher
Message:

Removing useless global variables
Interpolation only done once a month
Changing test orders in implicit coagulation routine as some compilers didn't like it

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ5/trunk/libf/phylmd/StratAer/cond_evap_tstep_mod.F90

    r2690 r2695  
    1818      USE aerophys
    1919      USE infotrac
    20       USE YOMCST
     20      USE YOMCST, ONLY : RPI
    2121
    2222      IMPLICIT NONE
     
    139139      USE aerophys
    140140      USE infotrac
     141      USE YOMCST, ONLY : RPI
    141142
    142143      IMPLICIT NONE
    143 
    144       include "YOMCST.h"
    145144
    146145      ! input variables
     
    175174        ff(IK,:)=0.0
    176175        DO k=1, nbtr_bin
    177           IF (k.LE.(nbtr_bin-1).AND.Vbin_wet(k).LE.Vnew.AND.Vnew.LT.Vbin_wet(k+1)) THEN
    178             ff(IK,k)= Vbin_wet(k)/Vnew*(Vbin_wet(k+1)-Vnew)/(Vbin_wet(k+1)-Vbin_wet(k))
     176          IF (k.LE.(nbtr_bin-1)) THEN
     177            IF (Vbin_wet(k).LE.Vnew.AND.Vnew.LT.Vbin_wet(k+1)) THEN
     178              ff(IK,k)= Vbin_wet(k)/Vnew*(Vbin_wet(k+1)-Vnew)/(Vbin_wet(k+1)-Vbin_wet(k))
     179            ENDIF
    179180          ENDIF
    180181          IF (k.EQ.1.AND.Vnew.LE.Vbin_wet(k)) THEN
    181182            ff(IK,k)= 1.
    182183          ENDIF
    183           IF (k.GT.1.AND.Vbin_wet(k-1).LT.Vnew.AND.Vnew.LT.Vbin_wet(k)) THEN
    184             ff(IK,k)= 1.-ff(IK,k-1)
     184          IF (k.GT.1) THEN
     185            IF (Vbin_wet(k-1).LT.Vnew.AND.Vnew.LT.Vbin_wet(k)) THEN
     186              ff(IK,k)= 1.-ff(IK,k-1)
     187            ENDIF
    185188          ENDIF
    186189          IF (k.EQ.nbtr_bin.AND.Vnew.GE.Vbin_wet(k)) THEN
Note: See TracChangeset for help on using the changeset viewer.