*/=======================================================================
*/ Correction set for LEM version 2.3.2
*/ Created: 6th May 2005
*/ Author: Ben Shipway
*/ Last Modified: 1st August 2005
*/ Modification history:
*/         6/5/2005 - Added Chris's modification to the time filter
*/         6/5/2005 - Added Adrian's correction to microphysics
*/         1/8/2005 - Added Ben's array bounds modification
*/=======================================================================
*/
*/
*/=======================================================================
*/=======================================================================
*/ LEM Version 2.3.2 Modset 
*/ ========================
*/
*/ Purpose:
*/ -------
*/ Allows forward stepping of scalars to be combined with
*/ leapfrog stepping of wind.
*/
*/ Method:
*/ -------
*/ Only apply time filter to scalars when they are integrated
*/ with the leapfrog scheme, irrespective of what timestepping
*/ scheme is used for wind components.
*/
*/ Chris Smith (APP) 25th April 2005
*IDENT FORSCAL
*DECLARE STEPFLDS
*D STEPFLDS.121
       IF (ISMTHP.EQ.1 .AND. IFORSCALP.EQ.0) THEN
*D STEPFLDS.133
       IF (ISMTHP.EQ.1 .AND. IFORSCALP.EQ.0) THEN
*D STEPFLDS.142
       IF (ISMTHP.EQ.1 .AND. IFORSCALP.EQ.0) THEN
*D STEPFLDS.156
       IF (ISMTHP.EQ.1 .AND. IFORSCALP.EQ.0) THEN
*D STEPFLDS.169
       IF (ISMTHP.EQ.1 .AND. IFORSCALP.EQ.0) THEN
*D STEPFLDS.178
       IF (ISMTHP.EQ.1 .AND. IFORSCALP.EQ.0) THEN
*DECLARE SWAPSMTH
*I SWAPSMTH.61
     & ,TSMTH_SCAL
*I SWAPSMTH.65
        IF (IFORSCALP.EQ.0) THEN
          TSMTH_SCAL=TSMTH
        ELSE
          TSMTH_SCAL=0.0
        ENDIF
*D SWAPSMTH.99
         ZTH(J,K)=TH(J,K)+TSMTH_SCAL*TH_NP1    ! finish off t-level n
*D SWAPSMTH.109
         ZQ(J,K,IQ)=Q(J,K,IQ)+TSMTH_SCAL*Q_NP1   ! finish off t-level n
*D SWAPSMTH.144
          OLZTHBAR(K)=OLTHBAR(K)+TSMTH_SCAL*TH_NP1
*D SWAPSMTH.153
          OLZQBAR(K,IQ)=OLQBAR(K,IQ)+TSMTH_SCAL*Q_NP1
*/
*/ End of modset.
*/=======================================================================
*/=======================================================================
*/
*/=======================================================================
*/=======================================================================
*IDENT MICFIX
*// Fix to stop overflows if QL has tiny values.
*// On the SX6 I've seen QL=1.E-150 which when raised to a power gives 
*// silly answers (eg. -1E+212)
*// 
*I NEWMICRO.837
          IF ( QLjk .LT. SMALLP) QLjk = 0.0
*/=======================================================================
*/=======================================================================
*/=======================================================================
*/=======================================================================
*//Prevents array bounds problems when compiling with ifc
*// Ben Shipway 6/7/2005
*ID PUDCORR
*D DIAG_NC.1881
            TEMP_JJP(1:((JJP-1)*IMICROP+1))
     &           =PUDDLE(IMICROP:JJP*IMICROP,NMET,IMICROP)
*/=======================================================================
*/=======================================================================